Neutral Variation Research Articles

A multiscale approach to understanding the shared blue-orange flower color polymorphism in two Lysimachia species

BackgroundPolymorphisms are common in nature, but they are rarely shared among closely related species. Polymorphisms could originate through convergence, ancestral polymorphism, or introgression. Although shared neutral genomic variation across species is commonplace, few examples of shared functional traits exist. The blue-orange petal color polymorphisms in two closely related species, Lysimachia monelli and L. arvensis were investigated with UV-vis reflectance spectra, flavonoid biochemistry, and transcriptome comparisons followed by climate niche analysis.ResultsSimilar color morphs between species have nearly identical reflectance spectra, flavonoid biochemistry, and ABP gene expression patterns. Transcriptome comparisons reveal two orange-specific genes directly involved in both blue-orange color polymorphisms: DFR-2 specificity redirects flux from the malvidin to the pelargonidin while BZ1-2 stabilizes the pelargonidin with glucose, producing the orange pelargonidin 3-glucoside. Moreover, a reduction of F3’5’H expression in orange petals also favors pelargonidin production. The climate niches for each color morph are the same between the two species for three temperature characteristics but differ for four precipitation variables.ConclusionsThe similarities in reflectance spectra, biochemistry, and ABP genes suggest that a single shift from blue-to-orange shared by both lineages is the most plausible explanation. Our evidence suggests that this persistent flower color polymorphism may represent an ancestrally polymorphic trait that has transcended speciation, yet future analyses are necessary to confidently reject the alternative hypotheses.

Deflection estimation of reinforced concrete beams using long-gauge optic sensors

Deflection is a crucial indicator in structural health monitoring, directly impacting the service life, safety, and economic viability of structures. However, current indirect methods for measuring deflection, such as multi-element FBG sensors or accelerometers, often have limitations. While they may be suitable for flexure-dominated reinforced concrete (RC) specimens in the elastic stage, they neglect the significant influence of nonlinearity and shear deformation. This paper introduces a new approach for estimating the deflection of reinforced concrete (RC) beams using distributed long-gauge optic sensors. The method tackles the limitations of current methods by incorporating the crucial factors of shear deformation, neutral axis position variation and shear stiffness degradation. Experimental studies are first conducted to validate the accuracy and stability of long-gauge optic sensors in small-scale RC members. Algorithms are then utilized the measurement data to estimate the deflection. The results demonstrate the excellent performance of long-gauge optic sensors in measuring shear deformation and vertical deflection of RC beams, as evidenced by comparisons with LVDT measurements. Additionally, a method for decoupling flexural and shear deformation is proposed. This method offers a preliminary assessment of the reinforced concrete (RC) beam type by quantifying the relative significance of shear deformation.

Where the "ruber" Meets the Road: Using the Genome of the Red Diamond Rattlesnake to Unravel the Evolutionary Processes Driving Venom Evolution.

Understanding the proximate and ultimate causes of phenotypic variation is fundamental in evolutionary research, as such variation provides the substrate for selection to act upon. Although trait variation can arise due to selection, the importance of neutral processes is sometimes understudied. We presented the first reference-quality genome of the Red Diamond Rattlesnake (Crotalus ruber) and used range-wide 'omic data to estimate the degree to which neutral and adaptive evolutionary processes shaped venom evolution. We characterized population structure and found substantial genetic differentiation across two populations, each with distinct demographic histories. We identified significant differentiation in venom expression across age classes with substantially reduced but discernible differentiation across populations. We then used conditional redundancy analysis to test whether venom expression variation was best predicted by neutral divergence patterns or geographically variable (a)biotic factors. Snake size was the most significant predictor of venom variation, with environment, prey availability, and neutral sequence variation also identified as significant factors, though to a lesser degree. By directly including neutrality in the model, our results confidently highlight the predominant, yet not singular, role of life history in shaping venom evolution.

The protein domains of vertebrate species in which selection is more effective have greater intrinsic structural disorder.

The nearly neutral theory of molecular evolution posits variation among species in the effectiveness of selection. In an idealized model, the census population size determines both this minimum magnitude of the selection coefficient required for deleterious variants to be reliably purged, and the amount of neutral diversity. Empirically, an 'effective population size' is often estimated from the amount of putatively neutral genetic diversity and is assumed to also capture a species' effectiveness of selection. A potentially more direct measure of the effectiveness of selection is the degree to which selection maintains preferred codons. However, past metrics that compare codon bias across species are confounded by among-species variation in %GC content and/or amino acid composition. Here, we propose a new Codon Adaptation Index of Species (CAIS), based on Kullback-Leibler divergence, that corrects for both confounders. We demonstrate the use of CAIS correlations, as well as the Effective Number of Codons, to show that the protein domains of more highly adapted vertebrate species evolve higher intrinsic structural disorder.

Population Genomic Scans for Natural Selection and Demography.

Uncovering the fundamental processes that shape genomic variation in natural populations is a primary objective of population genetics. These processes include demographic effects such as past changes in effective population size or gene flow between structured populations. Furthermore, genomic variation is affected by selection on nonneutral genetic variants, for example, through the adaptation of beneficial alleles or balancing selection that maintains genetic variation. In this article, we discuss the characterization of these processes using population genetic models, and we review methods developed on the basis of these models to unravel the underlying processes from modern population genomic data sets. We briefly discuss the conditions in which these approaches can be used to infer demography or identify specific nonneutral genetic variants and cases in which caution is warranted. Moreover, we summarize the challenges of jointly inferring demography and selective processes that affect neutral variation genome-wide.

Comparative Genomics Supports Ecologically Induced Selection as a Putative Driver of Banded Penguin Diversification.

The relative importance of genetic drift and local adaptation in facilitating speciation remains unclear. This is particularly true for seabirds, which can disperse over large geographic distances, providing opportunities for intermittent gene flow among distant colonies that span the temperature and salinity gradients of the oceans. Here, we delve into the genomic basis of adaptation and speciation of banded penguins, Galápagos (Spheniscus mendiculus), Humboldt (Spheniscus humboldti), Magellanic (Spheniscus magellanicus), and African penguins (Spheniscus demersus), by analyzing 114 genomes from the main 16 breeding colonies. We aim to identify the molecular mechanism and genomic adaptive traits that have facilitated their diversifications. Through positive selection and gene family expansion analyses, we identified candidate genes that may be related to reproductive isolation processes mediated by ecological thermal niche divergence. We recover signals of positive selection on key loci associated with spermatogenesis, especially during the recent peripatric divergence of the Galápagos penguin from the Humboldt penguin. High temperatures in tropical habitats may have favored selection on loci associated with spermatogenesis to maintain sperm viability, leading to reproductive isolation among young species. Our results suggest that genome-wide selection on loci associated with molecular pathways that underpin thermoregulation, osmoregulation, hypoxia, and social behavior appears to have been crucial in local adaptation of banded penguins. Overall, these results contribute to our understanding of how the complexity of biotic, but especially abiotic, factors, along with the high dispersal capabilities of these marine species, may promote both neutral and adaptive lineage divergence even in the presence of gene flow.

The relationship between neutral genetic diversity and performance in wild arthropod populations.

Larger effective populations (Ne) are characterized by higher genetic diversity, which is expected to predict population performance (average individual performance that influences fitness). Empirical studies of the relationship between neutral diversity and performance mostly represent species with small Ne, while there is limited data from the species-rich and ecologically important arthropods that are assumed to have large Ne but are threatened by massive declines. We performed a systematic literature search and used meta-analytical models to test the prediction of a positive association between neutral genetic diversity and performance in wild arthropods. From 14 relevant studies of 286 populations, we detected a weak (r = 0.15) but nonsignificant positive association both in the full data set (121 effect sizes) and a reduced data set accounting for dependency (14 effect sizes). Theory predicts that traits closely associated with fitness show a relatively stronger correlation with neutral diversity; this relationship was upheld for longevity and marginally for reproduction. Our analyses point to major knowledge gaps in our understanding of relationships between neutral diversity and performance. Future studies using genome-wide data sets across populations could guide more powerful designs to evaluate relationships between adaptive, deleterious and neutral diversity and performance.

Loss of genetic variation and ancestral sex determination system in North American northern pike characterized by whole-genome resequencing.

The northern pike Esox lucius is a freshwater fish with low genetic diversity but ecological success throughout the Northern Hemisphere. Here, we generate an annotated chromosome-level genome assembly of 941 Mbp in length with 25 chromosome-length scaffolds. We then genotype 47 northern pike from Alaska through New Jersey at a genome-wide scale and characterize a striking decrease in genetic diversity along the sampling range. Individuals west of the North American Continental Divide have substantially higher diversity than those to the east (e.g. Interior Alaska and St. Lawrence River have on average 181 and 64K heterozygous SNPs per individual, or a heterozygous SNP every 5.2 and 14.6 kbp, respectively). Individuals clustered within each population with strong support, with numerous private alleles observed within each population. Evidence for recent population expansion was observed for a Manitoba hatchery and the St. Lawrence population (Tajima's D = -1.07 and -1.30, respectively). Several chromosomes have large regions with elevated diversity, including LG24, which holds amhby, the ancestral sex determining gene. As expected amhby was largely male-specific in Alaska and the Yukon and absent southeast to these populations, but we document some amhby(-) males in Alaska and amhby(+) males in the Columbia River, providing evidence for a patchwork of presence of this system in the western region. These results support the theory that northern pike recolonized North America from refugia in Alaska and expanded following deglaciation from west to east, with probable founder effects resulting in loss of both neutral and functional diversity (e.g. amhby).

The impact of geographic isolation and host shifts on population divergence of the rare cicada Subpsaltria yangi

The contributions of divergent selection and spatial isolation to population divergence are among the main focuses of evolutionary biology. Here we employed integrated methods to explore genomic divergence, demographic history and calling-song differentiation in the cicada Subpsaltria yangi, and compared the genotype and calling-song phenotype of different populations occurring in distinct habitats. Our results indicate that this species comprises four main lineages with unique sets of haplotypes and calling-song structure, which are distinctly associated with geographic isolation and habitats. The populations occurring on the Loess Plateau underwent substantial expansion at ∼0.130–0.115 Ma during the Last Interglacial. Geographic distance and host shift between pairs of populations predict genomic divergence, with geographic distance and acoustical signal together explaining > 60% of the divergence among populations. Differences in calling songs could reflect adaptation of populations to novel environments with different host plants, habitats and predators, which may have resulted from neutral divergence at the molecular level followed by natural selection. Geomorphic barriers and climate oscillations associated with Pleistocene glaciation may have been primary factors in shaping the population genetic structure of this species. Ultimately this may couple with a host shift in leading toward allopatric speciation in S. yangi, i.e., isolation by distance. Our findings improve understanding of divergence in allopatry of herbivorous insects, and may inform future studies on the molecular mechanisms underlying the association between genetic/phenotypic changes and adaptation of insects to novel niches and host plants.

Artificial intelligence enables unified analysis of historical and landscape influences on genetic diversity

While genetic variation in any species is potentially shaped by a range of processes, phylogeography and landscape genetics are largely concerned with inferring how environmental conditions and landscape features impact neutral intraspecific diversity. However, even as both disciplines have come to utilize SNP data over the last decades, analytical approaches have remained for the most part focused on either broad-scale inferences of historical processes (phylogeography) or on more localized inferences about environmental and/or landscape features (landscape genetics). Here we demonstrate that an artificial intelligence model-based analytical framework can consider both deeper historical factors and landscape-level processes in an integrated analysis. We implement this framework using data collected from two Brazilian anurans, the Brazilian sibilator frog (Leptodactylus troglodytes) and granular toad (Rhinella granulosa). Our results indicate that historical demographic processes shape most the genetic variation in the sibulator frog, while landscape processes primarily influence variation in the granular toad. The machine learning framework used here allows both historical and landscape processes to be considered equally, rather than requiring researchers to make an a priori decision about which factors are important.

The substrate quality of CK2 target sites has a determinant role on their function and evolution

Most biological processes are regulated by signaling modules that bind to short linear motifs. For protein kinases, substrates may have full or only partial matches to the kinase recognition motif, a property known as "substrate quality." However, it is not clear whether differences in substrate quality represent neutral variation or if they have functional consequences. We examine this question for the kinase CK2, which has many fundamental functions. We show that optimal CK2 sites are phosphorylated at maximal stoichiometries and found in many conditions, whereas minimal substrates are more weakly phosphorylated and have regulatory functions. Optimal CK2 sites tend to be more conserved, and substrate quality is often tuned by selection. For intermediate sites, increases or decreases in substrate quality may be deleterious, as we demonstrate for a CK2 substrate at the kinetochore. The results together suggest a strong role for substrate quality in phosphosite function and evolution. A record of this paper's transparent peer review process is included in the supplemental information.

Inherent versus random variation in fitness of elephant seals: offspring quality and quantity

Variation in reproductive success is the basis of evolution and allows species to respond to the environment, but only when it is based on fixed individual variation that is heritable. Several recent studies suggest that observed variation in reproduction is due to chance, not inherent individual differences. Our aim was to quantify inherent versus neutral variation in fitness of northern elephant seal ( Mirounga angustirostris (Gill, 1866)) females, including both quality and quantity of their offspring. Using 44 years of observations at Año Nuevo in California, we assembled lifetime pup production of 1065 individual females and mass at weaning for 2120 of their pups. Females varied significantly in mean lifetime mass of their pups, with 28% of the variance due to fixed individual differences among mothers. Variation was repeatable over 6 years of a mother’s lifetime and heritable ( h = 0.48). Moreover, pup mass at weaning was associated with future lifetime fitness, since larger pups had a higher chance of surviving to breed. Larger pups, however, did not produce more offspring once breeding, and lifetime pup production was not heritable. Traits related to offspring quality in elephant seals were inherently different among females, but variation in pup production was neutral.

Temporally resolved relative krypton neutral density during breathing mode of a hall effect thruster recorded by TALIF

A key issue in the development of theory and models for plasma propulsion devices is to describe the instabilities and fluctuations of the devices. It has been widely recognized that many Hall effect thrusters (HETs) exhibit oscillations at frequencies in the range of ∼ 20 kHz. These ionization-related oscillations are generally referred to as Breathing Mode oscillations and have been the subject of considerable research. Here, for the first time, we report direct temporally resolved measurements of the ground state neutral density variation during the period of the oscillation. We used the laser-based Two-Photon Absorption Laser Induced Fluorescence (TALIF) technique to measure neutrals within the plume of a 1.5 kW HET operating on krypton (Kr). Our TALIF scheme employs a frequency-doubled, pulsed dye laser operating at ∼ 212 nm to probe ground state Kr atoms. A novel phase-binning approach is used to recover the time-dependent signal by assigning the timing of each collected TALIF signal (laser shot) relative to the phase of the discharge current. We find that the neutral density fluctuates quite strongly over the period of the oscillation, and that this fluctuation leads the current fluctuation as expected.

Selection drives divergence of eye morphology in sympatric Heliconius butterflies.

When populations experience different sensory conditions, natural selection may favor sensory system divergence, affecting peripheral structures and/or downstream neural pathways. We characterized the outer eye morphology of sympatric Heliconius butterflies from different forest types and their first-generation reciprocal hybrids to test for adaptive visual system divergence and hybrid disruption. In Panama, Heliconius cydno occurs in closed forests, whereas Heliconius melpomene resides at the forest edge. Among wild individuals, H. cydno has larger eyes than H. melpomene, and there are heritable, habitat-associated differences in the visual brain structures that exceed neutral divergence expectations. Notably, hybrids have intermediate neural phenotypes, suggesting disruption. To test for similar effects in the visual periphery, we reared both species and their hybrids in common garden conditions. We confirm that H. cydno has larger eyes and provide new evidence that this is driven by selection. Hybrid eye morphology is more H. melpomene-like despite body size being intermediate, contrasting with neural trait intermediacy. Overall, our results suggest that eye morphology differences between H. cydno and H. melpomene are adaptive and that hybrids may suffer fitness costs due to a mismatch between the peripheral visual structures and previously described neural traits that could affect visual performance.

Fine scale diversity in the lava: genetic and phenotypic diversity in small populations of Arctic charr Salvelinus alpinus.

A major goal in evolutionary biology is to understand the processes underlying phenotypic variation in nature. Commonly, studies have focused on large interconnected populations or populations found along strong environmental gradients. However, studies on small fragmented populations can give strong insight into evolutionary processes in relation to discrete ecological factors. Evolution in small populations is believed to be dominated by stochastic processes, but recent work shows that small populations can also display adaptive phenotypic variation, through for example plasticity and rapid adaptive evolution. Such evolution takes place even though there are strong signs of historical bottlenecks and genetic drift. Here we studied 24 small populations of the freshwater fish Arctic charr (Salvelinus alpinus) found in groundwater filled lava caves. Those populations were found within a few km2-area with no apparent water connections between them. We studied the relative contribution of neutral versus non-neutral evolutionary processes in shaping phenotypic divergence, by contrasting patterns of phenotypic and neutral genetic divergence across populations in relation to environmental measurements. This allowed us to model the proportion of phenotypic variance explained by the environment, taking in to account the observed neutral genetic structure. These populations originated from the nearby Lake Mývatn, and showed small population sizes with low genetic diversity. Phenotypic variation was mostly correlated with neutral genetic diversity with only a small environmental effect. Phenotypic diversity in these cave populations appears to be largely the product of neutral processes, fitting the classical evolutionary expectations. However, the fact that neutral processes did not explain fully the phenotypic patterns suggests that further studies can increase our understanding on how neutral evolutionary processes can interact with other forces of selection at early stages of divergence. The accessibility of these populations has provided the opportunity for long-term monitoring of individual fish, allowing tracking how the environment can influence phenotypic and genetic divergence for shaping and maintaining diversity in small populations. Such studies are important, especially in freshwater, as habitat alteration is commonly breaking populations into smaller units, which may or may not be viable.

Mutation Rate and Effective Population Size of the Model Cooperative Bacterium Myxococcus xanthus.

Intrinsic rates of genetic mutation have diverged greatly across taxa and exhibit statistical associations with several other parameters and features. These include effective population size (Ne), genome size, and gametic multicellularity, with the latter being associated with both increased mutation rates and decreased effective population sizes. However, data sufficient to test for possible relationships between microbial multicellularity and mutation rate (µ) are lacking. Here, we report estimates of two key population-genetic parameters, Ne and µ, for Myxococcus xanthus, a bacterial model organism for the study of aggregative multicellular development, predation, and social swarming. To estimate µ, we conducted an ∼400-day mutation accumulation experiment with 46 lineages subjected to regular single colony bottlenecks prior to clonal regrowth. Upon conclusion, we sequenced one clonal-isolate genome per lineage. Given collective evolution for 85,323 generations across all lines, we calculate a per base-pair mutation rate of ∼5.5 × 10-10 per site per generation, one of the highest mutation rates among free-living eubacteria. Given our estimate of µ, we derived Ne at ∼107 from neutral diversity at four-fold degenerate sites across two dozen M. xanthus natural isolates. This estimate is below average for eubacteria and strengthens an already clear negative correlation between µ and Ne in prokaryotes. The higher and lower than average mutation rate and Ne for M. xanthus, respectively, amplify the question of whether any features of its multicellular life cycle-such as group-size reduction during fruiting-body development-or its highly structured spatial distribution have significantly influenced how these parameters have evolved.

Neutral diversity in experimental metapopulations

New automated and high-throughput methods allow the manipulation and selection of numerous bacterial populations. In this manuscript we are interested in the neutral diversity patterns that emerge from such a setup in which many bacterial populations are grown in parallel serial transfers, in some cases with population-wide extinction and splitting events. We model bacterial growth by a birth–death process and use the theory of coalescent point processes. We show that there is a dilution factor that optimises the expected amount of neutral diversity for a given number of cycles, and study the power law behaviour of the mutation frequency spectrum for different experimental regimes. We also explore how neutral variation diverges between two recently split populations by establishing a new formula for the expected number of shared and private mutations. Finally, we show the interest of such a setup to select a phenotype of interest that requires multiple mutations.

Background Selection From Unlinked Sites Causes Nonindependent Evolution of Deleterious Mutations.

Background selection describes the reduction in neutral diversity caused by selection against deleterious alleles at other loci. It is typically assumed that the purging of deleterious alleles affects linked neutral variants, and indeed simulations typically only treat a genomic window. However, background selection at unlinked loci also depresses neutral diversity. In agreement with previous analytical approximations, in our simulations of a human-like genome with a realistically high genome-wide deleterious mutation rate, the effects of unlinked background selection exceed those of linked background selection. Background selection reduces neutral genetic diversity by a factor that is independent of census population size. Outside of genic regions, the strength of background selection increases with the mean selection coefficient, contradicting the linked theory but in agreement with the unlinked theory. Neutral diversity within genic regions is fairly independent of the strength of selection. Deleterious genetic load among haploid individuals is underdispersed, indicating nonindependent evolution of deleterious mutations. Empirical evidence for underdispersion was previously interpreted as evidence for global epistasis, but we recover it from a non-epistatic model.

Identification of low levels of neutral and functional genetic diversity in South African bontebok (Damaliscus pygargus pygargus).

Bontebok (Damaliscus pygargus pygargus) and blesbok (D. p. phillipsi) are classified as separate sub-species. The blesbok has a widespread distribution throughout South Africa and is listed as least concern by the International Union for Conservation of Nature (IUCN) Red List of Threatened Species. Bontebok on the other hand is endemic within the Cape Floristic Region of the Western Cape in South Africa and has been listed as near-threatened species on the IUCN Red List of Threatened Species. Bontebok populations experienced a severe bottleneck and were brought back from the brink of extinction in the 1830s. Currently, the subspecies is threatened by hybridisation with blesbok resulting in fertile offspring. To date, molecular investigations using neutral markers have determined that genetic diversity in pure South African bontebok was significantly lower than in pure blesbok. Here, we investigated genetic diversity in bontebok, blesbok and hybrid individuals using microsatellites and an adaptive marker (toll-like receptor two (TLR2)). The study of single nucleotide polymorphisms (SNPs) revealed five mutations in TLR2 in different individuals and subspecies of D. pygargus. This included three non-synonymous and two synonymous mutations. The three amino acid substitution mutations were predicted to have no effect on protein function. Two of the five mutations, one of which resulted in an amino acid substitution, were not present in bontebok. The other three mutations were present to varying frequencies in the three groups. We confirm low adaptive and neutral diversity in bontebok. These mutations provide insights into the genetic diversity and relationships among the two sub-species of D. pygargus and may have implications for their conservation and management.

Local adaptive variation in a highly migratory fish: The smooth hammerhead shark Sphyrna zygaena

Populations of highly mobile species that undertake long distance migrations are typically considered to be panmictic. Nonetheless, mechanisms related to behavior or local environmental conditions promote genetic isolation in the absence of physical barriers. Highly migratory shark species exhibit varying levels of fidelity to specific regions, shaping the genetic architecture of different populations and resulting in geographically based genetic variation with potential adaptive value. An understanding of the genetic variation of highly migratory species is needed to develop effective conservation strategies. This study aimed to assess the neutral and adaptive variation of the smooth hammerhead shark (Sphyrna zygaena) in the northern Mexican Pacific (NMP) via single nucleotide polymorphisms (SNPs). We analyzed 1480 SNPs in 92 individuals from four geographic regions in the NMP, of which 1469 SNPs were neutral loci (n-SNP), and 11 were putatively under selection (o-SNP) using four genoma scan methods. Genetic diversity was geographically similar among regions (Ho = 0.275). The neutral variation showed panmixia (n-SNPs; FST = 0.0012, p = 0.44), which may be associated with the high dispersal capacity of S. zygaena. A pattern of adaptive variation between individuals from the Gulf of California and Pacific coast was revealed using o-SNPs FST-based methods (24 oSNPs; FST = 0.061, p < 0.001), which may be promoted by individual preferences based on physiological limitations. The estimated effective population size (Ne) of S. zygaena was 1390 individuals, which is theoretically optimal for the population to persist over time.