The genome of eukaryotic species is a complex structure that experience many different interactions within itself and with the surrounding environment. The genetic architecture of a phenotype (that is, the set of genetic elements affecting a trait of the organism) plays a fundamental role in understanding the adaptation process of a species to, for example, different climate environments, or to its interaction with other species. Thus, it is fundamental to study the different aspects of the genetic architecture of the species and its relationship with its surronding environment. Aspects such as modularity (the number of genetic units and the degree to which each unit is affecting a trait of the organism), pleiotropy (the number of different effects that a genetic unit can have on an organism) or linkage (the degree of association between the different genetic units) are essential to understand the genetic architecture and to interpret the effects of selection on the genome. Indeed, the knowledge of the different aspects of the genetic architecture could clarify whether genes are affected by multiple aspects of the environment or, on the contrary, are affected by only specific aspects [1,2].

The work performed by Lotterhos et al. [3] sought to understand the genetic architecture of the adaptation to different environments in lodgepole pine (Pinus contorta), considering as candidate SNPs those previously detected as a result of its extreme association patterns to different environmental variables or to extreme population differentiation. This consideration is very important because the study is only relevant if the studied markers are under the effect of selection. Otherwise, the genetic architecture of the adaptation to different environments would be masked by other (neutral) kind of associations that would be difficult to interpret [4,5]. In order to understand the relationship between genetic architecture and adaptation, it is relevant to detect the association networks of the candidate SNPs with climate variables (a way to measure modularity) and if these SNPs (and loci) are affected by single or multiple environments (a way to measure pleiotropy).

The authors used co-association networks, an innovative approach in this field, to analyse the interaction between the environmental information and the genetic polymorphism of each individual. This methodology is more appropriate than other multivariate methods - such as analysis based on principal components - because it is possible to cluster SNPs based on associations with similar environmental variables. In this sense, the co-association networks allowed to both study the genetic and physical linkage between different co-associations modules but also to compare two different models of evolution: a Modular environmental response architecture (specific genes are affected by specific aspects of the environment) or a Universal pleiotropic environmental response architecture (all genes are affected by all aspects of the environment). The representation of different correlations between allelic frequency and environmental factors (named galaxy biplots) are especially informative to understand the effect of the different clusters on specific aspects of the environment (for example, the co-association network ‘Aridity’ shows strong associations with hot/wet versus cold/dry environments).

The analysis performed by Lotterhos et al. [3], although it has some unavoidable limitations (e.g., only extreme candidate SNPs are selected, limiting the results to the stronger effects; the genetic and physical map is incomplete in this species), includes relevant results and also implements new methodologies in the field. To highlight some of them: the preponderance of a Modular environmental response architecture (evolution in separated modules), the detection of physical linkage among SNPs that are co-associated with different aspects of the environment (which was unexpected a priori), the implementation of co-association networks and galaxy biplots to see the effect of modularity and pleiotropy on different aspects of environment. Finally, this work contains remarkable introductory Figures and Tables explaining unambiguously the main concepts [6] included in this study. This work can be treated as a starting point for many other future studies in the field.

References

[1] Hancock AM, Brachi B, Faure N, Horton MW, Jarymowycz LB, Sperone FG, Toomajian C, Roux F & Bergelson J. 2011. Adaptation to climate across the Arabidopsis thaliana genome. Science 334: 83–86. doi: 10.1126/science.1209244
[2] Wagner GP & Zhang J. The pleiotropic structure of the genotype­phenotype map: the evolvability of complex organisms. Nature Review Genetics 12: 204–213. doi: 10.1038/nrg2949
[3] Lotterhos KE, Yeaman S, Degner J, Aitken S, Hodgins K. 2018. Modularity of genes involved in local adaptation to climate despite physical linkage. bioRxiv 202481, ver. 4 peer-reviewed by Peer Community In Evolutionary Biology. doi: 10.1101/202481
[4] Lotterhos KE & Whitlock MC. 2014. Evaluation of demographic history and neutral parameterization on the performance of FST outlier tests. Molecular Ecology 23: 2178–2192. doi: 10.1111/mec.12725
[5] Lotterhos KE & Whitlock MC. 2015. The relative power of genome scans to detect local adaptation depends on sampling design and statistical method. Molecular Ecology 24: 1031–1046. doi: 10.1111/mec.13100
[6] Paaby AB & Rockman MV. 2013. The many faces of pleiotropy. Trends in Genetics 29: 66-73. doi: 10.1016/j.tig.2012.10.010

Inbreeding, or mating between relatives, generally lowers fitness [1]. Mating between genetically similar individuals can result in higher levels of homozygosity and consequently a higher frequency with which recessive disease alleles may be expressed within a population. Reduced fitness as a consequence of inbreeding, or inbreeding depression, can vary between individuals, sexes, populations and species [2], but remains a pervasive challenge for many organisms with small local population sizes, including humans [3]. But all is not lost for individuals within small populations, because an array of mechanisms can be employed to evade the negative effects of inbreeding [4], including sib-mating avoidance and dispersal [5, 6].

Despite thorough investigation of inbreeding and sib-mating avoidance in the laboratory, only very few studies have ventured into the field besides studies on vertebrates and eusocial insects. The study of Collet et al. [7] is a surprising exception, where the effect of male density and frequency of relatives on inbreeding avoidance was tested in the laboratory, after which robust field collections and microsatellite genotyping were used to infer relatedness and dispersal in natural populations. The parasitic wasp Venturia canescens is an excellent model system to study inbreeding, because mating success was previously found to decrease with increasing relatedness between mates in the laboratory [8] and this species thus suffers from inbreeding depression [9–11]. The authors used an elegant design combining population genetics and model simulations to estimate relatedness of mating partners in the field and compared that with a theoretical distribution of potential mate encounters when random mating is assumed. One of the most important findings of this study is that mating between siblings is not avoided in this species in the wild, despite negative fitness effects when inbreeding does occur. Similar findings were obtained for another insect species, the field cricket Gryllus campestris [12], which leaves us to wonder whether inbreeding tolerance could be more common in nature than currently appreciated.

The authors further looked into sex-specific dispersal patterns between two patches located a few hundred meters apart. Females were indeed shown to be more related within a patch, but no genetic differences were observed between males, suggesting that V. canescens males more readily disperse. Moreover, microsatellite data at 18 different loci did not reveal genetic differentiation between populations approximately 300 kilometers apart. Gene flow is thus occurring over considerable distances, which could play an important role in the ability of this species to avoid negative fitness consequences of inbreeding in nature.

Another interesting aspect of this work is that discrepancies were found between laboratory- and field-based data. What is the relevance of laboratory-based experiments if they cannot predict what is happening in the wild? Many, if not most, biologists (including us) bring our model system into the laboratory to control, at least to some extent, the plethora of environmental factors that could potentially affect our system (in ways that we do not want). Most behavioral studies on mating patterns and sexual selection are conducted in standardized laboratory conditions, but sexual selection is in essence social selection, because an individual’s fitness is partly determined by the phenotype of its social partners (i.e. the social environment) [13]. The social environment may actually dictate the expression of female mate choice and it is unclear how potential laboratory-induced social biases affect mating outcome. In V. canescens, findings using field-caught individuals paint a completely opposite picture of what was previously shown in the laboratory, i.e. sib-avoidance is not taking place in the field. It is likely that density, level of relatedness, sex ratio in the field, and/or the size of experimental arenas in the lab are all factors affecting mate selectivity, as we have previously shown in a butterfly [14–16]. If females, for example, typically only encounter a few males in sequence in the wild, it may be problematic for them to express choosiness when confronted simultaneously with two or more males in the laboratory. A recent study showed that, in the wild, female moths take advantage of staying in groups to blur male choosiness [17]. It is becoming more and more clear that what we observe in the laboratory may not actually reflect what is happening in nature [18]. Instead of ignoring the species-specific life history and ecological features of our favorite species when conducting lab experiments, we suggest that it is time to accept that we now have the theoretical foundations to tease apart what in this “environmental noise” actually shapes sexual selection in nature. Explicitly including ecology in studies on sexual selection will allow us to make more meaningful conclusions, i.e. rather than “this is what may happen in the wild”, we would be able to state “this is what often happens in nature”.

References

[1] Charlesworth D & Willis JH. 2009. The genetics of inbreeding depression. Nat. Rev. Genet. 10: 783–796. doi: 10.1038/nrg2664
[2] Hedrick PW & Garcia-dorado A. 2016. Understanding inbreeding depression, purging, and genetic rescue. Trends Ecol. Evol. 31: 940–952. doi: 10.1016/j.tree.2016.09.005
[3] Bittles AH & Black ML. 2010. Consanguinity, human evolution, and complex diseases. Proc. Natl. Acad. Sci. United States Am. 107: 1779–1786. doi: 10.1073/pnas.0906079106
[4] Pusey A & Wolf M. 1996. Inbreeding avoidance in animals. Trends Ecol. Evol. 11: 201–206. doi: 10.1016/0169-5347(96)10028-8
[5] Greenwood PJ & Harvey PH. 1978. Inbreeding and dispersal in the great tit. Nature 271: 52–54. doi: 10.1038/271052a0
[6] Szulkin M & Sheldon BC. 2008. Dispersal as a means of inbreeding avoidance in a wild bird population. Proc. R. Soc. B 275: 703–711. doi: 10.1098/rspb.2007.0989
[7] Collet M, Amat I, Sauzet S, Auguste A, Fauvergue X, Mouton L, Desouhant E. 2018. Insects and incest: sib-mating tolerance in natural populations of a parasitoid wasp. bioRxiv 169268, ver. 4 peer-reviewed by Peer Community In Evolutionary Biology. doi: 10.1101/169268
[8] Metzger M, Bernstein C, Hoffmeister TS & Desouhant E. 2010. Does kin recognition and sib-mating avoidance limit the risk of genetic incompatibility in a parasitic wasp ? PLoS One 5: e13505. doi: 10.1371/journal.pone.0013505
[9] Beukeboom LW. 2001. Single-locus complementary sex determination in the Ichneumonid Venturia canescens. Netherlands J. Zool. 51: 1–15. doi: 10.1163/156854201X00017
[10] Vayssade C, de Fazio C, Quaglietti B, Auguste A, Ris N, Fauvergue X. 2014. Inbreeding depression in a parasitoid wasp with single- locus complementary sex determination. PLoS One 9: 1–8. doi: 10.1371/journal.pone.0097733
[11] Chuine A, Sauzet S, Debias F & Desouhant E. 2015. Consequences of genetic incompatibility on fitness and mate choice: the male point of view. Biol. J. Linn. Soc. 114: 279–286. doi: 10.1111/bij.12421
[12] Bretman A, Rodri R & Tregenza T. 2011. Fine-scale population structure , inbreeding risk and avoidance in a wild insect population. Mol. Ecol. 20: 3045–3055. doi: 10.1111/j.1365-294X.2011.05140.x
[13] West-Eberhard MJ. 2014. Darwin’s forgotten idea: The social essence of sexual selection. Neurosci. Biobehav. Rev. 46: 501–508. doi: 10.1016/j.neubiorev.2014.06.015
[14] Holveck M-J, Gauthier A-L & Nieberding CM 2015. Dense, small and male-biased cages exacerbate male-male competition and reduce female choosiness in Bicyclus anynana. Anim. Behav. 104: 229–245. doi: 10.1016/j.anbehav.2015.03.025
[15] Nieberding, CM & Holveck M-J 2017. Laboratory social environment biases mating outcome: a first quantitative synthesis in a butterfly. Behav. Ecol. Sociobiol. 71: 117. doi: 10.1007/s00265-017-2346-9
[16] Nieberding CM & Holveck M-J. (In prep). Comentary on Kehl et al. 2018: "Young male mating success is associated with sperm number but not with male sex pheromone titres". Front. Ecol. Evol.
[17] Wijk M Van, Heath J, Lievers R, Schal C & Groot AT. 2017. Proximity of signallers can maintain sexual signal variation under stabilizing selection. Sci. Rep. 7: 18101. doi: 10.1038/s41598-017-17327-9
[18] Miller CW & Svensson EI. 2014. Sexual selection in complex environments. Annu. Rev. Entomol. 59: 427–445. doi: 10.1146/annurev-ento-011613-162044

Perhaps because of its deadliness, the 2013-2016 Ebola Virus (EBOV) epidemics in West-Africa has led to unprecedented publication and sharing of full virus genome sequences. This was both rapid (90 full genomes were shared within weeks [1]) and important (more than 1500 full genomes have been released overall [2]). Furthermore, the availability of the metadata (especially GPS location) has led to depth analyses of the geographical spread of the epidemics [3].
In their work, Dellicour et al. [4] pursue earlier phylogeographical investigations in an original and yet simple approach to address questions of key public health importance. The originality of the approach is dual. First, from a technical standpoint, they capture the spread of infectious diseases in a continuous framework using a novel model that allows for rare long-distance dispersal events. Second, in a more classical discrete meta-population framework, they simulate the effect of public health interventions by pruning the phylogenetic tree and assessing how this affects key parameters. For instance, to simulate the effect of closing borders they remove subsets of the phylogeny that involved dispersal between countries and to simulate the effect of protecting a region by quarantine they remove all the leaves (i.e. the infections sampled) from this region. This phylogeny pruning is both original and simple. It is however limited because it currently assumes that policies are 100% effective and earlier modelling work on human influenza showed that long distance travel bans had to be implemented with >99% efficiency in order to slow epidemic growth from a time scale of days to weeks [5].
From a biological standpoint, Dellicour et al. [4] corroborate earlier findings that highly populated locations (>1,000,000 inhabitants) were crucial in explaining the magnitude of the epidemics but also show the importance of the transmission between the three capital cities. They also show that rare long-distance dispersing events of the virus are not key to explaining the magnitude of the epidemics (even though they assume 100% efficiency of suppressing long-distance event). Finally, thanks to their continuous model they estimate the speed of spread of the epidemics and are able to detect the effect of border closing on this speed.
Overall, this study [4], which involves state-of-the-art Bayesian inference methods of infection phylogenies using MCMC, stands out because of its effort to simulate public health interventions. It stands as an encouragement for the development of intervention models with increased realism and for even faster and larger virus sequence data sharing.

References

[1] Gire et al. 2014. Genomic surveillance elucidates Ebola virus origin and transmission during the 2014 outbreak. Science 345: 1369–1372. doi: 10.1126/science.1259657.
[2] Holmes EC, Dudas G, Rambaut A and Andersen KG. 2016. The evolution of Ebola virus: insights from the 2013-2016 epidemic. Nature 538: 193–200. doi: 10.1038/nature19790.
[3] Dudas et al. 2017. Virus genomes reveal factors that spread and sustained the Ebola epidemic. Nature 544: 309–315 (2017). doi: 10.1038/nature22040.
[4] Dellicour S, Baele G, Dudas G, Faria NR, Pybus OG, Suchard MA, Rambaud A and Lemey P. 2018. Phylodynamic assessment of intervention strategies for the West African Ebola virus outbreak. bioRxiv, 163691, ver. 3 peer-reviewed by Peer Community In Evolutionary Biology. doi: 10.1101/163691.
[5] Hollingsworth TD, Ferguson NM and Anderson RM. 2006. Will travel restrictions control the international spread of pandemic influenza? Nature Medicine 12, 497–499. doi: 10.1038/nm0506-497.

The influence of environmental change over geological time on the tempo and mode of biological diversification is a hot topic in biogeography. Of central interest are questions about where, when, and how fast lineages proliferated, suffered extinction, and migrated in response to tectonic events, the waxing and waning of dominant biomes, etc. In this context, the dynamic conditions of the Miocene have received much attention, from studies of many clades and biogeographic regions. Here, Chazot et al. [1] present an exemplary analysis of butterflies (tribe Ithomiini) in the Neotropics, examining their diversification across the Andes and Amazon. They infer sharp contrasts between these regions in the late Miocene: accelerated diversification during orogeny of the Andes, and greater extinction in the Amazon associated during the Pebas system, with interchange and local diversification increasing following the Pebas during the Pliocene.
Two features of this study stand out. First is the impressive taxon sampling (340 out of 393 extant species). Second is the use of ancestral range reconstructions to compute per-lineage rates of colonization between regions, and rates of speciation within regions, through time. The latter allows for relatively fine-grained comparisons across the 2 fundamental dimensions of historical biogeography, space and time, and is key to the main results. The method resonated with me because I performed a similar analysis in a study showing evidence for uplift-driven diversification in the Hengduan Mountains of China [2]. This analysis is complemented by a variety of other comparative methods for inferring variable diversification across clades, through time, and in response to external factors. Overall, it represents a very nice contribution to our understanding of the effects of Miocene/Pliocene environmental change on the evolution of Neotropical biodiversity.

References

[1] Chazot N, Willmott KR, Lamas G, Freitas AVL, Piron-Prunier F, Arias CF, Mallet J, De-Silva DL and Elias M. 2017. Renewed diversification following Miocene landscape turnover in a Neotropical butterfly radiation. BioRxiv 148189, ver 4 of 19th December 2017. doi: 10.1101/148189

[2] Xing Y, and Ree RH. 2017. Uplift-driven diversification in the Hengduan Mountains, a temperate biodiversity hotspot. Proceedings of the National Academy of Sciences of the United States of America, 114: E3444-E3451. doi: 10.1073/pnas.1616063114

Recent methodological developments and increased genome sequencing efforts have introduced the tantalizing possibility of inferring ancestral phenotypes using DNA from contemporary species. One intriguing application of this idea is to exploit the apparent correlation between substitution rates and body size to infer ancestral species' body sizes using the inferred patterns of substitution rate variation among species lineages based on genomes of extant species [1].
The recommended paper by Figuet et al. [2] examines the utility of such approaches by analyzing the Cetartiodactyla, a clade of large mammals that have mostly well resolved phylogenetic relationships and a reasonably good fossil record. This combination of genomic data and fossils allows a direct comparison between body size predictions obtained from the genomic data and empirical evidence from the fossil record. If predictions seem good in groups such as the Cetartiodactyla, where there is independent evidence from the fossil record, this would increase the credibility of predictions made for species with less abundant fossils.
Figuet et al. [2] analyze transcriptome data for 41 species and report a significant effect of body mass on overall substitution rate, synonymous vs. non-synonymous rates, and the dynamics of GC-content, thus allowing a prediction of small ancestral body size in this group despite the fact that the extant species that were analyzed are nearly all large.
A comparative method based solely on morphology and phylogenetic relationships would be very unlikely to make such a prediction. There are many sources of uncertainty in the variables and parameters associated with these types of approaches: phylogenetic uncertainty (topology and branch lengths), uncertainty about inferred substitution rates, and so on. Although the authors do not account for all these sources of uncertainty the fact that their predicted body sizes appear sensible is encouraging and undoubtedly the methods will become more statistically sophisticated over time.

References

[1] Romiguier J, Ranwez V, Douzery EJP and Galtier N. 2013. Genomic evidence for large, long-lived ancestors to placental mammals. Molecular Biology and Evolution 30: 5–13. doi: 10.1093/molbev/mss211

[2] Figuet E, Ballenghien M, Lartillot N and Galtier N. 2017. Reconstruction of body mass evolution in the Cetartiodactyla and mammals using phylogenomic data. bioRxiv, ver. 3 of 4th December 2017. 139147. doi: 10.1101/139147