Identification of a Unique Genomic Region in Sweet Chestnut (Castanea sativa Mill.) That Controls Resistance to Asian Chestnut Gall Wasp Dryocosmus kuriphilus Yasumatsu.

The Asian chestnut gall wasp (ACGW) (Hymenoptera Dryocosmus kuriphilus Yasumatsu) is a severe pest of sweet chestnut (Castanea sativa Mill.) with a strong impact on growth and nut production. A comparative field trial in Central Italy, including provenances from Spain, Italy, and Greece, was screened for ACGW infestation over consecutive years. The Greek provenance Hortiatis expressed a high proportion of immune plants and was used to perform a genome-wide association study based on DNA pool sequencing (Pool-GWAS) by comparing two DNA pools from 25 susceptible and 25 resistant plants. DNA pools were sequenced with 50X coverage depth. Sequence reads were aligned to a C. mollissima reference genome and the pools were compared to identify SNPs associated with resistance. Twenty-one significant SNPs were identified and highlighted a small genomic region on pseudochromosome 3 (Chr 3), containing 12 candidate genes of three gene families: Cytochrome P450, UDP-glycosyltransferase, and Rac-like GTP-binding protein. Functional analyses revealed a putative metabolic gene cluster related to saccharide biosynthesis in the genomic regions associated with resistance that could be involved in the production of a toxic metabolite against parasites. The comparison with previous genetic studies confirmed the involvement of Chr 3 in the control of resistance to ACGW.

Rapid vertebrate speciation via isolation, bottlenecks, and drift.

Speciation is often driven by selective processes like those associated with viability, mate choice, or local adaptation, and "speciation genes" have been identified in many eukaryotic lineages. In contrast, neutral processes are rarely considered as the primary drivers of speciation, especially over short evolutionary timeframes. Here, we describe a rapid vertebrate speciation event driven primarily by genetic drift. The White Sands pupfish (Cyprinodon tularosa) is endemic to New Mexico's Tularosa Basin where the species is currently managed as two Evolutionarily significant units (ESUs) and is of international conservation concern (Endangered). Whole-genome resequencing data from each ESU showed remarkably high and uniform levels of differentiation across the entire genome (global FST ≈ 0.40). Despite inhabiting ecologically dissimilar springs and streams, our whole-genome analysis revealed no discrete islands of divergence indicative of strong selection, even when we focused on an array of candidate genes. Demographic modeling of the joint allele frequency spectrum indicates the two ESUs split only ~4 to 5 kya and that both ESUs have undergone major bottlenecks within the last 2.5 millennia. Our results indicate the genome-wide disparities between the two ESUs are not driven by divergent selection but by neutral drift due to small population sizes, geographic isolation, and repeated bottlenecks. While rapid speciation is often driven by natural or sexual selection, here we show that isolation and drift have led to speciation within a few thousand generations. We discuss these evolutionary insights in light of the conservation management challenges they pose.

Adaptation to seasonal reproduction and environment-associated factors drive temporal and spatial differentiation in northwest Atlantic herring despite gene flow.

Understanding how marine organisms adapt to local environments is crucial for predicting how populations will respond to global climate change. The genomic basis, environmental factors and evolutionary processes involved in local adaptation are however not well understood. Here we use Atlantic herring, an abundant, migratory and widely distributed marine fish with substantial genomic resources, as a model organism to evaluate local adaptation. We examined genomic variation and its correlation with environmental variables across a broad environmental gradient, for 15 spawning aggregations in Atlantic Canada and the United States. We then compared our results with available genomic data of northeast Atlantic populations. We confirmed that population structure lies in a fraction of the genome including likely adaptive genetic variants of functional importance. We discovered 10 highly differentiated genomic regions distributed across four chromosomes. Nine regions show strong association with seasonal reproduction. One region, corresponding to a known inversion on chromosome 12, underlies a latitudinal pattern discriminating populations north and south of a biogeographic transition zone on the Scotian Shelf. Genome-environment associations indicate that winter seawater temperature best correlates with the latitudinal pattern of this inversion. The variation at two so-called 'islands of divergence' related to seasonal reproduction appear to be private to the northwest Atlantic. Populations in the northwest and northeast Atlantic share variation at four of these divergent regions, simultaneously displaying significant diversity in haplotype composition at another four regions, which includes an undescribed structural variant approximately 7.7 Mb long on chromosome 8. Our results suggest that the timing and geographic location of spawning and early development may be under diverse selective pressures related to allelic fitness across environments. Our study highlights the role of genomic architecture, ancestral haplotypes and selection in maintaining adaptive divergence in species with large population sizes and presumably high gene flow.

Estimating microhaplotype allele frequencies from low-coverage or pooled sequencing data.

BACKGROUND: Microhaplotypes have the potential to be more cost-effective than SNPs for applications that require genetic panels of highly variable loci. However, development of microhaplotype panels is hindered by a lack of methods for estimating microhaplotype allele frequency from low-coverage whole genome sequencing or pooled sequencing (pool-seq) data. RESULTS: We developed new methods for estimating microhaplotype allele frequency from low-coverage whole genome sequence and pool-seq data. We validated these methods using datasets from three non-model organisms. These methods allowed estimation of allele frequency and expected heterozygosity at depths routinely achieved from pooled sequencing. CONCLUSIONS: These new methods will allow microhaplotype panels to be designed using low-coverage WGS and pool-seq data to discover and evaluate candidate loci. The python script implementing the two methods and documentation are available at https://www.github.com/delomast/mhFromLowDepSeq .

Integrating Pool-seq uncertainties into demographic inference.

Next-generation sequencing of pooled samples (Pool-seq) is a popular method to assess genome-wide diversity patterns in natural and experimental populations. However, Pool-seq is associated with specific sources of noise, such as unequal individual contributions. Consequently, using Pool-seq for the reconstruction of evolutionary history has remained underexplored. Here we describe a novel Approximate Bayesian Computation (ABC) method to infer demographic history, explicitly modelling Pool-seq sources of error. By jointly modelling Pool-seq data, demographic history and the effects of selection due to barrier loci, we obtain estimates of demographic history parameters accounting for technical errors associated with Pool-seq. Our ABC approach is computationally efficient as it relies on simulating subsets of loci (rather than the whole-genome) and on using relative summary statistics and relative model parameters. Our simulation study results indicate Pool-seq data allows distinction between general scenarios of ecotype formation (single versus parallel origin) and to infer relevant demographic parameters (e.g. effective sizes and split times). We exemplify the application of our method to Pool-seq data from the rocky-shore gastropod Littorina saxatilis, sampled on a narrow geographical scale at two Swedish locations where two ecotypes (Wave and Crab) are found. Our model choice and parameter estimates show that ecotypes formed before colonization of the two locations (i.e. single origin) and are maintained despite gene flow. These results indicate that demographic modelling and inference can be successful based on pool-sequencing using ABC, contributing to the development of suitable null models that allow for a better understanding of the genetic basis of divergent adaptation.

Characterization of the male-specific region containing the candidate sex-determining gene in Amur catfish (Silurus asotus) using third-generation- and pool-sequencing data.

Amur catfish (Silurus asotus) is an ecologically and economically important fish species in Asia. Here, we assembled the female and male Amur catfish genomes, with genome sizes of 757.15 and 755.44 Mb, respectively, at the chromosome level using nanopore and Hi-C technologies. Consistent with the known diploid chromosome count, both genomes contained 29 chromosome-size scaffolds covering 98.80 and 98.73 % of the complete haplotypic assembly with scaffold N50 of 28.87 and 27.29 Mb, respectively. The female (n = 40) and male (n = 40) pools were re-sequenced. Comparative analysis of sequencing and re-sequencing data from both sexes confirmed the presence of an XX/XY sex determination system in Amur catfish and revealed Chr5 as the sex chromosome containing an approximately 400 kb Y-specific region (MSY). Gene annotation revealed a male-specific duplicate of amhr2, namely amhr2y, in MSY, which is male-specific in different wild populations and expressed only in the testes. Amur catfish shared partially syntenic MSY and amhr2y genes with the southern catfish (S. meridionalis, Chr24), which were located on different chromosomes. High sequence divergence between amhr2y and amhr2 and high sequence similarity with amhr2y were observed in both species. These results indicate the common origin of the sex-determining (SD) gene and transition of amhr2y in the two Silurus species. Accumulation of repetitive elements in the MSY of both species may be the main driver of the transition of amhr2y. Overall, our study provides valuable catfish genomic resources. Moreover, determination of amhr2y as the candidate SD gene in Amur catfish provides another example of amhr2 as the SD gene in fish.

Genetic incompatibilities in reciprocal hybrids between populations of Tigriopus californicus with low to moderate mitochondrial sequence divergence.

All mitochondrial-encoded proteins and RNAs function through interactions with nuclear-encoded proteins, which are critical for mitochondrial performance and eukaryotic fitness. Coevolution maintains inter-genomic (i.e., mitonuclear) compatibility within a taxon, but hybridization can disrupt coevolved interactions, resulting in hybrid breakdown. Thus, mitonuclear incompatibilities may be important mechanisms underlying reproductive isolation and, potentially, speciation. Here we utilize Pool-seq to assess the effects of mitochondrial genotype on nuclear allele frequencies in fast- and slow-developing reciprocal inter-population F2 hybrids between relatively low-divergence populations of the intertidal copepod Tigriopus californicus. We show that mitonuclear interactions lead to elevated frequencies of coevolved (i.e., maternal) nuclear alleles on two chromosomes in crosses between populations with 1.5% or 9.6% fixed differences in mitochondrial DNA nucleotide sequence. However, we also find evidence of excess mismatched (i.e., noncoevolved) alleles on three or four chromosomes per cross, respectively, and of allele frequency differences consistent with effects involving only nuclear loci (i.e., unaffected by mitochondrial genotype). Thus, our results for low-divergence crosses suggest an underlying role for mitonuclear interactions in variation in hybrid developmental rate, but despite substantial effects of mitonuclear coevolution on individual chromosomes, no clear bias favoring coevolved interactions overall.

Ace and ace-like genes of invasive redlegged earth mite: copy number variation, target-site mutations, and their associations with organophosphate insensitivity.

BACKGROUND: Invasive Australian populations of redlegged earth mite, Halotydeus destructor (Tucker), are evolving increasing organophosphate resistance. In addition to the canonical ace gene, the target gene of organophosphates, the H. destructor genome contains many radiated ace-like genes that vary in copy number and amino acid sequence. In this work, we characterise copy number and target-site mutation variation at the canonical ace and ace-like genes and test for potential associations with organophosphate insensitivity. This was achieved through comparisons of whole-genome pool-seq data from alive and dead mites following organophosphate exposure. RESULTS: A combination of increased copy number and target-site mutations at the canonical ace was associated with organophosphate insensitivity in H. destructor. Resistant populations were segregating for G119S, A201S, F331Y at the canonical ace. A subset of populations also had copy numbers of canonical ace > 2, which potentially helps overexpress proteins carrying these target-site mutations. Haplotypes possessing different copy numbers and target-site mutations of the canonical ace gene may be under selection across H. destructor populations. We also detected some evidence that increases in copy number of radiated ace-like genes are associated with organophosphate insensitivity, which might suggest potential roles in sequestration or breakdown of organophosphates. CONCLUSION: Different combinations of target-site mutations and (or) copy number variation in the canonical ace and ace-like genes may provide non-convergent ways for H. destructor to respond to organophosphate selection. However, these changes may only play a partial role in organophosphate insensitivity, which appears to have a polygenic architecture. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Management and conservation implications of cryptic population substructure for two commercially exploited fishes (Merluccius spp.) in southern Africa.

Genomic information can aid in the establishment of sustainable management plans for commercially exploited marine fishes, aiding in the long-term conservation of these resources. The southern African hakes (Merluccius capensis and M. paradoxus) are commercially valuable demersal fishes with similar distribution ranges but exhibiting contrasting life histories. Using a comparative framework based on Pool-Seq genome-wide SNP data, we investigated whether the evolutionary processes that shaped extant patterns of diversity and divergence are shared among these two congeneric fishes, or unique to each one. Our findings revealed that M. capensis and M. paradoxus show similar levels of genome-wide diversity, despite different census sizes and life-history features. In addition, M. capensis shows three highly structured geographic populations across the Benguela Current region (one in the northern Benguela and two in the southern Benguela), with no consistent genome-environment associations detected. In contrast, although population structure and outlier analyses suggested panmixia for M. paradoxus, reconstruction of its demographic history suggested the presence of an Atlantic-Indian Ocean subtle substructuring pattern. Therefore, it appears that M. paradoxus might be composed by two highly connected populations, one in the Atlantic and one in the southwest Indian Ocean. The reported similar low levels of genomic diversity, as well as newly discovered genetically distinct populations in both hake species can thus assist in informing and improving conservation and management plans for the commercially important southern African Merluccius.

Genomic Architecture of Hybrid Male Sterility in a Species Without Sex Chromosomes (Tigriopus californicus, Copepoda: Harpacticoida).

Sterility among hybrids is one of the most prevalent forms of reproductive isolation delineating species boundaries and is expressed disproportionately in heterogametic XY males. While hybrid male sterility (HMS) due to the "large X effect" is a well-recognized mechanism of reproductive isolation, it is less clear how HMS manifests in species that lack heteromorphic sex chromosomes. We evaluated differences in allele frequencies at approximately 460,000 SNPs between fertile and sterile F2 interpopulation male hybrids to characterize the genomic architecture of HMS in a species without sex chromosomes (Tigriopus californicus). We tested associations between HMS and mitochondrial-nuclear and/or nuclear-nuclear signatures of incompatibility. Genomic regions associated with HMS were concentrated on a single chromosome with the same primary 2-Mbp regions identified in one pair of reciprocal crosses. Gene Ontology analysis revealed that annotations associated with spermatogenesis were the most overrepresented within the implicated region, with nine protein-coding genes connected with this process found in the quantitative trait locus of chromosome 2. Our results indicate that a narrow genomic region was associated with the sterility of male hybrids in T. californicus and suggest that incompatibilities among select nuclear loci may replace the large X effect when sex chromosomes are absent.

A major 6 Mb superlocus is involved in pyrethroid resistance in the common bed bug Cimex lectularius.

In the last few years, the bed bug Cimex lectularius has been an increasing problem worldwide, mainly due to the development of insecticide resistance to pyrethroids. The characterization of resistance alleles is a prerequisite to improve surveillance and resistance management. To identify genomic variants associated with pyrethroid resistance in Cimex lectularius, we compared the genetic composition of two recent and resistant populations with that of two ancient-susceptible strains using a genome-wide pool-seq design. We identified a large 6 Mb "superlocus" showing particularly high genetic differentiation and association with the resistance phenotype. This superlocus contained several clustered resistance genes and was also characterized by a high density of structural variants (inversions, duplications). The possibility that this superlocus constitutes a resistance "supergene" that evolved after the clustering of alleles adapted to insecticide and after reduction in recombination is discussed.

When One's Not Enough: Colony Pool-Seq Outperforms Individual-Based Methods for Assessing Introgression in Apis mellifera mellifera.

The human management of honey bees (Apis mellifera) has resulted in the widespread introduction of subspecies outside of their native ranges. One well known example of this is Apis mellifera mellifera, native to Northern Europe, which has now been significantly introgressed by the introduction of C lineage honey bees. Introgression has consequences for species in terms of future adaptive potential and long-term viability. However, estimating introgression in colony-living haplodiploid species is challenging. Previous studies have estimated introgression using individual workers, individual drones, multiple drones, and pooled workers. Here, we compare introgression estimates via three genetic approaches: SNP array, individual RAD-seq, and pooled colony RAD-seq. We also compare two statistical approaches: a maximum likelihood cluster program (ADMIXTURE) and an incomplete lineage sorting model (ABBA BABA). Overall, individual approaches resulted in lower introgression estimates than pooled colonies when using ADMIXTURE. However, the pooled colony ABBA BABA approach resulted in generally lower introgression estimates than all three ADMIXTURE estimates. These results highlight that sometimes one individual is not enough to assess colony-level introgression, and future studies that do use colony pools should not be solely dependent on clustering programs for introgression estimates.

Balancing selection and candidate loci for survival and growth during larval development in the Mediterranean mussel, Mytilus galloprovincialis.

Many marine bivalves have complex life histories with distinct developmental processes and genetic mechanisms. Larval development for most bivalves is often a prolonged and crucial physiological stage, where they suffer mass mortality due to early-acting genetic load. In this study, we describe genetic changes taking place within a single generation of families of the Mediterranean mussel Mytilus galloprovincialis over 23 days of larval development. Using replicated cultures and a pooled sequencing approach, we demonstrate that temporal balancing selection at the majority of loci preserve genetic variation in the early developmental stages of M. galloprovincialis. Balancing selection may be the mechanism which maintains standing genetic variation within the mussel genome and may improve the chances of survival and shield larvae from high levels of genetic load. Additionally, we used changes in allele frequencies to identify potential size-associated SNPs and viability-associated SNPs and found that patterns of genetic changes in directionally selected SNPs cannot be simply explained by traditional theories of genetic purging or directional selection without consideration of balancing selection. Finally, we observed a negative correlation between larval growth rates and survival, implying a potential trade-off relationship between the 2 commercially relevant phenotypes.

Whole Genome Sequencing Provides Information on the Genomic Architecture and Diversity of Cultivated Gilthead Seabream (Sparus aurata) Broodstock Nuclei.

The gilthead seabream (Sparus aurata) is a species of relevance for the Mediterranean aquaculture industry. Despite the advancement of genetic tools for the species, breeding programs still do not often include genomics. In this study, we designed a genomic strategy to identify signatures of selection and genomic regions of high differentiation among populations of farmed fish stocks. A comparative DNA pooling sequencing approach was applied to identify signatures of selection in gilthead seabream from the same hatchery and from different nuclei that had not been subjected to genetic selection. Identified genomic regions were further investigated to detect SNPs with predicted high impact. The analyses underlined major genomic differences in the proportion of fixed alleles among the investigated nuclei. Some of these differences highlighted genomic regions, including genes involved in general metabolism and development already detected in QTL for growth, size, skeletal deformity, and adaptation to variation of oxygen levels in other teleosts. The obtained results pointed out the need to control the genetic effect of breeding programs in this species to avoid the reduction of genetic variability within populations and the increase in inbreeding level that, in turn, might lead to an increased frequency of alleles with deleterious effects.

Partial reuse of circadian clock genes along parallel clines of diapause in two moth species.

Understanding the molecular basis of repeated evolution improves our ability to predict evolution across the tree of life. Only since the last decade has high-throughput sequencing enabled comparative genome scans to thoroughly examine the repeatability of genetic changes driving repeated phenotypic evolution. The Asian corn borer (ACB), Ostrinia furnacalis (Guenée), and the European corn borer (ECB), Ostrinia nubilalis (Hübner), are two closely related moths displaying repeatable phenological adaptation to a wide range of climates on two separate continents, largely manifesting as changes in the timing of diapause induction and termination across latitude. Candidate genes underlying diapause variation in North American ECB have been previously identified. Here, we sampled seven ACB populations across 23 degrees of latitude in China to elucidate the genetic basis of diapause variation and evolutionary mechanisms driving parallel clinal responses in the two species. Using pooled whole-genome sequencing (Pool-seq) data, population genomic analyses revealed hundreds of single nucleotide polymorphisms (SNP) whose allele frequencies covaried with mean diapause phenotypes along the cline. Genes involved in circadian rhythm were over-represented among candidate genes with strong signatures of spatially varying selection. Only one of two circadian clock genes associated with diapause evolution in ECB showed evidence of reuse in ACB (period [per]), but per alleles were not shared between species nor with their outgroup, implicating independent mutational paths. Nonetheless, evidence of adaptive introgression was discovered at putative diapause loci located elsewhere in the genome, suggesting that de novo mutations and introgression might both underlie the repeated phenological evolution.

Bait-ER: A Bayesian method to detect targets of selection in Evolve-and-Resequence experiments.

For over a decade, experimental evolution has been combined with high-throughput sequencing techniques. In so-called Evolve-and-Resequence (E&R) experiments, populations are kept in the laboratory under controlled experimental conditions where their genomes are sampled and allele frequencies monitored. However, identifying signatures of adaptation in E&R datasets is far from trivial, and it is still necessary to develop more efficient and statistically sound methods for detecting selection in genome-wide data. Here, we present Bait-ER - a fully Bayesian approach based on the Moran model of allele evolution to estimate selection coefficients from E&R experiments. The model has overlapping generations, a feature that describes several experimental designs found in the literature. We tested our method under several different demographic and experimental conditions to assess its accuracy and precision, and it performs well in most scenarios. Nevertheless, some care must be taken when analysing trajectories where drift largely dominates and starting frequencies are low. We compare our method with other available software and report that ours has generally high accuracy even for trajectories whose complexity goes beyond a classical sweep model. Furthermore, our approach avoids the computational burden of simulating an empirical null distribution, outperforming available software in terms of computational time and facilitating its use on genome-wide data. We implemented and released our method in a new open-source software package that can be accessed at https://doi.org/10.5281/zenodo.7351736.

The Dynamics of Adaptation to Stress from Standing Genetic Variation and de novo Mutations.

Adaptation from standing genetic variation is an important process underlying evolution in natural populations, but we rarely get the opportunity to observe the dynamics of fitness and genomic changes in real time. Here, we used experimental evolution and Pool-Seq to track the phenotypic and genomic changes of genetically diverse asexual populations of the yeast Saccharomyces cerevisiae in four environments with different fitness costs. We found that populations rapidly and in parallel increased in fitness in stressful environments. In contrast, allele frequencies showed a range of trajectories, with some populations fixing all their ancestral variation in <30 generations and others maintaining diversity across hundreds of generations. We detected parallelism at the genomic level (involving genes, pathways, and aneuploidies) within and between environments, with idiosyncratic changes recurring in the environments with higher stress. In particular, we observed a tendency of becoming haploid-like in one environment, whereas the populations of another environment showed low overall parallelism driven by standing genetic variation despite high selective pressure. This work highlights the interplay between standing genetic variation and the influx of de novo mutations in populations adapting to a range of selective pressures with different underlying trait architectures, advancing our understanding of the constraints and drivers of adaptation.

Genomic differentiation in Pacific cod using Pool-Seq.

Patterns of genetic differentiation across the genome can provide insight into selective forces driving adaptation. We used pooled whole genome sequencing, gene annotation, and environmental covariates to evaluate patterns of genomic differentiation and to investigate mechanisms responsible for divergence among proximate Pacific cod (Gadus macrocephalus) populations from the Bering Sea and Aleutian Islands and more distant Washington Coast cod. Samples were taken from eight spawning locations, three of which were replicated to estimate consistency in allele frequency estimation. A kernel smoothing moving weighted average of relative divergence (F ST) identified 11 genomic islands of differentiation between the Aleutian Islands and Bering Sea samples. In some islands of differentiation, there was also elevated absolute divergence (d XY) and evidence for selection, despite proximity and potential for gene flow. Similar levels of absolute divergence (d XY) but roughly double the relative divergence (F ST) were observed between the distant Bering Sea and Washington Coast samples. Islands of differentiation were much smaller than the four large inversions among Atlantic cod ecotypes. Islands of differentiation between the Bering Sea and Aleutian Island were associated with SNPs from five vision system genes, which can be associated with feeding, predator avoidance, orientation, and socialization. We hypothesize that islands of differentiation between Pacific cod from the Bering Sea and Aleutian Islands provide evidence for adaptive differentiation despite gene flow in this commercially important marine species.

Discovering candidate SNPs for resilience breeding of red clover.

Red clover is a highly valuable crop for the ruminant industry in the temperate regions worldwide. It also provides multiple environmental services, such as contribution to increased soil fertility and reduced soil erosion. This study used 661 single nucleotide polymorphism (SNP) markers via targeted sequencing using seqSNP, to describe genetic diversity and population structure in 382 red clover accessions. The accessions were selected from NordGen representing red clover germplasm from Norway, Sweden, Finland and Denmark as well as from Lantmännen, a Swedish seed company. Each accession was represented by 10 individuals, which was sequenced as a pool. The mean Nei's standard genetic distance between the accessions and genetic variation within accessions were 0.032 and 0.18, respectively. The majority of the accessions had negative Tajima's D, suggesting that they contain significant proportions of rare alleles. A pairwise FST revealed high genetic similarity between the different cultivated types, while the wild populations were divergent. Unlike wild populations, which exhibited genetic differentiation, there was no clear differentiation among all cultivated types. A principal coordinate analysis revealed that the first principal coordinate, distinguished most of the wild populations from the cultivated types, in agreement with the results obtained using a discriminant analysis of principal components and cluster analysis. Accessions of wild populations and landraces collected from southern and central Scandinavia showed a higher genetic similarity to Lantmännen accessios. It is therefore possible to link the diversity of the environments where wild populations were collected to the genetic diversity of the cultivated and wild gene pools. Additionally, least absolute shrinkage and selection operator (LASSO) models revealed associations between variation in temperature and precipitation and SNPs within genes controlling stomatal opening. Temperature was also related to kinase proteins, which are known to regulate plant response to temperature stress. Furthermore, the variation between wild populations and cultivars was correlated with SNPs within genes regulating root development. Overall, this study comprehensively investigated Nordic European red clover germplasm, and the results provide forage breeders with valuable information for further selection and development of red clover cultivars.

Diversity of sex chromosomes in Sulawesian medaka fishes.

Recent genetic and genomic studies have revealed tremendous diversity in sex chromosomes across diverse taxa. Closely related species with different sex chromosomes provide us excellent opportunities to investigate the driving forces and the consequences of sex chromosome turnover. In the present study, we investigated the diversity of sex chromosomes of 13 Oryzias species from Sulawesi, Indonesia, which diversified during the last 4.86 million years. Using pooled sequencing, we found sex chromosomes in nine species that all had XY systems, with a species being possibly modified by multiple loci. Seven species (O. woworae, O. asinua, O. wolasi, O. matanensis, O. celebensis, O. hadiatyae, and O. dopingdopingensis) share linkage group (LG) 24 as sex chromosomes; however, they differed in the length and magnitude of sequence divergence between the X and Y chromosomes. The sex chromosome of O. eversi was LG4, which has not been reported as a sex chromosome in any other medaka species. In O. sarasinorum, LG16 and LG22 are associated with sex. Although LG16 was found to be sex-linked in another medaka species previously examined, the sex-determining regions did not overlap. No significant signatures for sex chromosomes were identified in the other four species (O. marmoratus, O. nigrimas, O. nebulosus, and O. orthognathus). Frequent turnovers and the great diversity of the sex chromosomes will make Sulawesian medaka species a model system for investigating the driving forces and consequences of sex chromosome turnover.