Neutral Variation Research Articles

Effects of hybridization on pelvic morphology: A macaque model

Ancient DNA analyses have shown that interbreeding between hominin taxa occurred multiple times. Although admixture is often reflected in skeletal phenotype, the relationship between the two remains poorly understood, hampering interpretation of the hominin fossil record. Direct study of this relationship is often impossible due to the paucity of hominin fossils and difficulties retrieving ancient genetic material. Here, we use a sample of known ancestry hybrids between two closely related nonhuman primate taxa (Indian and Chinese Macaca mulatta) to investigate the effect of admixture on skeletal morphology. We focus on pelvic shape, which has potential fitness implications in hybrids, as mismatches between maternal pelvic and fetal cranial morphology are often fatal to mother and offspring. As the pelvis is also one of the skeletal regions that differs most between Homo sapiens and Neanderthals, investigating the pelvic consequences of interbreeding could be informative regarding the viability of their hybrids. We find that the effect of admixture in M. mulatta is small and proportional to the relatively small morphological difference between the parent taxa. Sexual dimorphism appears to be the main determinant of pelvic shape in M. mulatta. The lack of difference in pelvic shape between Chinese and Indian M. mulatta is in contrast to that between Neanderthals and H. sapiens, despite a similar split time (in generations) between the hybridizing pairs. Greater phenotypic divergence between hominins may relate to adaptations to disparate environments but may also highlight how the unique degree of cultural buffering in hominins allowed for greater neutral divergence. In contrast to some previous work identifying extreme morphologies in first- and second-generation hybrids, here the relationship between pelvic shape and admixture is linear. This linearity may be because most sampled animals have a multigenerational admixture history or because of relatively high constraints on the pelvis compared with other skeletal regions.

Genomic time-series data show that gene flow maintains high genetic diversity despite substantial genetic drift in a butterfly species.

Effective population size affects the efficacy of selection, rate of evolution by drift and neutral diversity levels. When species are subdivided into multiple populations connected by gene flow, evolutionary processes can depend on global or local effective population sizes. Theory predicts that high levels of diversity might be maintained by gene flow, even very low levels of gene flow, consistent with species long-term effective population size, but tests of this idea are mostly lacking. Here, we show that Lycaeides butterfly populations maintain low contemporary (variance) effective population sizes (e.g. ~200 individuals) and thus evolve rapidly by genetic drift. However, populations harboured high levels of genetic diversity consistent with an effective population size several orders of magnitude larger. We hypothesized that the differences in the magnitude and variability of contemporary versus long-term effective population sizes were caused by gene flow of sufficient magnitude to maintain diversity but only subtly affect evolution on generational timescales. Consistent with this hypothesis, we detected low but nontrivial gene flow among populations. Furthermore, using short-term population-genomic time-series data, we documented patterns consistent with predictions from this hypothesis, including a weak but detectable excess of evolutionary change in the direction of the mean (migrant gene pool) allele frequencies across populations and consistency in the direction of allele frequency change over time. The documented decoupling of diversity levels and short-term change by drift in Lycaeides has implications for our understanding of contemporary evolution and the maintenance of genetic variation in the wild.

The role of neutral and adaptive genomic variation in population diversification and speciation in two ground squirrel species of conservation concern

Understanding the neutral (demographic) and adaptive processes leading to the differentiation of species and populations is a critical component of evolutionary and conservation biology. In this context, recently diverged taxa represent a unique opportunity to study the process of genetic differentiation. Northern and southern Idaho ground squirrels (Urocitellus brunneus-NIDGS, and U.endemicus-SIDGS, respectively) are a recently diverged pair of sister species that have undergone dramatic declines in the last 50years and are currently found in metapopulations across restricted spatial areas with distinct environmental pressures. Here we genotyped single-nucleotide polymorphisms (SNPs) from buccal swabs with restriction site-associated DNA sequencing (RADseq). With these data we evaluated neutral genetic structure at both the inter- and intraspecific level, and identified putatively adaptive SNPs using population structure outlier detection and genotype-environment association (GEA) analyses. At the interspecific level, we detected a clear separation between NIDGS and SIDGS, and evidence for adaptive differentiation putatively linked to torpor patterns. At the intraspecific level, we found evidence of both neutral and adaptive differentiation. For NIDGS, elevation appears to be the main driver of adaptive differentiation, while neutral variation patterns match and expand information on the low connectivity between some populations identified in previous studies using microsatellite markers. For SIDGS, neutral substructure generally reflected natural geographical barriers, while adaptive variation reflected differences in land cover and temperature, as well as elevation. These results clearly highlight the roles of neutral and adaptive processes for understanding the complexity of the processes leading to species and population differentiation, which can have important conservation implications in susceptible and threatened species.

Proton loss of inner radiation belt during geomagnetic storm of 2018 based on CSES satellite observation* *Project supported by the Research Fund from the National Institute of Natural Hazards, Ministry of Emergency Management of China (Grant No. 2021-JBKY- 11), the National Natural Science Foundation of China (Grant Nos. 41904149 and 12173038), and the Stable Support Projects of Basic Scientific Research

The proton distribution in inner radiation belt is often affected by strong geomagnetic storm disturbance. Based on the data of the sun-synchronous CSES satellite, which carries with several high energy particle payloads and was launched in February 2018, we analyzed the extensive proton variations in the inner radiation belt in a wide energy range of 2 MeV– 220 MeV during 2018 major geomagnetic storm. The result indicates that the loss mechanism of protons was energy dependence which is consistent with some previous studies. For protons at low energy 2 MeV–20 MeV, the fluxes were decreased during main phase of the storm and did not come back quickly during the recovery phase, which is likely to be caused by Coulomb collision due to neutral atmosphere density variation. At higher energy 30 MeV–100 MeV, it was confirmed that the magnetic field line curvature scattering plays a significant role in the proton loss phenomenon during this storm. At highest energies > 100 MeV, the fluxes of protons kept a stable level and did not exhibit a significant loss during this storm.

Spatial population structure: patterns of adaptation in populations of the water hyacinth grasshopper Cornops aquaticum (Bruner 1906)

AbstractThe water hyacinth grasshopper Cornops aquaticum (Bruner 1906) constitutes an appropriate model to assess phenotypic and karyotypic variability in the context of differentiation and adaptation of insect populations because it occurs over a wide latitudinal range. This study represents a general analysis of phenotype, karyotype and molecular variation in native populations of C. aquaticum in South America. This is also relevant because this insect is considered a promising biological control agent of water hyacinth, a native South American aquatic plant but a pest in South Africa. Along Paraná and Uruguay River Basins, body size correlated negatively with latitude, and positively so with temperature and rainfall in both sexes. To test whether the chromosomal and phenotypic patterns were adaptive, we compared them with neutral microsatellite loci variation in populations from the medium and lower course of the Paraná River. Firstly, the lack of pairwise association between karyotype and phenotype distance matrixes with that of neutral loci suggested non-neutrality. Secondly, phenotypic differentiation for all morphometric traits (PST) was significantly larger than molecular differentiation (FST), indicating a prevailing divergence selection effect on the observed phenotypic patterns. Finally, the phenotypic and genotypic spatial structures – inferred from Bayesian approaches – were discordant: neutral genetic structure clustered together most populations except for the two southernmost, downstream ones, whereas phenotypic spatial structure groups together all the deltaic populations and singles out the two northernmost ones. The results suggest directional selection leading to higher centric fusion frequencies in the downstream populations and favouring morphometric optimal differences in relation to the environment.

Balancing selection, genetic drift, and human-mediated introgression interplay to shape MHC (functional) diversity in Mediterranean brown trout.

The extraordinary polymorphism of major histocompatibility complex (MHC) genes is considered a paradigm of pathogen‐mediated balancing selection, although empirical evidence is still scarce. Furthermore, the relative contribution of balancing selection to shape MHC population structure and diversity, compared to that of neutral forces, as well as its interaction with other evolutionary processes such as hybridization, remains largely unclear. To investigate these issues, we analyzed adaptive (MHC‐DAB gene) and neutral (11 microsatellite loci) variation in 156 brown trout (Salmo trutta complex) from six wild populations in central Italy exposed to introgression from domestic hatchery lineages (assessed with the LDH gene). MHC diversity and structuring correlated with those at microsatellites, indicating the substantial role of neutral forces. However, individuals carrying locally rare MHC alleles/supertypes were in better body condition (a proxy of individual fitness/parasite load) regardless of the zygosity status and degree of sequence dissimilarity of MHC, hence supporting balancing selection under rare allele advantage, but not heterozygote advantage or divergent allele advantage. The association between specific MHC supertypes and body condition confirmed in part this finding. Across populations, MHC allelic richness increased with increasing admixture between native and domestic lineages, indicating introgression as a source of MHC variation. Furthermore, introgression across populations appeared more pronounced for MHC than microsatellites, possibly because initially rare MHC variants are expected to introgress more readily under rare allele advantage. Providing evidence for the complex interplay among neutral evolutionary forces, balancing selection, and human‐mediated introgression in shaping the pattern of MHC (functional) variation, our findings contribute to a deeper understanding of the evolution of MHC genes in wild populations exposed to anthropogenic disturbance.

Sensory-based quantification of male colour patterns in Trinidadian guppies reveals no support for parallel phenotypic evolution in multivariate trait space.

Parallel evolution, in which independent populations evolve along similar phenotypic trajectories, offers insights into the repeatability of adaptive evolution. Here, we revisit a classic example of parallelism, that of repeated evolution of brighter males in the Trinidadian guppy (Poecilia reticulata). In guppies, colonisation of low predation habitats is associated with emergence of 'more colourful' phenotypes since predator-induced viability selection for crypsis weakens while sexual selection by female preference for conspicuousness remains strong. Our study differs from previous investigations in three respects. First, we adopted a multivariate phenotyping approach to characterise parallelism in multitrait space. Second, we used ecologically-relevant colour traits defined by the visual systems of the two selective agents (i.e., guppy, predatory cichlid). Third, we estimated population genetic structure to test for adaptive (parallel) evolution against a model of neutral phenotypic divergence. We find strong phenotypic differentiation that is inconsistent with a neutral model but very limited support for the predicted pattern of greater conspicuousness at low predation. Effects of predation regime on each trait were in the expected direction, but weak, largely nonsignificant, and explained little among-population variation. In multitrait space, phenotypic trajectories of lineages colonising low from high predation regimes were not parallel. Our results are consistent with reduced predation risk facilitating adaptive differentiation, potentially by female choice, but suggest that this proceeds in independent directions of multitrait space across lineages. Pool-sequencing data also revealed SNPs showing greater differentiation than expected under neutrality, among which some are found in genes contributing to colour pattern variation, presenting opportunities for future genetic study.

Chromosome Fusion Affects Genetic Diversity and Evolutionary Turnover of Functional Loci but Consistently Depends on Chromosome Size.

Major changes in chromosome number and structure are linked to a series of evolutionary phenomena, including intrinsic barriers to gene flow or suppression of recombination due to chromosomal rearrangements. However, chromosome rearrangements can also affect the fundamental dynamics of molecular evolution within populations by changing relationships between linked loci and altering rates of recombination. Here, we build chromosome-level assembly Eueides isabella and, together with a recent chromosome-level assembly of Dryas iulia, examine the evolutionary consequences of multiple chromosome fusions in Heliconius butterflies. These assemblies pinpoint fusion points on 10 of the 20 autosomal chromosomes and reveal striking differences in the characteristics of fused and unfused chromosomes. The ten smallest autosomes in D. iulia and E. isabella, which have each fused to a longer chromosome in Heliconius, have higher repeat and GC content, and longer introns than predicted by their chromosome length. When fused, these characteristics change to become more in line with chromosome length. The fusions also led to reduced diversity, which likely reflects increased background selection and selection against introgression between diverging populations, following a reduction in per-base recombination rate. We further show that chromosome size and fusion impact turnover rates of functional loci at a macroevolutionary scale. Together these results provide further evidence that chromosome fusion in Heliconius likely had dramatic effects on population level processes shaping rates of neutral and adaptive divergence. These effects may have impacted patterns of diversification in Heliconius, a classic example of an adaptive radiation.

Genomic data support management of anadromous Arctic Char fisheries in Nunavik by highlighting neutral and putatively adaptive genetic variation.

Distinguishing neutral and adaptive genetic variation is one of the main challenges in investigating processes shaping population structure in the wild, and landscape genomics can help identify signatures of adaptation to contrasting environments. Arctic Char (Salvelinus alpinus) is an anadromous salmonid and the most harvested fish species by Inuit people, including in Nunavik (Québec, Canada), one of the most recently deglaciated regions in the world. Unlike many other anadromous salmonids, Arctic Char occupy coastal habitats near their natal rivers during their short marine phase restricted to the summer ice‐free period. Our main objective was to document putatively neutral and adaptive genomic variation in anadromous Arctic Char populations from Nunavik and bordering regions to inform local fisheries management. We used genotyping by sequencing (GBS) to genotype 18,112 filtered single nucleotide polymorphisms (SNP) in 650 individuals from 23 sampling locations along >2000 km of coastline. Our results reveal a hierarchical genetic structure, whereby neighboring hydrographic systems harbor distinct populations grouped by major oceanographic basins: Hudson Bay, Hudson Strait, Ungava Bay, and Labrador Sea. We found genetic diversity and differentiation to be consistent both with the expected postglacial recolonization history and with patterns of isolation‐by‐distance reflecting contemporary gene flow. Results from three gene–environment association methods supported the hypothesis of local adaptation to both freshwater and marine environments (strongest associations with sea surface and air temperatures during summer and salinity). Our results support a fisheries management strategy at a regional scale, and other implications for hatchery projects and adaptation to climate change are discussed.

Broad Concordance in the Spatial Distribution of Adaptive and Neutral Genetic Variation across an Elevational Gradient in Deer Mice.

When species are continuously distributed across environmental gradients, the relative strength of selection and gene flow shape spatial patterns of genetic variation, potentially leading to variable levels of differentiation across loci. Determining whether adaptive genetic variation tends to be structured differently than neutral variation along environmental gradients is an open and important question in evolutionary genetics. We performed exome-wide population genomic analysis on deer mice sampled along an elevational gradient of nearly 4,000 m of vertical relief. Using a combination of selection scans, genotype−environment associations, and geographic cline analyses, we found that a large proportion of the exome has experienced a history of altitude-related selection. Elevational clines for nearly 30% of these putatively adaptive loci were shifted significantly up- or downslope of clines for loci that did not bear similar signatures of selection. Many of these selection targets can be plausibly linked to known phenotypic differences between highland and lowland deer mice, although the vast majority of these candidates have not been reported in other studies of highland taxa. Together, these results suggest new hypotheses about the genetic basis of physiological adaptation to high altitude, and the spatial distribution of adaptive genetic variation along environmental gradients.

Molecular Signatures of Reticulate Evolution within the Complex of European Pine Taxa

Speciation mechanisms, including the role of interspecific gene flow and introgression in the emergence of new species, are the major focus of evolutionary studies. Inference of taxonomic relationship between closely related species may be challenged by past hybridization events, but at the same time, it may provide new knowledge about mechanisms responsible for the maintenance of species integrity despite interspecific gene flow. Here, using nucleotide sequence variation and utilizing a coalescent modeling framework, we tested the role of hybridization and introgression in the evolutionary history of closely related pine taxa from the Pinus mugo complex and P. sylvestris. We compared the patterns of polymorphism and divergence between taxa and found a great overlap of neutral variation within the P. mugo complex. Our phylogeny reconstruction indicated multiple instances of reticulation events in the past, suggesting an important role of interspecific gene flow in the species divergence. The best-fitting model revealed P. mugo and P. uncinata as sister species with basal P. uliginosa and asymmetric migration between all investigated species after their divergence. The magnitude of interspecies gene flow differed greatly, and it was consistently stronger from representatives of P. mugo complex to P. sylvestris than in the opposite direction. The results indicate the prominent role of reticulation evolution in those forest trees and provide a genetic framework to study species integrity maintained by selection and local adaptation.

Characterizing the demographic history and prion protein variation to infer susceptibility to chronic wasting disease in a naïve population of white-tailed deer (Odocoileus virginianus).

Assessments of the adaptive potential in natural populations are essential for understanding and predicting responses to environmental stressors like climate change and infectious disease. Species face a range of stressors in human‐dominated landscapes, often with contrasting effects. White‐tailed deer (Odocoileus virginianus; deer) are expanding in the northern part of their range following decreasing winter severity and increasing forage availability. Chronic wasting disease (CWD), a prion disease affecting deer, is likewise expanding and represents a major threat to deer and other cervids. We obtained tissue samples from free‐ranging deer across their native range in Ontario, Canada, which has yet to detect CWD in wild populations. We used high‐throughput sequencing to assess neutral genomic variation and variation in the prion protein gene (PRNP) that is partly responsible for the protein misfolding when deer contract CWD. Neutral variation revealed a high number of rare alleles and no population structure, and demographic models suggested a rapid historical population expansion. Allele frequencies of PRNP variants associated with CWD susceptibility and disease progression were evenly distributed across the landscape and consistent with deer populations not infected with CWD. We estimated the selection coefficient of CWD, with simulations showing an observable and rapid shift in PRNP allele frequencies that coincides with the start of a novel CWD outbreak. Sustained surveillance of genomic and PRNP variation can be a useful tool for guiding management practices, which is especially important for CWD‐free regions where deer are managed for ecological and economic benefits.

Environmental adaptation is stronger for abundant rather than rare microorganisms in wetland soils from the Qinghai-Tibet Plateau.

Disentangling the biogeographic patterns of rare and abundant microbes is essential in order to understand the generation and maintenance of microbial diversity with respect to the functions they provide. However, little is known about ecological assembly processes and environmental adaptation of rare and abundant microbes across large spatial-scale wetlands. Using Illumina sequencing and multiple statistical analyses, we characterized the taxonomic and phylogenetic diversity of rare and abundant bacteria and fungi in Qinghai-Tibet Plateau wetland soils. Abundant microbial taxa exhibited broader environmental thresholds and stronger phylogenetic signals for ecological traits than rare ones. By contrast, rare taxa showed higher sensitivity to environmental changes and closer phylogenetic clustering than abundant ones. The null model analysis revealed that dispersal limitation belonging to stochastic process dominated community assemblies of abundant bacteria, and rare and abundant fungi, while variable selection belonging to deterministic process governed community assembly of rare bacteria. Neutral model analysis and variation partitioning analysis further confirmed that abundant microbes were less environmentally constrained. Soil ammonia nitrogen was the crucial factor in mediating the balance between stochasticity and determinism of both rare and abundant microbes. Abundant microbes may have better environmental adaptation potential and are less dispersed by environmental changes than rare ones. Our findings extend knowledge of the adaptation of rare and abundant microbes to ongoing environmental change and could facilitate prediction of biodiversity loss caused probably by climate change and human activity in the Qinghai-Tibet Plateau wetlands.

Role of the transverse ligament of the ulnar collateral ligament of the elbow: a biomechanical study.

BackgroundThe ulnar collateral ligament (UCL) complex of the elbow plays a primary role in valgus and posteromedial stability of the elbow. The anterior oblique ligament (AOL) of the UCL is believed to provide the majority of resistance to external forces on the medial elbow. The transverse ligament (TL) of the UCL is generally thought to have minimal contribution to the elbow’s overall stability. However, recent studies have suggested a more significant role for the TL. The primary aim of this study was to identify the TL’s contribution to the stability of the elbow joint in determining the joint stiffness and neutral zone variation in internal rotation.MethodsTwelve cadaveric elbows, set at a 90° flexion angle, were tested by applying an internal rotational force on the humerus to generate a medial opening torque at the level of the elbow. The specimens were preconditioned with 10 cycles of humeral internal rotation with sinusoidal torque ranging from 0 to 5 Nm. Elbow stiffness measures and joint neutral zone were first evaluated in its integrity during a final ramp loading. The test was subsequently repeated after cutting the TL at 33%, 66%, and 100% followed by the AOL in the same fashion.ResultsThe native UCL complex joint stiffness to internal rotation measured 1.52 ± 0.51 Nm/°. The first observable change occurred with 33% sectioning of the AOL, with further sectioning of the AOL minimizing the joint stiffness to 1.26 ± 0.32 Nm/° (P = .004). A 33% resection of the TL found an initial neutral zone variation of 0.376 ± 0.23° that increased to 0.771 ± 0.41° (P < .01) at full resection. These values were marginal when compared with the full resection of the AOL for which we have found 3.69 ± 1.65° (P < .01).ConclusionThe TL had no contribution to internal rotation elbow joint stiffness at a flexion angle of 90°. However, sequential sectioning of the TL was found to significantly increase the joint neutral zone when compared with the native cadaveric elbow at a flexion angle of 90°. This provides evidence toward the TL having some form of contribution to the elbow’s overall stability.

Genomic signatures of divergent selection are associated with social behaviour for spinner dolphin ecotypes.

Understanding the genomic basis of adaptation is critical for understanding evolutionary processes and predicting how species will respond to environmental change. Spinner dolphins in the eastern tropical Pacific (ETP) present a unique system for studying adaptation. Within this large geographical region are four spinner dolphin ecotypes with weak neutral genetic divergence and no obvious barriers to gene flow, but strong spatial variation in morphology, behaviour and habitat. These ecotypes have large population sizes, which could reduce the effects of drift and facilitate selection. To identify genomic regions putatively under divergent selective pressures between ecotypes, we used genome scans with 8994 RADseq single nucleotide polymorphisms (SNPs) to identify population differentiation outliers and genotype-environment association outliers. Gene ontology enrichment analyses indicated that outlier SNPs from both types of analyses were associated with multiple genes involved in social behaviour and hippocampus development, including 15 genes associated with the human social disorder autism. Evidence for divergent selection on social behaviour is supported by previous evidence that these spinner dolphin ecotypes differ in mating systems and associated social behaviours. In particular, three of the ETP ecotypes probably have a polygynous mating system characterized by strong premating competition among males, whereas the fourth ecotype probably has a polygynandrous mating system characterized by strong postmating competition such as sperm competition. Our results provide evidence that selection for social behaviour may be an evolutionary force driving diversification of spinner dolphins in the ETP, potentially as a result of divergent sexual selection associated with different mating systems. Future studies should further investigate the potential adaptive role of the candidate genes identified here, and could probably find further signatures of selection using whole genome sequence data.

Opportunities and challenges in assessing climate change vulnerability through genomics

By investigating how past selection has affected allele frequencies across space, genomic tools are providing new insights into adaptive evolutionary processes. Now researchers are considering how this genomic information can be used to predict the future vulnerability of species under climate change. Genomic vulnerability assessments show promise, but challenges remain.

Transcriptional activity of transposable elements along an elevational gradient in Arabidopsis arenosa

BackgroundPlant genomes can respond rapidly to environmental changes and transposable elements (TEs) arise as important drivers contributing to genome dynamics. Although some elements were reported to be induced by various abiotic or biotic factors, there is a lack of general understanding on how environment influences the activity and diversity of TEs. Here, we combined common garden experiment with short-read sequencing to investigate genomic abundance and expression of 2245 consensus TE sequences (containing retrotransposons and DNA transposons) in an alpine environment in Arabidopsis arenosa. To disentangle general trends from local differentiation, we leveraged four foothill-alpine population pairs from different mountain regions. Seeds of each of the eight populations were raised under four treatments that differed in temperature and irradiance, two factors varying with elevation. RNA-seq analysis was performed on leaves of young plants to test for the effect of elevation and subsequently of temperature and irradiance on expression of TE sequences.ResultsGenomic abundance of the 2245 consensus TE sequences varied greatly between the mountain regions in line with neutral divergence among the regions, representing distinct genetic lineages of A. arenosa. Accounting for intraspecific variation in abundance, we found consistent transcriptomic response for some TE sequences across the different pairs of foothill-alpine populations suggesting parallelism in TE expression. In particular expression of retrotransposon LTR Copia (e.g. Ivana and Ale clades) and LTR Gypsy (e.g. Athila and CRM clades) but also non-LTR LINE or DNA transposon TIR MuDR consistently varied with elevation of origin. TE sequences responding specifically to temperature and irradiance belonged to the same classes as well as additional TE clades containing potentially stress-responsive elements (e.g. LTR Copia Sire and Tar, LTR Gypsy Reina).ConclusionsOur study demonstrated that the A. arenosa genome harbours a considerable diversity of TE sequences whose abundance and expression response varies across its native range. Some TE clades may contain transcriptionally active elements responding to a natural environmental gradient. This may further contribute to genetic variation between populations and may ultimately provide new regulatory mechanisms to face environmental challenges.

The inflated significance of neutral genetic diversity in conservation genetics.

The current rate of species extinction is rapidly approaching unprecedented highs, and life on Earth presently faces a sixth mass extinction event driven by anthropogenic activity, climate change, and ecological collapse. The field of conservation genetics aims at preserving species by using their levels of genetic diversity, usually measured as neutral genome-wide diversity, as a barometer for evaluating population health and extinction risk. A fundamental assumption is that higher levels of genetic diversity lead to an increase in fitness and long-term survival of a species. Here, we argue against the perceived importance of neutral genetic diversity for the conservation of wild populations and species. We demonstrate that no simple general relationship exists between neutral genetic diversity and the risk of species extinction. Instead, a better understanding of the properties of functional genetic diversity, demographic history, and ecological relationships is necessary for developing and implementing effective conservation genetic strategies.

Neural divergence and hybrid disruption between ecologically isolated Heliconius butterflies

The importance of behavioral evolution during speciation is well established, but we know little about how this is manifest in sensory and neural systems. A handful of studies have linked specific neural changes to divergence in host or mate preferences associated with speciation. However, the degree to which brains are adapted to local environmental conditions, and whether this contributes to reproductive isolation between close relatives that have diverged in ecology, remains unknown. Here, we examine divergence in brain morphology and neural gene expression between closely related, but ecologically distinct, Heliconius butterflies. Despite ongoing gene flow, sympatric species pairs within the melpomene-cydno complex are consistently separated across a gradient of open to closed forest and decreasing light intensity. By generating quantitative neuroanatomical data for 107 butterflies, we show that Heliconius melpomene and Heliconius cydno clades have substantial shifts in brain morphology across their geographic range, with divergent structures clustered in the visual system. These neuroanatomical differences are mirrored by extensive divergence in neural gene expression. Differences in both neural morphology and gene expression are heritable, exceed expected rates of neutral divergence, and result in intermediate traits in first-generation hybrid offspring. Strong evidence of divergent selection implies local adaptation to distinct selective optima in each parental microhabitat, suggesting the intermediate traits of hybrids are poorly matched to either condition. Neural traits may therefore contribute to coincident barriers to gene flow, thereby helping to facilitate speciation.

Dispersal and Good Habitat Quality Promote Neutral Genetic Diversity in Metapopulations.

Dispersal is a fundamental and crucial ecological process for a metapopulation to survive in heterogeneous or changing habitats. In this paper, we investigate the effect of the habitat quality and the dispersal on the neutral genetics diversity of a metapopulation. We model the metapopulation dynamics on heterogeneous habitats using a deterministic system of ordinary differential equations. We decompose the metapopulation into several neutral genetic fractions seeing as they could be located in different habitats. By using a mathematical model which describes their temporal dynamics inside the metapopulation, we provide the analytical results of their transient dynamics, as well as their asymptotic proportion in the different habitats. The diversity indices show how the genetic diversity at a global metapopulation scale is preserved by the correlation of two factors: the dispersal of the population, as well as the existence of adequate and sufficiently large habitats. The diversity indices show how the genetic diversity at a global metapopulation scale is preserved by the correlation of two factors: the dispersal of the population as well as the existence of adequate and sufficiently large habitats. Moreover, they ensure genetic diversity at the local habitat scale. In a source-sink metapopulation, we demonstrate that the diversity of the sink can be rescued if the condition of the sink is not too deteriorated and the migration from the source is larger than the migration from the sink. Furthermore, our study provides an analytical insight into the dynamics of the solutions of the systems of ordinary differential equations.