Genetic, but Not Behavioral, Evidence Supports the Distinctiveness of the Mealy Amazon Parrot in the Brazilian Atlantic Forest

: The presence of unidentiﬁed cryptic species within a species complex can obscure demographic trends of vulnerable species, impacting potential species conservation and management decisions. Previous work identiﬁed a taxonomic split between Central and South American populations of the mealy amazon ( Amazona farinosa ) that subsequently resulted in the elevation of these two populations to full species status ( Amazona guatemalae and A. farinosa, respectively). In that study, however, a third, geographically disjunct population from the Brazilian Atlantic Forest was insufﬁ-ciently sampled, limiting the ability of researchers to fully evaluate its genetic distinctiveness. Given that signiﬁcant levels of biodiversity and endemism are found in this region, we aimed to use genetic and behavioral data to determine if the Atlantic Forest population of A. f. farinosa represents a third cryptic species within the complex. We sequenced 6 genes (4 mitochondrial and 2 nuclear introns) from the Atlantic Forest population of A. f. farinosa to measure the genetic relationships between this population and all other recognized species and subspecies of the mealy amazon. In addition, we use spectrographic cross-correlation and an analysis of 29 acoustic parameters to determine whether the taxa diverge in their learned contact call structure and if the degree of vocal differentiation correlates to genetic structure. We found that the Atlantic Forest population of A. f. farinosa was genetically distinct from that of the greater Amazon basin, but the degree of differentiation was less than that separating the Central and South American taxa. Acoustic analysis revealed substantial variation in contact call structure within each clade. This variation created substantial overlap in acoustic space between the clades. In all, the degree of call divergence between clades did not correspond to the degree of genetic divergence between the same clades. The results suggest that in taxa with substantial geographic variation in learned calls, such as the mealy amazon, vocalizations may not be a useful tool in the identiﬁcation of cryptic species that are lifelong vocal learners. While these results do not support the elevation of the Brazilian Atlantic Forest population of the mealy amazon to full species status, given current trends of habitat loss in the Atlantic Forest as well as the imperiled status of large parrot species globally, we argue that this population nonetheless warrants special conservation and management consideration as a pool of unique genetic diversity within the southern mealy amazon species.


Introduction
A primary challenge for conservationists and wildlife managers is understanding the demographic trends of vulnerable species. This task is complicated if the basic taxonomic relationships between species are not well understood. This issue commonly arises within species complexes that contain unidentified cryptic species: morphologically indistinguishable but genetically distinct species that are mistakenly classified under a single The Atlantic Forest is the second largest rainforest system in South America, stretching across a broad strip of Brazil's southeastern Atlantic coastline, and is a hotspot of endemic biodiversity in the Neotropics [40]. Currently, the Atlantic Forest is separated from the greater Amazon basin by the dry Caatinga region to the northeast, the upland, wooded savanna of the Cerrado region to the north/northwest, and from the Andean forests to the west by the lowland, arid Chaco region of central South America [41]. Historically, this forest has undergone several periods of connection and separation from these other major forest regions, leading to the isolation of some species over both short (e.g., 10,000-20,000 years ago) and long (e.g., >3 million years ago) time scales [41,42]. This dynamic evolutionary history has created a region rich with endemism and cryptic species complexes [43][44][45][46][47] that is threatened today by high levels of habitat loss and fragmentation [48]. For this reason, the Atlantic Forest is an area of upmost concern for biodiversity conservation [48,49].
We aimed to address two main questions. First, is the Atlantic Forest population of A. f. farinosa a cryptic species within the mealy amazon species complex? Given the possible timeframes when the Atlantic Forest population of A. f. farinosa could have been separated from the larger Amazon basin population of A. f. farinosa, we hypothesize that it is as genetically dissimilar from the other southern mealy amazon subspecies as the northern species is from the southern species. Second, can we use the vocalizations of an open-ended vocal learner to identify potential cryptic species? If so, we hypothesize that the degree of divergence in contact call structure between and within the mealy amazon species, subspecies, and the Atlantic Forest population of A. f. farinosa will correlate with the genetic distance between these groups. To measure genetic differentiation, we expanded on the original dataset from [36] by sequencing six gene regions from newly sampled individuals of the Atlantic Forest population of A. f. farinosa to create both a phylogeny and a haplotype network that includes this group. The similarity of contact calls within and between groups, as well as their relationship to the genetic distance between groups, was then assessed using 29 acoustic parameter measurements and spectrographic cross-correlation values. We predicted that our genetic and behavioral (i.e., contact call) measures will show related degrees of differentiation, with the Atlantic Forest population of A. f. farinosa falling out as divergent from the Amazon basin population of A. f. farinosa and the southern mealy amazon species as a whole. These results will inform the study and preservation of this geographically distinct population of the mealy amazon.

Study Species and Data Summary
The mealy amazon species complex includes two recognized species and five subspecies. The northern mealy amazon A. guatemalae (including the virenticeps and guatemalae subspecies) ranges from southern Mexico to central Panama, while the southern mealy amazon A. farinosa is distributed in three geographically distinct locations: (1) the inornata subspecies extends from southern Panama through western Ecuador and eastern Colombia, (2) the chapmani and farinosa subspecies range from eastern Ecuador, eastern Peru, and northern Bolivia through the greater Amazon Basin, and (3) in an isolated population of the farinosa subspecies along Brazil's Atlantic coast [38] (Figure 1). Seven samples from the Atlantic Forest population of A. f. farinosa were unique to this study as all other clades within the species complex were sampled adequately by [36]; see Table S1 for details on the origins of the genetic samples.
Recordings of contact calls from individuals of all subspecies and the Atlantic Forest population of A. f. farinosa were obtained from the Macaulay library (Cornell Laboratory of Ornithology, Ithaca, NY, U.S.) and Xeno-canto (Xeno-canto Foundation for Nature Sounds, Netherlands). We used 148 high-quality sound files from 1954 to 2019 in our analysis; sound file metadata can be found in Table S2. Spectrograms of representative calls from each clade are shown at their corresponding recording location. Spectrograms were created using the 'spec_param' and 'specreator' functions in warbler to optimize the visualization parameters (overlap = 90, window length = 475, color levels = (−50, 0, 5)). Genetic samples from the [36] study are indicated by grey-outlined triangles on the map. The location of the new genetic samples used in our study is indicated by the purple-outlined triangle.

Phylogenetic Analysis and Haplotype Network
Blood samples were collected, according to SISBIO permit #71978-1 protocols following Brazilian environmental legislation, on Whatman FTA elute micro cards from seven individuals being held in captivity at the ASM Cambaquara rescue center on the island of Ilhabela, São Paulo, Brazil, after seizure as illegal pets or after being found injured in the wild. These individuals represent a single Atlantic Forest A. f. farinosa population (for complete sample metadata see Table S1). For each sample, genomic DNA was isolated and then PCR amplification was performed for the four mitochondrial and two nuclear intron gene fragments previously used in [36]: 12S rDNA (12S), 16S rDNA (16S), Cytochrome oxidase subunit I (COI), cytochrome b (CytB), tropomyosin alpha-subunit intron 5 (TROP), and transforming growth factor β-2 intron I (TGFB2). An additional mitochondrial gene fragment for NADH dehydrogenase 2 (ND2) was amplified for both the [36] samples and the new A. f. farinosa Atlantic Forest samples. A complete description of the number of individuals sampled from each population, the genes sequences per individual, GenBank accession numbers, and primers can be found in Tables S1, S3 and S4. All Spectrograms of representative calls from each clade are shown at their corresponding recording location. Spectrograms were created using the 'spec_param' and 'specreator' functions in warbler to optimize the visualization parameters (overlap = 90, window length = 475, color levels = (−50, 0, 5)). Genetic samples from the [36] study are indicated by grey-outlined triangles on the map. The location of the new genetic samples used in our study is indicated by the purple-outlined triangle.

Phylogenetic Analysis and Haplotype Network
Blood samples were collected, according to SISBIO permit #71978-1 protocols following Brazilian environmental legislation, on Whatman FTA elute micro cards from seven individuals being held in captivity at the ASM Cambaquara rescue center on the island of Ilhabela, São Paulo, Brazil, after seizure as illegal pets or after being found injured in the wild. These individuals represent a single Atlantic Forest A. f. farinosa population (for complete sample metadata see Table S1). For each sample, genomic DNA was isolated and then PCR amplification was performed for the four mitochondrial and two nuclear intron gene fragments previously used in [36]: 12S rDNA (12S), 16S rDNA (16S), Cytochrome oxidase subunit I (COI), cytochrome b (CytB), tropomyosin alpha-subunit intron 5 (TROP), and transforming growth factor β-2 intron I (TGFB2). An additional mitochondrial gene fragment for NADH dehydrogenase 2 (ND2) was amplified for both the [36] samples and the new A. f. farinosa Atlantic Forest samples. A complete description of the number of individuals sampled from each population, the genes sequences per individual, GenBank accession numbers, and primers can be found in Tables S1, S3 and S4. All PCR products were directly sequenced at the University of Texas at El Paso's Border Biomedical Research Center genomic analysis core facility with the same primers used for PCR. Bidirectional reads were screened for quality, assembled to an Amazona reference sequence (Table S3), and trimmed to the appropriate length to enable gene alignment using Sequencher v5.4.6 (Gene Codes Corporation, Ann Arbor, MI, U.S.). The individual gene alignments were concatenated into a supermatrix using the Bio.Nexus module from the Biopython project v1.76 [50]. Three outgroups used in the [36] analysis-A. kawalli, A. auropalliata, and A. amazonica-were included to root our phylogenetic tree and to provide a measure of relative genetic distance with which to compare recognized Amazona species to the clades within the mealy amazon species complex.
Two median joining haplotype networks [54] were constructed and visualized using PopART v1.7 [55]. The first network includes all of the available mitochondrial genes (12S, 16S, CytB, ND2) because we were able to obtain sequence data for both recognized species, at least four of the five subspecies (all minus virenticeps), and Atlantic Forest population of A. f. farinosa. The second network was created using only CytB data, as this most closely resembles the data set used by [36].

Genetic Differentiation
To quantify inter-and intraspecific genetic variation we calculated the between group mean genetic distance (i.e., the average number of nucleotide substitutions between clades) using the concatenated gene alignment and a p-distance model in MEGA X [56,57].

Acoustic Analysis
High-quality (i.e., low background noise, no signal overlap) contact calls were manually selected from the sound files using Raven Pro v1.5.0 (Cornell Laboratory of Ornithology, Ithaca, NY, U.S.). The accurate documentation of caller ID was lacking in most sound files with multiple calling individuals; therefore, only one call from one individual was randomly selected per call file. In call files where the researcher was able to identify more than one individual from counter-calling (i.e., the repeated, predictable response of a second individual to the calls of a first), one call was randomly selected from up to two individuals. In all, only one randomly selected high-quality (i.e., determined by each database to be "A" quality out of an A-F rating scale) call per individual was used in all downstream analyses. After call selection, the warbleR package v1.1.19 [58] in R v3.6.2 (The R Foundation) was used to perform additional quality control and data processing. Individual spectrograms were created with the 'specreator' function using a window length of 510 samples, a 90% overlap between windows, and a frequency limit of 0 kHz to 9kHz to reassess call quality. Poorly selected calls (i.e., calls where the temporal coordinates being read by a function did not match the actual start and end times of the calls) were reselected using the 'seltailor' function.
Differences in contact call structure between all populations were assessed both quantitatively and qualitatively. First, 26 acoustic parameter measurements-including various frequency and duration characteristics, entropy, skew, and kurtosis-were calculated for each call using the warbleR 'specan' function with a 90% overlap between windows and a bandpass filter from 0 kHz to 7 kHz (Table S5). The extracted data were then used in a principal component analysis performed with the 'prcomp' function from the base stats R package. The first two principal components were plotted to visualize the clustering patterns of calls by population. Next, the 'xcorr' function of the warbleR package was used to create a spectrographic cross-correlation matrix of all the calls. This matrix was transformed into distance measurements so that a multidimensional scaling analysis using the stats base package in R could be used to visualize the clustering patterns of calls by clade.
A Mantel test was performed using the vegan package v2.5-6 in R to determine whether the degree of genetic differentiation (measured as the average number of sequence differences) between clades correlates with vocal differentiation (measured as the average cross-correlation dissimilarity values).

Phylogenetic Tree and Haplotype Network
A total of 3444 bp from 6 genes (393 bp of 12S, 525 bp of 16S, 868 bp of CytB, 510 bp of ND2, 625 bp of TGFB2, 523 bp of TROP) were included in the maximum-likelihood phylogenetic analysis, though not all genes were able to be sequenced for every individual; see Table S3 for a complete list genes per individual. Notably, the forward reads for COI failed to sequence for all of the Atlantic Forest A. f. farinosa samples, so this gene was removed from the analysis. With the addition of the ND2 mtDNA sequence and the exclusion of the COI mtDNA sequence, the topology of our maximum likelihood majority rule consensus tree is nearly identical to the phylogeny reported by the original [36] study for taxa included in both studies ( Figure 2). The Atlantic Forest population of A. f. farinosa is recovered as a monophyletic group sister to the assemblage that includes the three recognized subspecies of A. farinosa (A. f. inornata, A. f. chapmani, and A. f. farinosa). This group does not, however, show the same level of genetic differentiation as the two currently recognized species of mealy amazon (A. guatemalae and A. farinosa; see insert of Figure 2). Diversity 2020, 12, x FOR PEER REVIEW 6 of 15 using the stats base package in R could be used to visualize the clustering patterns of calls by clade. A Mantel test was performed using the vegan package v2.5-6 in R to determine whether the degree of genetic differentiation (measured as the average number of sequence differences) between clades correlates with vocal differentiation (measured as the average cross-correlation dissimilarity values).

Phylogenetic Tree and Haplotype Network
A total of 3,444 bp from 6 genes (393bp of 12S, 525bp of 16S, 868bp of CytB, 510bp of ND2, 625bp of TGFB2, 523bp of TROP) were included in the maximum-likelihood phylogenetic analysis, though not all genes were able to be sequenced for every individual; see Table S3 for a complete list genes per individual. Notably, the forward reads for COI failed to sequence for all of the Atlantic Forest A. f. farinosa samples, so this gene was removed from the analysis. With the addition of the ND2 mtDNA sequence and the exclusion of the COI mtDNA sequence, the topology of our maximum likelihood majority rule consensus tree is nearly identical to the phylogeny reported by the original [36] study for taxa included in both studies ( Figure 2). The Atlantic Forest population of A. f. farinosa is recovered as a monophyletic group sister to the assemblage that includes the three recog-   The median joining haplotype network created using the CytB data returned a total of 16 haplotypes, with 28 sequence differences between the northern and southern species and 5 sequence differences between the Atlantic Forest populations of A. f. farinosa and their closest Amazon basin A. f. farinosa relative (Figure 3). In addition, the Atlantic Forest populations form a distinct cluster within the greater southern mealy amazon cluster. The haplotype network created with the reduced dataset of the four mitochondrial genes shows a similar pattern. Thirteen haplotypes were recovered, with 67 sequence differences separating the northern and southern mealy amazon species haplotypes and 10 sequence differences separating the Atlantic Forest populations of A. f. farinosa from their closest relative from the Amazon basin population of A. f. farinosa (Figure 3). This network lacks a representative of the A. g. virenticeps northern mealy amazon subspecies (Table S3) due to an insufficient amount of sequence data for the analysis software. However, we assume any virenticeps individuals would be closely related to the other northern mealy amazon subspecies A. g. guatamalae based on the topology of the CytB network. Together, these results consistently show that the Atlantic Forest population is a distinct group when compared to the other A. farinosa subspecies, albeit one that lacks the genetic distinctiveness of previously recognized full species.
Diversity 2020, 12, x FOR PEER REVIEW 7 of 15 between the Atlantic Forest clade and the other southern mealy amazon clades. The relative branch lengths of the phylogram represent sequence divergence.
The median joining haplotype network created using the CytB data returned a total of 16 haplotypes, with 28 sequence differences between the northern and southern species and 5 sequence differences between the Atlantic Forest populations of A. f. farinosa and their closest Amazon basin A. f. farinosa relative (Figure 3). In addition, the Atlantic Forest populations form a distinct cluster within the greater southern mealy amazon cluster. The haplotype network created with the reduced dataset of the four mitochondrial genes shows a similar pattern. Thirteen haplotypes were recovered, with 67 sequence differences separating the northern and southern mealy amazon species haplotypes and 10 sequence differences separating the Atlantic Forest populations of A. f. farinosa from their closest relative from the Amazon basin population of A. f. farinosa (Figure 3). This network lacks a representative of the A. g. virenticeps northern mealy amazon subspecies (Table S3) due to an insufficient amount of sequence data for the analysis software. However, we assume any virenticeps individuals would be closely related to the other northern mealy amazon subspecies A. g. guatamalae based on the topology of the CytB network. Together, these results consistently show that the Atlantic Forest population is a distinct group when compared to the other A. farinosa subspecies, albeit one that lacks the genetic distinctiveness of previously recognized full species.

Genetic Differentiation
Between-group mean genetic distances based on data from all genes included in this study indicate 2.9% average sequence divergence between northern and southern mealy

Genetic Differentiation
Between-group mean genetic distances based on data from all genes included in this study indicate 2.9% average sequence divergence between northern and southern mealy amazon subspecies, compared to 0.4% average sequence divergence between southern mealy amazon subspecies (including the Atlantic Forest population of A. f. farinosa) and 0.3% average sequence divergence between northern mealy amazon subspecies.

Call Similarity between Subspecies
A total of 150 calls (110 from the Macaulay Library and 40 from Xeno-canto) were used in our analysis of call similarity between the 5 recognized mealy amazon subspecies and the Atlantic Forest population of A. f. farinosa. Sampling locations and spectrograms of representative calls from each subspecies are shown in Figure 1. An initial visual inspection of the call spectrograms suggested there is some variation in the structure of contact calls between subspecies (Figure 1), however, plots of the first two components from a principal component analysis of acoustic parameter measurements and a multidimensional scaling analysis of the spectrographic cross-correlation matrix did not reveal distinct clusters of calls among the subspecies. Instead, we saw significant overlap of each clades' calls in acoustic space (Figure 4) due to a large amount of variation in calls within each subspecies. This pattern suggests that (a) substantial structural differences exist across each subspecies' range (e.g., the three example spectrograms of the A. f. farinosa subspecies of southern mealy amazon in Figure 1), and (b) variation within a subspecies is as great as variation among subspecies (Figure 4). amazon subspecies, compared to 0.4% average sequence divergence between southern mealy amazon subspecies (including the Atlantic Forest population of A. f. farinosa) and 0.3% average sequence divergence between northern mealy amazon subspecies.

Call Similarity Between Subspecies
A total of 150 calls (110 from the Macaulay Library and 40 from Xeno-canto) were used in our analysis of call similarity between the 5 recognized mealy amazon subspecies and the Atlantic Forest population of A. f. farinosa. Sampling locations and spectrograms of representative calls from each subspecies are shown in Figure 1. An initial visual inspection of the call spectrograms suggested there is some variation in the structure of contact calls between subspecies (Figure 1), however, plots of the first two components from a principal component analysis of acoustic parameter measurements and a multidimensional scaling analysis of the spectrographic cross-correlation matrix did not reveal distinct clusters of calls among the subspecies. Instead, we saw significant overlap of each clades' calls in acoustic space (Figure 4) due to a large amount of variation in calls within each subspecies. This pattern suggests that (a) substantial structural differences exist across each subspecies' range (e.g., the three example spectrograms of the A. f. farinosa subspecies of southern mealy amazon in Figure 1), and (b) variation within a subspecies is as great as variation among subspecies (Figure 4).

Assessment of Genetic Distance and Vocal Divergence
A qualitative assessment of the relationship between group mean genetic distance and average cross-correlation vocal similarity among clades of A. farinosa and A. guatemalae ( Figure 5) showed no clear pattern of association. This apparent lack of a correlation between vocal and genetic differentiation was supported by the Mantel test, which did not detect any statistical association between the two matrices (Mantel test: r = −0.012, matrix size = 6 clades, p = 0.451).

Assessment of Genetic Distance and Vocal Divergence
A qualitative assessment of the relationship between group mean genetic distance and average cross-correlation vocal similarity among clades of A. farinosa and A. guatemalae ( Figure 5) showed no clear pattern of association. This apparent lack of a correlation between vocal and genetic differentiation was supported by the Mantel test, which did not detect any statistical association between the two matrices (Mantel test: r = −0.012, matrix size = 6 clades, p = 0.451).

Discussion
We expanded upon a genetic dataset representing all recognized subspecies of the mealy amazon species complex to examine their relationships to the disjunct Atlantic Forest population of A. f. farinosa. In addition, we examined contact call similarity within and between these groups to determine if variation in this behavioral trait was correlated with underlying genetic structure. We found that while the Atlantic Forest population is recovered as a distinct clade, the degree of genetic differentiation between this population and the Amazon basin population of A. f. farinosa is substantially less than the degree of separation between the recognized northern and southern mealy amazon species. We also found that variation in a learned behavioral trait did not correlate to underlying genetic population structure. We discuss these results and their conservation implications in more detail below.

Genetic Relationships of the Mealy Amazon Clades
Our phylogeny and haplotype networks suggest that while the Brazilian Atlantic Forest population of A. f. farinosa is genetically and geographically distinct from other populations and subspecies of A. farinosa, it is not genetically distant enough to be considered a third cryptic species within the mealy amazon species complex. Specifically, the

Discussion
We expanded upon a genetic dataset representing all recognized subspecies of the mealy amazon species complex to examine their relationships to the disjunct Atlantic Forest population of A. f. farinosa. In addition, we examined contact call similarity within and between these groups to determine if variation in this behavioral trait was correlated with underlying genetic structure. We found that while the Atlantic Forest population is recovered as a distinct clade, the degree of genetic differentiation between this population and the Amazon basin population of A. f. farinosa is substantially less than the degree of separation between the recognized northern and southern mealy amazon species. We also found that variation in a learned behavioral trait did not correlate to underlying genetic population structure. We discuss these results and their conservation implications in more detail below.

Genetic Relationships of the Mealy Amazon Clades
Our phylogeny and haplotype networks suggest that while the Brazilian Atlantic Forest population of A. f. farinosa is genetically and geographically distinct from other populations and subspecies of A. farinosa, it is not genetically distant enough to be considered a third cryptic species within the mealy amazon species complex. Specifically, the betweengroup genetic distances for the three recognized southern mealy amazon subspecies and Atlantic Forest population had a range of 0.13-0.75% compared to genetic distances with a range of 1.9-3.8% between the northern and southern mealy amazon subspecies. Based on the magnitude of difference between the recently reclassified northern and southern mealy amazon species, we conclude that the Atlantic Forest population is best considered a genetically distinct population warranting formal recognition as a subspecies within the A. farinosa species.
Other studies aiming to identify potential cryptic species in parrots have employed comparable thresholds at which a genetically distinct clade should be considered a new subspecies or full cryptic species. In studies of three populations of the Cuban parrot subspecies Amazona leucocephala bahamensis, researchers concluded that sequence differences of 1.7-2.2% in the CR1 region warranted a subspecies designation for the Abaco phylogenetic species [33]. In the mulga parrot Psephotellus varius, an average sequence difference of 1.9% in the majority of the mitogenome (~84%) between the eastern and western populations was also justification for an elevation of the populations to separate subspecies [34]. The same conclusion was reached for northern and southern subspecies of the scarlet macaw (Ara macao) showing a 1.8% average sequence difference in a combined dataset of 12S, 16S, COI, and CytB [59], and for three recognized species within the Amazona ochrocephala complex with only 2% sequence differences in a combined dataset of ATP synthase 6 and 8, COI, ND2, and CytB despite their varied morphology [60]. In contrast, the authors of a study finding 4.4-5.1% sequence divergence in CytB between eastern and western populations of the ground parrot (Pezoporus wallicus) concluded that the western population should be reclassified as a new, cryptic species [35]. While the results of these studies are not directly comparable due to the various genes and analysis models used, the relative sequence differences serve as useful benchmarks when we consider the threshold at which a genetically distinct population is assigned a new taxonomic status. We should note that the range of genetic distances between the Atlantic Forest population of A. f. farinosa and all other A. farinosa subspecies for the CytB gene alone (0.76-0.99%) is equivalent to the interspecific ranges previously calculated for CytB in 88 avian genera [61]. However, given our more extensive data set and the established level of sequence divergence between the northern and southern mealy amazon species, we do not believe the Atlantic Forest population merits classification as a full species.
While it is likely our results represent the true phylogenetic relationship between the Atlantic Forest and Amazon basin populations of A. f. farinosa, there is a possibility our data suffer from a lack of geographic sampling diversity from the Atlantic Forest population (i.e., all individuals sampled originated from the same location). Wider geographic sampling of the mealy amazon across the Atlantic Forest may reveal additional genetic structure between populations within the Atlantic Forest (e.g., between the north-eastern region of the Atlantic Forest, often referred to as the Pernambuco Center of Endemism, and the southern or central regions of the Atlantic Forest) as well as between A. f. farinosa populations in the Amazon basin and Atlantic Forest. Patterns of intraspecific genetic differentiation within the Atlantic Forest are seen in other avian [62][63][64] and non-avian taxa [65] within the Atlantic Forest, and may be driven by the historic expansion of species from the Amazon basin to the central Atlantic Forest through gallery forests of the Cerrado region [42,66,67] or the relative stability and connectivity of the central and western regions of the Atlantic Forest in relation to historical climate patterns [68,69].

Utility of Behavioral Variation for Detecting Cryptic Species
The second goal of our study was to evaluate the potential utility of vocalization data in the identification of cryptic species. Given that behavioral data can often be collected with relative ease and at a very low cost, we aimed to determine whether variation in call structure correlated with genetic variation among groups. Such a relationship has been observed for a number of songbird species such as the variable antshrike (Thamnophilus caertulescens), red crossbill (Loxia curvirostra), and greenish warbler (Phylloscopus trochiloides) [13][14][15][16]. However, parrots and other open-ended vocal learners are able to alter their calls into adulthood, a phenomenon that may potentially lead to rapid call divergence or convergence irrespective of the underlying genetic structure of a population. For example, individually distinct call signatures may promote the recognition of social group members in fissionfusion societies and lead to a wide distribution of calls across the available acoustic space of a species or population [70]. Alternatively, because convergence on local call types (i.e., vocal dialects) signals group membership, the increased flexibility provided by open-ended learning may facilitate social integration throughout an individual's lifetime, especially after dispersal to unrelated populations [24], leading to an overall reduction in the total number of call types in a species or population.
Our results do not indicate a high degree of differentiation in contact call structure amongst groups. Cluster plots of each subspecies calls indicate substantial overlap of call structure and call features in acoustic space, and no correlation of call similarity with genetic structure (e.g., average cross-correlation similarity value for the Atlantic Forest population of A. f. farinosa and closely related A. f. chapmani is nearly identical, 0.419, to the similarity value of the Atlantic Forest population of A. f. farinosa and A. g. virenticeps, 0.418, from the northern species). Similarly, in yellow-naped amazons (Amazona auropalliata), a species well known for its regionally distinct vocal dialects, call diversity is unlinked to the underlying population genetic structure, indicating a strong preference of individuals to conform to local call types after dispersal or movement across dialect boundaries [71]. In contrast, the crimson rosella (Platycercus elegans) is one parrot species that shows clinal variation in their vocalizations and microsatellite genetic data [72]. Songbirds such as the white-crowned sparrow (Zonotrichia leucophrys) and suboscine passerines such as the variable antshrike (Thamnophilus caerulescens), which use innate vocalizations, also show greater correlation between acoustic and genetic variation [15,73,74], though this is not a consistent trend amongst all species (e.g., [65]).
Potential methodological limitations to our current study include inconsistencies in the quality and quantity of the available data from citizen science databases (e.g., lack of standardized recording protocols, variation in the type of recording equipment used, the unequal spatial distribution of sampling across each population's range). In addition, the call data were collected over a substantial timeframe (1950s to 2010s). Both factors may mask patterns of variation at the regional versus the subspecies level if call structures have changed over the decades or if information was lost in poor-quality recordings. These issues might be remedied by more aggressive pre-analysis filtering of the available data; however, attempting to apply stricter quality standards to our study led to a severe reduction in sample size that would have inhibited our current analyses. Overall, we suggest that vocalizations may be less useful as tools in the identification of cryptic species when those species are open-ended vocal learners, though further study is warranted.

Conservation Implications
Typically, conservation action is focused at the species level, creating a critical need for understanding the true relationship of disjunct and vulnerable populations to the species as a whole. In the case of the Atlantic Forest population of A. f. farinosa, our results do not support an elevation in species status but do indicate this population represents a unique pool of genetic diversity within the mealy amazon species complex. Currently, it is estimated that 7-16% of historic Atlantic Forest cover remains, and rates of habitat degradation or destruction have increased significantly in the past three decades [48,49], leading to restrictions in the ranges of local mealy amazon populations and increased conflict with humans (A. Saidenberg, unpublished data). Shrinking refuges could further endanger a population that, like most large parrot species, is already facing substantial threats from poaching and harvest for the pet trade [26,27]. For these reasons, we advocate for the special consideration of the Atlantic Forest population of mealy amazons in conservation and management decisions for the species, including the formal recognition of this population as a subspecies of the southern mealy amazon. Finally, we strongly encourage the continued study of parrot species complexes to help identify additional taxonomic groups warranting classification and protection.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/d13060273/s1, Table S1: Metadata and gene regions sequenced of the samples used in this study. Table S2: Metadata for the sound files used in this study. Table S3: Primer sequences used to amplify and sequence the gene regions used in this study. Table S4: Genbank accession numbers for new Atlantic forest A. f. farinosa gene sequences used in this study. Table S5: Acoustic parameters generated by the 'specan' function and used in our principal component analysis.