A divergent haplotype with a large deletion at the berry color locus causes a white-skinned phenotype in grapevine

The current genetic model explaining berry skin color in Vitis vinifera is incomplete and fails to predict berry skin color phenotypes for one allele of VvMybA1, referred to as VvMybA1_SUB. Our study focuses on this specific allele, revealing that the haplotype containing VvMybA1_SUB (haplotype F) represents an ancient lineage of the berry color locus. Within haplotype F, we identified two functional subhaplotypes, HapF1 and HapF2, associated with black-skinned phenotype, and one non-functional subhaplotype, HapFDEL, responsible for white-skinned phenotype. HapF1 likely originated from wild populations domesticated in the Near East and subsequently spread globally with the expansion of viticulture. In contrast, HapF2 has a more restricted distribution and may have emerged from hybridization events between cultivated grapevines and local wild populations as viticulture migrated to the Italian peninsula. Furthermore, we found that in white-skinned berry cultivar, HapF has undergone a large deletion at the berry color locus, removing the majority of the VvMybA genes. Previous works suggested a single common origin for white-skinned varieties during grapevine domestication. Our results challenge this notion, instead proposing that white-skinned grape cultivars arose at least twice during grapevine domestication history. Alongside the major haplotype A, some white-skinned cultivars, such as cv. ‘Sultanina’ harbor HapFDEL. Since HapFDEL is present only in table grape varieties, we suggest that it likely arose from a recent mutational event and dispersed along the ancient Silk Road into East Asia. These findings enhance our understanding of the genetic diversity and evolutionary trajectory of grapevine cultivars, offering insights into their domestication and spread across different geographical regions.

The Genetics of Host Plant Acceptance in Pea Aphids

The evolution of host-associated sympatric populations in phytophagous insects (so called "host races") connects adaptive divergence to barriers to gene flow. Pea aphid (Acyrthosiphon pisum) host races specialise on legume species, and because host plant choice leads to assortative mating, the genetic basis of host plant acceptance is key to understanding speciation. Aphids use smell and taste in their host plant selection. While chemosensory genes frequently emerge as "outliers" in genome scans, their link to plant acceptance behaviour remains unclear. We examined the genetic basis of host-associated phenotypes using an F2 cross between two pea aphid host-associated races (specialised on alfalfa-Medicago sativa- and pea-Pisum sativum), assaying behaviour on both host plants and conducting QTL and regional heritability analyses based on a high-resolution linkage map. We identified five regions of moderate effect associated with acceptance of alfalfa, two with pea acceptance and two with survival on alfalfa. Two QTLs, one for alfalfa and one for pea acceptance, are located within a large rearranged region on chromosome 1, while other QTLs linked to alfalfa acceptance and survival are in the same region on chromosome 3-linking host plant choice to fitness. These findings highlight the polygenic basis of acceptance behaviour and the role of gene clustering and chromosomal rearrangements in promoting coupling among barrier loci. We identified 60 chemosensory genes within regions connected to acceptance, 24 of which were divergent among pea aphid races in previous genome scan or gene expression analyses. Evidence linking these genes to acceptance phenotypes supports their role in determining host plant specificity and as barrier loci contributing to pea aphid speciation.

Allele-based modeling to predict phenological stages of grapevine hybrids under future climatic conditions

KEY MESSAGE

The integration of genetic values in ecophysiological models for phenological stages allows us to predict the effect of loci in future conditions. Modern grapevine (Vitis vinifera L.) breeding programs aim to create new varieties resistant to biotic and abiotic stresses, simultaneously. Developmental stages may affect many physiological processes in grapevine, especially berry composition. The shifts of phenological stages observed in the context of climate change challenge the selection of new varieties. In this paper, we evaluate how genotypes derived from a breeding program aimed at developing disease-resistant varieties may adapt to future climatic conditions. Specifically, we examine the genetic variability of three key phenological periods (February 15 to budbreak, budbreak to flowering, and flowering to véraison) using an ecophysiological model based on thermal requirements. Using high-density genetic information, we identified more than 18 quantitative trait loci for three phenological periods. By combining relevant allelic effects, we virtually constructed both an early and late composite genotype and evaluated their potential adaptation to future climatic conditions, using the greenhouse gas IPCC emissions scenario RCP 8.5 and simulated meteorological data at a local scale. While the early composite genotype may not outperform V. vinifera cv. Chardonnay under these projected conditions, the late composite genotype appears to remain suitably adapted through at least 2060. Our approach enables the prediction of allele-specific advantages on phenological stages across a range of future climate scenarios.

Litchi40K v1.0: a cost-effective, flexible, and versatile liquid SNP chip for genetic analysis and digitalization of germplasm resources in litchi

Genetic breeding and molecular identification in varieties depend on high-performance genotyping tools. The high heterozygosity of the litchi genome contributes to increased resequencing costs and elevated error rates in hybridization-based genotyping methods. In this study, a liquid chip named Litchi40K v1.0 was developed with high-depth resequencing data from 875 litchi samples, and its efficacy was validated across three different populations. In the L. chinensis var. fulvosus population, three subpopulations characterized by spatial distribution, and a total of 1110 genes were identified in the genomic regions with subpopulation differentiation. Additionally, a total of 30 significant signals associated with diverse agronomic traits were identified. The H002 haplotype of LITCHI02696, dominant in the Sub2 subgroup, significantly increased the soluble solid content in the L. chinensis var. fulvosus population. In a hybrid F1 population, a high-density genetic map was constructed and 79 dwarfing-related QTLs were identified with the liquid chip. An NAC transcription factor was identified as a candidate gene with a heterozygous frameshift variant in the male parent. To facilitate the digitization of germplasm resources, 384 SNPs were selected, and the DNA fingerprint map revealed clear genetic relationships and a total of 10 potential synonym groups or instances of bud mutations were identified in 164 main cultivated litchi varieties. This study provides cost-effective, flexible, and versatile liquid chip for genetic analysis and digitalization of germplasm resources in litchi.

Major Facilitator Superfamily transporters balance sugar metabolism in peach

Sugar content is a key determinant of peach (Prunus persica) fruit quality, influencing taste, consumer preferences, and market value. However, the roles of Major Facilitator Superfamily (MFS) transporters in sugar metabolism and regulation remain largely unexplored. This study employed a combination of spatial metabolomics, quantitative genetics, transcriptomics, comparative genomics, and functional genomics to investigate the role of 67 MFS members in balancing sugar metabolism during peach fruit development. Spatial metabolomics revealed dynamic sugar distribution patterns, with ERD6-like transporters (PpERDL16-1) and tonoplastic sugar transporters 1 (PpTST1) promoting sucrose accumulation and polyol/monosaccharide transporters 5 (PpPMT5-1) and sucrose transporters 4 (PpSUT4) reducing sucrose transport during fruit ripening. Functional studies confirmed these roles: PpERDL16-1 overexpression enhanced sucrose transport, and PpPMT5-1 or PpSUT4 silencing reduced sugar levels in peach fruit. Quantitative trait locus (QTL) mapping identified a major locus on chromosome 5, upstream of PpTST1, forming distinct haplotypes (Hap1 and Hap2). Hap1 was associated with lower PpTST1 expression and higher sugar and soluble solids content (SSC), while Hap2 was linked to higher PpTST1 expression and lower sugar content. This inverse relationship suggests that upstream genetic variants fine-tune PpTST1 expression in a context-dependent manner, potentially through interactions with transcription factors or epigenetic modifiers. Notably, PpTST1 overexpression increased sugar content but did not alter SSC, indicating compensatory mechanisms such as changes in organic acid metabolism or water content. These results illuminate the molecular mechanisms regulating sugar homeostasis in peach fruits, providing valuable targets for the genetic improvement of fruit quality through breeding programs.

Discrete genetic modules underlie divergent reproductive strategies in three-spined stickleback

A central challenge in biology is to understand how complex behaviors evolve. Given their importance to fitness, complex behavioral traits often evolve as an integrated package, but it is unclear whether suites of traits evolve through a few pleiotropic genetic changes, each affecting many behaviors, or by accumulating several changes that, when combined, give rise to an entire package of correlated traits. Typically, three-spined stickleback exhibit paternal care, a behavior that characterizes the entire Gasterosteidae family. However, an unusual "white" three-spined stickleback ecotype exhibits a suite of traits associated with the evolutionary loss of paternal care. In the white ecotype, males disperse embryos from their nests rather than care for them, build loose nests, exhibit high rates of courtship, and are relatively small in body size. These behavioral differences are apparent in stickleback reared in a common garden environment, suggesting the differences have a heritable basis. In an F2 intercross, we show that these traits are genetically uncorrelated and map to different genomic regions, suggesting that components of the white reproductive strategy segregate independently and evolved through the addition of multiple genetic changes. These results contribute to the growing body of evidence that the behavioral diversity observed in nature may often evolve by accumulating and combining alleles, each with modular effects, and show that this principle applies to a suite of behavioral traits that together form an integrated and adaptive strategy.

On the origin of an insular hybrid butterfly lineage

A new species can form through hybridization between species. Hybrid speciation in animals has been intensely debated, partly because hard evidence for the process has been difficult to obtain. Here, we report the discovery of a European hybrid butterfly lineage, a finding that can be considered surprising given the intense and long-term study of European butterflies. The lineage we describe is mainly inhabiting an island in the Baltic Sea and was previously designated as a subspecies (horkei) of one of the parental species (Aricia artaxerxes). By analyzing whole-genome resequencing data and developing a novel cluster analysis based on historical recombination events (Fisher junctions), we determine that horkei originated by hybridization between the nonsister species A. artaxerxes and A. agestis. This hybridization event occurred approximately 54,000 years ago, predating the last glaciation of the current distribution range. Horkei must therefore have persisted long enough to be able to colonize its current range, despite that this area lies between the current distributions of the parental species. The hybrid origin, the maintenance of genomic integrity across times of dramatic climate change, and the expression of a combination of parental traits suggest that horkei could be in the process of hybrid speciation.

Independent genetic basis of meiotic crossover positioning and interference in domestic pigs

Meiotic crossover patterning shows huge variation within and between chromosomes, individuals, and species, yet the molecular and evolutionary causes and consequences of this variation remain poorly understood. A key step is to understand the genetic architecture of the crossover rate, positioning, and interference to determine if these factors are governed by common or distinct genetic processes. Here, we investigate individual variation in autosomal crossover count, crossover position (measured as both intra-chromosomal shuffling and distance to telomere), and crossover interference in a large breeding population of domestic pigs (N = 82,474 gametes). We show that all traits are heritable in females at the gamete (h2 = 0.07-0.11) and individual mean levels (h2 = 0.08-0.41). In females, crossover count, and interference are strongly associated with RNF212, but crossover positioning is associated with SYCP2, MEI4, and PRDM9. Our results show that crossover positioning and rate/interference are driven by distinct genetic processes in female pigs and have the capacity to evolve independently.

Meta-Analysis of QTL Mapping and GWAS Reveal Candidate Genes for Heat Tolerance in Small Yellow Croaker, Larimichthys polyactis

High temperatures present considerable challenges to global fish growth and production, yet the genetic basis of heat tolerance remains underexplored. This study combines quantitative trait locus (QTL) mapping and genome-wide association studies (GWAS) to examine the genetics of heat tolerance in Larimichthys polyactis. As a result, a genetic linkage map was constructed with 3237 bin markers spanning 24 linkage groups and totaling 1900.84 centimorgans, using genotyping-by-sequencing of a full-sib family comprising 120 progeny and their two parents. Based on this genetic linkage map, QTL mapping identified four QTLs associated with heat tolerance, which encompassed 18 single nucleotide polymorphisms and harbored 648 genes within the QTL intervals. The GWAS further disclosed 76 candidate genes related to heat tolerance, 56 of which overlapped with the QTL results. Enrichment analysis indicated that these genes are involved in immune response, development, lipid metabolism, and endocrine regulation. qPCR validation of 14 prioritized genes, which were simultaneously enriched in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways, confirmed significant upregulation of smpd5, polr3d, rab11fip2, and gfpt1, along with downregulation of gpat4 and grk5 after 6 h of heat stress. These findings demonstrate their responsiveness to elevated high temperatures. This meta-analysis of QTL mapping and GWAS has successfully identified functional genes related to heat tolerance, enhancing understanding of the genetic architecture underlying this critical trait in L. polyactis. It also provides a molecular breeding tool to improve genetic traits associated with heat tolerance in cultured L. polyactis.

Rapid and dynamic evolution of a giant Y chromosome in Silene latifolia

Some plants have massive sex-linked regions. To test hypotheses about their evolution, we sequenced the genome of Silene latifolia, in which giant heteromorphic sex chromosomes were first discovered in 1923. It has long been known that the Y chromosome consists mainly of a male-specific region that does not recombine with the X chromosome and carries the sex-determining genes and genes with other male functions. However, only with a whole Y chromosome assembly can candidate genes be validated experimentally and their locations determined and related to the suppression of recombination. We describe the genomic changes as the ancestral chromosome evolved into the current XY pair, testing ideas about the evolution of large nonrecombining regions and the mechanisms that created the present recombination pattern.

Masu salmon species complex relationships and sex chromosomes revealed from analyses of the masu salmon (Oncorhynchus masou masou) genome assembly

Masu salmon (Oncorhynchus masou) are the only Pacific salmon endemic to Asia. Some researchers prefer to categorize these salmon into 4 subspecies (masu-Oncorhynchus masou masou, amago-Oncorhynchus masou ishikawae, Biwa-Oncorhynchus masou subsp., and Formosan-Oncorhynchus masou formosanus), while others prefer individual species designations. Even though the masu salmon fishery is thousands of years old, classification of the diversity within the masu salmon species complex remains elusive. In this study, a genetic map and reference genome assembly were generated for 1 species/subspecies (masu) to provide resources for understanding the species complex. In O. m. masou, the sex chromosome was determined to be chromosome 7. Resequenced genomes from 2 other putative subspecies (amago and Biwa) provided evidence that they do not share the same sex chromosome. Principal component and admixture analyses clustered the amago and Biwa salmon close together. This supported previous findings of a close relationship between amago and Biwa salmon and a more distant relationship to masu salmon for both. Additional analyses of the masu salmon species complex will benefit from using the new reference genome assembly.

The genetic architecture of resistance to flubendiamide insecticide in Helicoverpa armigera (Hübner)

Insecticide resistance is a major problem in food production, environmental sustainability, and human health. The cotton bollworm Helicoverpa armigera is a globally distributed crop pest affecting over 300 crop species. H. armigera has rapidly evolved insecticide resistance, making it one of the most damaging pests worldwide. Understanding the genetic basis of insecticide resistance provides insights to develop tools, such as molecular markers, that can be used to slow or prevent the evolution of resistance. We explore the genetic architecture of H. armigera resistance to a widely used insecticide, flubendiamide, using two complementary approaches: genome-wide association studies (GWAS) in wild-caught samples and quantitative trait locus (QTL) mapping in a controlled cross of susceptible and resistant laboratory strains. Both approaches identified one locus on chromosome 2, revealing two SNPs within 976 bp that can be used to monitor field resistance to flubendiamide. This was the only region identified using linkage mapping, though GWAS revealed additional sites associated with resistance. Other loci identified by GWAS in field populations contained known insecticide detoxification genes from the ATP-binding cassette family, ABCA1, ABCA3, ABCF2 and MDR1. Our findings revealed an oligogenic genetic architecture, contrasting previous reports of monogenic resistance associated with the ryanodine receptor. This work elucidates the genetic basis of rapidly evolving insecticide resistance and will contribute to developing effective insecticide resistance management strategies.

Pan-genome analyses of 11 Fraxinus species provide insights into salt adaptation in ash trees

Ash trees (Fraxinus) exhibit rich genetic diversity and wide adaptation to various ecological environments, and several species are highly salt tolerant. Dissecting the genomic basis of salt adaptation in Fraxinus is vital for its resistance breeding. Here, we present 11 high-quality chromosome-level genome assemblies for Fraxinus species, which reveal two unequal subgenome compositions and two recent whole-genome triplication events in their evolutionary history. A Fraxinus pan-genome was constructed on the basis of structural variations and revealed that presence-absence variations (PAVs) of transmembrane transport genes have likely contributed to salt adaptation in Fraxinus. Through whole-genome resequencing of an F1 population from an interspecies cross of F. velutina 'Lula 3' (salt tolerant) with F. pennsylvanica 'Lula 5' (salt sensitive), we mapped salt-tolerance PAV-based quantitative trait loci (QTLs) and pinpointed two PAV-QTLs and candidate genes associated with Fraxinus salt tolerance. Mechanistically, FvbHLH85 enhances salt tolerance by mediating reactive oxygen species and Na+/K+ homeostasis, whereas FvSWEET5 enhances salt tolerance by mediating osmotic homeostasis. Collectively, these findings provide valuable genomic resources for Fraxinus salt-resistance breeding and the research community.

Recombination landscape and karyotypic variations revealed by linkage mapping in the grapevine downy mildew pathogen Plasmopara viticola

Plasmopara viticola, the causal agent of grapevine downy mildew, is a biotrophic oomycete engaged in a tight coevolutionary relationship with its host. Rapid adaptation of the pathogen is favored by annual sexual reproduction that generates genotypic diversity. With the aim of studying the recombination landscape across the P. viticola genome, we generated 2 half-sibling F1 progenies (N = 189 and 162). Using targeted SNP sequencing, between 1,405 and 1,894 markers were included in parental linkage maps, and a consensus map was obtained by integrating 4,509 markers. The reference genome could be assembled into 17 pseudochromosomes, anchoring 88% of its physical length. We observed a strong collinearity between parental genomes and extensive synteny with the downy mildew Peronospora effusa. In the consensus map, the median recombination rate was 13.8 cM/Mb. The local recombination rate was highly variable along chromosomes, and recombination was suppressed in putative centromeric regions. Recombination rate was found negatively correlated with repeats' coverage and positively correlated with gene coverage. However, genes encoding secreted proteins and putative effectors were underrepresented in highly recombining regions. In both progenies, about 5% of the individuals presented karyotypic anomalies. Aneuploidies and triploidies almost exclusively originated from the male-transmitted chromosomes. Triploids resulted from fertilization by diploid gametes, but also from dispermy. Obligatory sexual reproduction each year may explain the lower level of karyotypic variation in P. viticola compared to other oomycetes. The linkage maps will be useful to guide future de novo chromosome-scale assemblies of P. viticola genomes and to perform forward genetics.

The recombination landscape of the barn owl, from families to populations

Homologous recombination is a meiotic process that generates diversity along the genome and interacts with all evolutionary forces. Despite its importance, studies of recombination landscapes are lacking due to methodological limitations and limited data. Frequently used approaches include linkage mapping based on familial data that provides sex-specific broad-scale estimates of realized recombination and inferences based on population linkage disequilibrium that reveal a more fine-scale resolution of the recombination landscape, albeit dependent on the effective population size and the selective forces acting on the population. In this study, we use a combination of these 2 methods to elucidate the recombination landscape for the Afro-European barn owl (Tyto alba). We find subtle differences in crossover placement between sexes that lead to differential effective shuffling of alleles. Linkage disequilibrium-based estimates of recombination are concordant with family-based estimates and identify large variation in recombination rates within and among linkage groups. Larger chromosomes show variation in recombination rates, while smaller chromosomes have a universally high rate that shapes the diversity landscape. We find that recombination rates are correlated with gene content, genetic diversity, and GC content. We find no conclusive differences in the recombination landscapes between populations. Overall, this comprehensive analysis enhances our understanding of recombination dynamics, genomic architecture, and sex-specific variation in the barn owl, contributing valuable insights to the broader field of avian genomics.

A chromosome-anchored reference assembly for the gray snapper, Lutjanus griseus

BACKGROUND

The gray snapper (Lutjanus griseus) is a marine reef fish commonly found in coastal and shelf waters of the tropical and subtropical western Atlantic Ocean. In this work, a draft reference genome was developed to support population genomic studies of gray snapper needed to assist with conservation and fisheries management efforts.

METHODS AND RESULTS

Hybrid assembly of PacBio and Illumina sequencing reads yielded a 1,003,098,032 bp reference across 2039 scaffolds with N50 and L50 values of 1,691,591 bp and 163 scaffolds, respectively. A linkage map was generated by genotyping parents and 286 offspring of a single pair cross using the double digest Restriction Associated DNA (ddRAD) protocol. The map featured 10,965 informative markers that were assigned to 24 linkage groups and used to scaffold the assembly. The anchored assembly spanned 962,844,722 bp (N50 = 41,865,368 bp, L50 = 11 scaffolds) in 24 pseudo chromosomes and yielded a BUSCO score of 95.2%. Annotation of the final assembly in Augustus revealed 39,070 candidate genes. Gene ontology annotation was obtained for 48.8% of the predicted genes. Analysis of shared syntenic regions revealed that each gray snapper chromosome matched a unique Japanese medaka, Oryzias latipes, counterpart and the reference showed a high degree of synteny with the closely related Lutjanus erythropterus assembly.

CONCLUSIONS

This resource will greatly enhance genomic studies of conservation and management of natural populations as well as efforts to develop breeding programs for this species and other lutjanids.

Repositioning of centromere-associated repeats during karyotype evolution in Oryzias fishes

The karyotype, which is the number and shape of chromosomes, is a fundamental characteristic of all eukaryotes. Karyotypic changes play an important role in many aspects of evolutionary processes, including speciation. In organisms with monocentric chromosomes, it was previously thought that chromosome number changes were mainly caused by centric fusions and fissions, whereas chromosome shape changes, that is, changes in arm numbers, were mainly due to pericentric inversions. However, recent genomic and cytogenetic studies have revealed examples of alternative cases, such as tandem fusions and centromere repositioning, found in the karyotypic changes within and between species. Here, we employed comparative genomic approaches to investigate whether centromere repositioning occurred during karyotype evolution in medaka fishes. In the medaka family (Adrianichthyidae), the three phylogenetic groups differed substantially in their karyotypes. The Oryzias latipes species group has larger numbers of chromosome arms than the other groups, with most chromosomes being metacentric. The O. javanicus species group has similar numbers of chromosomes to the O. latipes species group, but smaller arm numbers, with most chromosomes being acrocentric. The O. celebensis species group has fewer chromosomes than the other two groups and several large metacentric chromosomes that were likely formed by chromosomal fusions. By comparing the genome assemblies of O. latipes, O. javanicus, and O. celebensis, we found that repositioning of centromere-associated repeats might be more common than simple pericentric inversion. Our results demonstrated that centromere repositioning may play a more important role in karyotype evolution than previously appreciated.

Genes controlling hydrolysate toxin tolerance identified by QTL analysis of the natural Saccharomyces cerevisiae BCC39850

Lignocellulosic material can be converted to valorized products such as fuels. Pretreatment is an essential step in conversion, which is needed to increase the digestibility of the raw material for microbial fermentation. However, pretreatment generates by-products (hydrolysate toxins) that are detrimental to microbial growth. In this study, natural Saccharomyces strains isolated from habitats in Thailand were screened for their tolerance to synthetic hydrolysate toxins (synHTs). The Saccharomyces cerevisiae natural strain BCC39850 (toxin-tolerant) was crossed with the laboratory strain CEN.PK2-1C (toxin-sensitive), and quantitative trait locus (QTL) analysis was performed on the segregants using phenotypic scores of growth (OD600) and glucose consumption. VMS1, DET1, KCS1, MRH1, YOS9, SYO1, and YDR042C were identified from QTLs as candidate genes associated with the tolerance trait. CEN.PK2-1C knockouts of the VMS1, YOS9, KCS1, and MRH1 genes exhibited significantly greater hydrolysate toxin sensitivity to growth, whereas CEN.PK2-1C knock-ins with replacement of VMS1 and MRH1 genes from the BCC39850 alleles showed significant increased ethanol production titers compared with the CEN.PK2-1C parental strain in the presence of synHTs. The discovery of VMS1, YOS9, MRH1, and KCS1 genes associated with hydrolysate toxin tolerance in S. cerevisiae indicates the roles of the endoplasmic-reticulum-associated protein degradation pathway, plasma membrane protein association, and the phosphatidylinositol signaling system in this trait. KEY POINTS: • QTL analysis was conducted using a hydrolysate toxin-tolerant S. cerevisiae natural strain • Deletion of VMS1, YOS9, MRH1, and KCS1 genes associated with hydrolysate toxin-sensitivity • Replacement of VMS1 and MRH1 with natural strain alleles increased ethanol production titers in the presence of hydrolysate toxins.

Chromosome-scale genome assembly of acerola (Malpighia emarginata DC.)

Acerola (Malpighia emarginata DC.) is a tropical evergreen shrub that produces vitamin C-rich fruits. Increasing fruit nutrition is one of the main targets of acerola breeding programs. Genomic tools have been shown to accelerate plant breeding even in fruiting tree species, which generally have a long-life cycle; however, the availability of genomic resources in acerola, so far, has been limited. In this study, as a first step toward developing an efficient breeding technology for acerola, we established a chromosome-scale genome assembly of acerola using high-fidelity long-read sequencing and genetic mapping. The resultant assembly comprises 10 chromosome-scale sequences that span a physical distance of 1,032.5 Mb and contain 35,892 predicted genes. Phylogenetic analysis of genome-wide SNPs in 60 acerola breeding materials revealed 3 distinct genetic groups. Overall, the genomic resource of acerola developed in this study, including its genome and gene sequences, genetic map, and phylogenetic relationship among breeding materials, will not only be useful for acerola breeding but will also facilitate genomic and genetic studies on acerola and related species.

Too big to purge: persistence of deleterious Mutations in Island populations of the European Barn Owl (Tyto alba)

A key aspect of assessing the risk of extinction/extirpation for a particular wild species or population is the status of inbreeding, but the origin of inbreeding and the current mutational load are also two crucial factors to consider when determining survival probability of a population. In this study, we used samples from 502 barn owls from continental and island populations across Europe, with the aim of quantifying and comparing the level of inbreeding between populations with differing demographic histories. In addition to comparing inbreeding status, we determined whether inbreeding is due to non-random mating or high co-ancestry within the population. We show that islands have higher levels of inbreeding than continental populations, and that this is mainly due to small effective population sizes rather than recent consanguineous mating. We assess the probability that a region is autozygous along the genome and show that this probability decreased as the number of genes present in that region increased. Finally, we looked for evidence of reduced selection efficiency and purging in island populations. Among island populations, we found an increase in numbers of both neutral and deleterious minor alleles, possibly as a result of drift and decreased selection efficiency but we found no evidence of purging.

Chromosome-level genome assemblies and genetic maps reveal heterochiasmy and macrosynteny in endangered Atlantic Acropora

BACKGROUND

Over their evolutionary history, corals have adapted to sea level rise and increasing ocean temperatures, however, it is unclear how quickly they may respond to rapid change. Genome structure and genetic diversity contained within may highlight their adaptive potential.

RESULTS

We present chromosome-scale genome assemblies and linkage maps of the critically endangered Atlantic acroporids, Acropora palmata and A. cervicornis. Both assemblies and linkage maps were resolved into 14 chromosomes with their gene content and colinearity. Repeats and chromosome arrangements were largely preserved between the species. The family Acroporidae and the genus Acropora exhibited many phylogenetically significant gene family expansions. Macrosynteny decreased with phylogenetic distance. Nevertheless, scleractinians shared six of the 21 cnidarian ancestral linkage groups as well as numerous fission and fusion events compared to other distantly related cnidarians. Genetic linkage maps were constructed from one A. palmata family and 16 A. cervicornis families using a genotyping array. The consensus maps span 1,013.42 cM and 927.36 cM for A. palmata and A. cervicornis, respectively. Both species exhibited high genome-wide recombination rates (3.04 to 3.53 cM/Mb) and pronounced sex-based differences, known as heterochiasmy, with 2 to 2.5X higher recombination rates estimated in the female maps.

CONCLUSIONS

Together, the chromosome-scale assemblies and genetic maps we present here are the first detailed look at the genomic landscapes of the critically endangered Atlantic acroporids. These data sets revealed that adaptive capacity of Atlantic acroporids is not limited by their recombination rates. The sister species maintain macrosynteny with few genes with high sequence divergence that may act as reproductive barriers between them. In the Atlantic Acropora, hybridization between the two sister species yields an F1 hybrid with limited fertility despite the high levels of macrosynteny and gene colinearity of their genomes. Together, these resources now enable genome-wide association studies and discovery of quantitative trait loci, two tools that can aid in the conservation of these species.

Genetic linkage analysis of stable QTLs in Gossypium hirsutum RIL population revealed function of GhCesA4 in fiber development

INTRODUCTION

Upland cotton is an important allotetrapolyploid crop providing natural fibers for textile industry. Under the present high-level breeding and production conditions, further simultaneous improvement of fiber quality and yield is facing unprecedented challenges due to their complex negative correlations.

OBJECTIVES

The study was to adequately identify quantitative trait loci (QTLs) and dissect how they orchestrate the formation of fiber quality and yield.

METHODS

A high-density genetic map (HDGM) based on an intraspecific recombinant inbred line (RIL) population consisting of 231 individuals was used to identify QTLs and QTL clusters of fiber quality and yield traits. The weighted gene correlation network analysis (WGCNA) package in R software was utilized to identify WGCNA network and hub genes related to fiber development. Gene functions were verified via virus-induced gene silencing (VIGS) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 strategies.

RESULTS

An HDGM consisting of 8045 markers was constructed spanning 4943.01 cM of cotton genome. A total of 295 QTLs were identified based on multi-environmental phenotypes. Among 139 stable QTLs, including 35 newly identified ones, seventy five were of fiber quality and 64 yield traits. A total of 33 QTL clusters harboring 74 QTLs were identified. Eleven candidate hub genes were identified via WGCNA using genes in all stable QTLs and QTL clusters. The relative expression profiles of these hub genes revealed their correlations with fiber development. VIGS and CRISPR/Cas9 edition revealed that the hub gene cellulose synthase 4 (GhCesA4, GH_D07G2262) positively regulate fiber length and fiber strength formation and negatively lint percentage.

CONCLUSION

Multiple analyses demonstrate that the hub genes harbored in the QTLs orchestrate the fiber development. The hub gene GhCesA4 has opposite pleiotropic effects in regulating trait formation of fiber quality and yield. The results facilitate understanding the genetic basis of negative correlation between cotton fiber quality and yield.

QTL Mapping-Based Identification of Visceral White-Nodules Disease Resistance Genes in Larimichthys polyactis

Disease outbreaks in aquaculture have recently intensified. In particular, visceral white-nodules disease, caused by Pseudomonas plecoglossicida, has severely hindered the small yellow croaker (Larimichthys polyactis) aquaculture industry. However, research on this disease is limited. To address this gap, the present study employed a 100K SNP chip to genotype individuals from an F1 full-sib family, identify single nucleotide polymorphisms (SNPs), and construct a genetic linkage map for this species. A high-density genetic linkage map spanning a total length of 1395.72 cM with an average interval of 0.08 cM distributed across 24 linkage groups was obtained. Employing post-infection survival time as an indicator of disease resistance, 13 disease resistance-related quantitative trait loci (QTLs) were detected, and these regions included 169 genes. Functional enrichment analyses pinpointed 11 candidate disease resistance-related genes. RT-qPCR analysis revealed that the genes of chmp1a and arg1 are significantly differentially expressed in response to P. plecoglossicida infection in spleen and liver tissues, indicating their pivotal functions in disease resistance. In summary, in addition to successfully constructing a high-density genetic linkage map, this study reports the first QTL mapping for visceral white-nodules disease resistance. These results provide insight into the intricate molecular mechanisms underlying disease resistance in the small yellow croaker.

The Hypolimnas misippus Genome Supports a Common Origin of the W Chromosome in Lepidoptera

Moths and butterflies (Lepidoptera) have a heterogametic sex chromosome system with females carrying ZW chromosomes and males ZZ. The lack of W chromosomes in early-diverging lepidopteran lineages has led to the suggestion of an ancestral Z0 system in this clade and a B chromosome origin of the W. This contrasts with the canonical model of W chromosome evolution in which the W would have originated from the same homologous autosomal pair as the Z chromosome. Despite the distinct models proposed, the rapid evolution of the W chromosome has hindered the elucidation of its origin. Here, we present high-quality, chromosome-level genome assemblies of 2 Hypolimnas species (Hypolimnas misippus and Hypolimnas bolina) and use the H. misippus assembly to explore the evolution of W chromosomes in butterflies and moths. We show that in H. misippus, the W chromosome has higher similarity to the Z chromosome than any other chromosome, which could suggest a possible origin from the same homologous autosome pair as the Z chromosome. However, using genome assemblies of closely related species (ditrysian lineages) containing assembled W chromosomes, we present contrasting evidence suggesting that the W chromosome might have evolved from a B chromosome instead. Crucially, by using a synteny analysis to infer homology, we show that W chromosomes are likely to share a common evolutionary origin in Lepidoptera. This study highlights the difficulty of studying the evolution of W chromosomes and contributes to better understanding its evolutionary origins.

A high-density linkage map reveals broad- and fine-scale sex differences in recombination in the hihi (stitchbird; Notiomystis cincta)

Recombination, the process of DNA exchange between homologous chromosomes during meiosis, plays a major role in genomic diversity and evolutionary change. Variation in recombination rate is widespread despite recombination often being essential for progression of meiosis. One such variation is heterochiasmy, where recombination rates differ between sexes. Heterochiasmy has been observed across broad taxonomic groups, yet it remains an evolutionary enigma. We used Lep-MAP3, a pedigree-based software that is efficient in handling large datasets, to generate linkage maps for the hihi or stitchbird (Notiomystis cincta), utilising information from >36 K SNPs and 36 families. We constructed 29 linkage maps, including for the previously unscaffolded Z chromosome. The hihi is an endangered passerine endemic to Aotearoa New Zealand that is sexually dimorphic and exhibits high levels of sexual conflict, including sperm competition. Patterns in recombination in the hihi are consistent with those in other birds, including higher recombination rates in micro-chromosomes. Heterochiasmy in the hihi is male-biased, in line with predictions of the Haldane-Huxley rule, with the male linkage map being 15% longer. Micro-chromosomes exhibit heterochiasmy to a greater extent, contrary to that reported in other birds. At the intra-chromosomal level, heterochiasmy is higher nearer to chromosome ends and in gene-rich regions. Regions of extreme heterochiasmy are enriched for genes implicated in cell structure. This study adds an important contribution in assessing evolutionary theories of heterochiasmy and provides a framework for future studies investigating fine-scale heterochiasmy.

An ultra-dense linkage map identified quantitative trait loci corresponding to fruit quality- and size-related traits in red goji berry

Goji berries are a small-fruited shrub with industrial importance whose fruit considered beneficial in both fresh and dried forms. Current germplasms of goji berries include small fruits with a short shelf life, less sweet and bitter taste, and a lack of appropriate genetic information. This study aimed to employ whole genome resequencing to generate an ultra-dense bin linkage map and to elucidate the genetic basis of goji fruit quality and size using quantitative trait loci (QTL) mapping analysis in a cross-pollinated hybrid population. To achieve this goal, human sensory tests were carried out to determine the bitter taste (BT) and sweet taste (ST), and to quantify the soluble solid content (SSC), fruit firmness (FF), and fruit size-related traits of fresh goji fruits over three or four years. The results revealed that the goji bin linkage map based on resequencing spanned a total length of 966.42 cM and an average bin interval of 0.03 cM. Subsequent variant calling and ordering resulted in 3,058 bins containing 35,331 polymorphic markers across 12 chromosomes. A total of 99 QTLs, with individual loci in different environments explaining a phenotypic variance of 1.21-16.95% were identified for the studied traits. Ten major effects, including colocalized QTLs corresponding to different traits, were identified on chromosomes 1, 3, 5, 6, 7, and 8, with a maximum Logarithm of Odds (LOD) of 29.25 and 16.95% of explained phenotypic variance (PVE). In addition, four stable loci, one for FF, one for fruit weight (FW), and two for fruit shape index (FSI), were mainly mapped on chromosomes 5, 6, and 7, elucidating 2.10-16.95% PVE. These findings offer valuable insights into the genetic architecture of goji fruit traits along with identified specific loci and markers to further improve and develop sweeter, less bitter and larger fruited goji berry cultivars with extended shelf life.

The Fate of a Polygenic Phenotype Within the Genomic Landscapes of Introgression in the European Seabass Hybrid Zone

Unraveling the evolutionary mechanisms and consequences of hybridization is a major concern in biology. Many studies have documented the interplay between recombination and selection in modulating the genomic landscape of introgression, but few have considered how associations with phenotype may affect this landscape. Here, we use the European seabass (Dicentrarchus labrax), a key species in marine aquaculture that undergoes natural hybridization, to determine how selection on phenotype modulates the introgression landscape between Atlantic and Mediterranean lineages. We use a high-density single nucleotide polymorphism array to assess individual local ancestry along the genome and improve the mapping of muscle fat content, a polygenic trait that is divergent between lineages. Taking into account variation in recombination rates, we reveal a purging of Atlantic ancestry in the admixed Mediterranean populations. While Atlantic individuals had higher muscle fat content, we observed that genomic regions associated with this trait in Mediterranean populations displayed reduced introgression of Atlantic ancestry. These results emphasize how selection against maladapted alleles shapes the genomic landscape of introgression.

The Genetic Architecture of Recombination Rates is Polygenic and Differs Between the Sexes in Wild House Sparrows (Passer domesticus)

Meiotic recombination through chromosomal crossing-over is a fundamental feature of sex and an important driver of genomic diversity. It ensures proper disjunction, allows increased selection responses, and prevents mutation accumulation; however, it is also mutagenic and can break up favorable haplotypes. This cost-benefit dynamic is likely to vary depending on mechanistic and evolutionary contexts, and indeed, recombination rates show huge variation in nature. Identifying the genetic architecture of this variation is key to understanding its causes and consequences. Here, we investigate individual recombination rate variation in wild house sparrows (Passer domesticus). We integrate genomic and pedigree data to identify autosomal crossover counts (ACCs) and intrachromosomal allelic shuffling (r¯intra) in 13,056 gametes transmitted from 2,653 individuals to their offspring. Females had 1.37 times higher ACC, and 1.55 times higher r¯intra than males. ACC and r¯intra were heritable in females and males (ACC h2 = 0.23 and 0.11; r¯intra  h2 = 0.12 and 0.14), but cross-sex additive genetic correlations were low (rA = 0.29 and 0.32 for ACC and r¯intra). Conditional bivariate analyses showed that all measures remained heritable after accounting for genetic values in the opposite sex, indicating that sex-specific ACC and r¯intra can evolve somewhat independently. Genome-wide models showed that ACC and r¯intra are polygenic and driven by many small-effect loci, many of which are likely to act in trans as global recombination modifiers. Our findings show that recombination rates of females and males can have different evolutionary potential in wild birds, providing a compelling mechanism for the evolution of sexual dimorphism in recombination.

Sex chromosome turnover in hybridizing stickleback lineages

Recent discoveries of sex chromosome diversity across the tree of life have challenged the canonical model of conserved sex chromosome evolution and evoked new theories on labile sex chromosomes that maintain less differentiation and undergo frequent turnover. However, theories of labile sex chromosome evolution lack direct empirical support due to the paucity of case studies demonstrating ongoing sex chromosome turnover in nature. Two divergent lineages (viz. WL & EL) of nine-spined sticklebacks (Pungitius pungitius) with different sex chromosomes (linkage group [LG] 12 in the EL, unknown in the WL) hybridize in a natural secondary contact zone in the Baltic Sea, providing an opportunity to study ongoing turnover between coexisting sex chromosomes. In this study, we first identify an 80 kbp genomic region on LG3 as the sex-determining region (SDR) using whole-genome resequencing data of family crosses of a WL population. We then verify this region as the SDR in most other WL populations and demonstrate a potentially ongoing sex chromosome turnover in admixed marine populations where the evolutionarily younger and homomorphic LG3 sex chromosome replaces the older and heteromorphic LG12 sex chromosome. The results provide a rare glimpse of sex chromosome turnover in the wild and indicate the possible existence of additional yet undiscovered sex chromosome diversity in Pungitius sticklebacks.

Heterogeneous genomic architecture of skeletal armour traits in sticklebacks

Whether populations adapt to similar selection pressures using the same underlying genetic variants depends on population history and the distribution of standing genetic variation at the metapopulation level. Studies of sticklebacks provide a case in point: when colonizing and adapting to freshwater habitats, three-spined sticklebacks (Gasterosteus aculeatus) with high gene flow tend to fix the same adaptive alleles in the same major loci, whereas nine-spined sticklebacks (Pungitius pungitius) with limited gene flow tend to utilize a more heterogeneous set of loci. In accordance with this, we report results of quantitative trait locus (QTL) analyses using a backcross design showing that lateral plate number variation in the western European nine-spined sticklebacks mapped to 3 moderate-effect QTL, contrary to the major-effect QTL in three-spined sticklebacks and different from the 4 QTL previously identified in the eastern European nine-spined sticklebacks. Furthermore, several QTL were identified associated with variation in lateral plate size, and 3 moderate-effect QTL with body size. Together, these findings indicate more heterogenous and polygenic genetic underpinnings of skeletal armour variation in nine-spined than three-spined sticklebacks, indicating limited genetic parallelism underlying armour trait evolution in the family Gasterostidae.

Chromosome-level genome assemblies for two quinoa inbred lines from northern and southern highlands of Altiplano where quinoa originated

Quinoa is emerging as a key seed crop for global food security due to its ability to grow in marginal environments and its excellent nutritional properties. Because quinoa is partially allogamous, we have developed quinoa inbred lines necessary for molecular genetic analysis. Our comprehensive genomic analysis showed that the quinoa inbred lines fall into three genetic subpopulations: northern highland, southern highland, and lowland. Lowland and highland quinoa are the same species, but have very different genotypes and phenotypes. Lowland quinoa has relatively small grains and a darker grain color, and is widely tested and grown around the world. In contrast, the white, large-grained highland quinoa is grown in the Andean highlands, including the region where quinoa originated, and is exported worldwide as high-quality quinoa. Recently, we have shown that viral vectors can be used to regulate endogenous genes in quinoa, paving the way for functional genomics to reveal the diversity of quinoa. However, although a high-quality assembly has recently been reported for a lowland quinoa line, genomic resources of the quality required for functional genomics are not available for highland quinoa lines. Here we present high-quality chromosome-level genome assemblies for two highland inbred quinoa lines, J075 representing the northern highland line and J100 representing the southern highland line, using PacBio HiFi sequencing and dpMIG-seq. In addition, we demonstrate the importance of verifying and correcting reference-based scaffold assembly with other approaches such as linkage maps. The assembled genome sizes of J075 and J100 are 1.29 and 1.32 Gb, with contigs N50 of 66.3 and 12.6 Mb, and scaffold N50 of 71.2 and 70.6 Mb, respectively, comprising 18 pseudochromosomes. The repetitive sequences of J075 and J100 represent 72.6% and 71.5% of the genome, the majority of which are long terminal repeats, representing 44.0% and 42.7% of the genome, respectively. The de novo assembled genomes of J075 and J100 were predicted to contain 65,303 and 64,945 protein-coding genes, respectively. The high quality genomes of these highland quinoa lines will facilitate quinoa functional genomics research on quinoa and contribute to the identification of key genes involved in environmental adaptation and quinoa domestication.

Genetic analysis of agronomic and physiological traits associated with latex yield revealed complex genetic bases in Hevea brasiliensis

Hevea brasiliensis, a natural rubber producing species, is widely cultivated due to its high rubber yield potential. Natural rubber is synthesised in the rubber particles of laticifers. Latex diagnosis (LD) was established to characterise the physiological state of the laticiferous system by measuring its physiological parameters, i.e., sucrose, inorganic phosphorous (Pi), thiols and total solid content (TSC). Rubber clones are often classified in three groups i.e., quick starters, medium starters and slow starters. To better understand the genetic bases of latex yield, a biparental population was generated from a cross between the quick-starter clone PB 260 and the medium-starter clone SP 217. LD was performed during the peak latex production season and used to calculate sucrose loading. The agronomic and physiological parameters associated with latex yield led to the classification of genotypes according to the rubber clonal typology and to the identification of quantitative trait loci (QTL) using a high-density map. Inorganic phosphorous content was positively associated with yield during the first year of production thus enabling identification of quick-starter clones. In addition, the LD-based clonal typology led to determine the long-term yield potential and the use of appropriate ethephon stimulation. QTL analysis successfully identified several QTLs related to yield, sucrose, Pi and TSC. One QTL related to sucrose loading was identified in the same position as the QTL for sucrose on linkage group 1. To our knowledge, this is the first study to report QTL analysis for this trait. The use of a high-density map enables the identification of genes underlying QTLs. Several putative genes underlying QTLs related to yield, sucrose and TSC were identified.

Fine Mapping of QTLs for Alkaline Tolerance in Crucian Carp (Carassius auratus) Using Genome-Wide SNP Markers

Crucian carp (Carassius auratus) is widely distributed in the world and has become an economically freshwater fish. The population in Lake Dali Nur can tolerate the extreme alkaline environment with alkalinity over 50 mmol/L (pH 9.6), thus providing a special model for exploring alkali-tolerant molecular markers in an extremely alkaline environment. In this study, we constructed a high-density and high-resolution linkage map with 16,224 SNP markers based on genotyping-by-sequencing (GBS) consisting of 152 progenies and conducted QTL studies for alkali-tolerant traits. The total length of the linkage map was 3918.893 cM, with an average distance of 0.241 cM. Two QTLs for the ammonia-N-tolerant trait were detected on LG27 and LG45. A QTL for the urea-N-tolerant trait was detected on LG27. Interestingly, mapping the two QTLs on LG27 revealed that the mapped genes were both located in the intron of CDC42. GO functional annotation and KEGG enrichment analysis results indicated that the biological functions might be involved in the cell cycle, cellular senescence, MAPK, and Ras signaling pathways. These findings suggest that CDC42 may play an important role in the process of dealing with extremely alkaline environments.

A high-density linkage map and sex-determination loci in Pacific white shrimp (Litopenaeus vannamei)

BACKGROUND

Expansion of genomic resources for the Pacific white shrimp (Litopenaeus vannamei), such as the construction of dense genetic linkage maps, is crucial for the application of genomic tools in order to improve economically relevant traits. Sexual dimorphism exists in Pacific white shrimp, and the mapping of the sex-determination region in this species may help in future reproductive applications. We have constructed male, female, and sex-averaged high-density genetic maps using a 50 K single-nucleotide polymorphism (SNP) array, followed by a genome-wide association study (GWAS) to identify genomic regions associated with sex in white shrimp.

RESULTS

The genetic map yielded 15,256 SNPs assigned to 44 linkage groups (LG). The lengths of the male, female, and sex-averaged maps were 5,741.36, 5,461.20 and 5,525.26 cM, respectively. LG18 was found to be the largest for both sexes, whereas LG44 was the shortest for males and LG31 for females. A sex-determining region was found in LG31 with 21 statistically significant SNPs. The most important SNP was previously identified as a sex-linked marker and was able to identify 99% of the males and 88% of the females. Although other significant markers had a lower ability to determine sex, putative genes were intercepted or close to them. The oplophorus-luciferin 2-monooxygenase, serine/arginine repetitive matrix protein and spermine oxidase genes were identified as candidates with possible participation in important processes of sexual differentiation in shrimp.

CONCLUSIONS

Our results provide novel genomic resources for shrimp, including a high-density linkage map and new insights into the sex-determining region in L. vannamei, which may be usefulfor future genetics and reproduction applications.

Degenerate oligonucleotide primer MIG-seq: an effective PCR-based method for high-throughput genotyping

Next-generation sequencing (NGS) library construction often involves using restriction enzymes to decrease genome complexity, enabling versatile polymorphism detection in plants. However, plant leaves frequently contain impurities, such as polyphenols, necessitating DNA purification before enzymatic reactions. To overcome this problem, we developed a PCR-based method for expeditious NGS library preparation, offering flexibility in number of detected polymorphisms. By substituting a segment of the simple sequence repeat sequence in the MIG-seq primer set (MIG-seq being a PCR method enabling library construction with low-quality DNA) with degenerate oligonucleotides, we introduced variability in detectable polymorphisms across various crops. This innovation, named degenerate oligonucleotide primer MIG-seq (dpMIG-seq), enabled a streamlined protocol for constructing dpMIG-seq libraries from unpurified DNA, which was implemented stably in several crop species, including fruit trees. Furthermore, dpMIG-seq facilitated efficient lineage selection in wheat and enabled linkage map construction and quantitative trait loci analysis in tomato, rice, and soybean without necessitating DNA concentration adjustments. These findings underscore the potential of the dpMIG-seq protocol for advancing genetic analyses across diverse plant species.

Complete genome assembly provides a high-quality skeleton for pan-NLRome construction in melon

Melon (Cucumis melo L.), being under intensive domestication and selective breeding, displays an abundant phenotypic diversity. Wild germplasm with tolerance to stress represents an untapped genetic resource for discovery of disease-resistance genes. To comprehensively characterize resistance genes in melon, we generate a telomere-to-telomere (T2T) and gap-free genome of wild melon accession PI511890 (C. melo var. chito) with a total length of 375.0 Mb and a contig N50 of 31.24 Mb. The complete genome allows us to dissect genome architecture and identify resistance gene analogs. We construct a pan-NLRome using seven melon genomes, which include 208 variable and 18 core nucleotide-binding leucine-rich repeat receptors (NLRs). Multiple disease-related transcriptome analyses indicate that most up-regulated NLRs induced by pathogens are shell or cloud NLRs. The T2T gap-free assembly and the pan-NLRome not only serve as essential resources for genomic studies and molecular breeding of melon but also provide insights into the genome architecture and NLR diversity.

High-density genetic map and quantitative trait loci map of skin color in hawthorn (Crataegus pinnatifida bge. Var. major N.E.Br.)

Fruit skin color is an important trait of the hawthorn tree, which has an important influence on fruit quality. Crataegus pinnatifida Bge. var. Major N.E.Br. Is one of the most widely cultivated varieties in China and has a long history of medicinal use. In recent years, it has attracted the attention of the world due to its nutritional and medicinal values. Skin color is the focus of breeders and food processors. At present, skin color-related genes have still not been mapped. In this study, "Shandong Da Mianqiu" (♀, red skin color), "Da Huang Mianzha" (♂, yellow skin color) and 131 F1 hybrids were used to construct genetic map of hawthorn by RAD-seq, and QTL mapping was performed by combining these features with the hue angle and the observed color. In this study, 13,260 SNP was assigned to 17 linkage groups, with an integrated map covering 2,297.75 cM was constructed. A total of 5 QTLs related to hawthorn skin color were detected on LG1, LG3 and LG15. Whether hue angle or pericarp color acts as phenotype for QTL mapping, the candidate genes include bHLH086, WD repeat regions and Myb-like. bHLH, WD and Myb play an important role in the color regulation of Hawthorn skin color. These results lay a solid foundation for QTL mapping and molecular marker-assisted breeding of hawthorn.

Mapping quantitative trait loci associated with self-(in)compatibility in goji berries (Lycium barbarum)

BACKGROUND

Goji (Lycium barbarum L.) is a perennial deciduous shrub widely distributed in arid and semiarid regions of Northwest China. It is highly valued for its medicinal and functional properties. Most goji varieties are naturally self-incompatible, posing challenges in breeding and cultivation. Self-incompatibility is a complex genetic trait, with ongoing debates regarding the number of self-incompatible loci. To date, no genetic mappings has been conducted for S loci or other loci related to self-incompatibility in goji.

RESULTS

We used genome resequencing to create a high-resolution map for detecting de novo single-nucleotide polymorphisms (SNP) in goji. We focused on 229 F1 individuals from self-compatible '13-19' and self-incompatible 'new 9' varieties. Subsequently, we conducted a quantitative trait locus (QTL) analysis on traits associated with self-compatibility in goji berries. The genetic map consisted of 249,327 SNPs distributed across 12 linkage groups (LGs), spanning a total distance of 1243.74 cM, with an average interval of 0.002 cM. Phenotypic data related to self-incompatibility, such as average fruit weight, fruit rate, compatibility index, and comparable compatibility index after self-pollination and geitonogamy, were collected for the years 2021-2022, as well as for an extra year representing the mean data from 2021 to 2022 (2021/22). A total of 43 significant QTL, corresponding to multiple traits were identified, accounting for more than 11% of the observed phenotypic variation. Notably, a specific QTL on chromosome 2 consistently appeared across different years, irrespective of the relationship between self-pollination and geitonogamy. Within the localization interval, 1180 genes were annotated, including Lba02g01102 (annotated as an S-RNase gene), which showed pistil-specific expression. Cloning of S-RNase genes revealed that the parents had two different S-RNase alleles, namely S1S11 and S2S8. S-genotype identification of the F1 population indicated segregation of the four S-alleles from the parents in the offspring, with the type of S-RNase gene significantly associated with self-compatibility.

CONCLUSIONS

In summary, our study provides valuable insights into the genetic mechanism underlying self-compatibility in goji berries. This highlights the importance of further positional cloning investigations and emphasizes the importance of integration of marker-assisted selection in goji breeding programs.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

Hybrid speciation driven by multilocus introgression of ecological traits

Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.

The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

A self-compatible pear mutant derived from γ-irradiated pollen carries an 11-Mb duplication in chromosome 17

Self-compatibility is a highly desirable trait for pear breeding programs. Our breeding program previously developed a novel self-compatible pollen-part Japanese pear mutant (Pyrus pyrifolia Nakai), '415-1', by using γ-irradiated pollen. '415-1' carries the S-genotype S4dS5S5, with "d" indicating a duplication of S 5 responsible for breakdown of self-incompatibility. Until now, the size and inheritance of the duplicated segment was undetermined, and a reliable detection method was lacking. Here, we examined genome duplications and their inheritance in 140 F1 seedlings resulting from a cross between '515-20' (S1S3) and '415-1'. Amplicon sequencing of S-RNase and SFBB18 clearly detected S-haplotype duplications in the seedlings. Intriguingly, 30 partially triploid seedlings including genotypes S1S4dS5, S3S4dS5, S1S5dS5, S3S5dS5, and S3S4dS4 were detected among the 140 seedlings. Depth-of-coverage analysis using ddRAD-seq showed that the duplications in those individuals were limited to chromosome 17. Further analysis through resequencing confirmed an 11-Mb chromosome duplication spanning the middle to the end of chromosome 17. The duplicated segment remained consistent in size across generations. The presence of an S3S4dS4 seedling provided evidence for recombination between the duplicated S5 segment and the original S4haplotype, suggesting that the duplicated segment can pair with other parts of chromosome 17. This research provides valuable insights for improving pear breeding programs using partially triploid individuals.

Turnovers of Sex-Determining Mutation in the Golden Pompano and Related Species Provide Insights into Microevolution of Undifferentiated Sex Chromosome

The suppression of recombination is considered a hallmark of sex chromosome evolution. However, previous research has identified undifferentiated sex chromosomes and sex determination by single SNP in the greater amberjack (Seriola dumerili). We observed the same phenomena in the golden pompano (Trachinotus ovatus) of the same family Carangidae and discovered a different sex-determining SNP within the same gene Hsd17b1. We propose an evolutionary model elucidating the turnover of sex-determining mutations by highlighting the contrasting dynamics between purifying selection, responsible for maintaining W-linked Hsd17b1, and neutral evolution, which drives Z-linked Hsd17b1. Additionally, sporadic loss-of-function mutations in W-linked Hsd17b1 contribute to the conversion of W chromosomes into Z chromosomes. This model was directly supported by simulations, closely related species, and indirectly by zebrafish mutants. These findings shed new light on the early stages of sex chromosome evolution.

QTL analysis of femaleness in monoecious spinach and fine mapping of a major QTL using an updated version of chromosome-scale pseudomolecules

Although spinach is predominantly dioecious, monoecious plants with varying proportions of female and male flowers are also present. Recently, monoecious inbred lines with highly female and male conditions have been preferentially used as parents for F1-hybrids, rather than dioecious lines. Accordingly, identifying the loci for monoecism is an important issue for spinach breeding. We here used long-read sequencing and Hi-C technology to construct SOL_r2.0_pseudomolecule, a set of six pseudomolecules of spinach chromosomes (total length: 879.2 Mb; BUSCO complete 97.0%) that are longer and more genetically complete than our previous version of pseudomolecules (688.0 Mb; 81.5%). Three QTLs, qFem2.1, qFem3.1, and qFem6.1, responsible for monoecism were mapped to SOL_r2.0_pseudomolecule. qFem3.1 had the highest LOD score and corresponded to the M locus, which was previously identified as a determinant of monoecious expression, by genetic analysis of progeny from female and monoecious plants. The other QTLs were shown to modulate the ratio of female to male flowers in monoecious plants harboring a dominant allele of the M gene. Our findings will enable breeders to efficiently produce highly female- and male-monoecious parental lines for F1-hybrids by pyramiding the three QTLs. Through fine-mapping, we narrowed the candidate region for the M locus to a 19.5 kb interval containing three protein-coding genes and one long non-coding RNA gene. Among them, only RADIALIS-like-2a showed a higher expression in the reproductive organs, suggesting that it might play a role in reproductive organogenesis. However, there is no evidence that it is involved in the regulation of stamen and pistil initiation, which are directly related to the floral sex differentiation system in spinach. Given that auxin is involved in reproductive organ formation in many plant species, genes related to auxin transport/response, in addition to floral organ formation, were identified as candidates for regulators of floral sex-differentiation from qFem2.1 and qFem6.1.

The chromosome-scale genome and the genetic resistance machinery against insect herbivores of the Mexican toloache, Datura stramonium

Plant resistance refers to the heritable ability of plants to reduce damage caused by natural enemies, such as herbivores and pathogens, either through constitutive or induced traits like chemical compounds or trichomes. However, the genetic architecture-the number and genome location of genes that affect plant defense and the magnitude of their effects-of plant resistance to arthropod herbivores in natural populations remains poorly understood. In this study, we aimed to unveil the genetic architecture of plant resistance to insect herbivores in the annual herb Datura stramonium (Solanaceae) through quantitative trait loci mapping. We achieved this by assembling the species' genome and constructing a linkage map using an F2 progeny transplanted into natural habitats. Furthermore, we conducted differential gene expression analysis between undamaged and damaged plants caused by the primary folivore, Lema daturaphila larvae. Our genome assembly resulted in 6,109 scaffolds distributed across 12 haploid chromosomes. A single quantitative trait loci region on chromosome 3 was associated with plant resistance, spanning 0 to 5.17 cM. The explained variance by the quantitative trait loci was 8.44%. Our findings imply that the resistance mechanisms of D. stramonium are shaped by the complex interplay of multiple genes with minor effects. Protein-protein interaction networks involving genes within the quantitative trait loci region and overexpressed genes uncovered the key role of receptor-like cytoplasmic kinases in signaling and regulating tropane alkaloids and terpenoids, which serve as powerful chemical defenses against D. stramonium herbivores. The data generated in our study constitute important resources for delving into the evolution and ecology of secondary compounds mediating plant-insect interactions.

Gossypium purpurascens genome provides insight into the origin and domestication of upland cotton

INTRODUCTION

Allotetraploid upland cotton (Gossypium hirsutum L.) is native to the Mesoamerican and Caribbean regions, had been improved in the southern United States by the mid-eighteenth century, was then dispersed worldwide. However, a Hainan Island Native Cotton (HIC) has long been grown extensively on Hainan Island, China.

OBJECTIVES

Explore HIC's evolutionary relationship and genomic diversity with other tetraploid cottons, its origin and whether it was used for YAZHOUBU (Yazhou cloth, World Intangible Cultural Heritage) weaving, and the role of structural variations (SVs) in upland cotton domestication.

METHODS

We assembled a high-quality genome of one HIC plant. We performed phylogenetic analysis, divergence time estimation, principal component analysis and population differentiation estimation using cotton assemblies and/or resequencing data. SVs were detected by whole-genome comparison. A F2 population was used for linkage analysis and to study effects of SVs. Buoyancy and salt water tolerance tests for seeds were conducted.

RESULTS

We found that the HIC belongs to G. purpurascens. G. purpurascens is best classified as a primitive race of G. hirsutum. The potential for long range transoceanic dispersal of G. purpurascens seeds was proved. A set of SVs, selective sweep regions between G. hirsutum races and cultivars, and quantitative trait loci (QTLs) of eleven agronomic traits were obtained. SVs, especially large-scale SVs, were found to have important effects on cotton domestication and improvement. Of them, eight large-scale inversions strongly associated with yield and fiber quality have probably undergone artificial selection in domestication.

CONCLUSION

G. purpurascens including HIC is a primitive race of G. hirsutum, probably disperse to Hainan from Central America by floating on ocean currents, may have been partly domesticated, planted and was likely used for YAZHOUBU weaving in Hainan much earlier than the Pre-Columbian period. SV plays an important role in cotton domestication and improvement.

Postzygotic barriers persist despite ongoing introgression in hybridizing Mimulus species

The evolution of postzygotic isolation is thought to be a key step in maintaining species boundaries upon secondary contact, yet the dynamics and persistence of hybrid incompatibilities in naturally hybridizing species are not well understood. Here, we explore these issues using genetic mapping in three independent populations of recombinant inbred lines between naturally hybridizing monkeyflowers, Mimulus guttatus and Mimulus nasutus, from the sympatric Catherine Creek population. We discover that the three M. guttatus founders differ dramatically in admixture history, with nearly a quarter of one founder's genome introgressed from M. nasutus. Comparative genetic mapping in the three RIL populations reveals three new putative inversions, each one segregating among the M. guttatus founders, two due to admixture. We find strong, genome-wide transmission ratio distortion in all RILs, but patterns are highly variable among the three populations. At least some of this distortion appears to be explained by epistatic selection favouring parental genotypes, but tests of inter-chromosomal linkage disequilibrium also reveal multiple candidate Dobzhansky-Muller incompatibilities. We also map several genetic loci for hybrid pollen viability, including two interacting pairs that coincide with peaks of distortion. Remarkably, even with this limited sample of three M. guttatus lines, we discover abundant segregating variation for hybrid incompatibilities with M. nasutus, suggesting this population harbours diverse contributors to postzygotic isolation. Moreover, even with substantial admixture, hybrid incompatibilities between Mimulus species persist, suggesting postzygotic isolation might be a potent force in maintaining species barriers in this system.

Genomic resources for the Yellowfin tuna Thunnus albacares

BACKGROUND

The Yellowfin tuna (Thunnus albacares) is a large tuna exploited by major fisheries in tropical and subtropical waters of all oceans except the Mediterranean Sea. Genomic studies of population structure, adaptive variation or of the genetic basis of phenotypic traits are needed to inform fisheries management but are currently limited by the lack of a reference genome for this species. Here we report a draft genome assembly and a linkage map for use in genomic studies of T. albacares.

METHODS AND RESULTS

Illumina and PacBio SMRT sequencing were used in combination to generate a hybrid assembly that comprises 743,073,847 base pairs contained in 2,661 scaffolds. The assembly has a N50 of 351,587 and complete and partial BUSCO scores of 86.47% and 3.63%, respectively. Double-digest restriction associated DNA (ddRAD) was used to genotype the 2 parents and 164 of their F1 offspring resulting from a controlled breeding cross, retaining 19,469 biallelic single nucleotide polymorphism (SNP) loci. The SNP loci were used to construct a linkage map that features 24 linkage groups that represent the 24 chromosomes of yellowfin tuna. The male and female maps span 1,243.8 cM and 1,222.9 cM, respectively. The map was used to anchor the assembly in 24 super-scaffolds that contain 79% of the yellowfin tuna genome. Gene prediction identified 46,992 putative genes 20,203 of which could be annotated via gene ontology.

CONCLUSIONS

The draft reference will be valuable to interpret studies of genome wide variation in T. albacares and other Scombroid species.