Population genomics of eusocial insects: the costs of a vertebrate‐like effective population size

The evolution of reproductive division of labour and social life in social insects has lead to the emergence of several life‐history traits and adaptations typical of larger organisms: social insect colonies can reach masses of several kilograms, they start reproducing only when they are several years old, and can live for decades. These features and the monopolization of reproduction by only one or few individuals in a colony should affect molecular evolution by reducing the effective population size. We tested this prediction by analysing genome‐wide patterns of coding sequence polymorphism and divergence in eusocial vs. noneusocial insects based on newly generated RNA‐seq data. We report very low amounts of genetic polymorphism and an elevated ratio of nonsynonymous to synonymous changes – a marker of the effective population size – in four distinct species of eusocial insects, which were more similar to vertebrates than to solitary insects regarding molecular evolutionary processes. Moreover, the ratio of nonsynonymous to synonymous substitutions was positively correlated with the level of social complexity across ant species. These results are fully consistent with the hypothesis of a reduced effective population size and an increased genetic load in eusocial insects, indicating that the evolution of social life has important consequences at both the genomic and population levels.     

Population genetic structure and adaptation of malaria parasites on the edge of endemic distribution

To determine whether the major human malaria parasite Plasmodium falciparum exhibits fragmented population structure or local adaptation at the northern limit of its African distribution where the dry Sahel zone meets the Sahara, samples were collected from diverse locations within Mauritania over a range of ~1000 km. Microsatellite genotypes were obtained for 203 clinical infection samples from eight locations, and Illumina paired‐end sequences were obtained to yield high coverage genomewide single nucleotide polymorphism (SNP) data for 65 clinical infection samples from four locations. Most infections contained single parasite genotypes, reflecting low rates of transmission and superinfection locally, in contrast to the situation seen in population samples from countries further south. A minority of infections shared related or identical genotypes locally, indicating some repeated transmission of parasite clones without recombination. This caused some multilocus linkage disequilibrium and local divergence, but aside from the effect of repeated genotypes there was minimal differentiation between locations. Several chromosomal regions had elevated integrated haplotype scores (|iHS|) indicating recent selection, including those containing drug resistance genes. A genomewide F_ST scan comparison with previous sequence data from an area in West Africa with higher infection endemicity indicates that regional gene flow prevents genetic isolation, but revealed allele frequency differentiation at three drug resistance loci and an erythrocyte invasion ligand gene. Contrast of extended haplotype signatures revealed none to be unique to Mauritania. Discrete foci of infection on the edge of the Sahara are genetically highly connected to the wider continental parasite population, and local elimination would be difficult to achieve without very substantial reduction in malaria throughout the region.     

Sex, drugs and recombination: the wild life of Aspergillus

Throughout the eukaryotes, sexual reproduction is an almost universal phenomenon. However, within the Kingdom Fungi, this relationship is not so clear‐cut. Fungi exhibit a spectrum of reproductive modes and life‐cycles; amongst the better known species, sexual reproduction is often facultative, can be rare, and in over half of the known Ascomycota (the moulds) is unknown ( Taylor et al. 1999 ). However, over the last decade, it has become apparent that many of these asexual mitosporic taxa undergo cryptic recombination via unobserved mechanisms and that wholly asexual fungi are, in fact, a rarity ( Taylor et al. 1999, 2001 ; Heitman 2010 ). This revolution in our understanding of fungal sexuality has come about in two ways: Firstly, sexual reproduction leaves an imprint on fungal genomes by maintaining genes required for mating and by generating patterns of mutation and recombination restricted to meiotic processes. Secondly, scientists have become better at catching fungi in flagrante delicto. The genus Aspergillus is one such fungus where a combination of population genetics, genomics and taxonomy has been able to intuit the existence of sex, then to catch the fungus in the act and formally describe their sexual stages. So, why are sexy moulds exciting? One species in particular, Aspergillus flavus, is notorious for its ability to produce a diverse array of secondary metabolites, of which the polyketide aflatoxins (AF) are carcinogenic and others (such as cyclopiazonic acid) are toxigenic. Because of the predilection of A. flavus to grow on crops, such as peanuts, corn and cotton, biocontrol is widely used to mitigate infection by pre‐applying nonaflatoxigenic (AF?) strains to competitively exclude the wild‐type AF+ strains. However, the eventual fate in nature of these biocontrol strains is not known. In this issue of Molecular Ecology, Olarte et al. (2012) make an important contribution by using laboratory crosses of A. flavus to show that not only is AF highly heritable, but AF? strains can become AF+ via crossing over during meiosis. This observation has raised the spectre of cross‐breeding and non‐mendelian inheritance of AF between native and biocontrol strains of the fungus, leading to an increase in the natural diversity of the fungus with perhaps unanticipated consequences.     

Fine‐scale temporal adaptation within a salmonid population: mechanism and consequences

We demonstrate a clear example of local adaptation of seasonal timing of spawning and embryo development. The consequence is a population of pink salmon that is segmented into spawning groups that use the same limited habitat. We synthesize published observations with results of new analyses to demonstrate that genetic variation of these traits results in survival differentials related to that variation, and that density‐dependent embryo mortality and seasonally variable juvenile mortality are a mechanism of selection. Most examples of local adaptation in natural systems depend on observed correlations between environments and fitness traits, but do not fully demonstrate local adaptation: that the trait is genetically determined, exhibits different fitness in common environments or across different environments, and its variation is mechanistically connected to fitness differences. The geographic or temporal scales of local adaptation often remain obscure. Here, we show that heritable, fine‐scale differences of timing of reproductive migration in a pink salmon (Oncorhynchus gorbuscha) resulted in temporal structure that persisted several generations; the differences enable a density‐dependent population to pack more spawners into limited spawning habitat, that is, enhance its fitness. A balanced trade‐off of survivals results because embryos from early‐migrating fish have a lower freshwater survival (harsh early physical conditions and disturbance by late spawners), but emigrant fry from late‐migrating fish have lower marine survivals (timing of their vernal emergence into the estuarine environment). Such fine‐scale local adaptations increase the genetic portfolio of the populations and may provide a buffer against the impacts of climate change.     

Frequency and inheritance of non‐male producing clones in Daphnia magna: evolution towards sex specialization in a cyclical parthenogen?

In Daphnia (Cladocera, Crustacea), parthenogenetic reproduction alternates with sexual reproduction. Individuals of both sexes that belong to the same parthenogenetic line are genetically identical, and their sex is determined by the environment. Previously, non-male producing (NMP) genotypes have been described in species of the Daphnia pulex group. Such genotypes can only persist through phases of sexual reproduction if they co-occur with normal (MP) genotypes that produce both males and females, and thus the breeding system polymorphism is similar to gynodioecy (coexistence of females with hermaphrodites), which is well known in plants. Here we show that the same breeding system polymorphism also occurs in Daphnia magna, a species that has diverged from D. pulex more than 100 MY ago. Depending on the population, between 0% and 40% of D. magna females do not produce males when experimentally exposed to a concentration of the putative sex hormone methyl farnesoate that normally leads to male-only clutches. Natural broods of these NMP females never contained males, contrasting with high proportions of male offspring in MP females from the same populations. The results from a series of crossing experiments suggest that NMP is determined by a dominant allele at a single nuclear locus (or a several closely linked loci): NMP × MP crosses always yielded 50% NMP and 50% MP offspring, whereas MP × MP crosses always yielded 100% MP offspring. Based on cytochrome c oxidase subunit I-sequences, we found that NMP genotypes from different populations belong to three highly divergent mitochondrial lineages, potentially representing three independent evolutionary origins of NMP in D. magna. Thus, the evolution of NMP genotypes in cyclical parthenogens may be more common than previously thought. Moreover, MP genotypes that coexist with NMP genotypes may have responded to the presence of the latter by partially specializing on male production. Hence, these populations of D. magna may be a model for an evolutionary transition from a purely environmental to a partially genetic sex determination system.     

Flowering time QTL in natural populations of Arabidopsis thaliana and implications for their adaptive value

The genetic basis of phenotypic traits is of great interest to evolutionary biologists, but their contribution to adaptation in nature is often unknown. To determine the genetic architecture of flowering time in ecologically relevant conditions, we used a recombinant inbred line population created from two locally adapted populations of Arabidopsis thaliana from Sweden and Italy. Using these RILs, we identified flowering time QTL in growth chambers that mimicked the natural temperature and photoperiod variation across the growing season in each native environment. We also compared the genomic locations of flowering time QTL to those of fitness (total fruit number) QTL from a previous three‐year field study. Ten total flowering time QTL were found, and in all cases, the Italy genotype caused early flowering regardless of the conditions. Two QTL were consistent across chamber environments, and these had the largest effects on flowering time. Five of the fitness QTL colocalized with flowering time QTL found in the Italy conditions, and in each case, the local genotype was favoured. In contrast, just two flowering time QTL found in the Sweden conditions colocalized with fitness QTL and in only one case was the local genotype favoured. This implies that flowering time may be more important for adaptation in Italy than Sweden. Two candidate genes (FLC and VIN3) underlying the major flowering time QTL found in the current study are implicated in local adaptation.     

Angiosperm to Gymnosperm host‐plant switch entails shifts in microbiota of the Welwitschia bug,Probergrothius angolensis (Distant, 1902)

The adaptation of herbivorous insects to new host plants is key to their evolutionary success in diverse environments. Many insects are associated with mutualistic gut bacteria that contribute to the host's nutrition and can thereby facilitate dietary switching in polyphagous insects. However, how gut microbial communities differ between populations of the same species that feed on different host plants remains poorly understood. Most species of Pyrrhocoridae (Hemiptera: Heteroptera) are specialist seed‐feeders on plants in the family Malvaceae, although populations of one species, Probergrothius angolensis, have switched to the very distantly related Welwitschia mirabilis plant in the Namib Desert. We first compared the development and survival of laboratory populations of Pr. angolensis with two other pyrrhocorids on seeds of Welwitschia and found only Pr. angolensis was capable of successfully completing its development. We then collected Pr. angolensis in Namibia from Malvaceae and Welwitschia host plants, respectively, to assess their bacterial and fungal community profiles using high‐throughput amplicon sequencing. Comparison with long‐term laboratory‐reared insects indicated stable associations of Pr. angolensis with core bacteria (Commensalibacter, Enterococcus, Bartonella and Klebsiella), but not with fungi or yeasts. Phylogenetic analyses of core bacteria revealed relationships to other insect‐associated bacteria, but also found new taxa indicating potential host‐specialized nutritional roles. Importantly, the microbial community profiles of bugs feeding on Welwitschia versus Malvaceae revealed stark and consistent differences in the relative abundance of core bacterial taxa that correlate with the host‐plant switch; we were able to reproduce this result through feeding experiments. Thus, a dynamic gut microbiota may provide a means for insect adaptation to new host plants in new environments when food plants are extremely divergent.     

Longevity is associated with relative brain size in birds

Brain size of vertebrates has long been recognized to evolve in close association with basic life‐history traits, including lifespan. According to the cognitive buffer hypothesis, large brains facilitate the construction of behavioral responses against novel socioecological challenges through general cognitive processes, which should reduce mortality and increase lifespan. While the occurrence of brain size–lifespan correlation has been well documented in mammals, much less evidence exists for a robust link between brain size and longevity in birds. The aim of this study was to use phylogenetically controlled comparative approach to test for the relationship between brain size and longevity among 384 avian species from 23 orders. We used maximum lifespan and maximum reproductive lifespan as the measures of longevity and accounted for a set of possible confounding effects, such as allometry, sampling effort, geographic patterns, and life‐history components (clutch size, incubation length, and mode of development). We found that both measures of longevity positively correlated with relative (residual) brain size. We also showed that major diversification of brain size preceded diversification of longevity in avian evolution. In contrast to previous findings, the effect of brain size on longevity was consistent across lineages with different development patterns, although the relatively low strength of this correlation could likely be attributed to the ubiquity of allomaternal care associated with the altricial mode of development. Our study indicates that the positive relationship between brain size and longevity in birds may be more general than previously thought.     

The evolution of reproductive strategies in dasyurid marsupials: implications of molecular phylogeny

Dasyurid marsupials show a remarkable diversity of reproductive patterns ranging from aseasonal polyoestry to restricted annual breeding in which males synchronously die after a brief mating season. Previous studies have categorized dasyurid reproduction into six strategies, defined on the basis of five life-history characters. We provide an up-to-date summary of reproductive traits in dasyurid species and examine the evolution of these characters on a phylogeny for the family recently obtained from DNA sequence data. Our results suggest that reproductive evolution in modern dasyurids is characterized by a basal separation of subfamily lineages employing Strategy II (monoestrous females, restricted breeding season, 11 months to maturity; Dasyurinae) and Strategy V (polyoestrous females, extended breeding season, 8–11 months to maturity; Sminthopsinae). Strategies I (male die-off) and III (facultative polyoestry) appear to have arisen several times from Strategy II or V ancestors, and Strategy IV appears to have arisen within Sminthopsis from a Strategy V ancestor. Strategy VI (aseasonal breeding) has arisen independently in each of the four major dasyurid lineages (tribes), and is highly (but not perfectly) correlated with New Guinean endemism. This scenario is not strongly affected if reproductive characters are optimized on an alternative phylogeny more consistent with morphology-based opinions on species relationships. When evaluated in light of current habitat associations and geographic distributions, the reproductive data suggest that the Miocene diversification of modern dasyurids may have been correlated with the invasion of dry forest or woodland habitats.     

Genetic basis of octanoic acid resistance in Drosophila sechellia: functional analysis of a fine‐mapped region

Drosophila sechellia is a species of fruit fly endemic to the Seychelles islands. Unlike its generalist sister species, D. sechellia has evolved to be a specialist on the host plant Morinda citrifolia. This specialization is interesting because the plant's fruit contains secondary defence compounds, primarily octanoic acid (OA), that are lethal to most other Drosophilids. Although ecological and behavioural adaptations to this toxic fruit are known, the genetic basis for evolutionary changes in OA resistance is not. Prior work showed that a genomic region on chromosome 3R containing 18 genes has the greatest contribution to differences in OA resistance between D. sechellia and D. simulans. To determine which gene(s) in this region might be involved in the evolutionary change in OA resistance, we knocked down expression of each gene in this region in D. melanogaster with RNA interference (RNAi) (i) ubiquitously throughout development, (ii) during only the adult stage and (iii) within specific tissues. We identified three neighbouring genes in the Osiris family, Osiris 6 (Osi6), Osi7 and Osi8, that led to decreased OA resistance when ubiquitously knocked down. Tissue‐specific RNAi, however, showed that decreasing expression of Osi6 and Osi7 specifically in the fat body and/or salivary glands increased OA resistance. Gene expression analyses of Osi6 and Osi7 revealed that while standing levels of expression are higher in D. sechellia, Osi6 expression is significantly downregulated in salivary glands in response to OA exposure, suggesting that evolved tissue‐specific environmental plasticity of Osi6 expression may be responsible for OA resistance in D. sechellia.     

A new multiplex SNP genotyping assay for detecting hybridization and introgression between the M and S molecular forms of Anopheles gambiae

The M and S forms of Anopheles gambiae have been the subject of intense study, but are morphologically indistinguishable and can only be identified using molecular techniques. PCR‐based assays to distinguish the two forms have been designed and applied widely. However, the application of these assays towards identifying hybrids between the two forms, and backcrossed hybrids in particular, has been problematic as the currently available diagnostic assays are based on single locus and/or are located within a multicopy gene. Here, we present an alternative genotyping method for detecting hybridization and introgression between M and S molecular forms based on a multilocus panel of single‐nucleotide polymorphisms (SNPs) fixed between the M and S forms. The panel of SNPs employed is located in so‐called islands of divergence leading us to describe this method as the ‘Divergence Island SNP’ (DIS) assay. We show this multilocus SNP genotyping approach can robustly and accurately detect F1 hybrids as well as backcrossed individuals.     

The costs of being dark: the genetic basis of melanism and its association with fitness‐related traits in the sand cricket

Melanism is an important component of insect cuticle and serves numerous functions that enhance fitness. Despite its importance, there is little information on its genetic basis or its phenotypic and genetic correlation with fitness‐related traits. Here, we examine the heritability of melanism in the wing dimorphic sand cricket and determine its phenotypic and genetic correlation with wing morphology, gonad mass and size of the dorso‐longitudinal muscles (the principle flight muscles). Previously demonstrated trade‐offs among these traits are significant factors in the evolution of life history variation. Using path analysis, we show that melanization is causally related to gonad mass, but not flight muscle mass. Averaged over the sexes, the heritability of melanism was 0.61, the genetic correlation with gonad mass was ?0.36 and with wing morph was 0.51. The path model correctly predicted the ranking of melanization score in lines selected for increased ovary mass, increased flight muscle mass, an index that increased both traits and an unselected control. Our results support the general hypothesis that melanization is costly for insects and negatively impacts investment in early reproduction.     

对我国分子生态学研究近期发展战略的一些思考

分子生态学是多学科交叉的整合性研究领域, 是运用进化生物学理论解决宏观生物学问题的科学。经过半个多世纪的发展, 本学科已日趋成熟, 它不仅已经广泛渗透到宏观生物学的众多学科领域, 而且已经成为连接和融合很多不同学科的桥梁, 是目前最具活力的研究领域之一。其研究的范畴, 从最基础的理论和方法技术, 到格局和模式的发现和描述, 到对过程和机制的深入探讨, 再到付诸于实践的行动和规划指导等各个层次。分子生态学的兴起给宏观生物学带来了若干飞跃性的变化, 使宏观生物学由传统的以观察、测量和推理为主的描述性研究转变为以从生物和种群的遗传构成的变化和历史演化背景上检验、证明科学假设及揭示机制和规律为主的机制性/解释性研究, 因而使得对具有普遍意义的科学规律、生态和进化过程及机制的探索成为可能。分子生态学已经进入组学研究时代, 这使得阐明复杂生态过程、生物地理过程和适应性演化过程的机制性研究由原来难以企及的梦想变成完全可以实现的探求; 它也带来了全新的挑战, 其中最有深远影响的将是对分子生态学研究至关重要的进化生物学基础理论方面的突破, 例如遗传变异理论、种群分化理论、表观遗传因素的作用, 乃至进化生物学的基本知识构架等等。这些方面的进展必将使宏观生物学迎来一场空前的革命, 并对生态学的所有分支学科产生重大影响, 甚至催生诸如生态表观组学这样的新分支学科。对于中国科学家来说, 分子生态学组学时代的开启, 更是一个千载难逢的机遇, 为提出和建立生命科学的新方法、新假说、新思想和新理论提供了莫大的探索空间——此前我们对宏观生物学方法、理论和思想的发展贡献很小。然而, 限制组学时代重大突破的关键因素是理论、概念、理念、实验方法或分析方法方面的创新和突破, 这正是我国分子生态学研究最薄弱的环节。我国教育部门应尽快调整生命科学本科生培养的理念和方法, 以培养具备突出创新潜力的年轻一代后备人才; 同时, 科研项目资助部门和研究人员不仅应清醒地认识本学科领域的发展态势, 更要及时调整思路, 树立新的项目管理理念和治学 理念。     

Predictive ecology: systems approaches

The world is experiencing significant, largely anthropogenically induced, environmental change. This will impact on the biological world and we need to be able to forecast its effects. In order to produce such forecasts, ecology needs to become more predictive--to develop the ability to understand how ecological systems will behave in future, changed, conditions. Further development of process-based models is required to allow such predictions to be made. Critical to the development of such models will be achieving a balance between the brute-force approach that naively attempts to include everything, and over simplification that throws out important heterogeneities at various levels. Central to this will be the recognition that individuals are the elementary particles of all ecological systems. As such it will be necessary to understand the effect of evolution on ecological systems, particularly when exposed to environmental change. However, insights from evolutionary biology will help the development of models even when data may be sparse. Process-based models are more common, and are used for forecasting, in other disciplines, e.g. climatology and molecular systems biology. Tools and techniques developed in these endeavours can be appropriated into ecological modelling, but it will also be necessary to develop the science of ecoinformatics along with approaches specific to ecological problems. The impetus for this effort should come from the demand coming from society to understand the effects of environmental change on the world and what might be performed to mitigate or adapt to them.     

Disease and colony establishment in the dampwood termite Zootermopsis angusticollis: survival and fitness consequences of infection in primary reproductives

Pathogens have likely infl uenced life-history evolution in social insects because their nesting ecology and sociality can exacerbate the risk of disease transmission and place demands on the immune system that ultimately can impact colony survival and growth. The costs of the maintenance and induction of immune function may be particularly significant in termites, which have a nitrogen-poor diet. We examined the effect of fungal exposure on survival and reproduction during colony foundation in the dampwood termite Zootermopsis angusticollis by experimentally pairing male and female primary reproductives and exposing them to single (‘acute’) and multiple (‘serial’) dosages of conidia of the fungus Metarhizium anisopliae and recording their survival and fitness over a 560 day period. The number of eggs laid 70 days post-pairing was significantly reduced relative to controls in the serial-exposure but not the acute-exposure treatment. Reproduction thus appeared to be more resilient to a single pathogen exposure than to serial challenge to the immune system. The impact of fungal exposure was transient: all surviving colonies had similar reproductive output after 300 days post-pairing. Our results suggest that disease can have significant survival and fitness costs during the critical phase of colony foundation but that infection at this time may not necessarily impact long-term colony growth. Received 25 February 2005; revised 27 September and 20 October 2005; accepted 20 December 2005.     

Evolutionarily accelerated invasions: the rate of dispersal evolves upwards during the range advance of cane toads

Human activities are changing habitats and climates and causing species' ranges to shift. Range expansion brings into play a set of powerful evolutionary forces at the expanding range edge that act to increase dispersal rates. One likely consequence of these forces is accelerating rates of range advance because of evolved increases in dispersal on the range edge. In northern Australia, cane toads have increased their rate of spread fivefold in the last 70 years. Our breeding trials with toads from populations spanning the species' invasion history in Australia suggest a genetic basis to dispersal rates and interpopulation genetic variation in such rates. Toads whose parents were from the expanding range front dispersed faster than toads whose parents were from the core of the range. This difference reflects patterns found in their field-collected mothers and fathers and points to heritable variance in the traits that have accelerated the toads' rate of invasion across tropical Australia over recent decades. Taken together with demonstrated spatial assortment by dispersal ability occurring on the expanding front, these results point firmly to ongoing evolution as a driving force in the accelerated expansion of toads across northern Australia.     

Life-history evolution in guppies viii: the demographics of density regulation in guppies (poecilia reticulata)

In prior research, we found the way guppy life histories evolve in response to living in environments with a high or low risk of predation is consistent with life-history theory that assumes no density dependence. We later found that guppies from high-predation environments experience higher mortality rates than those from low-predation environments, but the increased risk was evenly distributed across all age/size classes. Life-history theory that assumes density-independent population growth predicts that life histories will not evolve under such circumstances, yet we have shown with field introduction experiments that they do evolve. However, theory that incorporates density regulation predicts this pattern of mortality can result in the patterns of life-history evolution we had observed. Here we report on density manipulation experiments performed in populations of guppies from low-predation environments to ask whether natural populations normally experience density regulation and, if so, to characterize the short-term demographic changes that underlie density regulation. Our experiments reveal that these populations are density regulated. Decreased density resulted in higher juvenile growth, decreased juvenile mortality rates, and increased reproductive investment by adult females. Increased density causes reduced offspring size, decreased fat storage by adult females, and increased adult mortality.     

Habitat quality ranking depends on habitat-independent environmental factors: a model and results from Callosobruchus maculatus

1. Survival and reproduction in different habitats of a heterogeneous environment are affected by intrinsic habitat properties as well as by habitat-independent extrinsic factors. It is argued that the ranking of habitats according to quality may depend on extrinsic environmental factors that act with the same intensity in all habitats. This can happen if, across habitats, high reproductive success in the absence of extrinsic mortality is coupled with late reproduction.
2. In a model of a semelparous organism with overlapping generations that uses two habitats, the effect of different patterns of extrinsic mortality on habitat quality are analysed, as measured by the reproductive value of an egg, which estimates its expected contribution to the future gene pool.
3. Measurements of fitness components in Callosobruchus maculatus (F.) (Coleoptera: Bruchidae) on four host species demonstrate that, across hosts, long development can be associated with high reproductive success in the absence of extrinsic mortality. As a result, the quality ranking of three out of the four hosts would be different in a stable population regulated by intraspecific competition and egg parasitoids than in a growing population or a stable population regulated by age-independent larval mortality.
4. It is suggested that habitat-independent mortality may play a role in the evolution of ecological niches.