Metapopulation genetic management of the Macquarie perch

This repository contains data and R script script used in the manuscript

A shift to metapopulation genetic management for persistence of a species threatened by fragmentation: the case of an endangered Australian freshwater fish, by Pavlova A, Tonkin Z, Pearce L, Robledo-Ruiz D, Lintermans M, Ingram B, Lyon J, Beitzel M, Broadhurst B, Rourke ML, Sturgiss F, Lake E, Castrejón-Figueroa, Stocks JR, and Sunnucks P. Molecular Ecology MEC-24-1090.R1, accepted pending minor revision.

The content of this data repository:

MaccaGM_SNPs.R –– R script for analyses of DArT SNP genotypic data and plotting the results of analyses of JeDi pipeline

Input files for MaccaGM_SNPs.R:

Report_DMacq23-8576_14_moreOrders_SNP_mapping_1.csv –– DArT genotypes in original format

Covariate_MaccaGM_recaller25inds_then_max50per_pop.csv –– covariate file for individuals

Unbiased genetic diversity: outputs of JeDi pipeline

Estimates of unbiased heterozygosity (piawka_uHe) resulting from JeDi pipeline run for all sites, biallelic sites and tri- and tetra-allelic sites. Analyses were run with 9 different settings (A-H). Estimates of individual unbiased heterozygosity from JeDi pipeline are appended to the individual covariate file with results of analyses of SNP dataset (SNP heterozygosity, PHt, and membership in STRUCTURE clusters (X1of15 to X15of15).

ind.gen.div_nonaPHt_noNAuHe.with.PHt.residuals.noB_qual30_minDP10_SETA.csv –– Output of SetA: Piawka run with Settings: no (-B), QUAL30, MinDP10, No doubleton filter

ind.gen.div_nonaPHt_noNAuHe.with.PHt.residuals.B_qual30_minDP10_SETB.csv –– Output of SetB: Piawka run with Settings: (-B), QUAL30, MinDP10, No doubleton filter

ind.gen.div_nonaPHt_noNAuHe.with.PHt.residuals.B_qual40_minDP10_SETC.csv –– Output of SetC: Piawka run with Settings: (-B), QUAL40, MinDP10, No doubleton filter

ind.gen.div_nonaPHt_noNAuHe.with.PHt.residuals.B_qual30_minDP12_setD.csv –– Output of SetD: Piawka run with Settings: (-B), QUAL30, MinDP12, No doubleton filter

ind.gen.div_nonaPHt_noNAuHe.with.PHt.residuals.B_qual30_minDP10_doubleton_SETE.csv–– Output of SetE: Piawka run with Settings: (-B), QUAL30, MinDP10, doubleton filter

ind.gen.div_nonaPHt_noNAuHe.with.PHt.residuals.B_qual50_minDP15_SETF.csv–– Output of SetF: Piawka run with Settings : (-B), QUAL50, MinDP15, No doubleton filter

ind.gen.div_nonaPHt_noNAuHe.with.PHt.residuals.noB_qual30_minDP10_doub_SETG.csv–– Output of SetG: Piawka run with Settings: No (-B), QUAL30, MinDP10, doubleton filter

ind.gen.div_nonaPHt_noNAuHe.with.PHt.residuals.noB_qual30_minDP10_doub_sing_SETH.csv–– Output of SetH: Piawka run with Settings: No (-B), QUAL30, MinDP10, doubleton and singleton filter

ind.gen.div_nonaPHt_noNAuHe.with.PHt.residuals.noB_qual30_minDP15_doub_sing_SETI.csv–– Output of SetI: Piawka run with Settings: No (-B), QUAL30, MinDP15, doubleton and singleton filter

Estimates of population nucleotide diversity (pi): outputs of JeDi pipeline run under three settings:

genomic_pi_table_SETA.tsv—pi from JeDi pipeline run under Settings A: no (-B), QUAL30, MinDP10, No doubleton filter

genomic_pi_table SETG.tsv—pi from JeDi pipeline run under Settings G: No (-B), QUAL30, MinDP10, doubleton filter

genomic_pi_table - SETH.tsv—pi from JeDi pipeline run under Settings H: No (-B), QUAL30, MinDP10, doubleton and singleton filter

Population summaries of the genetic estimates from SNP dataset (Ho- observed heterozygosity, He-expected heterozygosity, private.alleles.dartr.one2rest and PA- number of private alleles per population including and excluding three with sample size <13, Mean_Allelic_Richness- allelic richness, Ne- effective population size estimated by LDNe) and sequencing dataset (pi_SetI and pi_SetH- nucleotide diversity estimated with settings I and H, mean_uHe_SetI and mean_uHe_SetH- mean unbiased heterozygosity estimated with setting I and H, respectively):

Pop.genetic.div.estimates.csv––collated population outputs of the script MaccaGM_SNPs.R and JeDi pipeline

Metapop2 files:

Metapop_res.csv –– results of the Metapop2 analyses to create a new population (Fig. 4 of the manuscript)

Metapop2_24sims.zip––input and output files for 24 realistic simulations (Table 3 of the main manuscript).

SRA_PRJNA1242510_accession_numbers.txt––GenBank SRA accession numbers for DArT sequences

Appendices for the manuscript:

Pavlova_etal_AppendixA_Genetic_augmentations.xlsx–– Genetic augmentation of Macquarie perch via translocations and stocking conducted from 2010 onwards.

Pavlova_et_al_AppendixB_JeDi_pipeline_revised.png–– Schematic representation of revised JeDi pipeline for estimating unbiased individual heterozygosity and population nucleotide diversity from reduced-representation sequencing data, in presence of reference genome.

Pavlova_et_al_MEC_Supplemental_Information_revised.docx––Supplementary Information for the revised version of the main manuscript: S1. Additional tables and figures, S2. Details and optimization of the JeDi pipeline.

We acknowledge the First Nations throughout Australia, recognise their continuing connection to land, waters and culture, and pay our respects to their Elders past, present and emerging. This research was conducted on Ngarigo, Ngambri and Ngunnawal, Taungurung, Wirajuri and Wurundjeri Woi-wurrung Countries.