DOI: 10.5281/zenodo.8132980
title: Exon capture museomics deciphers the continent-wide nine-banded armadillo species complex and identifies a new species endemic to the Guiana Shield

description: <strong>Exon capture museomics deciphers the continent-wide nine-banded armadillo species complex and identifies a new species endemic to the Guiana Shield</strong> Mathilde Barthe, Loïs Rancilhac, Maria C. Arteaga, Anderson Feijó, Marie-Ka Tilak, Fabienne Justy, W. J. Loughry, Colleen M. McDonough, Benoit de Thoisy, François Catzeflis<sup>✝</sup>, Guillaume Billet, Lionel Hautier, Benoit Nabholz, and Frédéric Delsuc <strong>Description of available files:</strong> <strong>01_Figures_&_tables_of_the_main_text.zip</strong> - <strong>Figure 1: </strong>Phylogenetic relationships reconstructed by maximum likelihood and maps representing the distribution of individuals according to their lineage. - <strong>Figure 2: </strong>Assignment of individuals to lineages according to phylogenetic analyses, admixture analysis and phylogenetic delimitation. - <strong>Figure 3: </strong>Principal Component Analysis of genetic variance. - <strong>Figure 4: </strong>Distribution map and genetic composition of individuals of the four recognized species. <strong>02_Supplementary_tables_&_figures.zip</strong> - <strong>Figure S1: </strong>Mitochondrial genome depth of coverage - <strong>Figure S2: </strong>Calculation of mitochondrial lineage support for detecting contamination. - <strong>Figure S3: </strong>Mitochondrial lineage support for each individual - <strong>Figure S4: </strong>a) Inbreeding coefficient and b) heterozygosity estimate for individuals according to cleaning steps. - <strong>Figure S5: </strong>Percentage of missing data by captured locus - <strong>Figure S6: </strong> Phylogenetic relationships of the 62 <em>Dasypus</em> individuals obtained using Astral on the 832 ML gene trees from the captured nuclear loci reconstructed with IQ-Tree and ModelFinder. - <strong>Figure S7: </strong>Results of analyses to detect introgression. - <strong>Figure S8: </strong>Cross validation errors according to the number of clusters (K) investigated. - <strong>Figure S9: </strong>Detailed analysis of the sub-structure within the newly recognized <em>D. novemcinctus </em>(Southern lineage)<em>.</em> - <strong>Figure S10: </strong>Species delimitation estimated by bPTP-h. - <strong>Figure S11: </strong>Species delimitation estimated using GMYC. - <strong>Figure S12: </strong>Heatmaps of pairwise genetic indexes between lineages. - <strong>Figure S13: </strong>Effect of filters on admixture results. - <strong>Figure S14: </strong>Updated map from Arteaga et al. (2020). - <strong>Figure S15: </strong>Maximum likelihood phylogenetic tree of 212 pb of the 16s ribosomal RNA of five individuals analyzed in Abba et al., (2018) and three from this study. - <strong>Table S1: </strong>List of biological samples with detailed information - <strong>Table S2a: </strong>Quality statistics by locus after filtering steps - <strong>Table S2b: </strong>Quality statistics by individuals after filtering steps - <strong>Table S3: </strong>Comparison of the lineage of the nineteen individuals in common with Arteaga et al. (2020) and our study. - <strong>Table S4: </strong>Adult cranial measurements (in millimeters) of the four <em>Dasypus</em> species recognized in this study following Feijó & Cordeiro-Estrela (2016). - <strong>Table S5:</strong> Adult external measurements (in millimeters) of the four <em>Dasypus</em> species recognized in this study. <strong>03_Mitogenomes.zip</strong> ● <strong>Mitogenome_reference_Dasypus_novemcinctus.fasta: </strong>Mitogenome reference used to mapped reads and extract mitochondrial DNA. ● <strong>Concatenated_mitochondrial_genes.fasta: </strong>Concatenated nucleotide sequences of 15 mitochondrial genes (13 protein-coding + 2 rRNAs). Sites with more than 50% missing data were excluded resulting in a total of 13,924 sites. ● <strong>Concatenated_mitochondrial_genes_partition.txt: </strong>Partition file of concatenated sequences of the 15 mitochondrial genes (13 protein-coding + 2 rRNAs). ● <strong>Concatenated_mitochondrial_genes_TESTNEW.treefile : </strong>Maximum likelihood phylogenetic tree inferred from the concatenated sequences of the 15 mitochondrial genes using IQ-TREE under a partitioned model applying ModelFinder on each partition. ● <strong>Depth_coverage_mitogenomes.csv: </strong>Table of mean depth of coverage and proportion of missing data (Ns) of the 72 reconstructed mitochondrial genomes sequenced for this study. <strong>04_Reanalyses.zip</strong> ● <strong>Dloop_alignment.fasta: </strong>Alignment of the D-loop sequences obtained in this study with those from Arteaga et al. (2020). ● <strong>Dloop_alignment.treefile: </strong>Maximum likelihood phylogenetic tree inferred from the D-loop alignment using IQ-TREE (GTR+G model). ● <strong>Abba_shotgun.fasta: </strong>Alignment of the 16S rRNA of individuals from this study and those from Abba et al. (2018). ● <strong>Abba_shotgun.fasta.treefile: </strong>Maximum likelihood phylogenetic tree inferred from the 16S rRNA alignment using IQ-TREE (GTR+G model). <strong>05_Contamination_exploration.zip</strong> ● <strong>Mitochondrial_diagnostic_positions.csv: </strong>Table of the 350 diagnostic mitochondrial positions used to estimate proportion of reads supporting each lineage. Position number refers to the Complete_mitogenome_alignment.fasta file. ● <strong>Read_support_to_diagnostic_positions.csv: </strong>For each individual, this table reports the Diagnostic Rate (proportion of diagnostic positions per lineage supported by at least 3 reads), the Read Proportion (mean read proportion supporting diagnostic positions per lineage), Index (proportion of synapomorphies per lineage normalized by average frequency of reads supporting these synapomorphies) and the type of tissue (museum or fresh tissue). ● <strong>Contamination_exploration.R:</strong> R script used to plot read support to lineages and the effect of tissue type (fresh or museum). <strong>06_Nuclear_dataset.zip</strong> ● <strong>TATU_1000exons4baits.fasta: </strong>Reference sequences of 1,000 exons and flanking regions used to define the probes for exon capture extracted from the <em>Dasypus novemcinctus</em> genome. ● <strong>Dasypus_capture_Final_Baits_Set.fas: </strong>Sequences of the 16,146 probes used to capture the 997 nuclear loci (exons and flanking regions). ● <strong>Diploid_837_nuclear_loci.fasta: </strong>Diploid sequences of the 837 nuclear loci for the 62 individuals in PopPhyl format (Locus|lineage|individual|Allele). ● <strong>Mean_coverage_by_individuals.csv: </strong>Table of mean depth of coverage, horizontal coverage, and number of loci per individual after filtering. ● <strong>Mean_coverage_by_loci.csv: </strong>Table of mean depth of coverage, horizontal coverage and number of loci per loci after filtering. <strong>07_Disentangle_genotyping_errors.zip</strong> ● <strong>Table_of_heterozygosity_and_inbreeging_coefficient.csv: </strong>Table of heterozygosity (He) and inbreeding coefficient (F) estimated for each cleaning steps: initial data, after correction of heterozygous positions (must be supported by a proportion of reads between 0.3 and 0.7), and after exclusion of 159 potentially paralogous loci. ● <strong>Plot_effect_of_cleaning_on_He&F.R: </strong>R script used to plot the effect of cleaning steps on heterozygosity (He) and inbreeding coefficient (F). <strong>08_Distribution_maps.zip</strong> ● <strong>Coordinates_according_mito_nuclear_lineages.csv: </strong>Table of GPS coordinates of individuals according to their mitochondrial and nuclear lineages. ● <strong>Plot_mito_nuclear_distribution.R: </strong>R script used to plot individuals on the Neotropical map according to their mitochondrial and nuclear lineages in Figure 1. ● <strong>Mitochondrial_distribution.pdf: </strong>Geographical distribution of the 75 individuals according to their mitochondrial lineage. ● <strong>Nuclear_distribution.pdf: </strong>Geographical distribution of the 58 individuals according to their nuclear lineage. <strong>09_Phylogenetic_inference.zip</strong> ● Phylogram_Tree ○ <strong>Concatenated_nuclear_loci.fasta: </strong>Concatenated sequences of the 837 nuclear loci representing a total of 506,355 sites. ○ <strong>Concatenated_nuclear_loci_partition.txt: </strong>Partition file for the 837 nuclear loci concatenation. ○ <strong>Concatenated_nuclear_loci_TESTNEW.treefile: </strong>Maximum likelihood phylogenetic tree inferred from the 837 nuclear loci concatenation using IQ-TREE under a partitioned model applying ModelFinder on each partition. ● Ultrametric_Tree ○ <strong>Ultrametric_tree_concatenated_nuclear_loci.treefile: </strong>Ultrametric tree inferred from the 837 nuclear loci concatenation (Concatenated_nuclear_loci.fasta in Phylogram_Tree folder) using a partitioned model applying ModelFinder on each partition (Concatenated_nuclear_loci_partition.txt in Phylogram_Tree folder). The ML phylogram (Concatenated_nuclear_loci_TESTNEW.treefile in Phylogram_Tree folder) was used as a guide tree. The root was dated at 6 Mya. ● Gene_Tree ○ <strong>Concatenate_gene_tree.treefile:</strong> File containing all gene trees reconstructed using IQ-TREE applying ModelFinder to each gene. ○ <strong>Astral_consensus_tree.txt: </strong>Summary species tree reconstructed with Astral using Concatenate_gene_tree_TESTNEW.treefile ● <strong>Introgression analyses:</strong> ○ <strong>Concordance_factors_Dasypus.csv </strong> ○ <strong>Topology_Weighting_Dasypus_plots.R</strong> ○ <strong>SnaQ_results_hmax0.out </strong> ○ <strong>SnaQ_results_hmax1.out </strong> ○ <strong>SnaQ_results_hmax2.out </strong> ○ <strong>SnaQ_results_hmax3.out </strong> ○ <strong>twisst_guianensis_spmap.txt </strong> ○ <strong>twisst_guianensis_Weights</strong> ○ <strong>twisst_mexico_spmap.txt</strong> ○ <strong>twisst_mexico_Weights</strong> <strong>10_Species_delimitation.zip</strong> ● BPP ○ <strong>input_for_bpp.phy: </strong>Sequence alignments of the 837 nuclear loci in phylip format. ○ <strong>lineage_for_BPP: </strong>Correspondence file between individuals and lineages. ○ <strong>r1 and r2: </strong>folders containing config files (bpp.ctl) and outputs of the BPP analysis. ● bPTP ○ <strong>PTPh_Support_Partition.txt: </strong>Details of the most supported species partition. ○ <strong>PTPh_tree_partition.png:</strong> Tree illustrating the most supported species partition. ● GMYC ○ <strong>Script_GMYC.R: </strong>R script used to run the GMYC delimitation method on the ultrametric tree (11_Phylogenetic_inference/Ultrametric_Tree/Ultrametric_tree_concatenated_nuclear_loci.treefile). ○ <strong>Figure_GMYC.png:</strong> Figure illustrating the results of the GMYC species delimitation analysis. <strong>11_Population_genetic_analyses.zip</strong> ● PCA ○ <strong>Input_for_PCA.fasta: </strong>Diploid sequences of the 57 individuals (DNO-MC21 and DPI-L29 excluded) in PopPhyl format (Locus|species|individual|allele). ○ <strong>PCA_Output: </strong>Output of the PopPhyl2PCA analysis using the Input_for_PCA.fasta file. ○ <strong>Script_to_plot_PCA.R:</strong> R script used to plot PCA according to the mitochondrial lineage and nuclear composition (Admixture results). ● ADMIXTURE ○ <strong>lineage_for_Admixture.list: </strong>Correspondence between individuals and lineages file. ○ <strong>Input_Admixture.*:</strong> 19,872 SNPs from nuclear data across the <em>Dasypus</em> complex. ○ <strong>Output_Admixture.k.*: </strong>Output from the Admixture analysis according to K values (from 1 to 7). ○ <strong>Output_Admixture.cv.error: </strong>Summary of the error value according to K. ○ <strong>Plot_Admixture.R:</strong> R script used to plot Admixture results reordered by phylogeny. ○ <strong>Plot_map_distribution_admixture.R:</strong> R script used to plot Admixture results on the Neotropical map. ● Stats_Da_Dxy_GDI ○ <strong>Pairwise_genetic_statistics.csv: </strong>Summary statistics computed using ABCstat_global.txt from the DILSmcsnp program for all pairwise combinations of individuals from the different lineages. ○ <strong>Pairwise_GDI.csv:</strong> Genetic Differentiation Index estimates for all pairwise combinations of individuals from the different lineages. ○ <strong>Plot_genetic_statistics.R:</strong> R script used to plot mean genetic statistics between lineages. ● Sublineage_structure ○ <strong>ADMIXTURE</strong> ■ <strong>Plot_map_distribution_sublineage_admixture.R: </strong>R script used to plot Admixture results on the Neotropical map. ○ <strong>PCA</strong> ■ <strong>Input_for_PCA_southern_lineage.fasta: </strong>Diploid sequences of the 24 individuals of the Southern lineage in PopPhyl format (Locus|species|individual|allele). ■ <strong>PCA_Output_southern_lineage: </strong>Output of the PopPhyl2PCA analysis using the Input_for_PCA_southern_lineage.fasta file. ■ <strong>Script_to_plot_PCA_sublineage.R:</strong> R script to plot PCA according to the mitochondrial lineage and nuclear composition (Admixture results) focussing on individuals from the Southern lineage. <strong>12_Morpho_molecular_distribution.zip</strong> ● Coordinates_according_morphogroup_lineages.csv: GPS coordinates of individuals used in Hautier et al. (2017) according to their morphogroup. ● Plot_map_distribution_morpho_admixture.R: R script used to plot Admixture results and the individuals from Hautier et al. (2017) on the Neotropical map in Figure 5. ● Skull_lateral_*.png: Illustration of the lateral view of the skull of four individuals representing each species. ● Skull_sinuses_*.png: Illustration of the skull and paranasal sinuses of four individuals representing each species.