Transcriptome analyses of primitively eusocial wasps reveal novel insights into the evolution of sociality and the origin of alternative phenotypes

BACKGROUND: Understanding how alternative phenotypes arise from the same genome is a major challenge in modern biology. Eusociality in insects requires the evolution of two alternative phenotypes - workers, who sacrifice personal reproduction, and queens, who realize that reproduction. Extensive work on honeybees and ants has revealed the molecular basis of derived queen and worker phenotypes in highly eusocial lineages, but we lack equivalent deep-level analyses of wasps and of primitively eusocial species, the latter of which can reveal how phenotypic decoupling first occurs in the early stages of eusocial evolution. RESULTS: We sequenced 20 Gbp of transcriptomes derived from brains of different behavioral castes of the primitively eusocial tropical paper wasp Polistes canadensis. Surprisingly, 75% of the 2,442 genes differentially expressed between phenotypes were novel, having no significant homology with described sequences. Moreover, 90% of these novel genes were significantly upregulated in workers relative to queens. Differential expression of novel genes in the early stages of sociality may be important in facilitating the evolution of worker behavioral complexity in eusocial evolution. We also found surprisingly low correlation in the identity and direction of expression of differentially expressed genes across similar phenotypes in different social lineages, supporting the idea that social evolution in different lineages requires substantial de novo rewiring of molecular pathways. CONCLUSIONS: These genomic resources for aculeate wasps and first transcriptome-wide insights into the origin of castes bring us closer to a more general understanding of eusocial evolution and how phenotypic diversity arises from the same genome.

Inferring the mode of origin of polyploid species from next-generation sequence data.

Many eukaryote organisms are polyploid. However, despite their importance, evolutionary inference of polyploid origins and modes of inheritance has been limited by a need for analyses of allele segregation at multiple loci using crosses. The increasing availability of sequence data for nonmodel species now allows the application of established approaches for the analysis of genomic data in polyploids. Here, we ask whether approximate Bayesian computation (ABC), applied to realistic traditional and next-generation sequence data, allows correct inference of the evolutionary and demographic history of polyploids. Using simulations, we evaluate the robustness of evolutionary inference by ABC for tetraploid species as a function of the number of individuals and loci sampled, and the presence or absence of an outgroup. We find that ABC adequately retrieves the recent evolutionary history of polyploid species on the basis of both old and new sequencing technologies. The application of ABC to sequence data from diploid and polyploid species of the plant genus Capsella confirms its utility. Our analysis strongly supports an allopolyploid origin of C. bursa-pastoris about 80 000 years ago. This conclusion runs contrary to previous findings based on the same data set but using an alternative approach and is in agreement with recent findings based on whole-genome sequencing. Our results indicate that ABC is a promising and powerful method for revealing the evolution of polyploid species, without the need to attribute alleles to a homeologous chromosome pair. The approach can readily be extended to more complex scenarios involving higher ploidy levels.

High-precision dating and origin of synsedimentary volcanism in the Late Carboniferous Salvan-Dorenaz basin (Aiguilles-Rouges Massif, Western Alps)

Late Variscan volcanic activity is documented in the Late Carboniferous Salvan-Dorenaz sedimentary basin and in the neighboring basement units of the Aiguilles-Rouges and Mont-Blanc crystalline massifs (Western Alps). Precise U/Pb isotopic dating, zircon morphology and geochemical analyses indicate that volcanism occurred during short-lived pulses and that coexisting crustal and mantle sources were involved in the production of melts. Volcanic and subvolcanic products were emplaced along major N-S to NNE-SSW transtensional fracture zones, similar to the ones that governed intense basement exhumation and that favored the formation and filling of the Late Carboniferous Salvan-Dorenaz continental basin. In the Aiguilles-Rouges massif, dacitic flows outcropping at the base of the Salvan-Dorenaz basin erupted at 308 +/- 3 Ma; they represent the surface equivalent of the nearby Vallorcine peraluminous granite and associated rhyolitic dykes (311 +/- 17 Ma). In the Mont Blanc massif, calc-alkaline rhyolitic dykes were emplaced simultaneously (307 +/- 2 Ma) at shallow crustal levels, but they derive from deeper magma sources denoting enhanced mantellic activity. Recently identified tuffs and volcaniclastic layers embedded at different levels of the Salvan-Dorenaz stratigraphic record testify a 295 +3/-4 Ma old episode of highly explosive volcanism from distant volcanic centers, possibly located in the Aar-Gotthard massifs (Central Alps). Their zircon typology is highly heterogeneous. documenting wall-rock contamination of the melts and/or admixture of crustal sediments, whereas consistent subpopulations point to high-temperature magmas of deep-seated origin and alkaline affinity. The dated volcanic layers from the Salvan-Dorenaz basin set the beginning of the detrital sedimentation at 308 +/- 3 Ma and constrain the deposition of 1.5-1.7 km thick of elastic sediments within a time span of 10-15 Ma. These results infer minimum, long-term subsidence rates during basin evolution in the order of >0.1 mm/a, while in the surrounding basement units estimated exhumation rates are in the range of 1 mm/a. All dated rocks contain inherited zircon populations about 350, 450 or 600 Ma old.

Chromosomal evolution and zoogeographic origin of southeast Asian shrews (genus Crocidura).

To the origins and evolution of Indomalayan shrews, we investigated the chromosomal variations of 14 species of Crocidura from SE Asia. Intraspecific polymorphism was mainly due to variation in the number of short chromosomal arms but C. lepidura and C. hutanis showed a polymorphism due to a centric fusion. The undifferentially stained karyotypes were similar in 9 species, all possessing 2n = 38 and FN = 54-56 (68); C. fuliginosa had 2n = 40 and FN = 54-58. These karyotypes are close to the presumed ancestral state for the genus Crocidura. Four species from Sulawesi had a reduced diploid number (2n = 30-34), a trend not observed among other SE Asian species but present in few Palaearctic taxa. Compared to the apparent stasis of karyotypic evolution observed among other SE Asian species, the high degree of interspecific differences reported among Sulawesian shrews is unusual and needs further investigation. Stasis and reduction in diploid number found in both Indomalayan and Palaeractic species suggest that these two groups share a common ancestry. This is in sharp contrast to most Afrotropical species which evolved towards higher diploid and fundamental numbers. The zoogeographical implications of these results are discussed.

Origin of cadmium enrichments in carbonate rocks deposited in the Alpine Tethys area during the Middle-Late Jurassic

ABSTRACTThe pollution of air, soil and water by heavy metals through anthropogenic activities is an object of numerous environmental studies since long times. A number of natural processes, such as volcanic activity, hydrothermal fluid circulation and weathering of metal-rich deposits may lead to an additional and potentially important input and accumulation of heavy metals in the environment. In the Swiss and French Jura Mountains, anomalous high cadmium (Cd) concentrations (up to 16 ppm) in certain soils are related to the presence of underlying Cd-enriched (up to 21 ppm) carbonate rocks of Middle to Late Jurassic age. The aim of this study is to understand the processes controlling Cd incorporation into carbonate rocks of Middle and Late Jurassic age and to reconstruct the sedimentary and environmental conditions, which have led to Cd enrichments in these sedimentary rocks.Cd concentrations in studied hemipelagic sections in France vary between 0.1 and 0.5 ppm (mean 0.15 ppm). Trace-element behavior and high Mn concentrations suggest that sediment accumulation occurred in a well-oxygenated environment. Increases in Cd contents in the bulk-rock carbonate sediments may be related to increases in surface-water productivity under oxic conditions and important remineralization of organic matter within the water column. In platform settings preserved in the Swiss Jura Mountains, no correlation is observed between Cd contents and evolution of environmental conditions. Cd concentrations in these platform sections are often below the detection limit, with isolated peaks of up to 21 ppm. These important Cd enrichments are associated with peaks in Zn concentrations and are present in carbonate rocks independently of facies and age. The high Cd contents in these shallow-water carbonate rocks are partly related to the presence of disseminated, Cd-rich (up to 1.8%), sphalerite (ZnS) mineralization. The basement rocks are considered to be the source of metals for sulfide mineralization in the overlying Jurassic strata, as the sphalerite Pb isotope pattern is comparable to that of granite rocks from the nearby southern Black Forest crystalline massif. The Rb-Sr ages of sphalerite samples indicate that a main phase of sphalerite formation occurred near the boundary between the late Middle and early Late Jurassic, at around 162 Ma, as a result of enhanced tectonic and hydrothermal activity in Europe, related to the opening of the Central Atlantic and to the tectonic/thermal subsidence during spreading of the Alpine Tethys. I therefore propose to use unusually high Cd concentrations in carbonates as a tracer of tectonic activity in the Jura Mountains area, especially in the case when important enrichments in Zn co-occur. Paleoproductivity reconstructions based on records of authigenic Cd may be compromised not only by post-depositional redistribution, but also by incorporation of additional Cd from hydrothermal solutions circulating in the rock.The circulation of metal-rich hydrothermal fluids through the sediment sequence, in addition to specific environmental conditions during sedimentation, contributes to the incorporation of Cd into the carbonate rocks. However, only hydrothermal activity has led to the unusually high concentrations of Cd in carbonate rocks of Bajocian-Oxfordian age, through the formation of sphalerite mineralization.

Loss of egg yolk genes in mammals and the origin of lactation and placentation

Embryonic development in nonmammalian vertebrates depends entirely on nutritional reserves that are predominantly derived from vitellogenin proteins and stored in egg yolk. Mammals have evolved new resources, such as lactation and placentation, to nourish their developing and early offspring. However, the evolutionary timing and molecular events associated with this major phenotypic transition are not known. By means of sensitive comparative genomics analyses and evolutionary simulations, we here show that the three ancestral vitellogenin-encoding genes were progressively lost during mammalian evolution (until around 30-70 million years ago, Mya) in all but the egg-laying monotremes, which have retained a functional vitellogenin gene. Our analyses also provide evidence that the major milk resource genes, caseins, which have similar functional properties as vitellogenins, appeared in the common mammalian ancestor approximately 200-310 Mya. Together, our data are compatible with the hypothesis that the emergence of lactation in the common mammalian ancestor and the development of placentation in eutherian and marsupial mammals allowed for the gradual loss of yolk-dependent nourishment during mammalian evolution

On the origin of the inertia: the modified Newtonian dynamics theory

The sameness between the inertial mass and the gravitational mass is an assumption and not a consequence of the equivalent principle is shown. In the context of the Sciama’s inertia theory, the sameness between the inertial mass and the gravitational mass is discussed and a certain condition which must be experimentally satisfied is given. The inertial force proposed by Sciama, in a simple case, is derived from the Assis’ inertia theory based in the introduction of a Weber type force. The origin of the inertial force is totally justified taking into account that the Weber force is, in fact, an approximation of a simple retarded potential, see [18, 19]. The way how the inertial forces are also derived from some solutions of the general relativistic equations is presented. We wonder if the theory of inertia of Assis is included in the framework of the General Relativity. In the context of the inertia developed in the present paper we establish the relation between the constant acceleration a0 , that appears in the classical Modified Newtonian Dynamics (M0ND) theory, with the Hubble constant H0 , i.e. a0 ≈ cH0 .

Metabolic origin of δ15N values in nitrogenous compounds from Brassica napus L. leaves

Nitrogen isotope composition (δ15N) in plant organic matter is currently used as a natural tracer of nitrogen acquisition efficiency. However, the δ15N value of whole leaf material does not properly reflect the way in which N is assimilated because isotope fractionations along metabolic reactions may cause substantial differences among leaf compounds. In other words, any change in metabolic composition or allocation pattern may cause undesirable variability in leaf δ15N. Here, we investigated the δ15N in different leaf fractions and individual metabolites from rapeseed (Brassica napus) leaves. We show that there were substantial differences in δ15N between nitrogenous compounds (up to 30 ) and the content in (15N enriched) nitrate had a clear influence on leaf δ15N. Using a simple steady-state model of day metabolism, we suggest that the δ15N value in major amino acids was mostly explained by isotope fractionation associated with isotope effects on enzyme-catalysed reactions in primary nitrogen metabolism. δ15N values were further influenced by light versus dark conditions and the probable occurrence of alternative biosynthetic pathways. We conclude that both biochemical pathways (that fractionate between isotopes) and nitrogen sources (used for amino acid production) should be considered when interpreting the δ15N value of leaf nitrogenous compounds

Conserved developmental expression of Fezf in chordates and Drosophila and the origin of the Zona Limitans Intrathalamica (ZLI) brain organizer

Background: The zona limitans intrathalamica (ZLI) and the isthmus organizer (IsO) are two major secondary organizers of vertebrate brain development. These organizers are located at the interface of the expression domains of key patterning genes (Fezf-Irx and Otx-Gbx, respectively). To gain insights into the evolutionary origin of the ZLI, we studied Fezf in bilaterians. Results: In this paper, we identified a conserved sequence motif (Fezf box) in all bilaterians. We report the expression pattern of Fezf in amphioxus and Drosophila and compare it with those of Gbx, Otx and Irx. We found that the relative expression patterns of these genes in vertebrates are fully conserved in amphioxus and flies, indicating that the genetic subdivisions defining the location of both secondary organizers in early vertebrate brain development were probably present in the last common ancestor of extant bilaterians. However, in contrast to vertebrates, we found that Irx-defective flies do not show an affected Fezf expression pattern. Conclusions: The absence of expression of the corresponding morphogens from cells at these conserved genetic boundaries in invertebrates suggests that the organizing properties might have evolved specifically in the vertebrate lineage by the recruitment of key morphogens to these conserved genetic locations.

Conserved developmental expression of Fezf in chordates and Drosophila and the origin of the Zona Limitans Intrathalamica (ZLI) brain organizer

Background: The zona limitans intrathalamica (ZLI) and the isthmus organizer (IsO) are two major secondary organizers of vertebrate brain development. These organizers are located at the interface of the expression domains of key patterning genes (Fezf-Irx and Otx-Gbx, respectively). To gain insights into the evolutionary origin of the ZLI, we studied Fezf in bilaterians. Results: In this paper, we identified a conserved sequence motif (Fezf box) in all bilaterians. We report the expression pattern of Fezf in amphioxus and Drosophila and compare it with those of Gbx, Otx and Irx. We found that the relative expression patterns of these genes in vertebrates are fully conserved in amphioxus and flies, indicating that the genetic subdivisions defining the location of both secondary organizers in early vertebrate brain development were probably present in the last common ancestor of extant bilaterians. However, in contrast to vertebrates, we found that Irx-defective flies do not show an affected Fezf expression pattern. Conclusions: The absence of expression of the corresponding morphogens from cells at these conserved genetic boundaries in invertebrates suggests that the organizing properties might have evolved specifically in the vertebrate lineage by the recruitment of key morphogens to these conserved genetic locations.

Metal dealing at the origin of the Chordata phylum: The metallothionein system and metal overload response in Amphioxus.

Non-vertebrate chordates, specifically amphioxus, are considered of the utmost interest for gaining insight into the evolutionary trends, i.e. differentiation and specialization, of gene/protein systems. In this work, MTs (metallothioneins), the most important metal binding proteins, are characterized for the first time in the cephalochordate subphylum at both gene and protein level, together with the main features defining the amphioxus response to cadmium and copper overload. Two MT genes (BfMT1 and BfMT2) have been identified in a contiguous region of the genome, as well as several ARE (antioxidant response element) and MRE (metal response element) located upstream the transcribed region. Their corresponding cDNAs exhibit identical sequence in the two lancelet species (B. floridae and B. lanceolatum), BfMT2 cDNA resulting from an alternative splicing event. BfMT1 is a polyvalent metal binding peptide that coordinates any of the studied metal ions (Zn, Cd or Cu) rendering complexes stable enough to last in physiological environments, which is fully concordant with the constitutive expression of its gene, and therefore, with a metal homeostasis housekeeping role. On the contrary, BfMT2 exhibits a clear ability to coordinate Cd(II) ions, while it is absolutely unable to fold into stable Cu (I) complexes, even as mixed species. This identifies it as an essential detoxification agent, which is consequently only induced in emergency situations. The cephalochordate MTs are not directly related to vertebrate MTs, neither by gene structure, protein similarity nor metal-binding behavior of the encoded peptides. The closest relative is the echinoderm MT, which confirm proposed phylogenetic relationships between these two groups. The current findings support the existence in most organisms of two types of MTs as for their metal binding preferences, devoted to different biological functions: multivalent MTs for housekeeping roles, and specialized MTs that evolve either as Cd-thioneins or Cu-thioneins, according to the ecophysiological needs of each kind of organisms.

On the origin of the selectivity of plasmidic H-NS towards horizontally acquired DNA: Linking H-NS oligomerization and cooperative DNA binding

The nucleoid-associated protein H-NS is a global modulator of the expression of genes associated with adaptation to environmental changes. A variant of H-NS expressed in the R27 plasmid was previously shown to selectively modulate the expression of horizontally acquired genes, with minimal effects on core genes that are repressed by the chromosomal form of H-NS. Both H-NS proteins are formed by an oligomerization domain and a DNA-binding domain, which are connected by a linker that is highly flexible in the absence of DNA. We studied DNA binding by means of oligomer-forming chimeric proteins in which domains of the chromosomal and plasmidic variants are exchanged, as well as in monomeric truncated forms containing the DNA-binding domain and variable portions of the linker. Point mutations in the linker were also examined in full-length and truncated H-NS constructs. These experiments show that the linker region contributes to DNA binding affinity and that it is a main component of the distinct DNA binding properties of chromosomal and plasmidic H-NS. We propose that interactions between the linker and DNA limit the flexibility of the connection between H- NS oligomerization and DNA binding and provide an allosteric indirect readout mechanism to detect long- range distortions of DNA, thus enabling discrimination between core and horizontally acquired DNA.

Protein loop compaction and the origin of the effect of arginine and glutamic acid mixtures on solubility, stability and transient oligomerization of proteins

Addition of a 50 mM mixture of l-arginine and l-glutamic acid (RE) is extensively used to improve protein solubility and stability, although the origin of the effect is not well understood. We present Small Angle X-ray Scattering (SAXS) and Nuclear Magnetic Resonance (NMR) results showing that RE induces protein compaction by collapsing flexible loops on the protein core. This is suggested to be a general mechanism preventing aggregation and improving resistance to proteases and to originate from the polyelectrolyte nature of RE. Molecular polyelectrolyte mixtures are expected to display long range correlation effects according to dressed interaction site theory. We hypothesize that perturbation of the RE solution by dissolved proteins is proportional to the volume occupied by the protein. As a consequence, loop collapse, minimizing the effective protein volume, is favored in the presence of RE.

Comparative genomics of Hedgehog loci in chordates and the origin of Shh regulatory novelties

The origin and evolution of the complex regulatory landscapes of some vertebrate developmental genes, often spanning hundreds of Kbp and including neighboring genes, remain poorly understood. The Sonic Hedgehog (Shh) genomic regulatory block (GRB) is one of the best functionally characterized examples, with several discrete enhancers reported within its introns, vast upstream gene-free region and neighboring genes (Lmbr1 and Rnf32). To investigate the origin and evolution of this GRB, we sequenced and characterized the Hedgehog (Hh) loci from three invertebrate chordate amphioxus species, which share several early expression domains with Shh. Using phylogenetic footprinting within and between chordate lineages, and reporter assays in zebrafish probing >30 Kbp of amphioxus Hh, we report large sequence and functional divergence between both groups. In addition, we show that the linkage of Shh to Lmbr1 and Rnf32, necessary for the unique gnatostomate-specific Shh limb expression, is a vertebrate novelty occurred between the two whole-genome duplications.