LECHYTIIDAE

CHTHONIIDAE

Figure 32. Maximum Likelihood of representative COI haplotypes of pseudoscorpionida from this study and a sub-set of reference specimens from WAM, Helix and GenBank. Proposed lineages from this study is indicated by the shaded, coloured boxes. Report by Western Australian Museum (Molecular Systematics Unit)

PSEUDOSCORPIONES

In this study a total of three specimens were extracted for molecular analysis, all of which produced good quality COI sequences. A summary of specimen identifications is depicted in Table 11. A total of three distinct lineages/species belonging to 2 different families, were detected using Maximum likelihood analysis (Figure 32), none of which matched known lineages. The two chthoniid lineages, identified by Bennelongia as belonging to the genus Tyrannochthonius, ranged between 11.4–15.4% genetic difference between Helix lineage GC1 (Koodaideri) and an unknown species from GenBank respectively. The lineage of Lechytia (LECHYTIIDAE) was also unique and differed from described species from GenBank (Lechytia sini; USA) by 19.2%. The family Lechytiidae includes only a single genus, Lechytia, which currently contains 24 described species worldwide, and only one species described from Australia (Lechytia libita; QLD; Harvey, 2006). Unfortunately there are no available sequences of Lechytia from Australia or from a number of undescribed species lodged with the WA Museum. The distribution of these lineages within the project site is depicted in Figure 33.

Table 11. Summary of lineages detected in the Project area. % Divergence from closest known lineage is shown with lineage name in parentheses.

Family Genus Species/Lineage Specimen Divergence RTIO

Chthoniidae Tyrannochthonius WAM-CHTH001 1 11.4% (GC1:Helix; Koodaideri) RT43 Chthoniidae Tyrannochthonius WAM-CHTH002 1 15.4% (Austrochthonius EU559513.1) RT56 Lechytiidae Lechytia WAM-LECYT001 1 19.2% (Lechytia sini KX263367.1) RT56

Figure 33. Map of pseudoscorpion lineages (CHTHONIIDAE & LECHYTIIDAE) detected within the project area.

Figure 34. Maximum Likelihood of representative COI haplotypes of from this study and a sub-set of reference specimens from WAM, Helix and GenBank. Proposed lineages from this study is indicated by the shaded, coloured box. Report by Western Australian Museum (Molecular Systematics Unit)

SCHIZOMIDA

A summary of the schizomid specimen identifications is depicted in Table 12. In this study, a total of three specimens were extracted for molecular analysis and produced good quality COI sequences. One new distinct lineage was detected from the project area using Maximum likelihood analysis (Figure 34). `WAM-DRAC001` is 13% different to the next closest lineage (EU272697.1) and around 14.2% different from described species Draculoides julianneae (EU272696.1) and (EU272684.1). The distribution of this lineage within the project site is depicted in Figure 29.

Table 12. Summary of the Schizomida lineage detected within the Project Area. % Divergence from closest known lineage is shown with lineage name in parentheses.

Family Genus Species/Lineage Specimen Divergence RTIO

Hubbardiidae Draculoides WAM-DRAC001 3 13% (EU272697.1 P sp) RT43

MYRIAPODA

CHILOPODA & DIPLOPODA

Specimens of Chilopoda or Diplopoda collected either failed to sequence or were not sub-sampled for molecular analysis and were identified taxonomically by Bennelongia (Appendix I)

PAUROPODA (TETRAMEROCERATA)

A summary of Pauropoda molecular identifications is depicted in Table 13. In this study 10 specimens likely belonging to the genus Decapauropus (Pauropodidae) were extracted for molecular analysis, and eight produced a good quality sequences. Five lineages were detected with large genetic divergences between species, ranging from 12.4–26.6% genetic divergence between the next closest lineage (Figure 35; Appendix II).

Very little molecular work has been undertaken for Pauropoda and there are no published papers on the molecular systematics of Pauropoda in Australia. GenBank sequences have been utilised from research undertaken by Qian et al. (unpublished), who examined Pauropoda from China and found deep barcode divergences between species. The distribution of this lineage within the project site is depicted in Figure 36.

Figure 35. Maximum Likelihood of representative COI haplotypes of Pauropoda from this study and a sub-set of reference specimens from GenBank. Proposed lineages from this study is indicated by the shaded, coloured boxes. Report by Western Australian Museum (Molecular Systematics Unit)

Table 13. Summary of Pauropoda lineages detected from the Project area. % Divergence from closest known lineage is shown with lineage name in parentheses.

Family Genus Species/Lineage Specimen Divergence RTIO Pauropodidae Decapauropus WAM-PAUD001 1 17.4% (WAM-PAUD002) RT43 Pauropodidae Decapauropus WAM-PAUD002 4 17.4% (WAM-PAUD001) RT43 Pauropodidae Decapauropus WAM-PAUD003 1 26.6% (Decapauropus crassescens JN974348.1) RT43 Pauropodidae Decapauropus WAM-PAUD004 1 12.4% (WAM-PAUD005) RT43 Pauropodidae Decapauropus WAM-PAUD005 1 12.4% (WAM-PAUD004) RT43

Figure 36. Map of Pauropoda lineages (PAUROPODIDAE) detected within the project area.

Report by Western Australian Museum (Molecular Systematics Unit)

SYMPHYLA (CEPHALOSTIGMATA)

A summary of the symphylan lineages detected from the project area is depicted in Table 14. In this study, seven specimens that were submitted for molecular analysis produced viable COI sequences after we optimised PCR protocols (unpublished; Figure 37). One additional specimen showed contamination and was eliminated from the analysis, one was not sub-sampled, and another was missing from the vial (WAMT146615). A distribution map of the lineages detected from the project is depicted in Figure 38, and images of WAMT144260 from Channar Mine is depicted in Figure 39.

Table 14. Summary of Symphylan lineages detected in this study. % Divergence from closest known lineage is shown with lineage name in parentheses.

Family Genus Species/Lineage Specimens Divergence from closest lineage RTIO Task

Scutigerellidae Scutigerella `WAM-SCUTI002` 1 25.6% (WAM-SYMPH001, Robe River) RTIO035 Scutigerellidae Scutigerella `WAM-SCUTI003` 2 16.1% (WAM-SCUTI005) RTIO043 Scutigerellidae Scutigerella `WAM-SCUTI004` 1 18.5% (WAM-SCUTI005) RTIO043 Scutigerellidae Scutigerella `WAM-SCUTI005` 1 16.1% (WAM-SCUTI003) RTIO035 Scutigerellidae Symphyella `WAM-SYMPH002` 1 25.6% (WAM-SYMPH001, Robe River) RTIO046

It is clear from the analysis that these specimens represent five distinct lineages from two genera and differed from the closest lineage from between 16.1–25.6% (Figure 37; Appendix II). Two specimens of the lineage WAM-SCUTI003 (WAMT145461 & T145462) differed from each other by 5.3% and differed from the next closed lineage (WAM-SCUTI005) by 16.1%. The closest lineage to Symphyella WAM-SYMPH002 (WAM144298) was Symphyella WAM-SYMPH001 (WAMT145913) from the Robe River Valley Mesa C, and differed by 25.6% (Figure 37; Appendix II).

Helix Molecular Solutions have also sequenced a number of symphylan specimens from the Robe River Valley and determined a number of distinct lineages (Unpublished). These lineages however cannot be compared to sequences produced from this molecular project as they used different primer pairs which do not overlap the Barcoding Gene region and were designed specifically for symphylans.

Very little molecular and morphological research has been undertaken for symphylans, and therefore results from this phylogenetic analysis should be used with caution. Very little inference can be made to the higher level classifications or relationships within this group until further molecular and morphological work has been undertaken, within the Pilbara region of WA. Generic level classifications of the Symphyla are based on morphological identifications undertaken by Bennelongia, and are supported at some level with the molecular analysis.

Symphyella

Scutigerella

Figure 37. Maximum Likelihood of all COI haplotypes of Symphylans (SCUTIGERELLIDAE) from the study area and a sub-set of reference specimens from WAM and GenBank. Proposed lineages from this study are indicated by the shaded, coloured boxes. Report by Western Australian Museum (Molecular Systematics Unit)

Figure 38. Distribution map of the Symphylan lineages (SCUTIGERELLIDAE) detected within the project area.

Figure 39. Automontage image a Symphylan lineage, `WAM-SCUTI002` (WAM T14460) collected from Channar Mine, approx. 14km SE Paraburdoo. Dorsal (top) and Lateral (Bottom)

Report by Western Australian Museum (Molecular Systematics Unit)

REFERENCES

Bradbury, J.H. & Williams, W.D. 1996. Freshwater amphipods from Barrow Island, Western Australia. Records of the Australian Museum 48: 33-74.

Bradbury, J.H. & Williams, W.D. 1997. The amphipod (Crustacea) stygofauna of Australia: description of new taxa (Melitidae, Neoniphargidae, Paramelitidae), and a synopsis of known species. Records of the Australian Museum 49 (3): 249-341.

Brown, L., Finston, T., Humphreys, G., Eberhard, S & Pinder, A. 2015. Groundwater oligochaetes show complex genetic patterns of distribution in the Pilbara region of Western Australia. Invertebrate Systematics 29: 405-420.

Castresana, J. 2000. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Molecular Biology and Evolution 17, 540-552.

Cooper, S.J.B., Bradbury, J.H., Saint, K.M., Remko Leys, S., Austin A.D. & Humphreys, W.F. 2007. Subterranean archipelago in the Australian arid zone: mitochondrial DNA phylogeography of amphipods from central Western Australia. Molecular Ecology 16: 1533-1544.

Cullen, K.L. 2018. Report on the molecular systematics of subfauna from Greater Paraburdoo, Western Australia – Initial Sampling. Unpublished report prepared for Rio Tinto (20 March, 2018).

Cullen, K.L. & Harvey, M.S., 2017. Report on the molecular systematics of subfauna from Mesa B, C, and Warramboo, Western Australia. Unpublished report prepared for Rio Tinto (11 December, 2017).

Eberhard, S.M., Halse, S.A. & Humphreys, W.F. 2005. Stygofauna in the Pilbara region, north-west Western Australia: a review. Journal of the Royal Society of Western Australia 88: 167-176.

Fagan-Jeffries, E. 2012. Hidden beneath the surface: exploring the evolution and systematics of Nedsia, a genus of subterranean amphipods from the Pilbara region of Western Australia. Unpublished Honours Thesis, The University of Adelaide, pp 1-49.

Finston, T., Hitchen, Y. & Berry, O. 2017. Report on the molecular systematics of Amphipoda from Mesa H– Phase 4.Unpublished report prepared for Biota Environmental Sciences (4 March 2017).

Finston, T., Hitchen, Y. & Berry, O. 2016a. Report on the molecular systematics of subfauna from the Robe River Valley – Phase 4. Unpublished report prepared for Biota Environmental Sciences (8 Dec 2016).

Report by Western Australian Museum (Molecular Systematics Unit)

Finston, T., Hitchen, Y. & Berry, O. 2016b. Report on the molecular systematics of subfauna from Mesa H and J. Unpublished report prepared for Biota Environmental Sciences (5 June 2016).

Finston, T., Hitchen, Y. & Berry, O. 2015. Report on the molecular systematics of subfauna from the Robe River Valley. Unpublished report prepared for Biota Environmental Sciences (6 Jan 2015).

Finston, T., Johnson, M.S., Eberhard, S.M., Cocking, J.S., McRae, J.M., Halse, S.A. & Knott, B. 2011. A new genus and two new species of groundwater paramelitid amphipods from the Pilbara, Western Australia: a combined molecular and morphological approach. Records of the Western Australian Museum 26: 154-178.

Finston, T., Johnson, M.S. & Knott, B. 2008. A new genus and species of stygobitic paramelitid amphipod from the Pilbara, Western Australia. Records of the Western Australian Museum 24: 295-410.

Finston, T., Johnson, M.S., Humphreys, W.F., Eberhard, S.M. & Halse, S.A. 2007. Cryptic speciation in two widespread subterranean amphipod genera reflects historical drainage patterns in an ancient landscape. Molecular Ecology 16: 355-365.

Finston, T. & Johnson, M.S. 2004. Geographic patterns of genetic diversity in subterranean amphipods of the Pilbara, Western Australia. Marine and Freshwater Research 55: 619-628.

Folmer, O., M. Black, W. Hoeh, R. Lutz, and R. Vrijenhoek (1994). DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Molecular Marine Biology and Biotechnology, 3: 294-299.

Giribet, G., McIntyre, E., Christian, E., Espinasa, L., Ferreira, R.L., Francke, O.F., Harvey, M.S., Isaia, M., Kovac, L., McCutchen, L., Souza, M.F.V.R., Zagmajster, M. 2014. The first phylogenetic analysis of Palpigradi (Arachnida) – the most enigmatic order. Invertebrate Systematics 28: 350-360.

Herbert, P.D.N., S. Ratnasingham., J.R. deWaard. (2003). Barcoding life: cytochrome c oxidase subunit I divergences among closely related species. Proceedings of the Royal Society of London B (supplement) 270:S96-S99.

Karanovic, I & Humphreys, W.F. 2014. Phylogeny and diversity of Timiriaseviinae ostracods (Podocopida, Cytheroidea) with description of one new species from arid Western Australia. Systematics and 12: 93-110.

Kearse, M., Moir, R., Wilson, A., Stones-Havas, S., Cheung, M., Sturrock, S., Buxton, S., Cooper, A., Markowitz, S., Duran, C., Thierer, T., Ashton, B., Mentjies, P., & Drummond, A. (2012). Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28(12): 1647-1649.

Report by Western Australian Museum (Molecular Systematics Unit)

Leijs, R., Bloechl, A. & Koenemann, S. 2011. Bogidiella veneris, A new species of subterranean amphipoda (Bogidiellidae) from Australia, with remarks on the systematics and Biogeography. Journal of Crustacean Biology 31(3): 566-575.

Perina, G., Camacho, A.I., Huey, J., Horwitz, P. & Koenders, A. 2018. Understanding subterranean variablility: the first genus of Bathynellidae (Bathynellacea, Crustacea) from Western Australia described through a morphological and multigene approach. Invertebrate Systematics 32: 423- 447.

Prendini, L., P.Weygoldt & W. C. Wheeler. 2005. Systematics of the Damon variegatus group of African whip (: Amblypygi): Evidence from behaviour, morphology and DNA. Organisms, Diversity & Evolution. 5: 203-236. https://doi.org/10.1016/j.ode.2004.12.004.

Ratnasingham, S. & Herbert, P.D.N. 2007. BOLD: The Barcode of Life Data System (www.barcodinglife.org). Molecular Ecology Notes. Pp. 1-10.

Talavera, G., and Castresana, J. 2007. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Systematic Biology 56, 564-577

Greater Paraburdoo Subterranean Fauna Survey

APPENDIX 4: Western Range 3D Habitat Modelling (AWT) – Long and cross sections

Appendix Page | 57

Greater Paraburdoo Subterranean Fauna Survey

Western Range 3D Habitat Modelling (AWT) – Long and cross sections

Appendix 4 - LS_1: Long section of Western Range (part 1), Current AWT habitat

Appendix Page | 58

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - LS_2: Long section of Western Range (part 1), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 59

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - LS_3: Long section of Western Range (part 2), Current AWT habitat

Appendix Page | 60

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - LS_4: Long section of Western Range (part 2), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 61

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - LS_5: Long section of Western Range (part 3), Current AWT habitat

Appendix Page | 62

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - LS_6: Long section of Western Range (part 3), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 63

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - LS_7: Long section of Western Range (part 4), Current AWT habitat

Appendix Page | 64

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - LS_8: Long section of Western Range (part 4), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 65

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - CS_1: Cross section of Western Range (part 1), Current AWT habitat

Appendix Page | 66

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - CS_2: Cross section of Western Range (part 1), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 67

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - CS_3: Cross section of Western Range (part 2), Current AWT habitat

Appendix Page | 68

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - CS_4: Cross section of Western Range (part 2), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 69

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - CS_5: Cross section of Western Range (part 3), Current AWT habitat

Appendix Page | 70

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - CS_6: Cross section of Western Range (part 3), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 71

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - CS_7: Cross section of Western Range (part 4), Current AWT habitat

Appendix Page | 72

Greater Paraburdoo Subterranean Fauna Survey

Appendix 4 - CS_8: Cross section of Western Range (part 4), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 73

Greater Paraburdoo Subterranean Fauna Survey

APPENDIX 5: Eastern Range 3D Habitat Modelling (AWT) – Long and cross sections

Appendix Page | 74

Greater Paraburdoo Subterranean Fauna Survey

Eastern Range 3D Habitat Modelling (AWT) – Long and cross sections

Appendix 5 - LS_1: Long section of Eastern Range (part 1), Current AWT habitat

Appendix Page | 75

Greater Paraburdoo Subterranean Fauna Survey

Appendix 5 - LS_2: Long section of Eastern Range (part 1), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 76

Greater Paraburdoo Subterranean Fauna Survey

Appendix 5 - LS_3: Long section of Eastern Range (part 2), Current AWT habitat

Appendix Page | 77

Greater Paraburdoo Subterranean Fauna Survey

Appendix 5 - LS_4: Long section of Eastern Range (part 2), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 78

Greater Paraburdoo Subterranean Fauna Survey

Appendix 5 - LS_5: Long section of Eastern Range (part 3), Current AWT habitat

Appendix Page | 79

Greater Paraburdoo Subterranean Fauna Survey

Appendix 5 - LS_6: Long section of Eastern Range (part 3), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 80

Greater Paraburdoo Subterranean Fauna Survey

Appendix 5 - CS_1: Cross section of Eastern Range (part 1), Current AWT habitat

Appendix Page | 81

Greater Paraburdoo Subterranean Fauna Survey

Appendix 5 - CS_2: Cross section of Eastern Range (part 1), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 82

Greater Paraburdoo Subterranean Fauna Survey

Appendix 5 - CS_3: Cross section of Eastern Range (part 2), Current AWT habitat

Appendix Page | 83

Greater Paraburdoo Subterranean Fauna Survey

Appendix 5 - CS_4: Cross section of Eastern Range (part 2), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 84

Greater Paraburdoo Subterranean Fauna Survey

Appendix 5 - CS_5: Cross section of Eastern Range (part 3), Current AWT habitat

Appendix Page | 85

Greater Paraburdoo Subterranean Fauna Survey

Appendix 5 - CS_6: Cross section of Eastern Range (part 3), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 86

Greater Paraburdoo Subterranean Fauna Survey

APPENDIX 6: Paraburdoo 3D Habitat Modelling (AWT/ BWT) – Long and cross sections

Appendix Page | 87

Greater Paraburdoo Subterranean Fauna Survey

Paraburdoo 3D Habitat Modelling (AWT/ BWT) – Long and cross sections

Appendix 6 - LS_1: Long section of Paraburdoo (part 1), Current AWT habitat

Appendix Page | 88

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - LS_2: Long section of Paraburdoo (part 1), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 89

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - LS_3: Long section of Paraburdoo (part 2), Current AWT habitat

Appendix Page | 90

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - LS_4: Long section of Paraburdoo (part 2), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 91

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - LS_5: Long section of Paraburdoo (part 3), Current AWT habitat

Appendix Page | 92

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - LS_6: Long section of Paraburdoo (part 3), Worst Case Scenario AWT habitat remaining after mining

Appendix Page | 93

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - LS_7: Long section of Paraburdoo (part 1), Current BWT habitat

Appendix Page | 94

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - LS_8: Long section of Paraburdoo (part 1), Worst Case Scenario BWT habitat remaining after mining

Appendix Page | 95

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - LS_9: Long section of Paraburdoo (part 2), Current BWT habitat

Appendix Page | 96

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - LS_10: Long section of Paraburdoo (part 2), Worst Case Scenario BWT habitat remaining after mining

Appendix Page | 97

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - LS_11: Long section of Paraburdoo (part 3), Current BWT habitat

Appendix Page | 98

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - LS_12: Long section of Paraburdoo (part 3), Worst Case Scenario BWT habitat remaining after mining

Appendix Page | 99

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - LS_13: Long section of Paraburdoo (part 3), Current BWT habitat

Appendix Page | 100

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - LS_14: Long section of Paraburdoo (part 4), Worst Case Scenario BWT habitat remaining after mining

Appendix Page | 101

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - CS_1: Cross section of Paraburdoo (part 1), Current BWT habitat

Appendix Page | 102

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - CS_2: Cross section of Paraburdoo (part 1), Worst Case Scenario BWT habitat remaining after mining

Appendix Page | 103

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - CS_3: Cross section of Paraburdoo (part 2), Current BWT habitat

Appendix Page | 104

Greater Paraburdoo Subterranean Fauna Survey

Appendix 6 - CS_4: Cross section of Paraburdoo (part 2), Worst Case Scenario BWT habitat remaining after mining

Appendix Page | 105

Greater Paraburdoo Subterranean Fauna Survey

APPENDIX 7a: Bores and drill holes visited and sampled during the survey

APPENDIX 7b: Database search data (WAM & PSS)

APPENDIX 7c: Specimen data – Current survey

APPENDIX 7d: Water physicochemical measurements observed during the current survey

Note: The above appendices 7a-7d have been provided to client in Excel Spreadsheet format. File name: “GPARA_SUBF_REV2_APPENDICES_7a-7d.xlsx”.

Appendix Page | 106