1
|
Miller JT, Clark BW, Reid NM, Karchner SI, Roach JL, Hahn ME, Nacci D, Whitehead A. Independently evolved pollution resistance in four killifish populations is largely explained by few variants of large effect. Evol Appl 2024; 17:e13648. [PMID: 38293268 PMCID: PMC10824703 DOI: 10.1111/eva.13648] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2023] [Revised: 01/02/2024] [Accepted: 01/09/2024] [Indexed: 02/01/2024] Open
Abstract
The genetic architecture of phenotypic traits can affect the mode and tempo of trait evolution. Human-altered environments can impose strong natural selection, where successful evolutionary adaptation requires swift and large phenotypic shifts. In these scenarios, theory predicts that adaptation is due to a few adaptive variants of large effect, but empirical studies that have revealed the genetic architecture of rapidly evolved phenotypes are rare, especially for populations inhabiting polluted environments. Fundulus killifish have repeatedly evolved adaptive resistance to extreme pollution in urban estuaries. Prior studies, including genome scans for signatures of natural selection, have revealed some of the genes and pathways important for evolved pollution resistance, and provide context for the genotype-phenotype association studies reported here. We created multiple quantitative trait locus (QTL) mapping families using progenitors from four different resistant populations, and using RAD-seq genetically mapped variation in sensitivity (developmental perturbations) following embryonic exposure to a model toxicant PCB-126. We found that one to two large-effect QTL loci accounted for resistance to PCB-mediated developmental toxicity. QTLs harbored candidate genes that govern the regulation of aryl hydrocarbon receptor (AHR) signaling. One QTL locus was shared across all populations and another was shared across three populations. One QTL locus showed strong signatures of recent natural selection in the corresponding wild population but another QTL locus did not. Some candidate genes for PCB resistance inferred from genome scans in wild populations were identified as QTL, but some key candidate genes were not. We conclude that rapidly evolved resistance to the developmental defects normally caused by PCB-126 is governed by few genes of large effect. However, other aspects of resistance beyond developmental phenotypes may be governed by additional loci, such that comprehensive resistance to PCB-126, and to the mixtures of chemicals that distinguish urban estuaries more broadly, may be more genetically complex.
Collapse
Affiliation(s)
- Jeffrey T. Miller
- Department of Environmental Toxicology, Center for Population Biology, Coastal and Marine Sciences InstituteUniversity of California, DavisDavisCaliforniaUSA
- Present address:
Molecular, Cellular, and Biomedical SciencesUniversity of New HampshireDurhamNew HampshireUSA
| | - Bryan W. Clark
- Office of Research and Development, Center for Environmental Measurement and Modeling, Atlantic Coastal Environmental Sciences DivisionUS Environmental Protection AgencyNarragansettRhode IslandUSA
| | - Noah M. Reid
- Department of Molecular & Cell BiologyUniversity of ConnecticutStorrsConnecticutUSA
| | - Sibel I. Karchner
- Biology DepartmentWoods Hole Oceanographic InstitutionWoods HoleMassachusettsUSA
| | - Jennifer L. Roach
- Department of Environmental Toxicology, Center for Population Biology, Coastal and Marine Sciences InstituteUniversity of California, DavisDavisCaliforniaUSA
| | - Mark E. Hahn
- Biology DepartmentWoods Hole Oceanographic InstitutionWoods HoleMassachusettsUSA
| | - Diane Nacci
- Office of Research and Development, Center for Environmental Measurement and Modeling, Atlantic Coastal Environmental Sciences DivisionUS Environmental Protection AgencyNarragansettRhode IslandUSA
| | - Andrew Whitehead
- Department of Environmental Toxicology, Center for Population Biology, Coastal and Marine Sciences InstituteUniversity of California, DavisDavisCaliforniaUSA
| |
Collapse
|
2
|
Miller JT, Clark BW, Reid NM, Karchner SI, Roach JL, Hahn ME, Nacci D, Whitehead A. Independently evolved pollution resistance in four killifish populations is largely explained by few variants of large effect. bioRxiv 2023:2023.04.07.536079. [PMID: 37066319 PMCID: PMC10104127 DOI: 10.1101/2023.04.07.536079] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/18/2023]
Abstract
The genetic architecture of phenotypic traits can affect the mode and tempo of trait evolution. Human-altered environments can impose strong natural selection, where successful evolutionary adaptation requires swift and large phenotypic shifts. In these scenarios, theory predicts the influence of few adaptive variants of large effect, but empirical studies that have revealed the genetic architecture of rapidly evolved phenotypes are rare, especially for populations inhabiting polluted environments. Fundulus killifish have repeatedly evolved adaptive resistance to extreme pollution in urban estuaries. Prior studies, including genome scans for signatures of natural selection, have revealed some of the genes and pathways important for evolved pollution resistance, and provide context for the genotype-phenotype association studies reported here. We created multiple quantitative trait locus (QTL) mapping families using progenitors from four different resistant populations, and genetically mapped variation in sensitivity (developmental perturbations) following embryonic exposure to a model toxicant PCB-126. We found that a few large-effect QTL loci accounted for resistance to PCB-mediated developmental toxicity. QTLs harbored candidate genes that govern the regulation of aryl hydrocarbon receptor (AHR) signaling, where some (but not all) of these QTL loci were shared across all populations, and some (but not all) of these loci showed signatures of recent natural selection in the corresponding wild population. Some strong candidate genes for PCB resistance inferred from genome scans in wild populations were identified as QTL, but some key candidate genes were not. We conclude that rapidly evolved resistance to the developmental defects normally caused by PCB-126 is governed by few genes of large effect. However, other aspects of resistance beyond developmental phenotypes may be governed by additional loci, such that comprehensive resistance to PCB-126, and to the mixtures of chemicals that distinguish urban estuaries more broadly, may be more genetically complex.
Collapse
Affiliation(s)
- Jeffrey T Miller
- Department of Environmental Toxicology, Center for Population Biology, Coastal and Marine Sciences Institute, University of California, Davis, CA
| | - Bryan W Clark
- US Environmental Protection Agency, Office of Research and Development, Center for Environmental Measurement and Modeling, Atlantic Coastal Environmental Sciences Division, Narragansett, RI
| | - Noah M Reid
- Department of Molecular & Cell Biology, University of Connecticut, Storrs, CT
| | - Sibel I Karchner
- Biology Department, Woods Hole Oceanographic Institution, Woods Hole, MA
| | - Jennifer L Roach
- Department of Environmental Toxicology, Center for Population Biology, Coastal and Marine Sciences Institute, University of California, Davis, CA
| | - Mark E Hahn
- Biology Department, Woods Hole Oceanographic Institution, Woods Hole, MA
| | - Diane Nacci
- US Environmental Protection Agency, Office of Research and Development, Center for Environmental Measurement and Modeling, Atlantic Coastal Environmental Sciences Division, Narragansett, RI
| | - Andrew Whitehead
- Department of Environmental Toxicology, Center for Population Biology, Coastal and Marine Sciences Institute, University of California, Davis, CA
| |
Collapse
|
3
|
LaFountain AM, McMahon HE, Reid NM, Yuan YW. To stripe or not to stripe: the origin of a novel foliar pigmentation pattern in monkeyflowers (Mimulus). New Phytol 2023; 237:310-322. [PMID: 36101514 PMCID: PMC10601762 DOI: 10.1111/nph.18486] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Accepted: 09/07/2022] [Indexed: 06/15/2023]
Abstract
The origin of phenotypic novelty is one of the most challenging problems in evolutionary biology. Although genetic regulatory network rewiring or co-option has been widely recognised as a major contributor, in most cases how such genetic rewiring/co-option happens is completely unknown. We have studied a novel foliar pigmentation pattern that evolved recently in the monkeyflower species Mimulus verbenaceus. Through genome-wide association tests using wild-collected samples, experimental crosses of laboratory inbred lines, gene expression analyses, and functional assays, we identified an anthocyanin-activating R2R3-MYB gene, STRIPY, as the causal gene triggering the emergence of the discrete, mediolateral anthocyanin stripe in the M. verbenaceus leaf. Chemical mutagenesis revealed the existence of upstream activators and repressors that form a 'hidden' prepattern along the leaf proximodistal axis, potentiating the unique expression pattern of STRIPY. Population genomics analyses did not reveal signatures of positive selection, indicating that nonadaptive processes may be responsible for the establishment of this novel trait in the wild. This study demonstrates that the origin of phenotypic novelty requires at least two separate phases, potentiation and actualisation. The foliar stripe pattern of M. verbenaceus provides an excellent platform to probe the molecular details of both processes in future studies.
Collapse
Affiliation(s)
- Amy M. LaFountain
- Department of Ecology and Evolutionary Biology, University of Connecticut, 75 North Eagleville Road, Storrs, CT, USA. 06269-3043
| | - Hayley E. McMahon
- Department of Ecology and Evolutionary Biology, University of Connecticut, 75 North Eagleville Road, Storrs, CT, USA. 06269-3043
| | - Noah M. Reid
- Institute for Systems Genomics, University of Connecticut, 67 North Eagleville Road, Storrs, CT, USA 06269-3197
| | - Yao-Wu Yuan
- Department of Ecology and Evolutionary Biology, University of Connecticut, 75 North Eagleville Road, Storrs, CT, USA. 06269-3043
- Institute for Systems Genomics, University of Connecticut, 67 North Eagleville Road, Storrs, CT, USA 06269-3197
| |
Collapse
|
4
|
Oziolor EM, Reid NM, Yair S, Lee KM, Guberman VerPloeg S, Bruns PC, Shaw JR, Whitehead A, Matson CW. Adaptive introgression enables evolutionary rescue from extreme environmental pollution. Science 2019; 364:455-457. [DOI: 10.1126/science.aav4155] [Citation(s) in RCA: 114] [Impact Index Per Article: 22.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2018] [Accepted: 03/25/2019] [Indexed: 12/19/2022]
Abstract
Radical environmental change that provokes population decline can impose constraints on the sources of genetic variation that may enable evolutionary rescue. Adaptive toxicant resistance has rapidly evolved in Gulf killifish (Fundulus grandis) that occupy polluted habitats. We show that resistance scales with pollution level and negatively correlates with inducibility of aryl hydrocarbon receptor (AHR) signaling. Loci with the strongest signatures of recent selection harbor genes regulating AHR signaling. Two of these loci introgressed recently (18 to 34 generations ago) from Atlantic killifish (F. heteroclitus). One introgressed locus contains a deletion in AHR that confers a large adaptive advantage [selection coefficient (s) = 0.8]. Given the limited migration of killifish, recent adaptive introgression was likely mediated by human-assisted transport. We suggest that interspecies connectivity may be an important source of adaptive variation during extreme environmental change.
Collapse
|
5
|
Rodgers R, Roach JL, Reid NM, Whitehead A, Duvernell DD. Phylogenomic analysis of Fundulidae (Teleostei: Cyprinodotiformes) using RNA-sequencing data. Mol Phylogenet Evol 2017; 121:150-157. [PMID: 29289545 DOI: 10.1016/j.ympev.2017.12.030] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2017] [Revised: 11/27/2017] [Accepted: 12/27/2017] [Indexed: 11/18/2022]
Abstract
Fishes of the New World cyprinodontiform family Fundulidae display a wide variety of tolerance to environmental conditions, making them a valuable model system for comparative, evolutionary, and environmental studies. Despite numerous attempts to resolve the phylogenetic relationships of family Fundulidae, the basal structure of the phylogeny remains unresolved. The lack of a robust and fully resolved phylogeny for family Fundulidae and its most speciose genus Fundulus is an impediment to future research. This study utilized novel RNA-sequencing data for phylogenetic inference among16 members of Fundulidae to better refine the basal nodes of the family and confront long-standing questions regarding (1) the monophyletic status of genus Fundulus, and validity of the Lucania and recently synonymized Adinia genera; (2) the relationship of the west coast endemic Fundulus parvipinnis and F. lima to other Fundulus species; and (3) the validity of subgeneric classifications. In addition, previously published nuclear gene sequences for 32 Fundulidae species were re-analyzed in combination with novel RNA-sequencing data. Maximum likelihood and Bayesian analyses generated identical phylogenies with strong statistical support at nearly all nodes, demonstrating the utility of RNA-sequencing data in constructing robust phylogenies not achievable by previous methods. While many past hypothesized evolutionary relationships for Fundulidae were reinforced, several alternative relationships are hypothesized at basal nodes resulting in a re-analysis of the deeper structure of family Fundulidae. These results reveal family Fundulidae as a paraphyletic grouping of members of genus Fundulus and Lucania and supports the previous synonymy of genus Adinia with genus Fundulus.
Collapse
Affiliation(s)
- Rachel Rodgers
- Biological Sciences, Southern Illinois University Edwardsville, Edwardsville, IL 62026, USA
| | - Jennifer L Roach
- Department of Environmental Toxicology, University of California, Davis, CA 95616, USA
| | - Noah M Reid
- Department of Environmental Toxicology, University of California, Davis, CA 95616, USA; Department of Molecular and Cell Biology, University of Connecticut, Storrs, CT 06269, USA
| | - Andrew Whitehead
- Department of Environmental Toxicology, University of California, Davis, CA 95616, USA
| | - David D Duvernell
- Biological Sciences, Southern Illinois University Edwardsville, Edwardsville, IL 62026, USA; Biological Sciences, Missouri University of Science and Technology, Rolla, MO 65401, USA.
| |
Collapse
|
6
|
Holland EB, Goldstone JV, Pessah IN, Whitehead A, Reid NM, Karchner SI, Hahn ME, Nacci DE, Clark BW, Stegeman JJ. Ryanodine receptor and FK506 binding protein 1 in the Atlantic killifish (Fundulus heteroclitus): A phylogenetic and population-based comparison. Aquat Toxicol 2017; 192:105-115. [PMID: 28942070 PMCID: PMC5662517 DOI: 10.1016/j.aquatox.2017.09.002] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/12/2017] [Revised: 09/01/2017] [Accepted: 09/02/2017] [Indexed: 05/12/2023]
Abstract
Non-dioxin-like polychlorinated biphenyls (NDL PCBs) activate ryanodine receptors (RyR), microsomal Ca2+ channels of broad significance. Teleost fish may be important models for NDL PCB neurotoxicity, and we used sequencing databases to characterize teleost RyR and FK506 binding protein 12 or 12.6kDa (genes FKBP1A; FKBP1B), which promote NDL PCB-triggered Ca2+ dysregulation. Particular focus was placed on describing genes in the Atlantic killifish (Fundulus heteroclitus) genome and searching available RNA-sequencing datasets for single nucleotide variants (SNV) between PCB tolerant killifish from New Bedford Harbor (NBH) versus sensitive killifish from Scorton Creek (SC), MA. Consistent with the teleost whole genome duplication (tWGD), killifish have six RyR genes, corresponding to a and b paralogs of mammalian RyR1, 2 and 3. The presence of six RyR genes was consistent in all teleosts investigated including zebrafish. Killifish have four FKBP1; one FKBP1b and three FKBP1a named FKBP1aa, FKBP1ab, likely from the tWGD and a single gene duplicate FKBP1a3 suggested to have arisen in Atherinomorphae. The RyR and FKBP1 genes displayed tissue and developmental stage-specific mRNA expression, and the previously uncharacterized RyR3, herein named RyR3b, and all FKBP1 genes were prominent in brain. We identified a SNV in RyR3b encoding missense mutation E1458D. In NBH killifish, 57% were heterozygous and 28% were homozygous for this SNV, whereas almost all SC killifish (94%) lacked the variant (n≥39 per population). The outlined sequence differences between mammalian and teleost RyR and FKBP1 together with outlined population differences in SNV frequency may contribute to our understanding of NDL PCB neurotoxicity.
Collapse
Affiliation(s)
- Erika B Holland
- Department of Biological Sciences, California State University of Long Beach, Long Beach, CA, USA; Department of Biology, Woods Hole Oceanographic Institution, Woods Hole MA, USA; Department of Molecular Biosciences, School of Veterinary Medicine, University of California Davis, Davis, CA, USA.
| | - Jared V Goldstone
- Department of Biology, Woods Hole Oceanographic Institution, Woods Hole MA, USA
| | - Isaac N Pessah
- Department of Molecular Biosciences, School of Veterinary Medicine, University of California Davis, Davis, CA, USA
| | - Andrew Whitehead
- Department of Environmental Toxicology, College of Agricultural and Environmental Sciences,University of California Davis, Davis, CA, USA
| | - Noah M Reid
- Department of Environmental Toxicology, College of Agricultural and Environmental Sciences,University of California Davis, Davis, CA, USA
| | - Sibel I Karchner
- Department of Biology, Woods Hole Oceanographic Institution, Woods Hole MA, USA
| | - Mark E Hahn
- Department of Biology, Woods Hole Oceanographic Institution, Woods Hole MA, USA
| | - Diane E Nacci
- Atlantic Ecology Division, Office of Research and Development, US Environmental Protection Agency, Narragansett, RI, USA
| | - Bryan W Clark
- Oak Ridge Institute for Science and Education at the United States Environmental Protection Agency, Office of Research and Development, Narragansett, RI, 02882, USA
| | - John J Stegeman
- Department of Biology, Woods Hole Oceanographic Institution, Woods Hole MA, USA
| |
Collapse
|
7
|
Whitehead A, Clark BW, Reid NM, Hahn ME, Nacci D. When evolution is the solution to pollution: Key principles, and lessons from rapid repeated adaptation of killifish ( Fundulus heteroclitus) populations. Evol Appl 2017; 10:762-783. [PMID: 29151869 PMCID: PMC5680427 DOI: 10.1111/eva.12470] [Citation(s) in RCA: 77] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2016] [Accepted: 02/10/2017] [Indexed: 12/18/2022] Open
Abstract
For most species, evolutionary adaptation is not expected to be sufficiently rapid to buffer the effects of human‐mediated environmental changes, including environmental pollution. Here we review how key features of populations, the characteristics of environmental pollution, and the genetic architecture underlying adaptive traits, may interact to shape the likelihood of evolutionary rescue from pollution. Large populations of Atlantic killifish (Fundulus heteroclitus) persist in some of the most contaminated estuaries of the United States, and killifish studies have provided some of the first insights into the types of genomic changes that enable rapid evolutionary rescue from complexly degraded environments. We describe how selection by industrial pollutants and other stressors has acted on multiple populations of killifish and posit that extreme nucleotide diversity uniquely positions this species for successful evolutionary adaptation. Mechanistic studies have identified some of the genetic underpinnings of adaptation to a well‐studied class of toxic pollutants; however, multiple genetic regions under selection in wild populations seem to reflect more complex responses to diverse native stressors and/or compensatory responses to primary adaptation. The discovery of these pollution‐adapted killifish populations suggests that the evolutionary influence of anthropogenic stressors as selective agents occurs widely. Yet adaptation to chemical pollution in terrestrial and aquatic vertebrate wildlife may rarely be a successful “solution to pollution” because potentially adaptive phenotypes may be complex and incur fitness costs, and therefore be unlikely to evolve quickly enough, especially in species with small population sizes.
Collapse
Affiliation(s)
- Andrew Whitehead
- Department of Environmental Toxicology University of California Davis Davis CA USA
| | - Bryan W Clark
- Atlantic Ecology Division National Health and Environmental Effects Research Laboratory Office of Research and Development Oak Ridge Institute for Science and Education US Environmental Protection Agency Narragansett RI USA
| | - Noah M Reid
- Department of Molecular and Cell Biology University of Connecticut Storrs CT USA
| | - Mark E Hahn
- Department of Biology Woods Hole Oceanographic Institution Woods Hole MA USA.,Superfund Research Program Boston University Boston MA USA
| | - Diane Nacci
- Atlantic Ecology Division National Health and Environmental Effects Research Laboratory Office of Research and Development US Environmental Protection Agency Narragansett RI USA
| |
Collapse
|
8
|
Reid NM, Jackson CE, Gilbert D, Minx P, Montague MJ, Hampton TH, Helfrich LW, King BL, Nacci DE, Aluru N, Karchner SI, Colbourne JK, Hahn ME, Shaw JR, Oleksiak MF, Crawford DL, Warren WC, Whitehead A. The landscape of extreme genomic variation in the highly adaptable Atlantic killifish. Genome Biol Evol 2017; 9:659-676. [PMID: 28201664 PMCID: PMC5381573 DOI: 10.1093/gbe/evx023] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2016] [Revised: 01/30/2017] [Accepted: 02/04/2017] [Indexed: 12/22/2022] Open
Abstract
Understanding and predicting the fate of populations in changing environments require knowledge about the mechanisms that support phenotypic plasticity and the adaptive value and evolutionary fate of genetic variation within populations. Atlantic killifish (Fundulus heteroclitus) exhibit extensive phenotypic plasticity that supports large population sizes in highly fluctuating estuarine environments. Populations have also evolved diverse local adaptations. To yield insights into the genomic variation that supports their adaptability, we sequenced a reference genome and 48 additional whole genomes from a wild population. Evolution of genes associated with cell cycle regulation and apoptosis is accelerated along the killifish lineage, which is likely tied to adaptations for life in highly variable estuarine environments. Genome-wide standing genetic variation, including nucleotide diversity and copy number variation, is extremely high. The highest diversity genes are those associated with immune function and olfaction, whereas genes under greatest evolutionary constraint are those associated with neurological, developmental, and cytoskeletal functions. Reduced genetic variation is detected for tight junction proteins, which in killifish regulate paracellular permeability that supports their extreme physiological flexibility. Low-diversity genes engage in more regulatory interactions than high-diversity genes, consistent with the influence of pleiotropic constraint on molecular evolution. High genetic variation is crucial for continued persistence of species given the pace of contemporary environmental change. Killifish populations harbor among the highest levels of nucleotide diversity yet reported for a vertebrate species, and thus may serve as a useful model system for studying evolutionary potential in variable and changing environments.
Collapse
Affiliation(s)
- Noah M Reid
- Department of Environmental Toxicology, University of California, Davis, CA 95616
| | - Craig E Jackson
- School of Public and Environmental Affairs, Indiana University, Bloomington, IN 47405
| | - Don Gilbert
- Biology Department, Indiana University, Bloomington, IN 47405
| | - Patrick Minx
- McDonnell Genome Institute, Washington University School of Medicine, St Louis, MO 63108
| | - Michael J Montague
- McDonnell Genome Institute, Washington University School of Medicine, St Louis, MO 63108
| | - Thomas H Hampton
- Department of Microbiology and Immunology, Dartmouth College Geisel School of Medicine, Hanover, NH 03755
| | - Lily W Helfrich
- Biology Department, Woods Hole Oceanographic Institution, Woods Hole, MA 02543
| | - Benjamin L King
- Mount Desert Island Biological Laboratory, Salisbury Cove, ME 04672
| | - Diane E Nacci
- US Environmental Protection Agency, Office of Research and Development, Narragansett, RI, 02882
| | - Neel Aluru
- Biology Department, Woods Hole Oceanographic Institution, Woods Hole, MA 02543
| | - Sibel I Karchner
- Biology Department, Woods Hole Oceanographic Institution, Woods Hole, MA 02543
| | - John K Colbourne
- School of Biosciences, University of Birmingham, United Kingdom, B15 2TT
| | - Mark E Hahn
- Biology Department, Woods Hole Oceanographic Institution, Woods Hole, MA 02543
| | - Joseph R Shaw
- School of Public and Environmental Affairs, Indiana University, Bloomington, IN 47405
| | - Marjorie F Oleksiak
- Department of Marine Biology and Ecology, Rosenstiel School of Marine and Atmospheric Science, University of Miami, Miami, FL 33149
| | - Douglas L Crawford
- Department of Marine Biology and Ecology, Rosenstiel School of Marine and Atmospheric Science, University of Miami, Miami, FL 33149
| | - Wesley C Warren
- McDonnell Genome Institute, Washington University School of Medicine, St Louis, MO 63108
| | - Andrew Whitehead
- Department of Environmental Toxicology, University of California, Davis, CA 95616
| |
Collapse
|
9
|
Reid NM, Proestou DA, Clark BW, Warren WC, Colbourne JK, Shaw JR, Karchner SI, Hahn ME, Nacci D, Oleksiak MF, Crawford DL, Whitehead A. The genomic landscape of rapid repeated evolutionary adaptation to toxic pollution in wild fish. Science 2016; 354:1305-1308. [PMID: 27940876 PMCID: PMC5206662 DOI: 10.1126/science.aah4993] [Citation(s) in RCA: 250] [Impact Index Per Article: 31.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2016] [Accepted: 10/31/2016] [Indexed: 01/20/2023]
Abstract
Atlantic killifish populations have rapidly adapted to normally lethal levels of pollution in four urban estuaries. Through analysis of 384 whole killifish genome sequences and comparative transcriptomics in four pairs of sensitive and tolerant populations, we identify the aryl hydrocarbon receptor-based signaling pathway as a shared target of selection. This suggests evolutionary constraint on adaptive solutions to complex toxicant mixtures at each site. However, distinct molecular variants apparently contribute to adaptive pathway modification among tolerant populations. Selection also targets other toxicity-mediating genes and genes of connected signaling pathways; this indicates complex tolerance phenotypes and potentially compensatory adaptations. Molecular changes are consistent with selection on standing genetic variation. In killifish, high nucleotide diversity has likely been a crucial substrate for selective sweeps to propel rapid adaptation.
Collapse
Affiliation(s)
- Noah M Reid
- Department of Environmental Toxicology, University of California, Davis, CA 95616, USA
| | - Dina A Proestou
- Agricultural Research Service, U.S. Department of Agriculture, Kingston, RI 02881, USA
| | - Bryan W Clark
- Oak Ridge Institute for Science and Education, Office of Research and Development, U.S. Environmental Protection Agency, Narragansett, RI 02882, USA
| | - Wesley C Warren
- McDonnell Genome Institute, Washington University School of Medicine, St. Louis, MO 63108, USA
| | - John K Colbourne
- School of Biosciences, University of Birmingham, Edgbaston B15 2TT, UK
| | - Joseph R Shaw
- School of Biosciences, University of Birmingham, Edgbaston B15 2TT, UK
- School of Public and Environmental Affairs, Indiana University, Bloomington, IN 47405, USA
| | - Sibel I Karchner
- Biology Department, Woods Hole Oceanographic Institution, Woods Hole, MA 02543, USA
- Boston University Superfund Research Program, Boston University, Boston, MA 02118, USA
| | - Mark E Hahn
- Biology Department, Woods Hole Oceanographic Institution, Woods Hole, MA 02543, USA
- Boston University Superfund Research Program, Boston University, Boston, MA 02118, USA
| | - Diane Nacci
- Office of Research and Development, U.S. Environmental Protection Agency, Narragansett, RI 02882, USA
| | - Marjorie F Oleksiak
- Department of Marine Biology and Ecology, Rosenstiel School of Marine and Atmospheric Science, University of Miami, Miami, FL 33149, USA
| | - Douglas L Crawford
- Department of Marine Biology and Ecology, Rosenstiel School of Marine and Atmospheric Science, University of Miami, Miami, FL 33149, USA
| | - Andrew Whitehead
- Department of Environmental Toxicology, University of California, Davis, CA 95616, USA.
| |
Collapse
|
10
|
Carstens BC, Gruenstaeudl M, Reid NM. Community trees: Identifying codiversification in the Páramo dipteran community. Evolution 2016; 70:1080-93. [PMID: 27061575 DOI: 10.1111/evo.12916] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2015] [Revised: 03/10/2016] [Accepted: 03/15/2016] [Indexed: 01/23/2023]
Abstract
Groups of codistributed species that responded in a concerted manner to environmental events are expected to share patterns of evolutionary diversification. However, the identification of such groups has largely been based on qualitative, post hoc analyses. We develop here two methods (posterior predictive simulation [PPS], Kuhner-Felsenstein [K-F] analysis of variance [ANOVA]) for the analysis of codistributed species that, given a group of species with a shared pattern of diversification, allow empiricists to identify those taxa that do not codiversify (i.e., "outlier" species). The identification of outlier species makes it possible to jointly estimate the evolutionary history of co-diversifying taxa. To evaluate the approaches presented here, we collected data from Páramo dipterans, identified outlier species, and estimated a "community tree" from species that are identified as having codiversified. Our results demonstrate that dipteran communities from different Páramo habitats in the same mountain range are more closely related than communities in other ranges. We also conduct simulation testing to evaluate this approach. Results suggest that our approach provides a useful addition to comparative phylogeographic methods, while identifying aspects of the analysis that require careful interpretation. In particular, both the PPS and K-F ANOVA perform acceptably when there are one or two outlier species, but less so as the number of outliers increases. This is likely a function of the corresponding degradation of the signal of community divergence; without a strong signal from a codiversifying community, there is no dominant pattern from which to detect an outlier species. For this reason, both the magnitude of K-F distance distribution and outside knowledge about the phylogeographic history of each putative member of the community should be considered when interpreting the results.
Collapse
Affiliation(s)
- Bryan C Carstens
- Department of Evolution, Ecology and Organismal Biology, The Ohio State University, 318 W 12th Avenue, Columbus, Ohio, 43210.
| | - Michael Gruenstaeudl
- Institut für Biologie-Botanik, Dahlem Centre of Plant Sciences, Freie Universität Berlin, Altensteinstraße 6, Berlin, 14195, Germany
| | - Noah M Reid
- Department of Environmental Toxicology, University of California-Davis, Davis, California, 95616
| |
Collapse
|
11
|
Reid NM, Whitehead A. Functional genomics to assess biological responses to marine pollution at physiological and evolutionary timescales: toward a vision of predictive ecotoxicology. Brief Funct Genomics 2015; 15:358-64. [PMID: 26700295 DOI: 10.1093/bfgp/elv060] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
Marine pollution is ubiquitous, and is one of the key factors influencing contemporary marine biodiversity worldwide. To protect marine biodiversity, how do we surveil, document and predict the short- and long-term impacts of pollutants on at-risk species? Modern genomics tools offer high-throughput, information-rich and increasingly cost-effective approaches for characterizing biological responses to environmental stress, and are important tools within an increasing sophisticated kit for surveiling and assessing impacts of pollutants on marine species. Through the lens of recent research in marine killifish, we illustrate how genomics tools may be useful for screening chemicals and pollutants for biological activity and to reveal specific mechanisms of action. The high dimensionality of transcriptomic responses enables their usage as highly specific fingerprints of exposure, and these fingerprints can be used to diagnose environmental problems. We also emphasize that molecular pathways recruited to respond at physiological timescales are the same pathways that may be targets for natural selection during chronic exposure to pollutants. Gene complement and sequence variation in those pathways can be related to variation in sensitivity to environmental pollutants within and among species. Furthermore, allelic variation associated with evolved tolerance in those pathways could be tracked to estimate the pace of environmental health decline and recovery. We finish by integrating these paradigms into a vision of how genomics approaches could anchor a modernized framework for advancing the predictive capacity of environmental and ecotoxicological science.
Collapse
|
12
|
Gruenstaeudl M, Reid NM, Wheeler GL, Carstens BC. Posterior predictive checks of coalescent models: P2C2M, an R package. Mol Ecol Resour 2015; 16:193-205. [DOI: 10.1111/1755-0998.12435] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2014] [Revised: 05/22/2015] [Accepted: 05/26/2015] [Indexed: 02/04/2023]
Affiliation(s)
- Michael Gruenstaeudl
- Department of Evolution, Ecology & Organismal Biology; Ohio State University; Columbus OH 43210 USA
| | - Noah M. Reid
- Department of Environmental Toxicology; University of California; Davis CA 95616 USA
| | - Gregory L. Wheeler
- Department of Evolution, Ecology & Organismal Biology; Ohio State University; Columbus OH 43210 USA
| | - Bryan C. Carstens
- Department of Evolution, Ecology & Organismal Biology; Ohio State University; Columbus OH 43210 USA
| |
Collapse
|
13
|
Carstens BC, Pelletier TA, Reid NM, Satler JD. How to fail at species delimitation. Mol Ecol 2013; 22:4369-83. [PMID: 23855767 DOI: 10.1111/mec.12413] [Citation(s) in RCA: 632] [Impact Index Per Article: 57.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2013] [Revised: 05/31/2013] [Accepted: 06/04/2013] [Indexed: 11/29/2022]
Abstract
Species delimitation is the act of identifying species-level biological diversity. In recent years, the field has witnessed a dramatic increase in the number of methods available for delimiting species. However, most recent investigations only utilize a handful (i.e. 2-3) of the available methods, often for unstated reasons. Because the parameter space that is potentially relevant to species delimitation far exceeds the parameterization of any existing method, a given method necessarily makes a number of simplifying assumptions, any one of which could be violated in a particular system. We suggest that researchers should apply a wide range of species delimitation analyses to their data and place their trust in delimitations that are congruent across methods. Incongruence across the results from different methods is evidence of either a difference in the power to detect cryptic lineages across one or more of the approaches used to delimit species and could indicate that assumptions of one or more of the methods have been violated. In either case, the inferences drawn from species delimitation studies should be conservative, for in most contexts it is better to fail to delimit species than it is to falsely delimit entities that do not represent actual evolutionary lineages.
Collapse
Affiliation(s)
- Bryan C Carstens
- Department of Evolution, Ecology and Organismal Biology, The Ohio State University, 318 W. 12th Avenue, Columbus, OH 43210-1293, USA.
| | | | | | | |
Collapse
|
14
|
Reid NM, Carstens BC. Phylogenetic estimation error can decrease the accuracy of species delimitation: a Bayesian implementation of the general mixed Yule-coalescent model. BMC Evol Biol 2012; 12:196. [PMID: 23031350 PMCID: PMC3503838 DOI: 10.1186/1471-2148-12-196] [Citation(s) in RCA: 324] [Impact Index Per Article: 27.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2012] [Accepted: 09/19/2012] [Indexed: 11/10/2022] Open
Abstract
Background Species are considered the fundamental unit in many ecological and evolutionary analyses, yet accurate, complete, accessible taxonomic frameworks with which to identify them are often unavailable to researchers. In such cases DNA sequence-based species delimitation has been proposed as a means of estimating species boundaries for further analysis. Several methods have been proposed to accomplish this. Here we present a Bayesian implementation of an evolutionary model-based method, the general mixed Yule-coalescent model (GMYC). Our implementation integrates over the parameters of the model and uncertainty in phylogenetic relationships using the output of widely available phylogenetic models and Markov-Chain Monte Carlo (MCMC) simulation in order to produce marginal probabilities of species identities. Results We conducted simulations testing the effects of species evolutionary history, levels of intraspecific sampling and number of nucleotides sequenced. We also re-analyze the dataset used to introduce the original GMYC model. We found that the model results are improved with addition of DNA sequence and increased sampling, although these improvements have limits. The most important factor in the success of the model is the underlying phylogenetic history of the species under consideration. Recent and rapid divergences result in higher amounts of uncertainty in the model and eventually cause the model to fail to accurately assess uncertainty in species limits. Conclusion Our results suggest that the GMYC model can be useful under a wide variety of circumstances, particularly in cases where divergences are deeper, or taxon sampling is incomplete, as in many studies of ecological communities, but that, in accordance with expectations from coalescent theory, rapid, recent radiations may yield inaccurate results. Our implementation differs from existing ones in two ways: it allows for the accounting for important sources of uncertainty in the model (phylogenetic and in parameters specific to the model) and in the specification of informative prior distributions that can increase the precision of the model. We have incorporated this model into a user-friendly R package available on the authors’ websites.
Collapse
Affiliation(s)
- Noah M Reid
- Department of Biological Science, Louisiana State University, Baton Rouge, LA 70803, USA.
| | | |
Collapse
|
15
|
Abstract
Data analysis in phylogeographic investigations is typically conducted in either a qualitative manner, or alternatively via the testing of null hypotheses. The former, where inferences about population processes are derived from geographical patterns of genetic variation, may be subject to confirmation bias and prone to overinterpretation. Testing the predictions of null hypotheses is arguably less prone to bias than qualitative approaches, but only if the tested hypotheses are biologically meaningful. As it is difficult to know a priori if this is the case, there is the general need for additional methodological approaches in phylogeographic research. Here, we explore an alternative method for analysing phylogeographic data that utilizes information theory to quantify the probability of multiple hypotheses given the data. We accomplish this by augmenting the model-selection procedure implemented in ima with calculations of Akaike Information Criterion scores and model probabilities. We generate a ranking of 17 models each representing a set of historical evolutionary processes that may have contributed to the evolution of Plethodon idahoensis, and then quantify the relative strength of support for each hypothesis given the data using metrics borrowed from information theory. Our results suggest that two models have high probability given the data. Each of these models includes population divergence and estimates of ancestral theta that differ from estimates of descendent theta, inferences consistent with prior work in this system. However, the models disagree in that one includes migration as a parameter and one does not, suggesting that there are two regions of parameter space that produce model likelihoods that are similar in magnitude given our data. Results of a simulation study suggest that when data are simulated with migration, most of the optimal models include migration as a parameter, and further that when all of the shared polymorphism results from incomplete lineage sorting, most of the optimal models do not. The results could also indicate a lack of precision, which may be a product of the amount of data that we have collected. In any case, the information-theoretic metrics that we have applied to the analysis of our data are statistically rigorous, as are hypothesis-testing approaches, but move beyond the 'reject/fail to reject' dichotomy of conventional hypothesis testing in a manner that provides considerably more flexibility to researchers.
Collapse
Affiliation(s)
- Bryan C Carstens
- Department of Biological Sciences, Louisiana State University, Baton Rouge, LA 70803, USA.
| | | | | |
Collapse
|
16
|
Campbell P, Reid NM, Zubaid A, Adnan AM, Kunz TH. Comparative Roosting Ecology of Cynopterus (Chiroptera: Pteropodidae) Fruit Bats in Peninsular Malaysia. Biotropica 2006. [DOI: 10.1111/j.1744-7429.2006.00203.x] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
17
|
Roberts N, Cruz-Orive LM, Reid NM, Brodie DA, Bourne M, Edwards RH. Unbiased estimation of human body composition by the Cavalieri method using magnetic resonance imaging. J Microsc 1993; 171:239-53. [PMID: 8246272 DOI: 10.1111/j.1365-2818.1993.tb03381.x] [Citation(s) in RCA: 119] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
The classical methods for estimating the volume of human body compartments in vivo (e.g. skin-fold thickness for fat, radioisotope counting for different compartments, etc.) are generally indirect and rely on essentially empirical relationships--hence they are biased to unknown degrees. The advent of modern non-invasive scanning techniques, such as X-ray computed tomography (CT) and magnetic resonance imaging (MRI) is now widening the scope of volume quantification, especially in combination with stereological methods. Apart from its superior soft tissue contrast, MRI enjoys the distinct advantage of not using ionizing radiations. By a proper landmarking and control of the scanner couch, an adult male volunteer was scanned exhaustively into parallel systematic MR 'sections'. Four compartments were defined, namely bone, muscle, organs and fat (which included the skin), and their corresponding volumes were easily and efficiently estimated by the Cavalieri method: the total section area of a compartment times the section interval estimates the volume of the compartment without bias. Formulae and nomograms are given to predict the errors and to optimize the design. To estimate an individual's muscle volume with a 5% coefficient of error, 10 sections and less than 10 min point counting (to estimate the relevant section areas) are required. Bone and fat require about twice as much work. To estimate the mean muscle volume of a population with the same error contribution, from a random sample of six subjects, the workload per subject can be divided by square root of 6, namely 4 min per subject. For a given number of sections planimetry would be as accurate but far more time consuming than point counting.
Collapse
Affiliation(s)
- N Roberts
- Magnetic Resonance Research Centre, University of Liverpool
| | | | | | | | | | | |
Collapse
|
18
|
Reid NM, Dickson DP, Greenwood C, Thompson A, Kadir FH, Moore GR. Evidence for Mössbauer spectroscopy for different forms of iron core in Pseudomonas aeruginosa bacterial ferritin. Biochem J 1990; 272:263-4. [PMID: 2124804 PMCID: PMC1149688 DOI: 10.1042/bj2720263] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Mössbauer spectroscopic studies of whole cells of Pseudomonas aeruginosa, grown under different conditions, indicate that the predominant form of iron in the cells varies significantly. These differences are interpreted in terms of differences in the nature of the iron cores of the bacterial ferritin, which result from different growth conditions.
Collapse
Affiliation(s)
- N M Reid
- Department of Physics, University of Liverpool, U.K
| | | | | | | | | | | |
Collapse
|
19
|
Ward RJ, O'Connell MJ, Dickson DP, Reid NM, Wade VJ, Mann S, Bomford A, Peters TJ. Biochemical studies of the iron cores and polypeptide shells of haemosiderin isolated from patients with primary or secondary haemochromatosis. Biochim Biophys Acta 1989; 993:131-3. [PMID: 2804120 DOI: 10.1016/0304-4165(89)90153-0] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Haemosiderin isolated from different iron-loading syndromes, primary haemochromatosis (PHC) and secondary haemochromatosis (SHC) biochemically exhibited differences in both their iron core and peptide composition. The rate of release of iron from PHC haemosiderin to oxalate was 3-fold greater than that from SHC haemosiderin. The major peptides separated by SDS-PAGE showed a major band at Mr 20,000 for PHC haemosiderin and at Mr 15,000 for SHC haemosiderin.
Collapse
Affiliation(s)
- R J Ward
- Department of Clinical Biochemistry, King's College School of Medicine and Dentistry, London, U.K
| | | | | | | | | | | | | | | |
Collapse
|
20
|
Dickson DP, Reid NM, Mann S, Wade VJ, Ward RJ, Peters TJ. Mössbauer spectroscopy, electron microscopy and electron diffraction studies of the iron cores in various human and animal haemosiderins. Biochim Biophys Acta 1988; 957:81-90. [PMID: 3179322 DOI: 10.1016/0167-4838(88)90159-8] [Citation(s) in RCA: 52] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Mössbauer spectroscopy has indicated significant differences in the iron-containing cores of various haemosiderins. In the present study, haemosiderin was isolated from a number of animal species including man. In addition, haemosiderin was isolated from patients with primary idiopathic haemochromatosis or with secondary (transfusional) iron-overload. The iron cores of the animal and normal human haemosiderin appear to be very similar by Mössbauer spectroscopy, and the electron diffraction data indicate a ferrihydrite structure similar to that of ferritin cores. The haemosiderin isolated from secondary iron-overload shows anomalous behaviour in its temperature-dependent Mössbauer spectra. This can be understood in terms of the microcrystalline goethite structure of the cores as indicated by electron diffraction. The haemosiderin cores obtained in the case of primary haemochromatosis have an amorphous Fe(III) oxide structure and show Mössbauer spectra characteristic of a magnetically disordered material, which only orders at very low temperatures.
Collapse
Affiliation(s)
- D P Dickson
- Department of Physics, University of Liverpool, U.K
| | | | | | | | | | | |
Collapse
|
21
|
Abstract
Haemosiderin iron cores isolated from patients with secondary haemochromatosis have a goethite-like (alpha-FeOOH) crystal structure whereas those from patients with primary haemochromatosis are amorphous Fe (III) oxide. Haemosiderin cores isolated from normal human spleen are crystalline ferrihydrite (5Fe2O3.9H2O). The disease-specific structures are significantly different from the ferrihydrite structure of associated ferritin cores. The results are important in understanding the biological processing of iron in pathological states and in the clinical treatment of iron-overload diseases.
Collapse
Affiliation(s)
- S Mann
- School of Chemistry, University of Bath, England
| | | | | | | | | | | | | |
Collapse
|
22
|
Lovett AM, Reid NM, Buckley JA, French JB, Cameron DM. Real-time analysis of breath using an atmospheric pressure ionization mass spectrometer. Biomed Mass Spectrom 1979; 6:91-7. [PMID: 369632 DOI: 10.1002/bms.1200060302] [Citation(s) in RCA: 40] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
A new technique is reported resulting in the direct, instantaneous analyses of trace compounds in breath. The analyses were performed using a commercial atmospheric pressure chemical ionization mass spectrometer (TAGA TM2000 APCI mass spectrometer). A known flow of breath sample is introduced into the ionization region of the mass spectrometer. The study includes the measurement and monitoring in real-time, of breath ammonia during a 24 hour and a 48 hour period. The ammonia profiles indicate a personalized daily pattern associated with each subject. This method appears to be of potential value in routine detection and treatment of hyperammonemia patients. Results also show that it is possible to obtain instantaneous analyses of several naturally occurring metabolites and other substances on breath in the ppm to ppt range, suggesting a number of diagnostic research applications.
Collapse
|