Patterns Of Genomic Variation Research Articles

Origin of resistance to pyrethroids in the redlegged earth mite (Halotydeus destructor) in Australia: repeated local evolution and migration.

Halotydeus destructor is a major pest of crops and pastures across southern parts of Australia. This invasive mite has been chemically controlled for over 50 years, but resistance to synthetic pyrethroids and organophosphates is developing. Understanding processes behind the emerging resistance is important for effective management efforts. We undertook a ddRAD pool-sequencing approach to analyse genome-wide single nucleotide polymorphism variation in H. destructor population samples at two scales: local resistance across a set of fields, and regional resistance across their Australian range, along with toxicology bioassays to screen for pyrethroid resistance. Spatial patterns of genomic variation and resistance at a local scale indicated that genetic similarity among samples was more closely correlated with distance along roads and fence-lines than with straight-line geographic distance. This pattern was particularly strong in resistant samples, which were also more related than susceptible samples, suggesting local spread of resistance within an area after it emerged. By contrast, regional data suggest resistance has emerged repeatedly within parts of Australia. Our de novo annotation of the H. destructor draft genome sequence and Bayesian analysis identified several candidate loci strongly associated with population-level resistance to pyrethroids, located in genomic regions that code for transmembrane transport and signalling proteins that have previously been linked to insecticide resistance in other arthropods. Our findings highlight multiple independent evolutionary events leading to resistance in H. destructor, and demonstrate the utility and cost-effectiveness of a cross-population, genome-wide association study to reveal processes underlying adaptive evolution in a non-model invasive species. © 2019 Society of Chemical Industry.

Hybridization and a loss of sex shape genome-wide diversity and the origin of species in the evening primroses (Oenothera, Onagraceae).

Hybridization is thought to promote speciation in at least two ways - by fixation of heterozygosity from diploid progenitors in allopolyploids, and by generation of transgressive phenotypes and shifting fitness optima during homoploid hybrid speciation. While recent studies support a hybrid origin for a growing number of species, the extent to which hybrid origins shape patterns of diversity in asexual species remains underexplored. Here we employed transcriptome sequencing and population genomic analysis to describe patterns of genomic variation in the 13 species belonging to Oenothera subsection Oenothera. Eight of these species are functionally asexual and arose by hybrid speciation from parents spanning a range of phylogenetic divergence. We showed that genomic divergence between parents has been retained as heterozygosity in functionally asexual species, and that genome-wide levels of heterozygosity in these asexuals largely reflects the divergence of parental haplotypes coupled with a breakdown in recombination and segregation across the genome. These results show that divergence between parental species and loss of sex in hybrids shape patterns of whole-genome diversity and the origin of asexual species.

Attenuation of Infectious Bronchitis Virus in Eggs Results in Different Patterns of Genomic Variation across Multiple Replicates.

The gammacoronavirus infectious bronchitis virus (IBV) causes an acute, highly contagious respiratory disease of poultry. Live attenuated vaccines are traditionally generated by serial passage of a virulent strain in embryonated chicken eggs; however, the molecular mechanism of attenuation is unknown. M41-CK, a virulent lab-adapted strain of IBV, was egg passaged over 100 times in four parallel independent replicates. All four final egg-passaged viruses were attenuated in vivo and exhibited similar growth phenotypes in adult chicken kidney cells and ex vivo tracheal organ cultures. The virus populations were sequenced by 454 pyrosequencing at the end of passaging, and the results showed that overall sequence diversity in the IBV population increased but the four replicates only had between 11 and 17 consensus-level single nucleotide polymorphisms (SNPs). Although hot spots of variation were identified in spike and nucleocapsid structural proteins as well as the 3' untranslated region, each attenuated virus possessed a different pattern of genomic variation. Overall, only a small number of consensus-level SNPs were acquired during egg passage, leaving a potentially short route back to virulence. These results highlight the unpredictable nature of attenuation by serial egg passage and the need to develop mechanisms to rationally attenuate IBV for the next generation of effective vaccines.IMPORTANCE Infectious bronchitis remains a major problem in the global poultry industry, despite the existence of many different vaccines. IBV vaccines are currently developed by serial passage of a virulent strain on embryonated hen's eggs until attenuation; however, little is known about the evolution of the viral population during the process of attenuation. High-throughput sequencing of four replicates of a serially egg-passaged IBV revealed a different pattern of genomic variation in each attenuated replicate and few consensus-level SNPs. This raises concerns that only a small number of genomic mutations are required to revert to a virulent phenotype, which may result in vaccine breakdown in the field. The observed hot spots of variation in the attenuated viruses have the potential to be used in the rational attenuation of virulent IBV for next-generation vaccine design.

Consequences of breed formation on patterns of genomic diversity and differentiation: the case of highly diverse peripheral Iberian cattle

BackgroundIberian primitive breeds exhibit a remarkable phenotypic diversity over a very limited geographical space. While genomic data are accumulating for most commercial cattle, it is still lacking for these primitive breeds. Whole genome data is key to understand the consequences of historic breed formation and the putative role of earlier admixture events in the observed diversity patterns.ResultsWe sequenced 48 genomes belonging to eight Iberian native breeds and found that the individual breeds are genetically very distinct with FST values ranging from 4 to 16% and have levels of nucleotide diversity similar or larger than those of their European counterparts, namely Jersey and Holstein. All eight breeds display significant gene flow or admixture from African taurine cattle and include mtDNA and Y-chromosome haplotypes from multiple origins. Furthermore, we detected a very low differentiation of chromosome X relative to autosomes within all analyzed taurine breeds, potentially reflecting male-biased gene flow.ConclusionsOur results show that an overall complex history of admixture resulted in unexpectedly high levels of genomic diversity for breeds with seemingly limited geographic ranges that are distantly located from the main domestication center for taurine cattle in the Near East. This is likely to result from a combination of trading traditions and breeding practices in Mediterranean countries. We also found that the levels of differentiation of autosomes vs sex chromosomes across all studied taurine and indicine breeds are likely to have been affected by widespread breeding practices associated with male-biased gene flow.

Spatiotemporally explicit demographic modelling supports a joint effect of historical barriers to dispersal and contemporary landscape composition on structuring genomic variation in a red‐listed grasshopper

Inferring the processes underlying spatial patterns of genomic variation is fundamental to understand how organisms interact with landscape heterogeneity and to identify the factors determining species distributional shifts. Here, we use genomic data (restriction site-associated DNA sequencing) to test biologically informed models representing historical and contemporary demographic scenarios of population connectivity for the Iberian cross-backed grasshopper Dociostaurus hispanicus, a species with a narrow distribution that currently forms highly fragmented populations. All models incorporated biological aspects of the focal taxon that could hypothetically impact its geographical patterns of genomic variation, including (a) spatial configuration of impassable barriers to dispersal defined by topographic landscapes not occupied by the species; (b) distributional shifts resulting from the interaction between the species bioclimatic envelope and Pleistocene glacial cycles; and (c) contemporary distribution of suitable habitats after extensive land clearing for agriculture. Spatiotemporally explicit simulations under different scenarios considering these aspects and statistical evaluation of competing models within an Approximate Bayesian Computation framework supported spatial configuration of topographic barriers to dispersal and human-driven habitat fragmentation as the main factors explaining the geographical distribution of genomic variation in the species, with no apparent impact of hypothetical distributional shifts linked to Pleistocene climatic oscillations. Collectively, this study supports that both historical (i.e., topographic barriers) and contemporary (i.e., anthropogenic habitat fragmentation) aspects of landscape composition have shaped major axes of genomic variation in the studied species and emphasizes the potential of model-based approaches to gain insights into the temporal scale at which different processes impact the demography of natural populations.

Genome-wide analyses of single nucleotide polymorphisms reveal the consequences of traditional mass-rearing on genetic variation in Aphytis melinus (Hymenoptera: Aphelinidae): the danger of putting all eggs in one basket.

Aphytis melinus DeBach (Hymenoptera: Aphelinidae) is a highly effective biocontrol agent of the California red scale Aonidiella aurantii (Maskell) (Hemiptera: Diaspididae). It is commercially reared and used for augmentative releases within integrated pest management programs. However, mass rearing of biocontrol agents can result in population bottlenecks and high levels of inbreeding and/or adaptation to the factitious rearing conditions. Although these factors can all negatively impact field performance of biocontrol agents, few empirical studies have examined the genetic consequences of mass rearing. We used double-digest RAD sequencing (ddRADseq) to investigate the effect of traditional mass rearing on genetic variation among insectary colonies of A. melinus relative to wild populations in native (Pakistan) and introduced (California) ranges. Analyses of up to 9700 single nucleotide polymorphisms (SNPs) revealed that insectary populations had less genomic variation than introduced populations. This was evidenced by fewer private alleles, reduced heterozygosity, and greater missing data in the insectary populations. Further, California insectaries formed a distinct genomic cluster relative to the other samples, a surprising result given that the insectary colonies were putatively established at different times and from different source populations. These differences were evident across most data sets also after we filtered out contaminant DNA from the most common host species (Aspidiotus nerii Bouché and A. aurantii). We hypothesize that this pattern would only result if: (i) directional selection for 'captive' phenotypes produces convergent patterns of genomic variation across insectaries; or (ii) the California insectary colonies were all founded from a unifying source population and/or that the insectaries regularly exchange 'genetic' stocks. We show that RADseq is an effective method to investigate the effects of mass rearing on genetics of biocontrol agents. © 2019 Society of Chemical Industry.

A population genomics insight by 2b‐RAD reveals populations' uniqueness along the Italian coastline in Leptopsammia pruvoti (Scleractinia, Dendrophylliidae)

AbstractAimMarine bioconstructions such as coralligenous formations are hotspot of biodiversity and play a relevant ecological role in the preservation of biodiversity by providing carbon regulation, protection and nursery areas for several marine species. For this reason, the European Union Habitat Directive included them among priority habitats to be preserved. Although their ecological role is well established, connectivity patterns are still poorly investigated, representing a limit in conservation planning. The present study pioneers a novel approach for the analysis of connectivity in marine bioconstructor species, which often lack suitable genetic markers, by taking advantage of next‐generation sequencing techniques. We assess the geographical patterns of genomic variation of the sunset cup coral Leptopsammia pruvoti Lacaze‐Duthiers, 1897, an ahermatypic, non‐zooxanthellate and solitary scleractinian coral species common in coralligenous habitats and distributed across the Mediterranean Sea.LocationThe Italian coastline (Western and Central Mediterranean).MethodsWe applied the restriction site‐associated 2b‐RAD approach to genotype over 1,000 high‐quality and filtered single nucleotide polymorphisms in 10 population samples.ResultsThe results revealed the existence of a strongly supported genetic structure, with highly significant pairwise FST values between all the population samples, including those collected about 5 km apart from each other. Moreover, genomic data indicate that the strongest barriers to gene flow are between the western (Ligurian–Tyrrhenian Sea) and the eastern side (Adriatic Sea) of the Italian peninsula.Main conclusionsThe strong differentiation found in L. pruvoti is similar to that found in other species of marine bioconstructors investigated in this area, but it strongly contrasts with the small differences found in many fish and invertebrates at the same geographical scale. All in one, our results highlight the importance of assessing connectivity in species belonging to coralligenous habitats as, due to their limited dispersal ability, they might require specific spatial conservation measures.

Biodiversity genomics of North American Dryobates woodpeckers reveals little gene flow across the D. nuttallii x D. scalaris contact zone

AbstractEvolutionary biologists have long used behavioral, ecological, and genetic data from contact zones between closely related species to study various phases of the speciation continuum. North America has several concentrations of avian contact zones, where multiple pairs of sister lineages meet, with or without hybridization. In a southern California contact zone, 2 species of woodpeckers, Nuttall’s Woodpecker (Dryobates nuttallii) and the Ladder-backed Woodpecker (D. scalaris), occasionally hybridize. We sampled these 2 species in a transect across this contact zone and included samples of their closest relative, the Downy Woodpecker (D. pubescens), to obtain large single nucleotide polymorphism panels using restriction-site associated DNA sequencing (RAD-seq). Furthermore, we used whole-genome resequencing data for 2 individuals per species to identify whether patterns of diversity inferred from RAD-seq were representative of whole-genome diversity. We found that these 3 woodpecker species are genomically distinct. Although low levels of gene flow occur between D. nuttallii and D. scalaris across the contact zone, there was no evidence for widespread genomic introgression between these 2 species. Overall patterns of genomic diversity from the RAD-seq and wholegenome datasets appear to be related to distributional range size and, by extension, are likely related to effective population sizes for each species.

Population genomics of rapid evolution in natural populations: polygenic selection in response to power station thermal effluents

BackgroundExamples of rapid evolution are common in nature but difficult to account for with the standard population genetic model of adaptation. Instead, selection from the standing genetic variation permits rapid adaptation via soft sweeps or polygenic adaptation. Empirical evidence of this process in nature is currently limited but accumulating.ResultsWe provide genome-wide analyses of rapid evolution in Fundulus heteroclitus populations subjected to recently elevated temperatures due to coastal power station thermal effluents using 5449 SNPs across two effluent-affected and four reference populations. Bayesian and multivariate analyses of population genomic structure reveal a substantial portion of genetic variation that is most parsimoniously explained by selection at the site of thermal effluents. An FST outlier approach in conjunction with additional conservative requirements identify significant allele frequency differentiation that exceeds neutral expectations among exposed and closely related reference populations. Genomic variation patterns near these candidate loci reveal that individuals living near thermal effluents have rapidly evolved from the standing genetic variation through small allele frequency changes at many loci in a pattern consistent with polygenic selection on the standing genetic variation.ConclusionsWhile the ultimate trajectory of selection in these populations is unknown and we survey only a minority of genomic loci, our findings suggest that polygenic models of adaptation may play important roles in large, natural populations experiencing recent selection due to environmental changes that cause broad physiological impacts.

Genomic epidemiology of the commercially important pathogen Renibacterium salmoninarum within the Chilean salmon industry.

Renibacterium salmoninarum is the causative agent of bacterial kidney disease (BKD), which is a commercially important disease of farmed salmonids. Typing by conventional methods provides limited information on the evolution and spread of this pathogen, as there is a low level of standing variation within the R. salmoninarum population. Here, we apply whole-genome sequencing to 42 R. salmoninarum isolates from Chile, primarily from salmon farms, in order to understand the epidemiology of BKD in this country. The patterns of genomic variation are consistent with multiple introductions to Chile, followed by rapid dissemination over a 30 year period. The estimated dates of introduction broadly coincide with major events in the development of the Chilean aquaculture industry. We find evidence for significant barriers to transmission of BKD in the Chilean salmon production chain that may also be explained by previously undescribed signals of host tropism in R. salmoninarum. Understanding the genomic epidemiology of BKD can inform disease intervention and improve sustainability of the economically important salmon industry. This article contains data hosted by Microreact.

A flicker of hope: Genomic data distinguish Northern Flicker taxa despite low levels of divergence

ABSTRACT Next-generation sequencing technologies are increasingly being employed to explore patterns of genomic variation in avian taxa previously characterized using morphology and/or traditional genetic markers. The hybridization dynamics of the Northern Flicker complex have received considerable attention, primarily due to the conspicuous plumage differences among these birds and the geographically extensive hybrid zone between the Red-shafted (Colaptes auratus cafer) and Yellow-shafted (Colaptes auratus auratus) flickers in the Great Plains region of North America. However, no traditional molecular techniques have been able to differentiate these 2 morphologically well-defined taxa from one another, or conclusively from the closely related Gilded Flicker (Colaptes chrysoides). Here, we use a next-generation sequencing approach to assess the genetic diversity and evolutionary history of these 3 taxa. We confirm the overall low levels of differentiation found using traditional molecular markers, but are...

Complex signatures of genomic variation of two non-model marine species in a homogeneous environment

BackgroundGenomic tools are increasingly being used on non-model organisms to provide insights into population structure and variability, including signals of selection. However, most studies are carried out in regions with distinct environmental gradients or across large geographical areas, in which local adaptation is expected to occur. Therefore, the focus of this study is to characterize genomic variation and selective signals over short geographic areas within a largely homogeneous region. To assess adaptive signals between microhabitats within the rocky shore, we compared genomic variation between the Cape urchin (Parechinus angulosus), which is a low to mid-shore species, and the Granular limpet (Scutellastra granularis), a high shore specialist.ResultsUsing pooled restriction site associated DNA (RAD) sequencing, we described patterns of genomic variation and identified outlier loci in both species. We found relatively low numbers of outlier SNPs within each species, and identified outlier genes associated with different selective pressures than those previously identified in studies conducted over larger environmental gradients. The number of population-specific outlier loci differed between species, likely owing to differential selective pressures within the intertidal environment. Interestingly, the outlier loci were highly differentiated within the two northernmost populations for both species, suggesting that unique evolutionary forces are acting on marine invertebrates within this region.ConclusionsOur study provides a background for comparative genomic studies focused on non-model species, as well as a baseline for the adaptive potential of marine invertebrates along the South African west coast. We also discuss the caveats associated with Pool-seq and potential biases of sequencing coverage on downstream genomic metrics. The findings provide evidence of species-specific selective pressures within a homogeneous environment, and suggest that selective forces acting on small scales are just as crucial to acknowledge as those acting on larger scales. As a whole, our findings imply that future population genomic studies should expand from focusing on model organisms and/or studying heterogeneous regions to better understand the evolutionary processes shaping current and future biodiversity patterns, particularly when used in a comparative phylogeographic context.

Patterns of geographic variation of thermal adapted candidate genes in Drosophila subobscura sex chromosome arrangements

BackgroundThe role of chromosomal arrangements in adaptation is supported by the repeatable clinal variation in inversion frequencies across continents in colonizing species such as Drosophila subobscura. However, there is a lack of knowledge on the genetic variation in genes within inversions, possibly targets of climatic selection, across a geographic latitudinal gradient. In the present study we analysed four candidate loci for thermal adaptation, located close to the breakpoints, in two chromosomal arrangements of the sex (A) chromosome of Drosophila subobscura with different thermal preferences. Individual chromosomes with A2 (the inverted arrangement considered warm adapted) or AST (the standard ancestral arrangement considered cold adapted) were sequenced across four European localities at varying latitudes, up to ~ 2500 Kms apart.ResultsImportantly, we found very low differentiation for each specific arrangement across populations as well as no clinal patterns of genomic variation. This suggests wide gene exchange along the cline. Differentiation between the sex chromosome arrangements was significant in the two more proximal regions relative to the AST orientation but not in the distal ones, independently of their location inside or outside the inversion. This can be possibly due to variation in the levels of gene flux and/or selection acting in these regions.ConclusionsGene flow appears to have homogenized the genetic content within-arrangement at a wide geographical scale, despite the expected diverse selective pressures in the specific natural environments of the different populations sampled. It is thus likely that the inversion frequency clines in this species are being maintained by local adaptation in face of gene flow. The differences between arrangements at non-coding regions might be associated with the previously observed differential gene expression in different thermal regimes. Higher resolution genomic scans for individual chromosomal arrangements performed over a large environmental gradient are needed to find the targets of selection and further elucidate the adaptive mechanisms maintaining chromosomal inversion polymorphisms.

Patterns of genomic variation in Chinese maize inbred lines and implications for genetic improvement.

Genetic relationships among Chinese maize germplasms reveal historical trends in heterotic patterns from Chinese breeding programs and identify line Dan340 as a potential genome donor for elite inbred line Zheng58. The characterization of the genetic relationships, heterotic patterns and breeding history of lines in maize breeding programs allows breeders to efficiently use maize germplasm for line improvement over time. In this study, 269 temperate inbred lines, most of which have been widely used in Chinese maize breeding programs since the 1970s, were genotyped using the Illumina MaizeSNP50 BeadChip, which contains 56,110 single-nucleotide polymorphisms. The STRUCTURE analysis, cluster analysis and principal coordinate analysis results consistently revealed seven groups, of which five were consistent with known heterotic groups within the Chinese maize germplasm-Domestic Reid, Lancaster, Zi330, Tang SPT and Tem-tropic I (also known as "P"). These genetic relationships also allowed us to determine the historical trends in heterotic patterns during the three decades from 1970 to 2000, represented by Mo17 from Lancaster, HuangZaoSi (HZS) from Tang SPT, Ye478 from Domestic Reid and P178 from Tem-tropic I heterotic groups. Mo17-related commercial hybrids were widely used in the 1970s and 1980s, followed by the release of HZS- and Ye478-related commercial hybrids in the 1980s and 1990s, and the introduction of Tem-tropic I group in the 1990s and 2000s. Additionally, we identified inbred line Dan340 as a potential genome donor for Zheng58, which is the female parent of the most widely grown commercial hybrid ZhengDan958 in China. We also reconstructed the recombination events of elite line HZS and its 14 derived lines. These findings provide useful information to direct future maize breeding efforts.

Demography and selection shape transcriptomic divergence in field crickets.

Gene flow, demography, and selection can result in similar patterns of genomic variation and disentangling their effects is key to understanding speciation. Here, we assess transcriptomic variation to unravel the evolutionary history of Gryllus rubens and Gryllus texensis, cryptic field cricket species with highly divergent mating behavior. We infer their demographic history and screen their transcriptomes for footprints of selection in the context of the inferred demography. We find strong support for a long history of bidirectional gene flow, which ceased during the late Pleistocene, and a bottleneck in G. rubens consistent with a peripatric origin of this species. Importantly, the demographic history has likely strongly shaped patterns of genetic differentiation (empirical FST distribution). Concordantly, FST -based selection detection uncovers a large number of outliers, likely comprising many false positives, echoing recent theoretical insights. Alternative genetic signatures of positive selection, informed by the demographic history of the sibling species, highlighted a smaller set of loci; many of these are candidates for controlling variation in mating behavior. Our results underscore the importance of demography in shaping overall patterns of genetic divergence and highlight that examining both demography and selection facilitates a more complete understanding of genetic divergence during speciation.

Patterns of Genomic Variation in the Opportunistic Pathogen Candida glabrata Suggest the Existence of Mating and a Secondary Association with Humans

SummaryCandida glabrata is an opportunistic fungal pathogen that ranks as the second most common cause of systemic candidiasis. Despite its genus name, this yeast is more closely related to the model yeast Saccharomyces cerevisiae than to other Candida pathogens, and hence its ability to infect humans is thought to have emerged independently. Moreover, C. glabrata has all the necessary genes to undergo a sexual cycle but is considered an asexual organism due to the lack of direct evidence of sexual reproduction. To reconstruct the recent evolution of this pathogen and find footprints of sexual reproduction, we assessed genomic and phenotypic variation across 33 globally distributed C. glabrata isolates. We cataloged extensive copy-number variation, which particularly affects genes encoding cell-wall-associated proteins, including adhesins. The observed level of genetic variation in C. glabrata is significantly higher than that found in Candida albicans. This variation is structured into seven deeply divergent clades, which show recent geographical dispersion and large within-clade genomic and phenotypic differences. We show compelling evidence of recent admixture between differentiated lineages and of purifying selection on mating genes, which provides the first evidence for the existence of an active sexual cycle in this yeast. Altogether, our data point to a recent global spread of previously genetically isolated populations and suggest that humans are only a secondary niche for this yeast.

A Large Panel of Drosophila simulans Reveals an Abundance of Common Variants.

The rapidly expanding availability of large NGS data sets provides an opportunity to investigate population genetics at an unprecedented scale. Drosophila simulans is the sister species of the model organism Drosophila melanogaster, and is often presumed to share similar demographic history. However, previous population genetic and ecological work suggests very different signatures of selection and demography. Here, we sequence a new panel of 170 inbred genotypes of a North American population of D. simulans, a valuable complement to the DGRP and other D. melanogaster panels. We find some unexpected signatures of demography, in the form of excess intermediate frequency polymorphisms. Simulations suggest that this is possibly due to a recent population contraction and selection. We examine the outliers in the D. simulans genome determined by a haplotype test to attempt to parse the contribution of demography and selection to the patterns observed in this population. Untangling the relative contribution of demography and selection to genomic patterns of variation is challenging, however, it is clear that although D. melanogaster was thought to share demographic history with D. simulans different forces are at work in shaping genomic variation in this population of D. simulans.

Genomic footprints of dryland stress adaptation in Egyptian fat-tail sheep and their divergence from East African and western Asia cohorts

African indigenous sheep are classified as fat-tail, thin-tail and fat-rump hair sheep. The fat-tail are well adapted to dryland environments, but little is known on their genome profiles. We analyzed patterns of genomic variation by genotyping, with the Ovine SNP50K microarray, 394 individuals from five populations of fat-tail sheep from a desert environment in Egypt. Comparative inferences with other East African and western Asia fat-tail and European sheep, reveal at least two phylogeographically distinct genepools of fat-tail sheep in Africa that differ from the European genepool, suggesting separate evolutionary and breeding history. We identified 24 candidate selection sweep regions, spanning 172 potentially novel and known genes, which are enriched with genes underpinning dryland adaptation physiology. In particular, we found selection sweeps spanning genes and/or pathways associated with metabolism; response to stress, ultraviolet radiation, oxidative stress and DNA damage repair; activation of immune response; regulation of reproduction, organ function and development, body size and morphology, skin and hair pigmentation, and keratinization. Our findings provide insights on the complexity of genome architecture regarding dryland stress adaptation in the fat-tail sheep and showcase the indigenous stocks as appropriate genotypes for adaptation planning to sustain livestock production and human livelihoods, under future climates.

Patterns of genomic variation in Coho salmon following reintroduction to the interior Columbia River.

Coho salmon were extirpated in the mid‐20th century from the interior reaches of the Columbia River but were reintroduced with relatively abundant source stocks from the lower Columbia River near the Pacific coast. Reintroduction of Coho salmon to the interior Columbia River (Wenatchee River) using lower river stocks placed selective pressures on the new colonizers due to substantial differences with their original habitat such as migration distance and navigation of six additional hydropower dams. We used restriction site‐associated DNA sequencing (RAD‐seq) to genotype 5,392 SNPs in reintroduced Coho salmon in the Wenatchee River over four generations to test for signals of temporal structure and adaptive variation. Temporal genetic structure among the three broodlines of reintroduced fish was evident among the initial return years (2000, 2001, and 2002) and their descendants, which indicated levels of reproductive isolation among broodlines. Signals of adaptive variation were detected from multiple outlier tests and identified candidate genes for further study. This study illustrated that genetic variation and structure of reintroduced populations are likely to reflect source stocks for multiple generations but may shift over time once established in nature.

Genomic variation in microbial populations inhabiting the marine subseafloor at deep-sea hydrothermal vents

Little is known about evolutionary drivers of microbial populations in the warm subseafloor of deep-sea hydrothermal vents. Here we reconstruct 73 metagenome-assembled genomes (MAGs) from two geochemically distinct vent fields in the Mid-Cayman Rise to investigate patterns of genomic variation within subseafloor populations. Low-abundance populations with high intra-population diversity coexist alongside high-abundance populations with low genomic diversity, with taxonomic differences in patterns of genomic variation between the mafic Piccard and ultramafic Von Damm vent fields. Populations from Piccard are significantly enriched in nonsynonymous mutations, suggesting stronger purifying selection in Von Damm relative to Piccard. Comparison of nine Sulfurovum MAGs reveals two high-coverage, low-diversity MAGs from Piccard enriched in unique genes related to the cellular membrane, suggesting these populations were subject to distinct evolutionary pressures that may correlate with genes related to nutrient uptake, biofilm formation, or viral invasion. These results are consistent with distinct evolutionary histories between geochemically different vent fields, with implications for understanding evolutionary processes in subseafloor microbial populations.