Benthic foraminifera in Tertiary sediments of DSDP Leg 90 holes

Eocene through Pliocene benthic foraminifers were examined from seven sites located at middle and lower bathyal depths on the Lord Howe Rise in the Tasman Sea, from another site at lower bathyal depths in the Coral Sea, and from a site in the intermediate-depth, hemipelagic province of the Chatham Rise, east of southern New Zealand. Age-related, depth-related, and bioprovincial faunal variations are documented in this chapter. One new species, Rectuvigerina tasmana, is named. The paleoecologic indications of several key groups, including the miliolids, uvigerinids, nuttallitids, and cibicidids, are combined with sedimentologic and stable isotopic tracers to interpret paleoceanographic changes in the Tasman Sea. Because the total stratigraphic ranges of many bathyal benthic foraminifers are not yet known, most endpoints in the Tasman Sea are considered ecologically controlled events. The disappearances of Uvigerina rippensis and Cibicidoidesparki and the first appearances of U. pigmaea, Sphaeroidina bulloides, and Rotaliatina sulcigera at the Eocene/Oligocene boundary can be considered evolutionary events, as also can the first appearance of Cibicides wuellerstorfi in Zone NN5. Species which are restricted to the lower bathyal zone except during discrete pulses, most of which are related to the development of glacial conditions, include Melonis pompilioides, M. sphaeroides, Pullenia quinqueloba, Nuttallides umbonifera, and U. hispido-costata. Middle bathyal indigenes include U. spinulosa, U. gemmaeformis, Ehrenbergina marwicki, R. sulcigera, and all rectuvigerinids except Rectuvigerina spinea. Although the miliolids first occurred at lower bathyal depths, they were more common in the middle bathyal zone. Although the Neogene hispido-costate uvigerinids first developed at lower bathyal depths and at higher middle latitude sites, in the later Neogene this group migrated to shallower depths and became predominant also in the middle bathyal zone. Despite the relatively similar sedimentologic settings at the six middle bathyal Tasman sites, there was extensive intrageneric and intraspecific geographic variation. Mililiolids, strongly ornamented brizalinids, bolivinitids, Bulimina aculeata, Osangularia culter, and strongly porous morphotypes were more common at higher latitudes. Osangularia bengalensis, striate brizalinids such as Brizalina subaenariensis, Gaudryina solida, osangularids in general, and finely porous morphotypes were more common in the subtropics. There was strong covariance between faunas at lower middle latitude, lower bathyal Site 591, and higher middle latitude, middle bathyal Site 593. The following oceanographic history of the Tasman Sea is proposed; using the stable isotopic record as evidence for glacials and examining the ecologic correlations between (1) miliolids and carbonate saturation, (2) nuttallitids and undersaturated, cooled, or "new" water masses, (3) uvigerinids with high organic carbon in the sediment and high rates of sediment accumulation, and (4) cibicidids and terrestrial organic carbon. The glacial located near the Eocene/Oligocene boundary is characterized by the penetration of cooler, more corrosive waters at intermediate depths in high southern latitudes. This may have caused overturn, upwelling pulses, in other Tasman areas. The development of Neogenelike conditions began in the late Oligocene (Zone NP24/NP25) with the evolution of several common Neogene species. A large number of Paleogene benthics disappeared gradually through the course of the early Miocene, which was not well preserved at any Tasman site. Corrosive conditions shallowed into the middle bathyal zone in several pulses during the early Miocene. The development of glacial conditions in the middle Miocene was accompanied by major changes throughout the Tasman Sea. Sediment accumulation rates increased and high-productivity faunas and corrosive conditions developed at all but the lowest-latitude Site 588. This increase in productivity and accumulation rate is attributed to the eutrophication of Antarctic water masses feeding Tasman current systems, as well as to invigorated circulation in general. It overlaps with the beginning of the Pacific High-productivity Episode (10-5 Ma). During the latest Miocene glacial episode, corrosive conditions developed at lower bathyal depths, while cooler water and lower nutrient levels shallowed to middle bathyal depths. Lower input of terrestrial organic carbon may be related to the lower nutrient levels of this time and to the termination of the Pacific High-productivity Episode. The moderate glacial episode during the mid-Pliocene (Zone NN15/NN16, ~3.2 Ma) corresponds to a decline in sediment accumulation rates and a reorganization of faunas unlike that of all other times. New genera proliferate and indices for cool, noncorrosive conditions and high organic carbon expand throughout the middle bathyal zone coeval with the sedimentation rate decreases. By the latest Pliocene (about 2.5 Ma), however, during another glacial episode, faunal patterns typical of this and later glacials develop throughout the Tasman Sea. Benthic foraminiferal patterns suggest increased input of terrestrial organic matter to Tasman Sea sediments during this episode and during later glacials.

Complexities in the genetic structure of Anopheles gambiae populations in west Africa as revealed by microsatellite DNA analysis

Chromosomal forms of Anopheles gambiae, given the informal designations Bamako, Mopti, and Savannah, have been recognized by the presence or absence of four paracentric inversions on chromosome 2. Studies of karyotype frequencies at sites where the forms occur in sympatry have led to the suggestion that these forms represent species. We conducted a study of the genetic structure of populations of An. gambiae from two villages in Mali, west Africa. Populations at each site were composed of the Bamako and Mopti forms and the sibling species, Anopheles arabiensis. Karyotypes were determined for each individual mosquito and genotypes at 21 microsatellite loci determined. A number of the microsatellites have been physically mapped to polytene chromosomes, making it possible to select loci based on their position relative to the inversions used to define forms. We found that the chromosomal forms differ at all loci on chromosome 2, but there were few differences for loci on other chromosomes. Geographic variation was small. Gene flow appears to vary among different regions within the genome, being lowest on chromosome 2, probably due to hitchhiking with the inversions. We conclude that the majority of observed genetic divergence between chromosomal forms can be explained by forces that need not involve reproductive isolation, although reproductive isolation is not ruled out. We found low levels of gene flow between the sibling species Anopheles gambiae and Anopheles arabiensis, similar to estimates based on observed frequencies of hybrid karyotypes in natural populations.

Ethnohistory, genetics, and cancer mortality in Europeans

Geographic variation in cancer rates is thought to be the result of two major factors: environmental agents varying spatially and the attributes, genetic or cultural, of the populations inhabiting the areas studied. These attributes in turn result from the history of the populations in question. We had previously constructed an ethnohistorical database for Europe since 2200 B.C., permitting estimates of the ethnic composition of modern European populations. We were able to show that these estimates correlate with genetic distances. In this study, we wanted to see whether they also correlate with cancer rates. We employed two data sets of cancer mortalities from 42 types of cancer for the European Economic Community and for Central Europe. We subjected spatial differences in cancer mortalities, genetic, ethnohistorical, and geographic distances to matrix permutation tests to determine the magnitude and significance of their association. Our findings are that distances in cancer mortalities are correlated more with ethnohistorical distances than with genetic distances. Possibly the cancer rates may be affected by loci other than the genetic systems available to us, and/or by cultural factors mediated by the ethnohistorical differences. We find it remarkable that patterns of frequently ancient ethnic admixture are still reflected in modern cancer mortalities. Partial correlations with geography suggest that local environmental factors affect the mortalities as well.

Genetics and the origin of European languages.

A new set of European genetic data has been analyzed to dissect independent patterns of geographic variation. The most important cause of European genetic variation has been confirmed to correspond to the migration of Neolithic farmers from the area of origin of agriculture in the Middle East. The next most important component of genetic variation is apparently associated with a north-south gradient possibly due to adaptation to cold climates but also to the differentiation of the Uralic and the Indo-European language-speaking people; however, the relevant correlations are not significantly different from zero after elimination of the spatial autocorrelation. The third component is highly correlated with the infiltration of the Yamna ("Kurgan") people, nomadic pastoralists who domesticated the horse and who have been claimed to have spread Indo-European languages to Europe; this association, which is statistically significant even when taking spatial autocorrelations into account, does not completely exclude the hypothesis of Indo-European as the language of Neolithic farmers. It is possible that both expansions were responsible for the spread of different subfamilies of Indo-European languages, but our genetic data cannot resolve their relative importance.

Descrição de um fóssil de Eira barbara (Carnivora, Mustelidae) do Pleistoceno final (Bacia do Acre) e morfologia comparada do sincrânio de E. barbara recente: implicações paleobiogeográficas, geográficas e taxonômicas

A espécie recente Eira barbara (Mustelidae, Carnivora) possui uma distribuição geográfica desde o México até o norte da Argentina. É um importante táxon a ser estudado como modelo anatômico dentre os mustelídeos, assim como um importante modelo para uma melhor compreensão e entendimento sobre a diversificação dos Mustelidae. Atualmente, os registros fósseis de E. barbara na América do Sul são bastante escassos e restritos às idades pleistocênicas, sendo que os estudos destes fósseis são frequentemente desprovidos de maiores esforços para realização de descrições morfológicas detalhadas e de estudos paleobiogeográficos. Assim como os estudos dos fósseis de E. barbara são limitados, constatou-se que o mesmo cenário é observado quanto aos estudos sobre a morfologia e biogeografia da espécie. Desta forma, o presente trabalho se propôs a: realizar uma revisão de todos os registros da espécie e de uma forma geral, contribuir para um melhor conhecimento sobre a morfologia sincraniana e sobre a história biogeográfica e paleobiogeográfica de E. barbara. Para tanto, os seguintes objetivos foram propostos: estudo e redescrição detalhada do fóssil UFAC-PV 036, proveniente do Pleistoceno final do Alto Rio Juruá do sudoeste da Amazônia Brasileira; descrição sincraniana comparada de estruturas morfológicas externas e internas, analisando caracteres intraespecíficos da espécie E. barbara; realização de análises multivariadas a fim de investigar variações geográficas sob o uso de caracteres craniométricos de E. barbara entre os diferentes biomas brasileiros. A revisão dos registros fósseis foi de grande importância para o estabelecimento dos verdadeiros registros de Eira na América do Sul e a redescrição de UFAC-PV 36 contribui para o melhor conhecimento morfológico e paleobiogeográfico da espécie. A descrição morfológica comparada do sincrânio de E. barbara contribui de forma significativa para o conhecimento sobre a morfologia da espécie bem como, a descrição de caracteres intraespecíficos proporcionam caracteres mais apropriados em matrizes morfológicas, fornecendo maior robustez nas análises filogenéticas futuras. Este trabalho propõe que E. barbara não possui diferenças craniométricas estatisticamente significativas entre os biomas brasileiros, porém, E. barbara caracteriza-se aqui como uma espécie dimórfica, na qual os machos possuem estruturas cranianas relativamente maiores do que as fêmeas.

Contrasting quantitative traits and neutral genetic markers for genetic resource assessment of Mesoamerican Cedrela odorata

We compared within-population variability and degree of population differentiation for neutral genetic markers (RAPDS) and eight quantitative traits in Central American populations of the endangered tree, Cedrela odorata. Whilst population genetic diversity for neutral markers (Shannon index) and quantitative traits (heritability, coefficient of additive genetic variation) were uncorrelated, both marker types revealed strong differentiation between populations from the Atlantic coast of Costa Rica and the rest of the species' distribution. The degree of interpopulation differentiation was higher for RAPD markers (F-ST 0.67 for the sampled Mesoamerican range) than for quantitative traits (Q(ST) = 0.30). Hence, the divergence in quantitative traits was lower than could have been achieved by genetic drift alone, suggesting that balancing selection for similar phenotypes in different populations of this species. Nevertheless, a comparison of pair-wise estimates of population differentiation in neutral genetic markers and quantitative traits revealed a strong positive correlation (r = 0.66) suggesting that, for C. odorata, neutral marker divergence could be used as a surrogate for adaptive gene divergence for conservation planning. The utility of this finding and suggested further work are discussed.

Divergent selection and the evolution of signal traits and mating preferences

Mating preferences are common in natural populations, and their divergence among populations is considered an important source of reproductive isolation during speciation. Although mechanisms for the divergence of mating preferences have received substantial theoretical treatment, complementary experimental tests are lacking. We conducted a laboratory evolution experiment, using the fruit fly Drosophila serrata, to explore the role of divergent selection between environments in the evolution of female mating preferences. Replicate populations of D. serrata were derived from a common ancestor and propagated in one of three resource environments: two novel environments and the ancestral laboratory environment. Adaptation to both novel environments involved changes in cuticular hydrocarbons, traits that predict mating success in these populations. Furthermore, female mating preferences for these cuticular hydrocarbons also diverged among populations. A component of this divergence occurred among treatment environments, accounting for at least 17.4% of the among- population divergence in linear mating preferences and 17.2% of the among-population divergence in nonlinear mating preferences. The divergence of mating preferences in correlation with environment is consistent with the classic by- product model of speciation in which premating isolation evolves as a side effect of divergent selection adapting populations to their different environments.

Imitating the neighbours: vocal dialect matching in a mimic-model system

Vocal mimicry provides a unique system for investigating song learning and cultural evolution in birds. Male lyrebirds produce complex vocal displays that include extensive and accurate mimicry of many other bird species. We recorded and analysed the songs of the Albert's lyrebird (Menura alberti) and its most commonly imitated model species, the satin bowerbird (Ptilonorhynchus violaceus), at six sites in southeast Queensland, Australia. We show that each population of lyrebirds faithfully reproduces the song of the local population of bowerbirds. Within a population, lyrebirds show less variation in song structure than the available variation in the songs of the models. These results provide the first quantitative evidence for dialect matching in the songs of two species that have no direct ecological relationship.

Identification, distribution, and probability of infection in sympatric gastropods by digenean parasites at Pa-hay-okee, Everglades National Park using molecules and morphology

In this study, I determined the identity, taxonomic placement, and distribution of digenetic trematodes parasitizing the snails Pomacea paludosa and Planorbella duryi at Pa-hay-okee, Everglades National Park. I also characterized temporal and geographic variation in the probability of parasite infection for these snails based on two years of sampling. Although studies indicate that digenean parasites may have important effects both on individual species and the structure of communities, there have been no studies of digenean parasitism on snails within the Everglades ecosystem. For example, the endangered Everglade Snail Kite, a specialist that feeds almost exclusively on Pomacea paludosa, and is known to be a definitive host of digenean parasites, may suffer direct and indirect effects from consumption of parasitized apple snails. Therefore, information on the diversity and abundance of parasites harbored in snail populations in the Everglades should be of considerable interest for management and conservation of wildlife. Juvenile digeneans (cercariae) representing 20 species were isolated from these two snails, representing a quadrupling of the number of species known. Species were characterized based on morphological, morphometric, and sequence data (18S rDNA, COI, and ITS). Species richness of shed cercariae from P. duryi was greater than P. paludosa, with 13 and 7 species respectively. These species represented 14 families. P. paludosa and P. duryi had no digenean species in common. Probability of digenean infection was higher for P. duryi than P. paludosa and adults showed a greater risk of infection than juveniles for both of these snails. Planorbella duryi showed variation in probability of infection between sampling sites and hydrological seasons. The number of unique combinations of multi-species infections was greatest among P. duryi individuals, while the overall percentage of multi-species infections was greatest in P. paludosa. Analyses of six frequently-observed multiple infections from P. duryi suggest the presence of negative interactions, positive interactions, and neutral associations between larval digeneans. These results should contribute to an understanding of the factors controlling the abundance and distribution of key species in the Everglades ecosystem and may in particular help in the management and recovery planning for the Everglade Snail Kite.

(Tables S1,2) Stable isotopes in muscle tissue of polar bears (Ursus maritimus) from various regions of the Arctic

The relative contribution of regional contamination versus dietary differences to geographic variation in polar bear (Ursus maritimus) contaminant levels is unknown. Dietary variation between Alaska, Canada, East Greenland, and Svalbard subpopulations was assessed by muscle nitrogen and carbon stable isotope (d15N, d13C) and adipose fatty acid (FA) signatures relative to their main prey (ringed seals). Western and southern Hudson Bay signatures were characterized by depleted d15N and d13C, lower proportions of C20 and C22 monounsaturated FAs and higher proportions of C18 and longer chain polyunsaturated FAs. East Greenland and Svalbard signatures were reversed relative to Hudson Bay. Alaskan and Canadian Arctic signatures were intermediate. Between-subpopulation dietary differences predominated over interannual, seasonal, sex, or age variation. Among various brominated and chlorinated contaminants, diet signatures significantly explained variation in adipose levels of polybrominated diphenyl ether (PBDE) flame retardants (14-15%) and legacy PCBs (18-21%). However, dietary influence was contaminant class-specific, since only low or nonsignificant proportions of variation in organochlorine pesticide (e.g., chlordane) levels were explained by diet. Hudson Bay diet signatures were associated with lower PCB and PBDE levels, whereas East Greenland and Svalbard signatures were associated with higher levels. Understanding diet/food web factors is important to accurately interpret contaminant trends, particularly in a changing Arctic.

Seawater carbonate chemistry and elemental composition of purple sea urchin (Strongylocentrotus purpuratus) calcite in a laboratory experiment

Ocean acidification will likely have negative impacts on invertebrates producing skeletons composed of calcium carbonate. Skeletal solubility is partly controlled by the incorporation of "foreign" ions (e.g. magnesium) into the crystal lattice of these skeletal structures, a process that is sensitive to a variety of biological and environmental factors. Here we explore effects of life stage, oceanographic region of origin, and changes in the partial pressure of carbon dioxide in seawater (pCO2) on trace elemental composition in the purple sea urchin (Strongylocentrotus purpuratus). We show that, similar to other urchin taxa, adult purple sea urchins have the ability to precipitate skeleton composed of a range of biominerals spanning low- to high-Mg calcites. Mg / Ca and Sr / Ca ratios were substantially lower in adult spines compared to adult tests. On the other hand, trace elemental composition was invariant among adults collected from four oceanographically distinct regions spanning a range of carbonate chemistry conditions (Oregon, Northern California, Central California, and Southern California). Skeletons of newly settled juvenile urchins that originated from adults from the four regions exhibited intermediate Mg / Ca and Sr / Ca between adult spine and test endmembers, indicating that skeleton precipitated during early life stages is more soluble than adult spines and less soluble than adult tests. Mean skeletal Mg / Ca or Sr / Ca of juvenile skeleton did not vary with source region when larvae were reared under present-day, global-average seawater carbonate conditions (400 µatm; pHT = 8.02 ± 0.03 1 SD; Omega calcite = 3.3 ± 0.2 1 SD). However, when reared under elevated pCO2 (900 µatm; pHT = 7.73 ± 0.03; Omega calcite = 1.8 ± 0.1), skeletal Sr / Ca in juveniles exhibited increased variance across the four regions. Although larvae from the northern populations (Oregon, Northern California, Central California) did not exhibit differences in Mg or Sr incorporation under elevated pCO2 (Sr / Ca = 2.10 ± 0.06 mmol/mol; Mg / Ca = 67.4 ± 3.9 mmol/mol), juveniles of Southern California origin partitioned ~8% more Sr into their skeletons when exposed to higher pCO2 (Sr / Ca = 2.26 ± 0.08 vs. 2.09 ± 0.005 mmol/mol 1 SD). Together these results suggest that the diversity of carbonate minerologies present across different skeletal structures and life stages in purple sea urchins does not translate into an equivalent geochemical plasticity of response associated with geographic variation or temporal shifts in seawater properties. Rather, composition of S. purpuratus skeleton precipitated during both early and adult life history stages appears relatively robust to spatial gradients and predicted future changes in carbonate chemistry. An exception to this trend may arise during early life stages, where certain populations of purple sea urchins may alter skeletal mineral precipitation rates and composition beyond a given pCO2 threshold. This potential for geochemical plasticity during early development in contrast to adult stage geochemical resilience adds to the growing body of evidence that ocean acidification can have differing effects across organismal life stages.

Impact of selection and demography on the diffusion of lactase persistence.

The lactase enzyme allows lactose digestion in fresh milk. Its activity strongly decreases after the weaning phase in most humans, but persists at a high frequency in Europe and some nomadic populations. Two hypotheses are usually proposed to explain the particular distribution of the lactase persistence phenotype. The gene-culture coevolution hypothesis supposes a nutritional advantage of lactose digestion in pastoral populations. The calcium assimilation hypothesis suggests that carriers of the lactase persistence allele(s) (LCT*P) are favoured in high-latitude regions, where sunshine is insufficient to allow accurate vitamin-D synthesis. In this work, we test the validity of these two hypotheses on a large worldwide dataset of lactase persistence frequencies by using several complementary approaches. Methodology We first analyse the distribution of lactase persistence in various continents in relation to geographic variation, pastoralism levels, and the genetic patterns observed for other independent polymorphisms. Then we use computer simulations and a large database of archaeological dates for the introduction of domestication to explore the evolution of these frequencies in Europe according to different demographic scenarios and selection intensities. Conclusions Our results show that gene-culture coevolution is a likely hypothesis in Africa as high LCT*P frequencies are preferentially found in pastoral populations. In Europe, we show that population history played an important role in the diffusion of lactase persistence over the continent. Moreover, selection pressure on lactase persistence has been very high in the North-western part of the continent, by contrast to the South-eastern part where genetic drift alone can explain the observed frequencies. This selection pressure increasing with latitude is highly compatible with the calcium assimilation hypothesis while the gene-culture coevolution hypothesis cannot be ruled out if a positively selected lactase gene was carried at the front of the expansion wave during the Neolithic transition in Europe.

Ecophysiology of Adonis distorta, a high-mountain species endemic of the Central Apennines

Morphological, anatomical and physiological plant and leaf traits of A. distorta, an endemic species of the Central Apennines on the Majella Massif, growing at 2,675 m a.s.l, were analyzed. The length of the phenological cycle starts immediately after the snowmelt at the end of May, lasting 128 ± 10 days. The low A. distorta height  (Hmax= 64 ± 4 mm) and total leaf area (TLA= 38 ± 9 cm2) associated to a high leaf mass area (LMA =11.8±0.6 mg cm−2) and a relatively high leaf tissue density (LTD = 124.6±14.3 mg cm−3) seem to be adaptive traits to the stress factors of the environment where it grows. From a physiological point of view, the high A. distorta photosynthetic rates (PN =19.6 ± 2.3 µmol m−2 s−1) and total chlorophyll content (Chla+b = 0.88 ± 0.13 mg g−1) in July are justified by the favorable temperature. PN decreases by 87% in September at the beginning of plant senescence. Photosynthesis and leaf respiration (RD) variations allow A. distorta to maintain a positive carbon balance during the growing season becoming indicative of the efficiency of plant carbon use. The results could be an important tool for conservation programmes of the A. distorta wild populations.

A Quantitative analysis of pulsed signals emitted by wild bottlenose dolphins

Common bottlenose dolphins (Tursiops truncatus), produce a wide variety of vocal emissions for communication and echolocation, of which the pulsed repertoire has been the most difficult to categorize. Packets of high repetition, broadband pulses are still largely reported under a general designation of burst-pulses, and traditional attempts to classify these emissions rely mainly in their aural characteristics and in graphical aspects of spectrograms. Here, we present a quantitative analysis of pulsed signals emitted by wild bottlenose dolphins, in the Sado estuary, Portugal (2011-2014), and test the reliability of a traditional classification approach. Acoustic parameters (minimum frequency, maximum frequency, peak frequency, duration, repetition rate and inter-click-interval) were extracted from 930 pulsed signals, previously categorized using a traditional approach. Discriminant function analysis revealed a high reliability of the traditional classification approach (93.5% of pulsed signals were consistently assigned to their aurally based categories). According to the discriminant function analysis (Wilk's Λ = 0.11, F3, 2.41 = 282.75, P < 0.001), repetition rate is the feature that best enables the discrimination of different pulsed signals (structure coefficient = 0.98). Classification using hierarchical cluster analysis led to a similar categorization pattern: two main signal types with distinct magnitudes of repetition rate were clustered into five groups. The pulsed signals, here described, present significant differences in their time-frequency features, especially repetition rate (P < 0.001), inter-click-interval (P < 0.001) and duration (P < 0.001). We document the occurrence of a distinct signal type-short burst-pulses, and highlight the existence of a diverse repertoire of pulsed vocalizations emitted in graded sequences. The use of quantitative analysis of pulsed signals is essential to improve classifications and to better assess the contexts of emission, geographic variation and the functional significance of pulsed signals.

Genetic Variation among Major Human Geographic Groups Supports a Peculiar Evolutionary Trend in PAX9

A total of 172 persons from nine South Amerindian, three African and one Eskimo populations were studied in relation to the Paired box gene 9 (PAX9) exon 3 (138 base pairs) as well as its 5' and 3' flanking intronic segments (232 bp and 220 bp, respectively) and integrated with the information available for the same genetic region from individuals of different geographical origins. Nine mutations were scored in exon 3 and six in its flanking regions; four of them are new South American tribe-specific singletons. Exon3 nucleotide diversity is several orders of magnitude higher than its intronic regions. Additionally, a set of variants in the PAX9 and 101 other genes related with dentition can define at least some dental morphological differences between Sub-Saharan Africans and non-Africans, probably associated with adaptations after the modern human exodus from Africa. Exon 3 of PAX9 could be a good molecular example of how evolvability works.