Adaptive radiations have helped shape how we view animal speciation, particularly classic examples such as Darwin's finches, Hawaiian fruitflies and African Great Lakes cichlids. These ‘island’ radiations are comparatively recent, making them particularly interesting because the mechanisms that caused diversification are still in motion. Here, we identify a new case of a recent bird radiation within a continentally distributed species group; the capuchino seedeaters comprise 11 Sporophila species originally described on the basis of differences in plumage colour and pattern in adult males. We use molecular data together with analyses of male plumage and vocalizations to understand species limits of the group. We find marked phenotypic variation despite lack of mitochondrial DNA monophyly and few differences in other putatively neutral nuclear markers. This finding is consistent with the group having undergone a recent radiation beginning in the Pleistocene, leaving genetic signatures of incomplete lineage sorting, introgressive hybridization and demographic expansions. We argue that this apparent uncoupling between neutral DNA homogeneity and phenotypic diversity is expected for a recent group within the framework of coalescent theory. Finally, we discuss how the ecology of open habitats in South America during the Pleistocene could have helped promote this unique and ongoing radiation.
Bird adaptive radiations have contributed greatly to our understanding of the process of speciation. Morphological innovations under specific ecological conditions have fuelled radiations like those of the drepanidine Hawaiian honeycreepers , while sexual selection may have contributed to forming speciose groups such as the Dendroica wood warblers . These, as well as many other celebrated examples (e.g. Phylloscopus warblers ; Diglossa flowerpiercers ) represent radiations that have accrued species over millions of years. There are comparatively far fewer known monophyletic assemblages of bird species originating recently over a short time span; perhaps, the most notable cases include Darwin's finches , the red crossbills of North America , African indigobirds  and Neotropical cardueline finches . These examples, together with other recent vertebrate radiations such as that of the cichlid fishes from the African Great Lakes , do not fully comply with the Biological Species Concept and thus provide key insight as to how permeable mechanisms of reproductive isolation and ultimately introgressive hybridization might contribute to the generation of biodiversity [10–13]. Clearly, additional examples of recent radiations in different ecological and historical contexts will provide further insights into the mechanisms that underpin diversification and speciation.
In this study, we identify a novel case of a recent continental bird radiation within the taxonomically complex Neotropical genus Sporophila. The capuchinos, as they are known colloquially, comprise a group of 11 seedeaters that are smaller than their congeners and characterized by cinnamon-based plumage colour patterns in reproductive males . Capuchinos are sexually dimorphic, with divergent plumage patterns among males, while females are mostly brown and olive with little apparent distinction among them. Capuchino species are thus diagnosed based on male nuptial plumage, because they otherwise show little differentiation in size and shape . Typically capuchino species are broadly sympatric (or even syntopic) with several other species of the group . Many are now rare, having limited geographical ranges with populations in decline owing to habitat loss and trapping for the pet trade . Most capuchinos are seasonal migrants, but little is known of the location of their wintering grounds. When not breeding, they are commonly seen in mixed flocks showing similar foraging behaviour and in some cases eclipse plumage . Phylogenetic analyses of mitochondrial DNA (mtDNA) sequence data suggest that the group is monophyletic and indicate that two of the 11 species (Sporophila minuta and Sporophila castaneiventris) are clearly diverged from the other members of the clade [16,17]. However, the remaining nine species (Sporophila bouvreuil, Sporophila cinnamomea, Sporophila hypochroma, Sporophila hypoxantha, Sporophila melanogaster, Sporophila nigrorufa, Sporophila palustris, Sporophila ruficollis and Sporophila zelichi) exhibit extremely low interspecific sequence divergence and apparent lack of reciprocal monophyly. These species are found predominantly south of the Amazon River and hereafter we refer to them as the southern capuchinos (see figure 1 and electronic supplementary material, figure S1 for representative illustrations and range maps).
Funk & Omland  and more recently McKay & Zink  examined the most commonly cited causes of lack of reciprocal monophyly at mitochondrial loci in animals, identifying imperfect taxonomy as the most frequent, a problem that may be even more severe in groups such as invertebrates that have received less attention. Other major causes are inadequate phylogenetic resolution (caused by too few phylogenetically informative characters or homoplasy), incomplete lineage sorting and interspecific hybridization. We use this conceptual framework to understand the processes responsible for the mitochondrial genetic pattern observed among the southern capuchinos [16,17], augmenting previous genetic analyses to include additional mitochondrial regions and both nuclear DNA microsatellites and sequences. We expand the analyses to include objective quantification of plumage and song traits among males, two key components of bird mate recognition systems . We find marked differences in colour and song among southern capuchino species, suggesting that their current designations are not artefacts of taxonomic practice. This result contrasts with the lack of neutral genetic differentiation among species. Overall, our genetic data imply that the southern capuchinos began to radiate within the Pleistocene but that the process is ongoing, consistent with demographic expansions from small ancestral populations and showing evidence of both incomplete lineage sorting and gene flow between species.
2. Material and methods
(a) Coloration analyses
The distinct differences in reproductive male plumage patterns among southern capuchinos provide the main suite of characters used to diagnose the nine currently recognized species. Because signals can only be understood with reference to the natural receiver, the differences in human and avian visual systems and colour perception can make human description of plumage colour inadequate for the study of many biological questions (see Benites et al.  and references therein). Thus, using species designations as a working hypothesis, we objectively quantified relative plumage differences (within the southern capuchinos and with respect to three outgroup species) using an avian visual perspective as modelled by Vorobyev & Osorio . Data were collected from adult male museum study skins: 104 individuals from seven southern capuchino species and 61 specimens from three non-capuchino Sporophila species (see electronic supplementary material, table S1). Plumage coloration reflectances were measured with a spectrometer on seven patches: back, belly, chest, crown, nape, throat and rump. A given colour stimulus can be described by the quantum catches of the light entering the eye by each receptor type present in the avian retina (Q1 = ultraviolet wavelength-sensitive; Q2 = short wavelength-sensitive; Q3 = medium wavelength-sensitive; and Q4 = long wavelength-sensitive). Colour differences were evaluated for each plumage region by estimating a distance in avian perceptual colour space (ΔS), representing the value by which any two plumage colour patches differ based on the quantum catches of each cone receptor and their respective noise-to-signal ratio . Since neither spectral sensitivity data nor cone cell type proportions are available for any of the Sporophila species, the blue tit (Parus caeruleus) was used as a representative passerine visual system (for justifications see the electronic supplementary material, S1 materials and methods). We calculated ΔS between identical patches of each pair of species and values were averaged across the seven measured patches to obtain a global ΔS representing overall plumage differences between taxon pairs. A neighbour-joining (NJ) tree was built from average pairwise ΔS distances, a strategy used throughout this study to clearly display pairwise distance matrices. A ΔS value of 1.0 ‘just noticeable difference’ (jnd) is taken as a general threshold value for the discrimination of two colours [21,22]. We also assessed relative plumage differences among taxa as well as the contribution of different spectral regions through a stepwise discriminant function analysis (DFA). A pairwise distance matrix between species was constructed using Euclidean distances between centroids in the space defined by the first two DFA functions. Statistical differences in DFA1 scores between pairs of taxa were assessed using multiple Mann–Whitney U-tests with sequential Bonferroni corrections . Finally, DFA results were compared with those obtained from a principal components analysis (PCA), which requires no a priori designations of species based on human perception. For details see the electronic supplementary material, S1 materials and methods.
(b) Song analyses
Our song dataset (electronic supplementary material, table S2) included vocalizations belonging to 131 individuals from 10 capuchino species and 30 individuals from three non-capuchino Sporophila species. We measured 14 variables for all songs and differences between species were assessed through a stepwise DFA, predicting group membership for each individual through a jackknifed classification procedure. We also classified syllables into six categories and calculated the overall proportion of each category used in the song of every species. Individual variation was assessed using a modified Jaccard similarity index . We calculated the mean among and within species values for the similarity index as well as 95% confidence intervals (CI). Transforming the index by 1 − x, we generated a species pairwise distance matrix (0 now representing identical syllable use). For details see the electronic supplementary material, S1 materials and methods.
(c) Genetic analyses
All 11 capuchinos were included in our genetic analysis and we used four other sympatric species as outgroups. When available, we included multiple individuals per species and from as many localities of their geographical distribution as possible. Our initial dataset included samples for which we previously obtained cytochrome c oxidase I (COI) sequences  to which we added cytochrome b (Cyt b) and mitochondrial control region (CR) data. For analyses based on the frequency of microsatellite DNA alleles, we increased sample sizes by including newly collected specimens and also museum study skins. Several nuclear markers were amplified with the objective of finding fixed single nucleotide polymorphisms between southern capuchino species. Details are provided in the electronic supplementary material, S1 materials and methods and table S3.
We constructed phylogenetic trees with MrBayes v. 3.1.2 [25,26] using concatenated mtDNA sequence data and used haplotype frequencies between southern capuchinos to estimate pairwise FST values for each gene with Arlequin v. 126.96.36.199 . Significance was tested through 1000 random permutations with sequential Bonferroni corrections . Pairwise FST matrices were displayed using NJ trees. Node ages were estimated using time to most recent common ancestor (TMRCA) with the BEAUti/BEAST v. 1.4.8 package  and data from Cyt b. For DNA microsatellite loci, genetic differentiation between species was assessed using FST, Nei's standard genetic distance  and Structure v. 2.3.3 . We also performed a principal coordinates analysis (PCoA) with GenAlEx v. 6  using distances calculated between all individuals genotyped. The southern capuchinos did not show fixed sites in CHD1Z, Fib5, MUSK, Numt2 or Numt3; however, Numt2 provided good resolution with respect to the outgroups and was used further to build phylogenetic trees. We tested whether southern capuchinos differed in the haplotype frequencies of the nuclear loci sequenced. A combination of the number of segregating sites and our sample sizes made recovery of useful haplotypes from all markers low. To increase the recovery, we arbitrarily divided the sequence of each gene in half and inferred the haplotypes separately for both the 5′ and 3′ portions of each locus using DnaSP v. 5.10 . Differences in haplotype frequencies between species were assessed using FST calculations in both the 5′ and 3′ fraction of each locus separately.
Mitochondrial and nuclear sequence data were used to explore the demographic history of the species in the southern capuchino clade. Our microsatellite markers showed evidence of linkage disequilibrium and deviation from Hardy–Weinberg expectations (see the electronic supplementary material, S1 materials and methods) and were excluded from these analyses. We tested for demographic expansions or contractions by using Fu's F-test  and by calculating the exponential growth parameter g . Splitting times between pairs of species, migration (introgression), and ancestral and current effective population sizes were estimated using IMa2 . More information is provided in the electronic supplementary material, S1 material and methods.
(a) Male southern capuchinos show significant phenotypic differences in coloration patterns and song
Seven southern capuchino and three other Sporophila species analysed differ markedly in plumage coloration. The average ΔS values (metric of colour differences in avian perceptual colour space) between pairs of species were in all cases above the discernable threshold according to the avian visual model used (figure 2a), suggesting that there are large overall differences in plumage among southern capuchino males. Moreover, males differed significantly in virtually every patch considered (electronic supplementary material, table S4). Interestingly, some southern capuchino males are more similar in plumage attributes to outgroup species than to members of their own clade (figure 2a). Statistically significant differences in plumage coloration among the species surveyed were also found through DFA (Wilk's λ = 1.26 × 10−6, χ1262 = 1888.53, p < 0.0001; figure 2b). Factor loadings for the 14 out of 28 variables retained by the stepwise DFA procedure (electronic supplementary material, table S5) suggest that every portion of the light spectrum and every body part measured (except the rump) accounts for differences among species, with the belly and crown contributing most strongly. Statistically significant differences between southern capuchino species were found in average DFA1 scores (figure 2c). Finally, we obtained similar results when plumage coloration differences were assessed at the individual level through PCA (electronic supplementary material, figure S2).
Songs of southern capuchino and other Sporophila species comprise distinct series of complex syllable types (electronic supplementary material, figure S3). Stepwise DFA (figure 3a) retained 13 out of 14 variables (see the electronic supplementary material, table S6 for factor loadings) and showed statistically significant differences among species (Wilk's λ = 3.36 × 10−4, χ1562 = 1175.82, p < 0.001). Our analysis distinguishes southern capuchino songs from those of S. castaneiventris and S. minuta (figure 3a) and more markedly from non-capuchino species. In comparison, song differences within the southern capuchinos are smaller with 66.5 per cent of cases reclassified to the correct species (41 of 54 misclassifications involving southern capuchinos). Despite this, significant differences are also found when outgroup species are excluded (Wilk's λ = 0.192, χ235 = 154.18, p<0.001; figure 3b). In this analysis performed only with the southern capuchinos, DFA1 represents variation in syllables that span a small range of frequencies, while DFA2 represents longer, more complex syllables that span a wide range of frequencies (see the electronic supplementary material, table S7). The cross-validation procedure had low success in reclassifying songs to the correct species (51.5% with a total of 49 errors). Although there is an overlap of attributes in southern capuchino vocalizations, clustering can be observed between songs from the same species (figure 3b), and significant differences in DFA1 scores exist among various species (figure 3c).
Sporophila species differ in the use of the six syllable categories in their songs (figure 3d), with larger differences seen between southern capuchinos and outgroup species than within the former group. A modified Jaccard similarity index  calculated between every pair of songs allowed us to assess variation in syllable-type usage within and between species. The average intraspecific value was 0.87 (range: 0.74–1), with an average lower 95% CI limit of 0.84. Interspecific values ranged from 0.95 (S. palustris versus S. zelichi) to 0.67 (S. melanogaster versus S. ruficollis). Only the species pair for which our sample sizes were the smallest (S. collaris and S. zelichi) had a 95% CI that included 1, suggesting that all other comparisons differed significantly. Approximately half of the 95% CIs of pairwise inter-southern capuchino comparisons exceeded 0.84, the threshold for intraspecific variation in syllable usage (figure 3e).
(b) Low genetic differentiation between southern capuchino species
Both the overall capuchinos and southern capuchinos are clearly monophyletic with high support. Sporophila minuta is the sister species to the southern capuchino radiation. By contrast, none of the southern capuchino species is monophyletic for mtDNA, despite including the rapidly evolving CR in our analyses (figure 4a). Indeed, all analyses using a suite of neutral genetic markers reveal extremely low genetic differentiation among southern capuchinos. All southern capuchinos do show statistically significant pairwise FST values with respect to S. minuta (figure 4a; average of 0.89 for COI; 0.76 for Cyt b; and 0.74 for CR; p < 0.01). Sporophila melanogaster has low but significant FST values with two other southern capuchinos consistently across mitochondrial loci: Sporophila hypoxantha (COI: 0.24; and CR: 0.15) and Sporophila ruficollis (COI: 0.47; Cyt b: 0.24; and CR: 0.26).
For DNA microsatellites, PCoA reveal lack of genetic structure within the southern capuchinos (electronic supplementary material, figure S4). We found evidence of linkage disequilibrium between some loci (see the electronic supplementary material, S1 materials and methods), contra previous studies [36,37]. This observation might result from undetected allele dropout or null alleles. Despite these potential issues, our Structure Bayesian assignment test results corroborate those from the PCoA, implying that the southern capuchino samples most probably comprise a single genetic cluster (K = 1). FST values between southern capuchino species were always below 0.01 and not significantly different from zero (figure 4b). Even though sample sizes for outgroup species were low and homoplasy could have overwritten some genetic signal, we found nearly all FST and Nei's genetic distance values for comparisons involving southern capuchinos and outgroups, or between outgroups, to be statistically significant and higher than intra-southern capuchino comparisons (figure 4b and the electronic supplementary material, figure S4).
Southern capuchino species did not exhibit fixed diagnostic sites in the five nuclear markers used in this study and the Bayesian tree produced using Numt2 DNA sequence data did not resolve below the level of the capuchinos (electronic supplementary material, figure S5). However, frequency-based haplotype differences were evident and most comparisons involving southern capuchinos and outgroups showed high and statistically significant FST values. Some significant FST comparisons were also found between S. hypoxantha and S. ruficollis (5′Numt2: 0.29; 3′Numt2: 0.15; and 5′MUSK: 0.36), S. hypoxantha and S. cinnamomea (3′Numt3: 0.51), and S. palustris and S. ruficollis (5′Numt2: 0.40) (figure 4c).
(c) Demographic history of the southern capuchino clade
Origins of the southern capuchinos date to the Pleistocene. The calculated TMRCA using the coalescent program BEAST (applying a 2.1 per cent per million year divergence time for Cyt b [28,38]) suggests that the mitochondrial ancestor of all southern capuchinos is approximately 422 000 years old (95% high posterior density interval: 211 000–690 000). The age of the ancestor between the southern capuchino clade and its sister species S. minuta was estimated in about 681 000 years before present (95% high posterior density interval: 409 000–989 000). Our analyses using IMa2  suggest that there is considerable gene flow among southern capuchinos (table 1). By contrast, there is practically no gene flow between the southern capuchino clade and its sister species S. minuta. The average splitting time between southern capuchinos calculated with IMa2 is approximately 110 000 years before present. The estimated splitting time between the southern capuchinos and S. minuta is approximately 140 000 years before present. These results must be interpreted with caution as parameters estimated using IMa2 relied on rough molecular clock calibrations and an approach that did not consider all southern capuchinos interacting at the same time in the model. As a consequence of these caveats, we consider that the most conservative interpretation of these results is to say that, overall, our data suggest that the southern capuchino radiation dates approximately to the period stretching from the middle Pleistocene to present. Average current effective population size estimates for southern capuchino species were approximately ninefold larger than the average ancestral population size (table 1), suggesting the possibility of demographic expansions. This is consistent with Fu's F-values, which were negative and statistically significant across markers for the southern capuchino clade itself, and for species for which sample sizes allowed us to perform the test (see the electronic supplementary material, table S8). The exponential growth parameter g was positive for the southern capuchino clade, S. hypoxantha and S. palustris, providing further evidence of demographic expansion during the southern capuchino radiation (see the electronic supplementary material, table S8).
We find extraordinary variation in phenotype (coloration and song) among southern capuchinos despite extremely low levels of neutral genetic differentiation and no species-level monophyly in any genetic marker assayed among any of the designated capuchino species. Our analyses are also consistent with demographic expansions in the clade as well as gene flow between some of its member species. Taken together, we suggest that this is a compelling example of an extremely rapid, recent and ongoing continental radiation, with species diverging in male plumage coloration patterns and song. Both incomplete lineage-sorting and introgressive hybridization appear to be responsible for the blurry genetic identity of these species, at least at neutral loci. Insufficient data as an explanation for lack of reciprocal monophyly between species  are not plausible as we see similar patterns across genetic markers with different modes of inheritance, mutation mechanisms and rates of divergence. Interestingly, the southern capuchinos add another example to the list of known rapid seedeater radiations [5–8], providing further support to the idea that finch-like birds may be particularly inclined to undergo such rapid bursts of speciation . Below, we discuss our findings, arguing that the current practice of recognizing multiple southern capuchino species is correct, and make inferences about the origin of the group, as well as its future in relation to the main conservation threats it faces.
(a) Are the southern capuchinos incipient species?
The southern capuchinos could not be diagnosed using neutral markers, a situation similar to that of Darwin's ground finches [39,40]. This result is expected if speciation events were recent and close together in time and does not necessarily imply panmixia. This genetic pattern has been termed hemiplasy  and is the result of rapid bursts of speciation causing ancestral polymorphisms to persist across nodes in a species tree . This imposes limitations on our ability to date speciation events that go beyond the technical difficulties related to the programs and molecular clock calibrations used. In this case, the age of the ancestor between the southern capuchinos and the sister species S. minuta would set the upper bound for speciation times (dating approximately to the middle Pleistocene). The lack of monophyly in neutral markers within species need not imply that these species are not valid biological entities. Differences in vocalizations, plumage or other morphological characters subjected to selection pressures (particularly in isolated populations) may evolve rapidly, creating a lag before the stochastic sorting of neutral molecular markers produces diagnosable lineages . Thus, differences in genes that underpin these diagnostic phenotypic characters would correctly delimit species, while neutral markers show the signature of incomplete lineage sorting or evidence of introgression . Genomic approaches such as amplified fragment length polymorphisms or next-generation sequencing will aid us in finding the crucial genes implicated in the evolution of the phenotypic differences among southern capuchinos.
Contrasting with lack of differentiation in neutral markers, we find striking diversity in male coloration pattern and song, suggesting that the existing species in the group are not a taxonomic artefact but a compelling example of an ongoing but incomplete radiation. This is further substantiated by cryptic plumage differences found between females in the group . Playback experiments are required to test if the observed phenotypic differences result in assortative mating. Thus far, an analysis of aggressive response to playback between S. hypoxantha and S. palustris males suggests that at least some southern capuchinos distinguish between conspecific and heterospecific song, virtually ignoring the latter (P. Benites, L. Campagna & P. L. Tubaro 2010, unpublished data). Future effort should focus on the phenotypic differences among species and their possible role as mechanisms of reproductive isolation.
(b) Inferences from the past and predictions for the future
Our data are consistent with the southern capuchinos originating from small isolated populations, expanding later to occupy their current ranges. Range fragmentation, as a consequence of Late Pleistocene marine ingressions and egressions, accelerated by sexual selection has been proposed as a possible cause for the radiation of the southern capuchinos [16,44]. Alternatively, Fjeldså & Rahbek  propose that specific habitat requirements could have kept small local populations isolated in different patches of suitable grassland habitat. Palaeoclimatic reconstructions indicate that open habitats (savannah, Caatinga and Cerrado) dominated South America during the last glacial maximum (approx. 27 000 years before present), accompanied by a marked retraction of the rainforest . There is evidence of south polar air incursions affecting temperatures and precipitation regimes during the various glacial cycles of the Quaternary, causing changes in the length of dry seasons and cycling between overall drier climate and humid conditions. These conditions are thought to have caused expansion and regression of open habitats and rainforest that continued into the Holocene [46–48]. The fluctuation in the prevalence of rainforest over open habitats and vice versa and the inter-digitation of these two biomes could have contributed to isolating small populations in islands of suitable habitat, or grassland refugia, in a scenario analogous to that proposed in the forest refugia hypothesis . Indeed, the contraction of open habitats during the Holocene owing to the changes in forest cover is thought to have affected the composition of communities of savannah specialists, leading to extinctions of several large mammals . Both marine ingressions and the existence of open habitat refugia could have contributed to isolating small ancestral populations of southern capuchinos, making it hard to distinguish these alternatives.
The future of the southern capuchinos is without doubt tied to human use of land. Expanding industrial agriculture is a great threat to grasslands and only small areas remain undisturbed , posing a risk to the highly sensitive southern capuchino species . The fact that southern capuchinos have been forced into the remaining patches of suitable habitat carries potential consequences for the viability of small populations with putatively incomplete mechanisms of reproductive isolation , and the very real possibility of hybridization. Forced breeding in syntopy could increase opportunities for misimprinting on heterospecifics that would lead to eventual hybridization [54–56]. Moreover, theory predicts that hybridization will increase when species are forced together into small patches of habitat, with habitat disturbance and when potential conspecific mates are rare [57–59], as is now the case for many species of this group. The latter possibility also seems plausible because southern capuchinos differ in their abundance, with taxa like S. hypoxantha and S. ruficollis greatly outnumbering the rare S. palustris, S. melanogaster, S. nigrorufa or S. zelichi. Ultimately, the fate of this unique group of birds may not be as enigmatic as its origin.
Field work for this study was authorized by the National Fauna authorities of Argentina, the Argentine National Parks Administration and the Office of Fauna of the Province of Corrientes, Argentina.
For tissue loans, we are indebted to A. T. Peterson and M. B. Robbins (KUNHM); R. T. Brumfield and D. L. Dittmann (LSUMZ); K. Cook and M. Adams (BMNH); C. S. Fontana (MCP); J. Fjeldså and J. B. Kristensen (ZMUC); E. Bermingham (STRI); M. J. Braun and J. Dean (USNM). For authorizing permits in Argentina, we thank the National Fauna authorities, the National Parks Administration and the Office of Fauna of Corrientes. We thank T. Webber (FLMNH), T. Bishop (Macaulay Library, Cornell Laboratory of Ornithology), J. Mazar Barnett and R. Straneck for providing recordings. We thank the Willi Hennig Society for a free edition of TNT 1.1. This project was supported by grants to PLT from ANPCyT, Argentina (PICT 2004-16-25171 and 2010-0805); CONICET, Argentina (PIP 112-200801-00741); UBA, Argentina (UBACyT 2010-2012); IDRC, Canada, and Richard Lounsbery Foundation; and by a grant from NSERC, Canada to SCL. Photographs in figure 1 by R. Güller (S. bouvreuil pileata and S. ruficollis) and P. Fenalti (S. melanogaster). We thank T. Price and two anonymous reviewers for their helpful comments on previous versions of this manuscript.
- Received October 14, 2011.
- Accepted November 11, 2011.
- This journal is © 2011 The Royal Society