The Evolution of Modularity in the Mammalian Skull II: Evolutionary Consequences

Changes in patterns and magnitudes of integration may influence the ability of a species to respond to selection. Consequently, modularity has often been linked to the concept of evolvability, but their relationship has rarely been tested empirically. One possible explanation is the lack of analytical tools to compare patterns and magnitudes of integration among diverse groups that explicitly relate these aspects to the quantitative genetics framework. We apply such framework here using the multivariate response to selection equation to simulate the evolutionary behavior of several mammalian orders in terms of their flexibility, evolvability and constraints in the skull. We interpreted these simulation results in light of the integration patterns and magnitudes of the same mammalian groups, described in a companion paper. We found that larger magnitudes of integration were associated with a blur of the modules in the skull and to larger portions of the total variation explained by size variation, which in turn can exert a strong evolutionary constraint, thus decreasing the evolutionary flexibility. Conversely, lower overall magnitudes of integration were associated with distinct modules in the skull, to smaller fraction of the total variation associated with size and, consequently, to weaker constraints and more evolutionary flexibility. Flexibility and constraints are, therefore, two sides of the same coin and we found them to be quite variable among mammals. Neither the overall magnitude of morphological integration, the modularity itself, nor its consequences in terms of constraints and flexibility, were associated with absolute size of the organisms, but were strongly associated with the proportion of the total variation in skull morphology captured by size. Therefore, the history of the mammalian skull is marked by a trade-off between modularity and evolvability. Our data provide evidence that, despite the stasis in integration patterns, the plasticity in the magnitude of integration in the skull had important consequences in terms of evolutionary flexibility of the mammalian lineages.

Craniometric Similarities Within and Between Human Populations in Comparison with Neutral Genetic Data

The statement that pairs of individuals from different populations are often more genetically similar than pairs from the same population is a widespread idea inside and outside the scientific community. Witherspoon et al. [""Genetic similarities within and between human populations,"" Genetics 176:351-359 (2007)] proposed an index called the dissimilarity fraction (omega) to access in a quantitative way the validity of this statement for genetic systems. Witherspoon demonstrated that, as the number of loci increases, omega decreases to a point where, when enough sampling is available, the statement is false. In this study, we applied the dissimilarity fraction to Howells`s craniometric database to establish whether or not similar results are obtained for cranial morphological traits. Although in genetic studies thousands of loci are available, Howells`s database provides no more than 55 metric traits, making the contribution of each variable important. To cope with this limitation, we developed a routine that takes this effect into consideration when calculating. omega Contrary to what was observed for the genetic data, our results show that cranial morphology asymptotically approaches a mean omega of 0.3 and therefore supports the initial statement-that is, that individuals from the same geographic region do not form clear and discrete clusters-further questioning the idea of the existence of discrete biological clusters in the human species. Finally, by assuming that cranial morphology is under an additive polygenetic model, we can say that the population history signal of human craniometric traits presents the same resolution as a neutral genetic system dependent on no more than 20 loci.

Evolution of Dendrocolaptes platyrostris (Aves: Furnariidae) between the South American open vegetation corridor and the Atlantic forest

The open vegetation corridor of South America is a region dominated by savanna biomes. It contains forests (i.e. riverine forests) that may act as corridors for rainforest specialists between the open vegetation corridor and its neighbouring biomes (i.e. the Amazonian and Atlantic forests). A prediction for this scenario is that populations of rainforest specialists in the open vegetation corridor and in the forested biomes show no significant genetic divergence. We addressed this hypothesis by studying plumage and genetic variation of the Planalto woodcreeper Dendrocolaptes platyrostris Spix (1824) (Aves: Furnariidae), a forest specialist that occurs in both open habitat and in the Atlantic forest. The study questions were: (1) is there any evidence of genetic continuity between populations of the open habitat and the Atlantic forest and (2) is plumage variation congruent with patterns of neutral genetic structure or with ecological factors related to habitat type? We used cytochrome b and mitochondrial DNA control region sequences to show that D. platyrostris is monophyletic and presents substantial intraspecific differentiation. We found two areas of plumage stability: one associated with Cerrado and the other associated with southern Atlantic Forest. Multiple Mantel tests showed that most of the plumage variation followed the transition of habitats but not phylogeographical gaps, suggesting that selection may be related to the evolution of the plumage of the species. The results were not compatible with the idea that forest specialists in the open vegetation corridor and in the Atlantic forest are linked at the population level because birds from each region were not part of the same genetic unit. Divergence in the presence of gene flow across the ecotone between both regions might explain our results. Also, our findings indicate that the southern Atlantic forest may have been significantly affected by Pleistocene climatic alteration, although such events did not cause local extinction of most taxa, as occurred in other regions of the globe where forests were significantly affected by global glaciations. Finally, our results neither support plumage stability areas, nor subspecies as full species. (C) 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 103, 801-820.

Worldwide sequence conservation of transmission-blocking vaccine candidate Pvs230 in Plasmodium vivax

Pfs230, surface protein of gametocyte/gamete of the human malaria parasite, Plasmodium falciparum, is a prime candidate of malaria transmission-blocking vaccine. Plasmodium vivax has an ortholog of Pfs230 (Pvs230), however, there has been no study in any aspects on Pvs230 to date. To investigate whether Pvs230 can be a vivax malaria transmission-blocking vaccine, we performed evolutionary and population genetic analysis of the Pvs230 gene (pvs230: PVX_003905). Our analysis of Pvs230 and its orthologs in eight Plasmodium species revealed two distinctive parts: an interspecies variable part (IVP) containing species-specific oligopeptide repeats at the N-terminus and a 7.5 kb interspecies conserved part (ICP) containing 14 cysteine-rich domains. Pvs230 was closely related to its orthologs, Pks230 and Pcys230, in monkey malaria parasites. Analysis of 113 pvs230 sequences obtained from worldwide, showed that nucleotide diversity is remarkably low in the non-repeat 8-kb region of pvs230 (theta pi = 0.00118) with 77 polymorphic nucleotide sites, 40 of which results in amino acid replacements. A signature of purifying selection but not of balancing selection was seen on pvs230. Functional and/or structural constraints may limit the level of polymorphism in pvs230. The observed limited polymorphism in pvs230 should ground for utilization of Pvs230 as an effective transmission-blocking vaccine. (C) 2011 Elsevier Ltd. All rights reserved.

Evaluation of petrogenetic models for intermediate and silicic plutonic rocks from the Sierra de Valle Fertil-La Huerta, Argentina: Petrologic constraints on the origin of igneous rocks in the Ordovician Famatinian-Puna paleoarc

The whole Valle Fertil-La Huerta section appears as a calc-alkaline plutonic suite typical of a destructive plate margin. New Sr and Nd isotopic whole-rock data and published whole-rock geochemistry suggest that the less-evolved intermediate (dioritic) rocks can be derived by magmatic differentiation, mainly by hornblende + plagioclase +/- Fe-Ti oxide fractional crystallization, from mafic (gabbroic) igneous precursors. Closed-system differentiation, however, cannot produce the typical intermediate (tonalitic) and silicic (granodioritic) plutonic rocks, which requires a preponderant contribution of crustal components. Intermediate and silicic plutonic rocks from Valle Fertil-La Huerta section have formed in a plate subduction setting where the thermal and material input of mantle-derived magmas promoted fusion of fertile metasedimentary rocks and favored mixing of gabbroic or dioritic magmas with crustal granitic melts. Magma mixing is observable in the field and evident in variations of chemical elemental parameters and isotopic ratios, revealing that hybridization coupled with fractionation of magmas took place in the crust. Consideration of the whole-rock geochemical and isotopic data in the context of the Famatinian-Puna magmatic belt as a whole demonstrates that the petrologic model postulated for the Sierra Valle Fertil-La Huerta section has the potential to explain the generation of plutonic and volcanic rocks across the Early Ordovician paleoarc from central and northwestern Argentina. As the petrologic model does not require the intervention of old Precambrian continental crust, the nature of the basement on which thick accretionary turbiditic sequences were deposited remains a puzzling aspect. Discussion in this paper provides insights into the nature of magmatic source rocks and mechanisms of magma generation in Cordilleran-type volcano-plutonic arcs of destructive plate margins. (C) 2010 Elsevier Ltd. All rights reserved.

A review of the tectonic evolution of the Sunsas belt, SW Amazonian Craton

The Sunsas-Aguapei province (1.20-0.95 Ga), SW Amazonian Craton, is a key area to study the heterogeneous effects of collisional events with Laurentia, which shows evidence of the Grenvillian and Sunsas orogens. The Sunsas orogen, characterized by an allochthonous collisional-type belt (1.11-1.00 Ga), is the youngest and southwestern most of the events recorded along the cratonic fringe. Its evolution occurred after a period of long quiescence and erosion of the already cratonized provinces (>1.30 Ga), that led to sedimentation of the Sunsas and Vibosi groups in a passive margin setting. The passive margin stage was roughly contemporary with intraplate tectonics that produced the Nova Brasilandia proto-oceanic basin (<1.21 Ga), the reactivation of the Ji-Parana shear zone network (1.18-1.12 Ga) and a system of aborted rifts that evolved to the Huanchaca-Aguapei basin (1.17-1.15 Ga). The Sunsas belt is comprised by the metamorphosed Sunsas and Vibosi sequences, the Rincon del Tigre mafic-ultramafic sill and granitic intrusive suites. The latter rocks yield epsilon(Nd(t)) signatures (-0.5 to -4.5) and geochemistry (S,1, A-types) suggesting their origin associated with a continental arc setting. The Sunsas belt evolution is marked by ""tectonic fronts"" with sinistral offsets that was active from c. 1.08 to 1.05 Ga, along the southern edge of the Paragua microcontinent where K/Ar ages (1.27-1.34 Ga) and the Huanchaca-Aguapei flat-lying cover attest to the earliest tectonic stability at the time of the orogen. The Sunsas dynamics is coeval with inboard crustal shortening, transpression and magmatism in the Nova Brasilandia belt (1.13-1.00 Ga). Conversely, the Aguapei aulacogen (0.96-0.91 Ga) and nearby shear zones (0.93-0.91 Ga) are the late tectonic offshoots over the cratonic margin. The post-tectonic to anorogenic stages took place after ca. 1.00 Ga, evidenced by the occurrences of intra-plate A-type granites, pegmatites, mafic dikes and sills, as well as of graben basins. Integrated interpretation of the available data related to the Sunsas orogen supports the idea that the main nucleus of Rodinia incorporated the terrains forming the SW corner of Amazonia and most of the Grenvillian margin, as a result of two independent collisional events, as indicated in the Amazon region by the Ji-Parana shear zone event and the Sunsas belt, respectively. (C) 2009 Elsevier Ltd. All rights reserved.

Age and geochemistry of mantle peridotites and diorite dykes from the Baldissero body: Insights into the Paleozoic-Mesozoic evolution of the Southern Alps

Trace element and isotopic data obtained for mantle spinel Iherzolites and diorite dykes from the Baldissero massif (Ivrea-Verbano Zone, Western Italy) provide new, valuable constraints on the petrologic and geodynamic evolution of the Southern Alps in Paleozoic to Mesozoic times. Whole rock and mineral chemistry indicates that Baldissero Iherzolites can be regarded as refractory mantle residues following limited melt extraction. In particular, the Light Rare Earth Elements (LREE)-depleted and fractionated compositions of whole rock and clinopyroxene closely match modelling results for refractory residues after low degrees (similar to 4-5%) of near-fractional melting of depleted mantle, possibly under garnet-facies conditions. Following this, the peridotite sequence experienced subsolidus re-equilibration at lithospheric spinel-facies conditions and intrusion of several generations of dykes. However, Iherzolites far from dykes show very modest metasomatic changes, as evidenced by the crystallisation of accessory titanian pargasite and the occurrence of very slight enrichments in highly incompatible trace elements (e.g. Nb). The Re-Os data for Iherzolites far from the dykes yield a 376 Ma (Upper Devonian) model age that is considered to record a partial melting event related to the Variscan orogenic cycle s.l. Dioritic dykes cutting the mantle sequence have whole rock, clinopyroxene and plagioclase characterised by high radiogenic Nd and low radiogenic Sr, which point to a depleted to slightly enriched mantle source. Whole rock and mafic phases of diorites have high Mg# values that positively correlate with the incompatible trace element concentrations. The peridotite at the dyke contact is enriched in orthopyroxene, iron and incompatible trace elements with respect to the Iherzolites far from dykes. Numerical simulations indicate that the geochemical characteristics of the diorites can be explained by flow of a hydrous, silica-saturated melt accompanied by reaction with the ambient peridotite and fractional crystallisation. The composition of the more primitive melts calculated in equilibrium with the diorite minerals show tholeiitic to transitional affinity. Internal Sm-Nd, three-point isochrons obtained for two dykes suggest an Upper Triassic-Lower Jurassic emplacement age (from 204 31 to 198 29 Ma). Mesozoic igneous events are unknown in the southern Ivrea-Verbano Zone (IVZ), but the intrusion of hydrous melts, mostly silica-saturated, have been well documented in the Finero region, i.e. the northernmost part of IVZ and Triassic magmatism with calc-alkaline to shoshonitic affinity is abundant throughout the Central-Eastern Alps. The geochemical and chronological features of the Baldissero diorites shed new light on the geodynamic evolution of the Southern Alps before the opening of the Jurassic Tethys. (C) 2010 Elsevier B.V. All rights reserved.

ENO adaptive method for solving one-dimensional conservation laws

In this work an efficient third order non-linear finite difference scheme for solving adaptively hyperbolic systems of one-dimensional conservation laws is developed. The method is based oil applying to the solution of the differential equation an interpolating wavelet transform at each time step, generating a multilevel representation for the solution, which is thresholded and a sparse point representation is generated. The numerical fluxes obtained by a Lax-Friedrichs flux splitting are evaluated oil the sparse grid by an essentially non-oscillatory (ENO) approximation, which chooses the locally smoothest stencil among all the possibilities for each point of the sparse grid. The time evolution of the differential operator is done on this sparse representation by a total variation diminishing (TVD) Runge-Kutta method. Four classical examples of initial value problems for the Euler equations of gas dynamics are accurately solved and their sparse solutions are analyzed with respect to the threshold parameters, confirming the efficiency of the wavelet transform as an adaptive grid generation technique. (C) 2008 IMACS. Published by Elsevier B.V. All rights reserved.

Does liquidity constraints matter for Ricardian equivalence?

The paper analysis a general equilibrium model with two periods, several households and a government that has to finance some expenditures in the first period. Households may have some private information either about their type (adverse selection) or about some action levei chosen in the first period that affects the probability of certain states of nature in the second period (moral hazard). Trade of financiai assets are intermediated by a finite collection of banks. Banks objective functions are determined in equilibrium by shareholders. Due to private information it may be optimal for the banks to introduce constraints in the set of available portfolios for each household as wellas household specific asset prices. In particular, households may face distinct interest rates for holding the risk-free asset. The government finances its expenditures either by taxing households in the first period or by issuing bonds in the first period and taxing households in the second period. Taxes may be state-dependent. Suppose government policies are neutml: i) government policies do not affect the distribution of wealth across households; and ii) if the government decides to tax a household in the second period there is a portfolio available for the banks that generates the Mme payoff in each state of nature as the household taxes. Tben, Ricardian equivalence holds if and only if an appropriate boundary condition is satisfied. Moreover, at every free-entry equilibrium the boundary condition is satisfied and thus Ricardian equivalence holds. These results do not require any particular assumption on the banks' objective function. In particular, we do not assume banks to be risk neutral.

Elucidating the Evolution of the Red Brocket Deer Mazama americana Complex (Artiodactyla; Cervidae)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

ORIGIN and EVOLUTION of DOMESTIC CATTLE

The origin and evolution of domestic cattle have recently moved to the forefront of the scientific literature in consideration of their links to human history and to decisions on Genetic Resources conservation strategies. DNA from modern and ancient Bos samples is being analysed to reconstruct, in cooperation with archaezoology, the main events and forces that shaped nowadays cattle genetic diversity. Still, a number of open questions remain, that hopefully will be answered with the help of new technologies and the combined analysis of worldwide data.

Reproductive phenology of a northeast Brazilian mangrove community: Environmental and biotic constraints

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

From Amazonia to the Atlantic forest: Molecular phylogeny of Phyzelaphryninae frogs reveals unexpected diversity and a striking biogeographic pattern emphasizing conservation challenges

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Anomaly-free constraints in neutrino seesaw models

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)