Module 7

Conservation Biochemistry

Biomarkers, eDNA, pollution biochemistry, bioremediation, and endocrine disruption

Conservation biology increasingly relies on biochemical tools to assess ecosystem health, detect hidden biodiversity, and remediate contaminated environments. From stress hormones in wildlife hair to environmental DNA in stream water, from mercury biomagnification through food webs to fungal enzymes that degrade petroleum, the molecular level reveals what traditional surveys cannot. This module explores how biochemistry serves as both diagnostic (detecting problems) and therapeutic (implementing solutions) in conservation.

7.1 Biomarkers of Ecosystem Health

A biomarker is any measurable biochemical, physiological, or morphological change in an organism that indicates exposure to or effects of environmental stressors. Unlike chemical monitoring of water or soil, biomarkers integrate cumulative exposure and reveal the biological significance of contamination.

Cortisol: The Stress Hormone

Cortisol, a glucocorticoid synthesized from cholesterol via the hypothalamic-pituitary-adrenal (HPA) axis, is elevated in chronically stressed wildlife. It can be measured non-invasively:

Hair/fur cortisol: Integrates HPA activity over weeks-months as cortisol is incorporated during hair growth. Elevated in bears near roads, wolves in fragmented habitat
Feather corticosterone: In birds, corticosterone (the avian equivalent) deposits in feathers during growth. Reflects breeding season stress
Fecal glucocorticoid metabolites: Non-invasive sampling of scat; metabolites reflect circulating levels with ~12–24 hr lag

The biosynthetic pathway: cholesterol is cleaved by CYP11A1 (cytochrome P450 side-chain cleavage enzyme) to pregnenolone, then sequentially hydroxylated by CYP17A1, CYP21A2, and CYP11B1 to produce cortisol. Chronic elevation suppresses immune function (via NF-\(\kappa\)B inhibition) and reproduction (via GnRH suppression).

Vitellogenin in Male Fish: Endocrine Disruption Marker

Vitellogenin (VTG) is an egg yolk precursor protein normally produced only by females under estrogen stimulation. Detection of VTG in male fish indicates exposure to estrogenic compounds in the water — a classic biomarker of endocrine disruption. Male trout downstream of sewage treatment plants can have VTG levels \(>10^6\)-fold higher than unexposed males.

Chlorophyll Fluorescence: Plant Stress Indicator

When photosystem II (PSII) absorbs light, the energy can follow three pathways: photochemistry, heat dissipation, or fluorescence re-emission. The ratio of variable to maximum fluorescence (\(F_v/F_m\)) measures PSII quantum yield:

PSII Maximum Quantum Yield

\[ \frac{F_v}{F_m} = \frac{F_m - F_0}{F_m} \]

where \(F_0\) is minimum fluorescence (all PSII reaction centers open, dark-adapted) and \(F_m\) is maximum fluorescence (all centers closed by a saturating pulse).

In healthy plants, \(F_v/F_m \approx 0.83\) (remarkably conserved across species). This value reflects the intrinsic efficiency of PSII photochemistry. The derivation:

\[ F_0 = \frac{k_f}{k_f + k_d + k_p} \cdot I_{\text{abs}} \]

\[ F_m = \frac{k_f}{k_f + k_d} \cdot I_{\text{abs}} \]

\[ \frac{F_v}{F_m} = 1 - \frac{F_0}{F_m} = 1 - \frac{k_f + k_d}{k_f + k_d + k_p} \cdot \frac{k_f + k_d}{k_f} \cdot \frac{k_f}{k_f + k_d + k_p} \]

Simplifying with \(k_f\) (fluorescence rate constant), \(k_d\)(thermal dissipation), and \(k_p\) (photochemistry):

\[ \frac{F_v}{F_m} = \frac{k_p}{k_f + k_d + k_p} \]

Stress (drought, heavy metals, ozone) damages the D1 protein of PSII, reducing \(k_p\)and causing \(F_v/F_m\) to decline. Values below 0.75 indicate significant photoinhibition; below 0.60 indicates severe stress.

7.2 Environmental DNA (eDNA)

Organisms constantly shed DNA into their environment — through skin cells, mucus, feces, urine, gametes, and decomposing tissue. This environmental DNA(eDNA) persists in water and soil for days to weeks and can be collected, amplified by PCR, and sequenced to identify species present without ever observing them directly.

eDNA Decay Kinetics

After an organism leaves an area, its eDNA degrades through nuclease activity, UV photolysis, microbial consumption, and hydrolysis. The decay follows first-order kinetics:

eDNA Decay Model

\[ C(t) = C_0 \cdot e^{-\lambda t} \]

where \(C_0\) is initial concentration, \(\lambda\) is the decay rate constant (d\(^{-1}\)), and the half-life is \(t_{1/2} = \ln 2 / \lambda\).

Empirical measurements show:

Freshwater: \(t_{1/2} \approx 1\text{--}7\) days. UV exposure and warm temperature accelerate decay
Marine water: \(t_{1/2} \approx 0.5\text{--}7\) days. Salinity and dilution are additional factors
Soil: \(t_{1/2} \approx 7\text{--}120\) days. DNA binds to clay minerals, slowing degradation
Sediment: \(t_{1/2} \approx 14\text{--}300\) days. Anoxic conditions dramatically slow nuclease activity

The decay rate depends on temperature via the Arrhenius equation:

\[ \lambda(T) = \lambda_0 \cdot \exp\!\left(-\frac{E_a}{R}\left(\frac{1}{T} - \frac{1}{T_0}\right)\right) \]

where \(E_a \approx 50\text{--}80\) kJ/mol for DNA hydrolysis. At 25\(^\circ\)C, decay is roughly 2\(\times\) faster than at 15\(^\circ\)C.

Detection Probability: Multiple Sampling

A single water sample may fail to detect a rare species even when eDNA is present, due to low concentration and stochastic sampling. The probability of detecting a species with \(n\) independent samples, each with per-sample detection probability\(p\), is:

Multi-Sample Detection Probability

\[ P_{\text{detect}} = 1 - (1 - p)^n \]

For a species with per-sample \(p = 0.3\), taking \(n = 5\) samples gives\(P_{\text{detect}} = 1 - 0.7^5 = 0.832\), and \(n = 10\) gives 0.972.

The minimum number of samples needed to achieve detection probability \(P^*\):

\[ n \geq \frac{\ln(1 - P^*)}{\ln(1 - p)} \]

For \(P^* = 0.95\) and \(p = 0.2\): \(n \geq \ln(0.05)/\ln(0.8) = 13.4\), so 14 samples are needed.

Metabarcoding Pipeline

eDNA metabarcoding uses universal primers (e.g., 12S rRNA for vertebrates, COI for invertebrates, ITS for fungi, 16S for bacteria) to amplify DNA from all species simultaneously. The workflow: water filtration (0.22–0.45 \(\mu\)m) \(\to\) DNA extraction \(\to\) PCR amplification \(\to\) high-throughput sequencing \(\to\) bioinformatics (ASV/OTU clustering, taxonomic assignment against reference databases like BOLD or GenBank).

7.3 Pollution Biochemistry

Heavy Metal Toxicity Mechanisms

Heavy metals disrupt biochemistry primarily through ionic mimicry — toxic metals substitute for essential ones in enzyme active sites:

Cadmium (Cd\(^{2+}\)): Replaces Zn\(^{2+}\) in metalloenzymes (carbonic anhydrase, carboxypeptidase, alcohol dehydrogenase). Same charge, similar ionic radius (0.95 vs 0.74 \(\text{\AA}\)), but Cd cannot perform the catalytic function. Also binds thiol groups in metallothioneins.
Lead (Pb\(^{2+}\)): Replaces Ca\(^{2+}\) in calmodulin and bone. Inhibits \(\delta\)-aminolevulinic acid dehydratase (ALAD) in heme biosynthesis, causing anemia. Also inhibits NMDA receptors in neurons.
Arsenic (As): As(V) mimics phosphate, entering cells via phosphate transporters and uncoupling oxidative phosphorylation. As(III) binds vicinal thiol groups in lipoic acid, inhibiting pyruvate dehydrogenase.

Mercury Methylation and Biomagnification

Inorganic mercury (Hg\(^{2+}\)) is converted to methylmercury (CH₃Hg\(^+\)) by sulfate-reducing bacteria in anoxic sediments. The methylation reaction:

\[ \text{Hg}^{2+} + \text{CH}_3\text{-Co(III)-corrinoid} \to \text{CH}_3\text{Hg}^+ + \text{Co(I)-corrinoid} \]

The methyl group is transferred from methylcobalamin (vitamin B₁₂ derivative) to mercury via the hgcAB gene cluster in anaerobic microbes.

Methylmercury biomagnifies through food webs because it is lipophilic and binds cysteine residues in proteins, creating a long biological half-life (~70 days in fish). The bioconcentration factor (BCF) quantifies accumulation:

Bioconcentration Factor

\[ \text{BCF} = \frac{C_{\text{organism}}}{C_{\text{water}}} = \frac{k_{\text{uptake}}}{k_{\text{elimination}}} \]

where \(k_{\text{uptake}}\) and \(k_{\text{elimination}}\) are the first-order rate constants for absorption and excretion, respectively.

Biomagnification occurs when BCF increases with trophic level. For methylmercury, the trophic magnification factor (TMF) is approximately 3–10\(\times\) per trophic level. A simple model for concentration at trophic level \(n\):

\[ C_n = C_{\text{water}} \cdot \text{BCF}_1 \cdot \text{TMF}^{(n-1)} \]

With \(C_{\text{water}} = 1\) ng/L, BCF\(_1 = 10^4\) for phytoplankton, and TMF = 5: trophic level 4 (large predatory fish) reaches \(10^4 \times 5^3 = 1.25 \times 10^6\) ng/kg = 1.25 mg/kg, exceeding the 0.3 mg/kg FDA action level.

7.4 Bioremediation Pathways

Bioremediation uses living organisms or their enzymes to degrade, transform, or sequester environmental contaminants. Three major approaches:

Phytoremediation: Hyperaccumulator Plants

Approximately 700 plant species are known to hyperaccumulate heavy metals, concentrating them at levels 100–1000\(\times\) higher than normal plants:

Thlaspi caerulescens: Zn/Cd hyperaccumulator. HMA4 transporter (P-type ATPase) loads metals into xylem; vacuolar sequestration via phytochelatins
Pteris vittata: As hyperaccumulator (brake fern). Reduces As(V) to As(III) and stores in vacuoles complexed with thiol peptides
Alyssum bertolonii: Ni hyperaccumulator. Chelation by histidine and organic acids (citrate, malate) in xylem sap

Microbial Petroleum Degradation

Hydrocarbon-degrading bacteria (e.g., Alcanivorax borkumensis, Pseudomonas putida) oxidize alkanes using alkane hydroxylase (AlkB) and cytochrome P450 (CYP153):

\[ \text{R-CH}_3 + \text{O}_2 + \text{NAD(P)H} + \text{H}^+ \xrightarrow{\text{AlkB}} \text{R-CH}_2\text{OH} + \text{NAD(P)}^+ + \text{H}_2\text{O} \]

The terminal hydroxylation is followed by sequential oxidation to aldehyde, then carboxylic acid, which enters \(\beta\)-oxidation for energy production.

Mycoremediation: Fungal Laccase

White-rot fungi (e.g., Phanerochaete chrysosporium, Trametes versicolor) produce extracellular laccase and lignin peroxidase that oxidize a remarkable range of pollutants, including PAHs, PCBs, dioxins, and pesticides. These enzymes evolved to degrade lignin — the most recalcitrant natural polymer — and their non-specific radical mechanism allows them to attack structurally similar synthetic pollutants.

Biodegradation Kinetics: Monod Model

Pollutant degradation by microbial communities follows Monod kinetics, coupling substrate consumption to microbial growth:

Monod Biodegradation Model

\[ \frac{dS}{dt} = -\frac{\mu_{\max} \cdot S}{K_s + S} \cdot \frac{X}{Y} \]

\[ \frac{dX}{dt} = \frac{\mu_{\max} \cdot S}{K_s + S} \cdot X - b \cdot X \]

where \(S\) = substrate concentration, \(X\) = microbial biomass,\(\mu_{\max}\) = maximum specific growth rate, \(K_s\) = half-saturation constant,\(Y\) = yield coefficient (biomass/substrate), and \(b\) = endogenous decay rate.

At low substrate concentrations (\(S \ll K_s\)), the kinetics become pseudo-first-order:

\[ \frac{dS}{dt} \approx -\frac{\mu_{\max}}{K_s \cdot Y} \cdot X \cdot S = -k_1 \cdot S \]

where \(k_1 = \mu_{\max} X / (K_s Y)\) is the apparent first-order rate constant, giving exponential decay: \(S(t) = S_0 e^{-k_1 t}\).

7.5 Endocrine Disruptors

Endocrine-disrupting chemicals (EDCs) are exogenous substances that interfere with hormonal signaling, even at extremely low concentrations (parts per trillion). The most studied class are xenoestrogens — synthetic chemicals that mimic 17\(\beta\)-estradiol (E2) at the estrogen receptor (ER).

Major Xenoestrogens

Bisphenol A (BPA): Plasticizer in polycarbonate and epoxy resins. Binds ER\(\alpha\) with RBA ~0.01–0.1% relative to E2
Atrazine: Herbicide. Induces aromatase (CYP19) expression, converting testosterone to estrogen in male amphibians, causing feminization
Ethinyl estradiol (EE2): Synthetic estrogen in contraceptive pills, the most potent EDC in wastewater. RBA ~100% (equal to E2)
PCBs: Persistent organic pollutants. Some congeners are estrogenic, others anti-estrogenic. Bioaccumulate in adipose tissue
Nonylphenol: Surfactant degradation product. Binds ER with RBA ~0.01%. Ubiquitous in wastewater

Competitive Binding: IC50 and Relative Binding Affinity

The binding of a xenoestrogen \(X\) to the estrogen receptor competes with the natural ligand E2. In a competitive binding assay, increasing concentrations of \(X\)displace radiolabeled E2 from the receptor:

Competitive Binding Model

\[ \frac{[\text{ER} \cdot \text{E2}]}{[\text{ER}]_{\text{total}}} = \frac{1}{1 + \frac{[\text{X}]}{K_i}} \cdot \frac{[\text{E2}]}{K_d + [\text{E2}]} \]

The \(\text{IC}_{50}\) is the concentration of competitor \(X\) that displaces 50% of bound E2. The relative binding affinity (RBA) normalizes this to E2:

Relative Binding Affinity

\[ \text{RBA} = \frac{\text{IC}_{50}(\text{E2})}{\text{IC}_{50}(\text{X})} \times 100\% \]

E2 has RBA = 100% by definition. BPA has RBA \(\approx\) 0.01–0.1%, meaning it binds ER 1000–10,000\(\times\) more weakly than E2. However, environmental concentrations of BPA can be \(10^4\text{--}10^6\times\) higher than E2, partially compensating for weak affinity.

The relationship to the Cheng-Prusoff equation for deriving \(K_i\) from \(\text{IC}_{50}\):

\[ K_i = \frac{\text{IC}_{50}}{1 + \frac{[\text{E2}]}{K_d}} \]

where \(K_d\) is the dissociation constant for E2-ER binding (~0.1–1 nM).

7.6 Mercury Biomagnification Pyramid

The diagram below illustrates how methylmercury concentration increases approximately 10-fold at each trophic level, from dissolved mercury in water through phytoplankton, zooplankton, small fish, large predatory fish, and finally to apex predators like bald eagles.

7.7 Computational Simulations

Four simulations exploring conservation biochemistry: (1) eDNA decay in different environments, (2) bioconcentration factors across taxa, (3) Monod-model biodegradation kinetics, and (4) mercury biomagnification through a food web.

Conservation Biochemistry: eDNA, BCF, Biodegradation & Biomagnification

Python

script.py163 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

np.random.seed(42)
fig, axes = plt.subplots(2, 2, figsize=(14, 11))

# ====== Panel 1: eDNA Decay in Different Environments ======
ax = axes[0, 0]
t = np.linspace(0, 60, 200)  # days

environments = {
    'Warm stream (25\u00b0C)': {'lambda': 0.7, 'color': '#ef4444'},
    'Cool lake (15\u00b0C)': {'lambda': 0.2, 'color': '#f59e0b'},
    'Marine (20\u00b0C)': {'lambda': 0.35, 'color': '#3b82f6'},
    'Soil': {'lambda': 0.05, 'color': '#a78bfa'},
    'Sediment': {'lambda': 0.015, 'color': '#2dd4bf'},
}

C0 = 1.0  # normalized initial concentration
for env, params in environments.items():
    C = C0 * np.exp(-params['lambda'] * t)
    half_life = np.log(2) / params['lambda']
    ax.plot(t, C, color=params['color'], linewidth=2.5,
            label=f"{env} (t\u00bd={half_life:.1f}d)")

# Detection threshold
ax.axhline(y=0.01, color='#f97316', linewidth=1, linestyle='--', alpha=0.7)
ax.text(45, 0.015, 'Detection threshold', color='#f97316', fontsize=10)

ax.set_xlabel('Time (days)', fontsize=12, color='#99f6e4')
ax.set_ylabel('Relative eDNA Concentration', fontsize=12, color='#99f6e4')
ax.set_title('eDNA Decay Kinetics', fontsize=14, fontweight='bold', color='#5eead4')
ax.set_yscale('log')
ax.set_ylim(0.001, 2)
ax.legend(fontsize=8.5, facecolor='#0a0a0a', edgecolor='#5eead4',
          labelcolor='#d1d5db', loc='upper right')
ax.set_facecolor('#0a0a0a')
ax.tick_params(colors='#94a3b8')
ax.spines['bottom'].set_color('#5eead4')
ax.spines['left'].set_color('#5eead4')
ax.spines['top'].set_visible(False)
ax.spines['right'].set_visible(False)
ax.grid(True, alpha=0.15, color='#5eead4')

# ====== Panel 2: BCF Across Species ======
ax = axes[0, 1]
organisms = ['Algae', 'Daphnia', 'Snail', 'Minnow', 'Trout', 'Pike', 'Osprey']
bcf_hg = [1e4, 5e4, 8e4, 2e5, 5e5, 1e6, 3e6]
bcf_pcb = [5e3, 2e4, 3e4, 1e5, 3e5, 8e5, 2e6]
bcf_ddt = [8e3, 3e4, 5e4, 1.5e5, 4e5, 9e5, 2.5e6]

x = np.arange(len(organisms))
width = 0.25
ax.bar(x - width, bcf_hg, width, color='#ef4444', label='MeHg', edgecolor='#0a0a0a')
ax.bar(x, bcf_pcb, width, color='#3b82f6', label='PCBs', edgecolor='#0a0a0a')
ax.bar(x + width, bcf_ddt, width, color='#2dd4bf', label='DDT', edgecolor='#0a0a0a')

ax.set_xticks(x)
ax.set_xticklabels(organisms, rotation=30, ha='right')
ax.set_yscale('log')
ax.set_xlabel('Organism', fontsize=12, color='#99f6e4')
ax.set_ylabel('Bioconcentration Factor', fontsize=12, color='#99f6e4')
ax.set_title('BCF Across Trophic Levels', fontsize=14, fontweight='bold', color='#5eead4')
ax.legend(fontsize=10, facecolor='#0a0a0a', edgecolor='#5eead4', labelcolor='#d1d5db')
ax.set_facecolor('#0a0a0a')
ax.tick_params(colors='#94a3b8')
ax.spines['bottom'].set_color('#5eead4')
ax.spines['left'].set_color('#5eead4')
ax.spines['top'].set_visible(False)
ax.spines['right'].set_visible(False)
ax.grid(True, alpha=0.15, color='#5eead4', axis='y')

# ====== Panel 3: Monod Biodegradation Kinetics ======
ax = axes[1, 0]
# Solve Monod equations numerically
dt_sim = 0.01  # days
t_total = 30
n_steps = int(t_total / dt_sim)

# Parameters
mu_max = 0.5   # per day
Ks = 50.0      # mg/L
Y = 0.4        # biomass yield
b = 0.05       # decay rate

# Initial conditions for different scenarios
scenarios = [
    {'S0': 200, 'X0': 5, 'label': 'High substrate (200 mg/L)', 'color': '#2dd4bf'},
    {'S0': 100, 'X0': 5, 'label': 'Medium substrate (100 mg/L)', 'color': '#f59e0b'},
    {'S0': 50, 'X0': 5, 'label': 'Low substrate (50 mg/L)', 'color': '#a78bfa'},
]

for sc in scenarios:
    S = np.zeros(n_steps)
    X = np.zeros(n_steps)
    t_arr = np.zeros(n_steps)
    S[0] = sc['S0']
    X[0] = sc['X0']

for i in range(1, n_steps):
        mu = mu_max * S[i-1] / (Ks + S[i-1])
        dS = -(mu / Y) * X[i-1] * dt_sim
        dX = (mu - b) * X[i-1] * dt_sim
        S[i] = max(0, S[i-1] + dS)
        X[i] = max(0, X[i-1] + dX)
        t_arr[i] = i * dt_sim

ax.plot(t_arr, S, color=sc['color'], linewidth=2.5, label=sc['label'])
    ax.plot(t_arr, X, color=sc['color'], linewidth=1.5, linestyle='--', alpha=0.6)

ax.set_xlabel('Time (days)', fontsize=12, color='#99f6e4')
ax.set_ylabel('Concentration (mg/L)', fontsize=12, color='#99f6e4')
ax.set_title('Monod Biodegradation (solid=S, dashed=X)', fontsize=13, fontweight='bold', color='#5eead4')
ax.legend(fontsize=9, facecolor='#0a0a0a', edgecolor='#5eead4', labelcolor='#d1d5db')
ax.set_facecolor('#0a0a0a')
ax.tick_params(colors='#94a3b8')
ax.spines['bottom'].set_color('#5eead4')
ax.spines['left'].set_color('#5eead4')
ax.spines['top'].set_visible(False)
ax.spines['right'].set_visible(False)
ax.grid(True, alpha=0.15, color='#5eead4')

# ====== Panel 4: Hg Biomagnification Through Food Web ======
ax = axes[1, 1]
trophic_levels = ['Water', 'Phytopl.', 'Zoopl.', 'Sm. Fish', 'Lg. Fish', 'Eagle']
hg_conc = [0.001, 0.01, 0.05, 0.2, 1.0, 10.0]  # ug/g or ug/L
tl_numbers = [0, 1, 2, 3, 4, 5]

colors_bar = ['#164e63', '#115e59', '#134e4a', '#0f766e', '#0d9488', '#14b8a6']
bars = ax.bar(trophic_levels, hg_conc, color=colors_bar, edgecolor='#5eead4', linewidth=1)

# FDA limit line
ax.axhline(y=0.3, color='#ef4444', linewidth=2, linestyle='--', alpha=0.8)
ax.text(4.3, 0.4, 'FDA action\nlevel 0.3 \u00b5g/g', color='#ef4444', fontsize=10,
        ha='center')

ax.set_yscale('log')
ax.set_xlabel('Trophic Level', fontsize=12, color='#99f6e4')
ax.set_ylabel('[MeHg] (\u00b5g/g or \u00b5g/L)', fontsize=12, color='#99f6e4')
ax.set_title('Mercury Biomagnification', fontsize=14, fontweight='bold', color='#5eead4')
ax.set_facecolor('#0a0a0a')
ax.tick_params(colors='#94a3b8')
ax.spines['bottom'].set_color('#5eead4')
ax.spines['left'].set_color('#5eead4')
ax.spines['top'].set_visible(False)
ax.spines['right'].set_visible(False)
ax.grid(True, alpha=0.15, color='#5eead4', axis='y')

# Add concentration labels on bars
for bar, conc in zip(bars, hg_conc):
    ax.text(bar.get_x() + bar.get_width()/2, conc * 1.5,
            f'{conc}', ha='center', va='bottom', color='#f97316',
            fontsize=10, fontweight='bold')

fig.patch.set_facecolor('#0a0a0a')
plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight',
            facecolor='#0a0a0a', edgecolor='none')
plt.close()
print("Conservation biochemistry simulations complete: eDNA decay, BCF, biodegradation, biomagnification.")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

References

Maceda-Veiga, A., Figuerola, J., Martínez-Silvestre, A. & Green, A.J. (2015). Inside the redbox: applications of haematology in wildlife monitoring and ecosystem health. Science of The Total Environment, 514, 322–332.
Sumpter, J.P. & Jobling, S. (1995). Vitellogenesis as a biomarker for estrogenic contamination of the aquatic environment. Environmental Health Perspectives, 103(Suppl 7), 173–178.
Baker, N.R. (2008). Chlorophyll fluorescence: a probe of photosynthesis in vivo. Annual Review of Plant Biology, 59, 89–113.
Thomsen, P.F. & Willerslev, E. (2015). Environmental DNA — an emerging tool in conservation for monitoring past and present biodiversity. Biological Conservation, 183, 4–18.
Deiner, K., Bik, H.M., Mächler, E., Seymour, M., Lacoursière-Roussel, A., Altermatt, F., ... & Bernatchez, L. (2017). Environmental DNA metabarcoding: transforming how we survey animal and plant communities. Molecular Ecology, 26(21), 5872–5895.
Tchounwou, P.B., Yedjou, C.G., Patlolla, A.K. & Sutton, D.J. (2012). Heavy metal toxicity and the environment. Experientia Supplementum, 101, 133–164.
Hsu-Kim, H., Kucharzyk, K.H., Zhang, T. & Deshusses, M.A. (2013). Mechanisms regulating mercury bioavailability for methylating microorganisms in the aquatic environment. Environmental Science & Technology, 47(6), 2441–2456.
Pilon-Smits, E. (2005). Phytoremediation. Annual Review of Plant Biology, 56, 15–39.
Pointing, S.B. (2001). Feasibility of bioremediation by white-rot fungi. Applied Microbiology and Biotechnology, 57(1), 20–33.
Gore, A.C., Chappell, V.A., Fenton, S.E., Flaws, J.A., Nadal, A., Prins, G.S., ... & Zoeller, R.T. (2015). EDC-2: The Endocrine Society’s second scientific statement on endocrine-disrupting chemicals. Endocrine Reviews, 36(6), E1–E150.

M6. Biodiversity-Stability M8. Climate Change & Adaptation