Figures and data in Non-invasive real-time genomic monitoring of the critically endangered kākāpō

Figures
Tables
Additional files

4 figures, 2 tables and 4 additional files

Figures

Figure 1

Download asset Open asset

Vertebrate biodiversity in New Zealand from soil eDNA.

(a) Cladogram of all species and genera detected by 12 S rRNA metabarcoding of soil samples from Whenua Hou and from parrot aviaries in the Dunedin Botanic Garden. (b) Relative taxon abundances of sampled locations averaged across replicates (*from left to right*: kākāpō display sites, feeding stations, abandoned nests, and *Nestor* parrot aviaries). Two different sites per location were sampled (*top and bottom*) at three different distances, and two aviaries of the *Nestor* species kea and kākā. For feeding station 2 (4m), both replicates resulted in dropouts.

Figure 2

Download asset Open asset

Resulting read length distribution (log10 scale) of nanopore sequencing of three exemplary soil samples (samples 3, 11 and 35; Supplementary file 1).

*Left*: Distribution of all passed (Q-score >7) reads; *right*: Distribution of all passed (Q-score >7) reads that map to the kākāpō reference genome using minimap2. The subset of mapped reads that have been accepted by selective sequencing (not ‘unblock’ reads; Methods) is highlighted in orange. The selective sequencing results are shown by (a) and (b) (Sample 3), (c) and (d) (Sample 11), (e) and (f) (Sample 35). The non-selective nanopore sequencing data is shown by (g) and (h) (Sample 3) and (i) and (j) (Sample 35). The selective runs result in many reads of ~500 bp length, which is the average sequencing length at which reads are long enough to be taken a decision upon and to be rejected.

Figure 3 with 1 supplement

Download asset Open asset

Individual identification from soil nanopore sequencing data.

(**a–c**) Distribution of haplotype agreement scores between all Whenua Hou kākāpō and (a) soil sample 3 (Moss’ display site), (b) soil sample 11 (Merv’s display site), and (c) soil sample 35 (Nora’s feeding station). (d) Mixing proportions [%; log10 scale] and (e) posterior means of individual assignment per sample (*y-axis*) assessed through Bayesian inference of individual assignments (see Materials and methods). The heatmaps show Sinbad’s omnipresent signal in the first column, the best hit when disregarding Sinbad in the second column, the second-best hit in the third column, and the mean values of all remaining Whenua Hou kākāpō in the last column.

Figure 3—figure supplement 1

Download asset Open asset

Stamen terrain map of the sampling sites of the three soil samples 3 (Moss’ site), 11 (Merv’s site), and 35 (Nora’s site), and of the radio transmitter signal receiver that recorded Sinbad’s presence far away from his home range (‘Sinbad’s bowl’) and in the middle of the sampling sites, on February 24, 2019, three days before our sampling efforts.

A terrain map with blurred contour lines is used to not disclose the exact home ranges of the individuals of this critically endangered species.

Figure 4

Download asset Open asset

The critically endangered kakapo and its New Zealand habitat.

(a) A kākāpō (*picture credit*: Lydia Uddstrom). (b) Map enhancement of sampling locations on Whenua Hou, New Zealand (*service layer credit*: Esri, Maxar, GeoEye, Earthstar, Geographics, CNES/Airbus DS, USDA, USGS, AeroGRID, IGN, and the GIS User Community).

Tables

Table 1

Number of passed nanopore reads and bases (Q-score >7), number of reads and bases mapping to the kākāpō reference genome, and relative amount [%] of mapped reads and bases per soil sample.

For samples 3 and 35, the results of the selective and the non-selective nanopore sequencing runs are shown.

Sample	Run	# reads	# bases	# mapped reads	# mapped bases	% mapped reads	% mapped bases
Sample 3	selective	903,250	313,408,229	137	67,480	0.015	0.022
	non-selective	3,942,038	6,134,919,092	2,739	7,983,165	0.069	0.130
Sample 11	selective	4,914,329	2,062,605,132	769	435,963	0.016	0.021
Sample 35	selective	2,151,935	802,633,273	563	324,114	0.026	0.040
	non-selective	4,133,117	7,708,401,775	3,854	10,720,441	0.093	0.139

Table 2

Details of the three soil samples subjected to nanopore sequencing.

DNA concentration after bead clean-up [ng/ul]; volume used as input for library preparation [ul] to achieve a DNA input amount of 1 µg per library preparation; amount of DNA in the final library [ng] used as input for sequencing; number of active pores per nanopore flow cell; and metadata of each sample.

Sample number	Concentration after clean-up [ng/ul]	Volume library preparation [ul]	Final amount of DNA in library [ng]	# pores	Metadata
3	91.2	12.3	584	1547	Display site of kākāpō individual Moss
11	202.0	5.0	257	1374	Display site of kākāpō individual Merv
35	119.0	12.0	615	1712	Feeding station of kākāpō individual Nora

Additional files

Supplementary file 1 Information about soil samples including details about sampling and DNA extraction. Sample ID, Sample Type, Distance from the sampling site [m], Amount of soil sampled [ug], DNA concentration after DNA extraction [ug/ml], Date of extraction (i.e. extraction batch), and Sample label, which is used throughout the manuscript.: https://cdn.elifesciences.org/articles/84553/elife-84553-supp1-v1.xlsx
Download elife-84553-supp1-v1.xlsx
Supplementary file 2 Number of 12 S rRNA reads per soil sample. All reads, reads that could be classified, reads that map to any Metazoan taxa, and reads that mapped to the taxonomic species level.: https://cdn.elifesciences.org/articles/84553/elife-84553-supp2-v1.xlsx
Download elife-84553-supp2-v1.xlsx
Supplementary file 3 Number of 12 S rRNA reads mapping to each detected species and genus (rows) across all soil samples (columns). For each taxon, the taxonomic rank, scientific and common name, and the NCBI taxonomic ID are indicated.: https://cdn.elifesciences.org/articles/84553/elife-84553-supp3-v1.xlsx
Download elife-84553-supp3-v1.xlsx
MDAR checklist: https://cdn.elifesciences.org/articles/84553/elife-84553-mdarchecklist1-v1.pdf
Download elife-84553-mdarchecklist1-v1.pdf