Breeding system and reproductive skew in a highly polygynous ant population

Abstract.: Factors affecting relatedness among nest members in ant colonies with high queen number are still poorly understood. In order to identify the major determinants of nest kin structure, we conducted a detailed analysis of the breeding system of the ant Formica exsecta. We estimated the number of mature queens by mark-release-recapture in 29 nests and dissected a sub-sample of queens to assess their reproductive status. We also used microsatellites to estimate relatedness within and between all classes of nestmates (queens, their mates, worker brood, queen brood and male brood). Queen number was very high, with an arithmetic mean of 253 per nest. Most queens (90%) were reproductively active, consistent with the genetic analyses revealing that there was only a minimal reproductive skew among nestmate queens. Despite the high queen number and low reproductive skew, almost all classes of individuals were significantly related to each other. Interestingly, the number of resident queens was a poor predictor of kin structure at the nest level, consistent with the observation that new queens are produced in bursts leading to highly fluctuating queen number across years. Queen number also varied tremendously across nests, with estimates ranging from five to several hundred queens. Accordingly, the harmonic mean queen number (40.5) was six times lower than the arithmetic mean. The variation in queen number was the most important factor of the breeding system contributing to a significant relatedness between almost all classes of nestmates despite a high average number of queens per nes

Superconductivity in CuxBi2Se3 and its implications for pairing in the undoped topological insulator

Bi2Se3 is one of a handful of known topological insulators. Here we show that copper intercalation in the van der Waals gaps between the Bi2Se3 layers, yielding an electron concentration of ~ 2 x 10^20cm-3, results in superconductivity at 3.8 K in CuxBi2Se3 for x between 0.12 and 0.15. This demonstrates that Cooper pairing is possible in Bi2Se3 at accessible temperatures, with implications for study of the physics of topological insulators and potential devices.Comment: 6 pages, 4 figure

How do species barriers decay? Concordance and local introgression in mosaic hybrid zones of mussels.

The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study 'replicated' instances of secondary contact between closely related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry-informative panel of such SNPs. We then compared their frequencies in newly sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bi-stable variants (Dobzhansky-Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact.ANR, France; Russian Science Foundatio

Maternal age effect and severe germ-line bottleneck in the inheritance of human mitochondrial DNA

The manifestation of mitochondrial DNA (mtDNA) diseases depends on the frequency of heteroplasmy (the presence of several alleles in an individual), yet its transmission across generations cannot be readily predicted owing to a lack of data on the size of the mtDNA bottleneck during oogenesis. For deleterious heteroplasmies, a severe bottleneck may abruptly transform a benign (low) frequency in a mother into a disease-causing (high) frequency in her child. Here we present a high-resolution study of heteroplasmy transmission conducted on blood and buccal mtDNA of 39 healthy mother–child pairs of European ancestry (a total of 156 samples, each sequenced at ∼20,000× per site). On average, each individual carried one heteroplasmy, and one in eight individuals carried a disease-associated heteroplasmy, with minor allele frequency ≥1%. We observed frequent drastic heteroplasmy frequency shifts between generations and estimated the effective size of the germ-line mtDNA bottleneck at only ∼30–35 (interquartile range from 9 to 141). Accounting for heteroplasmies, we estimated the mtDNA germ-line mutation rate at 1.3 × 10−8 (interquartile range from 4.2 × 10−9 to 4.1 × 10−8) mutations per site per year, an order of magnitude higher than for nuclear DNA. Notably, we found a positive association between the number of heteroplasmies in a child and maternal age at fertilization, likely attributable to oocyte aging. This study also took advantage of droplet digital PCR (ddPCR) to validate heteroplasmies and confirm a de novo mutation. Our results can be used to predict the transmission of disease-causing mtDNA variants and illuminate evolutionary dynamics of the mitochondrial genome

Postglacial expansion of the arctic keystone copepod calanus glacialis

Calanus glacialis, a major contributor to zooplankton biomass in the Arctic shelf seas, is a key link between primary production and higher trophic levels that may be sensitive to climate warming. The aim of this study was to explore genetic variation in contemporary populations of this species to infer possible changes during the Quaternary period, and to assess its population structure in both space and time. Calanus glacialis was sampled in the fjords of Spitsbergen (Hornsund and Kongsfjorden) in 2003, 2004, 2006, 2009 and 2012. The sequence of a mitochondrial marker, belonging to the ND5 gene, selected for the study was 1249 base pairs long and distinguished 75 unique haplotypes among 140 individuals that formed three main clades. There was no detectable pattern in the distribution of haplotypes by geographic distance or over time. Interestingly, a Bayesian skyline plot suggested that a 1000-fold increase in population size occurred approximately 10,000 years before present, suggesting a species expansion after the Last Glacial Maximum.GAME from the National Science Centre, the Polish Ministry of Science and Higher Education Iuventus Plus [IP2014 050573]; FCT-PT [CCMAR/Multi/04326/2013]; [2011/03/B/NZ8/02876

Evidence for Pervasive Adaptive Protein Evolution in Wild Mice

The relative contributions of neutral and adaptive substitutions to molecular evolution has been one of the most controversial issues in evolutionary biology for more than 40 years. The analysis of within-species nucleotide polymorphism and between-species divergence data supports a widespread role for adaptive protein evolution in certain taxa. For example, estimates of the proportion of adaptive amino acid substitutions (alpha) are 50% or more in enteric bacteria and Drosophila. In contrast, recent estimates of alpha for hominids have been at most 13%. Here, we estimate alpha for protein sequences of murid rodents based on nucleotide polymorphism data from multiple genes in a population of the house mouse subspecies Mus musculus castaneus, which inhabits the ancestral range of the Mus species complex and nucleotide divergence between M. m. castaneus and M. famulus or the rat. We estimate that 57% of amino acid substitutions in murids have been driven by positive selection. Hominids, therefore, are exceptional in having low apparent levels of adaptive protein evolution. The high frequency of adaptive amino acid substitutions in wild mice is consistent with their large effective population size, leading to effective natural selection at the molecular level. Effective natural selection also manifests itself as a paucity of effectively neutral nonsynonymous mutations in M. m. castaneus compared to humans

Convergent development of low-relatedness supercolonies in Myrmica ants.

Many ant species have independently evolved colony structures with multiple queens and very low relatedness among nestmate workers, but it has remained unclear whether low-relatedness kin structures can repeatedly arise in populations of the same species. Here we report a study of Danish island populations of the red ant Myrmica sulcinodis and show that it is likely that such repeated developments occur. Two microsatellite loci were used to estimate genetic differentiation (F(ST)) among three populations and nestmate relatedness within these populations. The F(ST) values were highly significant due to very different allele frequencies among the three populations with relatively few common alleles and relatively many rare alleles, possibly caused by single queen foundation and rare subsequent immigration. Given the isolation of the islands and the low investment in reproduction, we infer that each of the populations was most likely established by a single queen, even though all three extant populations now have within-colony relatedness 95%), and the genetic differentiation of nests showed a significantly positive correlation with the distance between them. Both male-biased sex-ratio and genetic viscosity are expected characteristics of populations where queens have very local dispersal and where new colonies are initiated through nest-budding. Based on a comparison with other M. sulcinodis populations we hypothesise a distinct succession of population types and suggest that this may be a possible pathway to unicoloniality, ie, development towards a complete lack of colony kin structure and unrelated nestmate workers

Variation in Base-Substitution Mutation in Experimental and Natural Lineages of Caenorhabditis Nematodes

Variation among lineages in the mutation process has the potential to impact diverse biological processes ranging from susceptibilities to genetic disease to the mode and tempo of molecular evolution. The combination of high-throughput DNA sequencing (HTS) with mutation-accumulation (MA) experiments has provided a powerful approach to genome-wide mutation analysis, though insights into mutational variation have been limited by the vast evolutionary distances among the few species analyzed. We performed a HTS analysis of MA lines derived from four Caenorhabditis nematode natural genotypes: C. elegans N2 and PB306 and C. briggsae HK104 and PB800. Total mutation rates did not differ among the four sets of MA lines. A mutational bias toward G:C→A:T transitions and G:C→T:A transversions was observed in all four sets of MA lines. Chromosome-specific rates were mostly stable, though there was some evidence for a slightly elevated X chromosome mutation rate in PB306. Rates were homogeneous among functional coding sequence types and across autosomal cores, arms, and tips. Mutation spectra were similar among the four MA line sets but differed significantly when compared with patterns of natural base-substitution polymorphism for 13/14 comparisons performed. Our findings show that base-substitution mutation processes in these closely related animal lineages are mostly stable but differ from natural polymorphism patterns in these two species

New methods for inferring the distribution of fitness effects for INDELs and SNPs

Small insertions and deletions (INDELs; ≤50bp) are the most common type of variability after SNPs. However, compared to SNPs, we know little about the distribution of fitness effects (DFE) of new INDEL mutations and how prevalent adaptive INDEL substitutions are. Studying INDELs has been difficult partly because identifying ancestral states at these sites is error-prone and misidentification can lead to severely biased estimates of the strength of selection. To solve these problems, we develop new maximum likelihood methods, which use polymorphism data to simultaneously estimate the DFE, the mutation rate, and the misidentification rate. These methods are applicable to both INDELs and SNPs. Simulations show that they can provide highly accurate results. We applied the methods to an INDEL polymorphism dataset in Drosophila melanogaster. We found that the DFE for polymorphic INDELs in protein-coding regions is bimodal, with the variants being either nearly neutral or strongly deleterious. Based on the DFE, we estimated that 71.5% - 83.7% of the INDEL substitutions that took place along the D. melanogaster lineage were fixed by positive selection, which is comparable to the prevalence of adaptive substitutions at non-synonymous sites. The new methods have been implemented in the software package anavar

Mutation Accumulation in a Selfing Population: Consequences of Different Mutation Rates between Selfers and Outcrossers

Currently existing theories predict that because deleterious mutations accumulate at a higher rate, selfing populations suffer from more intense genetic degradation relative to outcrossing populations. This prediction may not always be true when we consider a potential difference in deleterious mutation rate between selfers and outcrossers. By analyzing the evolutionary stability of selfing and outcrossing in an infinite population, we found that the genome-wide deleterious mutation rate would be lower in selfing than in outcrossing organisms. When this difference in mutation rate was included in simulations, we found that in a small population, mutations accumulated more slowly under selfing rather than outcrossing. This result suggests that under frequent and intense bottlenecks, a selfing population may have a lower risk of genetic extinction than an outcrossing population