Human promoter directionality is determined by transcriptional initiation and the opposing activities of INTS11 and CDK9.

RNA polymerase II (RNAPII) transcription initiates bidirectionally at many human protein-coding genes. Sense transcription usually dominates and leads to messenger RNA production, whereas antisense transcription rapidly terminates. The basis for this directionality is not fully understood. Here, we show that sense transcriptional initiation is more efficient than in the antisense direction, which establishes initial promoter directionality. After transcription begins, the opposing functions of the endonucleolytic subunit of Integrator, INTS11, and cyclin-dependent kinase 9 (CDK9) maintain directionality. Specifically, INTS11 terminates antisense transcription, whereas sense transcription is protected from INTS11-dependent attenuation by CDK9 activity. Strikingly, INTS11 attenuates transcription in both directions upon CDK9 inhibition, and the engineered recruitment of CDK9 desensitises transcription to INTS11. Therefore, the preferential initiation of sense transcription and the opposing activities of CDK9 and INTS11 explain mammalian promoter directionality.

Comparative transcriptomics reveal a novel tardigrade-specific DNA-binding protein induced in response to ionizing radiation.

Tardigrades are microscopic animals renowned for their ability to withstand extreme conditions, including high doses of ionizing radiation (IR). To better understand their radio-resistance, we first characterized induction and repair of DNA double- and single-strand breaks after exposure to IR in the model species Hypsibius exemplaris. Importantly, we found that the rate of single-strand breaks induced was roughly equivalent to that in human cells, suggesting that DNA repair plays a predominant role in tardigrades' radio-resistance. To identify novel tardigrade-specific genes involved, we next conducted a comparative transcriptomics analysis across three different species. In all three species, many DNA repair genes were among the most strongly overexpressed genes alongside a novel tardigrade-specific gene, which we named Tardigrade DNA damage Response 1 (TDR1). We found that TDR1 protein interacts with DNA and forms aggregates at high concentration suggesting it may condensate DNA and preserve chromosome organization until DNA repair is accomplished. Remarkably, when expressed in human cells, TDR1 improved resistance to Bleomycin, a radiomimetic drug. Based on these findings, we propose that TDR1 is a novel tardigrade-specific gene conferring resistance to IR. Our study sheds light on mechanisms of DNA repair helping cope with high levels of DNA damage inflicted by IR.

Surviving extreme radiation.

Tiny animals known as tardigrades use a combination of DNA repair machinery and a novel protein to mend their genome after intense ionizing radiation.

Comparative Cytogenetics of the Malagasy Ground Geckos of the Paroedura bastardi and Paroedura picta Species Groups.

We present a comparative chromosome study of several taxa of the Malagasy ground geckos of the Paroedura bastardi and P. picta species groups. We employed a preliminary molecular analysis using a trait of the mitochondrial 16S rRNA gene (of about 570 bp) to assess the taxonomic status of the samples studied and a cytogenetic analysis with standard karyotyping (5% Giemsa solution), silver staining (Ag-NOR staining) and sequential C-banding (C-banding + Giemsa and + fluorochromes). Our results show that all the taxa studied of the P. bastardi group (P. ibityensis, P. rennerae and P. cf. guibeae) have a similar karyotype composed of 2n = 34 chromosomes, with two metacentric pairs (1 and 3) and all other pairs being acrocentric. Chromosome diversification in the P. bastardi group was mainly linked to the diversification of heteromorphic sex chromosome systems (ZZ/ZW) in P. ibityensis and P. rennerae, while no heteromorphic sex chromosome pair was found in P. cf. guibeae. The two taxa investigated of the P. picta species group (here named P. picta and P. cf. picta based on molecular data) showed the same chromosome number of 2n = 36, mostly acrocentric elements, but differed in the number of metacentric elements, probably as a result of an inversion at chromosome pair 2. We highlight that the genus Paroedura is characterized by the independent diversification of heterogametic sex chromosomes in different evolutionary lineages and, similarly to other phylogenetically related gecko genera, by a progressive formation of a biarmed element by means of tandem fusions and inversions of distinct pairs.

Testosterone deficiency promotes arterial stiffening independent of sex chromosome complement.

BACKGROUND: Sex hormones and sex chromosomes play a vital role in cardiovascular disease. Testosterone plays a crucial role in men's health. Lower testosterone level is associated with cardiovascular and cardiometabolic diseases, including inflammation, atherosclerosis, and type 2 diabetes. Testosterone replacement is beneficial or neutral to men's cardiovascular health. Testosterone deficiency is associated with cardiovascular events. Testosterone supplementation to hypogonadal men improves libido, increases muscle strength, and enhances mood. We hypothesized that sex chromosomes (XX and XY) interaction with testosterone plays a role in arterial stiffening. METHODS: We used four core genotype male mice to understand the inherent contribution of sex hormones and sex chromosome complement in arterial stiffening. Age-matched mice were either gonadal intact or castrated at eight weeks plus an additional eight weeks to clear endogenous sex hormones. This was followed by assessing blood pressure, pulse wave velocity, echocardiography, and ex vivo passive vascular mechanics. RESULTS: Arterial stiffening but not blood pressure was more significant in castrated than testes-intact mice independent of sex chromosome complement. Castrated mice showed a leftward shift in stress-strain curves and carotid wall thinning. Sex chromosome complement (XX) in the absence of testosterone increased collagen deposition in the aorta and Kdm6a gene expression. CONCLUSION: Testosterone deprivation increases arterial stiffening and vascular wall remodeling. Castration increases Col1α1 in male mice with XX sex chromosome complement. Our study shows decreased aortic contractile genes in castrated mice with XX than XY sex chromosomes.

Cardiovascular disease is the leading cause of death worldwide. Cardiovascular disease presents differently in men and women. While men develop plaque buildup in large arteries, women develop buildup in the microvessels in the heart. Arterial stiffening, which is the hardening of arteries, increases with age in both men and women. Aging, coupled with the decline in sex hormones, exacerbates cardiovascular disease in women compared to men. Men with XY sex chromosomes have higher circulating testosterone, while women with XX sex chromosomes have increased circulating estradiol. The potential benefits of sex hormone replacement therapy are shown in men and women. Indeed, testosterone replacement deficiency is associated with adverse cardiovascular outcomes in men. Whether adverse events are dependent or independent of sex hormones' interaction with sex chromosomes is unknown. This study used the four core genotype mice comprising males with either XX or XY sex chromosome complement. We show castration increases arterial stiffening and collagen deposition on the arterial wall. We also identified the escapee and smooth muscle contractile genes that may play a role in arterial stiffening. Our data suggests that testosterone deprivation mediates arterial stiffening and remodeling.

Why not XY? Male monoecious sexual phenotypes challenge the female monoecious paradigm in Cannabis sativa L.

Monoecy in Cannabis sativa L. has long been considered an industrially important trait due to the increased uniformity it offers and was thought to be exclusively associated with XX females. The isolation and characterisation of a monoecious individual with XY chromosomes sourced from non-proprietary germplasm is reported for the first time. The chromosomal make up of this trait was confirmed through inflorescence structure, growth habit, PCR analysis and sexual phenotypes of progeny from a series of targeted crosses. The identification of an XY monoecious phenotype widens our understanding of monoecy in Cannabis and has important implications for breeding, particularly for producing F1-hybrid seed.

Transcriptional immune suppression and up-regulation of double-stranded DNA damage and repair repertoires in ecDNA-containing tumors.

Extrachromosomal DNA is a common cause of oncogene amplification in cancer. The non-chromosomal inheritance of ecDNA enables tumors to rapidly evolve, contributing to treatment resistance and poor outcome for patients. The transcriptional context in which ecDNAs arise and progress, including chromosomally-driven transcription, is incompletely understood. We examined gene expression patterns of 870 tumors of varied histological types, to identify transcriptional correlates of ecDNA. Here, we show that ecDNA-containing tumors impact four major biological processes. Specifically, ecDNA-containing tumors up-regulate DNA damage and repair, cell cycle control, and mitotic processes, but down-regulate global immune regulation pathways. Taken together, these results suggest profound alterations in gene regulation in ecDNA-containing tumors, shedding light on molecular processes that give rise to their development and progression.

Associations between mosaic loss of sex chromosomes and incident hospitalization for atrial fibrillation in the United Kingdom.

Background: Mosaic loss of chromosome Y (mLOY) in leukocytes of men reflects genomic instability from aging, smoking, and environmental exposures. A similar mosaic loss of chromosome X (mLOX) occurs among women. However, the associations between mLOY, mLOX, and risk of incident heart diseases are unclear. Methods: We estimated associations between mLOY, mLOX, and risk of incident heart diseases requiring hospitalization, including atrial fibrillation, myocardial infarction, ischemic heart disease, cardiomyopathy, and heart failure. We analyzed 190,613 men and 224,853 women with genotyping data from the UK Biobank. Among these participants, we analyzed 37,037 men with mLOY and 13,978 women with mLOX detected using Mosaic Chromosomal Alterations caller. Multivariable Cox regression was used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) of each incident heart disease in relation to mLOY in men and mLOX in women. Additionally, Mendelian randomization (MR) was conducted to estimate causal associations. Results: Among men, detectable mLOY was associated with elevated risk of atrial fibrillation (HR=1.06, 95%CI:1.03-1.11). The associations were apparent in both never-smokers (HR=1.07, 95%:1.01-1.14) and ever-smokers (HR=1.05, 95%CI:1.01-1.11) as well as men > and ≤60 years of age. MR analyses supported causal associations between mLOY and atrial fibrillation (HRMR-PRESSO=1.15, 95%CI:1.13-1.18). Among post-menopausal women, we found a suggestive inverse association between detectable mLOX and atrial fibrillation risk (HR=0.90, 95%CI:0.83-0.98). However, associations with mLOY and mLOX were not found for other heart diseases. Conclusions: Our findings suggest that mLOY and mLOX reflect sex-specific biological processes or exposure profiles related to incident atrial fibrillation requiring hospitalization.

ASAR lncRNAs control DNA replication timing through interactions with multiple hnRNP/RNA binding proteins.

ASARs are a family of very-long noncoding RNAs that control replication timing on individual human autosomes, and are essential for chromosome stability. The eight known ASAR lncRNAs remain closely associated with their parent chromosomes. Analysis of RNA-protein interaction data (from ENCODE) revealed numerous RBPs with significant interactions with multiple ASAR lncRNAs, with several hnRNPs as abundant interactors. An ~7 kb domain within the ASAR6-141 lncRNA shows a striking density of RBP interaction sites. Genetic deletion and ectopic integration assays indicate that this ~7 kb RNA binding protein domain contains functional sequences for controlling replication timing of entire chromosomes in cis. shRNA-mediated depletion of 10 different RNA binding proteins, including HNRNPA1, HNRNPC, HNRNPL, HNRNPM, HNRNPU, or HNRNPUL1, results in dissociation of ASAR lncRNAs from their chromosome territories, and disrupts the synchronous replication that occurs on all autosome pairs, recapitulating the effect of individual ASAR knockouts on a genome-wide scale. Our results further demonstrate the role that ASARs play during the temporal order of genome-wide replication, and we propose that ASARs function as essential RNA scaffolds for the assembly of hnRNP complexes that help maintain the structural integrity of each mammalian chromosome.

The chromatin-associated 53BP1 ortholog, HSR-9, regulates recombinational repair and X chromosome segregation in the Caenorhabditis elegans germ line.

53BP1 plays a crucial role in regulating DNA damage repair pathway choice and checkpoint signaling in somatic cells; however, its role in meiosis has remained enigmatic. In this study, we demonstrate that the Caenorhabditis elegans ortholog of 53BP1, HSR-9, associates with chromatin in both proliferating and meiotic germ cells. Notably, HSR-9 is enriched on the X chromosome pair in pachytene oogenic germ cells. HSR-9 is also present at kinetochores during both mitotic and meiotic divisions but does not appear to be essential for monitoring microtubule-kinetochore attachments or tension. Using cytological markers of different steps in recombinational repair, we found that HSR-9 influences the processing of a subset of meiotic double strand breaks into COSA-1-marked crossovers. Additionally, HSR-9 plays a role in meiotic X chromosome segregation under conditions where X chromosomes fail to pair, synapse, and recombine. Together, these results highlight that chromatin-associated HSR-9 has both conserved and unique functions in the regulation of meiotic chromosome behavior.

The First Identification of Homomorphic XY Sex Chromosomes by Integrating Cytogenetic and Transcriptomic Approaches in Plestiodon elegans (Scincidae).

The sex chromosomes of skinks are usually poorly differentiated and hardly distinguished by cytogenetic methods. Therefore, identifying sex chromosomes in species lacking easily recognizable heteromorphic sex chromosomes is necessary to fully understand sex chromosome diversity. In this paper, we employed cytogenetics, sex quantification of genes, and transcriptomic approaches to characterize the sex chromosomes in Plestiodon elegans. Cytogenetic examination of metaphases revealed a diploid number of 2n = 26, consisting of 12 macrochromosomes and 14 microchromosomes, with no significant heteromorphic chromosome pairs, speculating that the sex chromosomes may be homomorphic or poorly differentiated. The results of the sex quantification of genes showed that Calumenin (calu), COPI coat complex subunit γ 2 (copg2), and Smoothened (smo) were at half the dose in males as in females, suggesting that they are on the X chromosome. Transcriptomic data analysis from the gonads yielded the excess expression male-specific genes (n = 16), in which five PCR molecular markers were developed. Restricting the observed heterozygosity to males suggests the presence of homomorphic sex chromosomes in P. elegans, XX/XY. This is the first breakthrough in the study of the sex chromosomes of Plestiodon.

Consequences of partially recessive deleterious genetic variation for the evolution of inversions suppressing recombination between sex chromosomes1.

The evolution of suppressed recombination between sex chromosomes is widely hypothesized to be driven by sexually antagonistic selection (SA), where tighter linkage between the sex-determining gene(s) and nearby SA loci is favored when it couples male-beneficial alleles to the proto-Y chromosome, and female-beneficial alleles to the proto-X. Although difficult to test empirically, the SA selection hypothesis overshadows several alternatives, including an incomplete but often-repeated "sheltering" hypothesis which suggests that expansion of the sex-linked region (SLR) reduces the homozygous expression of deleterious mutations at selected loci. Here, we use population genetic models to evaluate the consequences of partially recessive deleterious mutational variation for the evolution of otherwise neutral chromosomal inversions expanding the SLR on proto-Y chromosomes. Both autosomal and SLR-expanding inversions face a race against time: lightly-loaded inversions are initially beneficial, but eventually become deleterious as they accumulate new mutations, after which their chances of fixing become negligible. In contrast, initially unloaded inversions eventually become neutral as their deleterious load reaches the same equilibrium as non-inverted haplotypes. Despite the differences in inheritance and indirect selection, SLR-expanding inversions exhibit similar evolutionary dynamics to autosomal inversions over many biologically plausible parameter conditions. Differences emerge when the population average mutation load is quite high; in this case large autosomal inversions that are lucky enough to be mutation-free can rise to intermediate to high frequencies where selection in homozygotes becomes important (Y-linked inversions never appear as homozygous karyotypes); conditions requiring either high mutation rates, highly recessive deleterious mutations, weak selection, or a combination thereof.

The role of conflict in the formation and maintenance of variant sex chromosome systems in mammals.

The XX/XY sex chromosome system is deeply conserved in therian mammals, as is the role of Sry in testis determination, giving the impression of stasis relative to other taxa. However, the long tradition of cytogenetic studies in mammals documents sex chromosome karyotypes that break this norm in myriad ways, ranging from fusions between sex chromosomes and autosomes to Y chromosome loss. Evolutionary conflict, in the form of sexual antagonism or meiotic drive, is the primary predicted driver of sex chromosome transformation and turnover. Yet conflict-based hypotheses are less considered in mammals, perhaps because of the perceived stability of the sex chromosome system. To address this gap, we catalogue and characterize all described sex chromosome variants in mammals, test for family-specific rates of accumulation, and consider the role of conflict between the sexes or within the genome in the evolution of these systems. We identify 152 species with sex chromosomes that differ from the ancestral state and find evidence for different rates of ancestral to derived transitions among families. Sex chromosome-autosome fusions account for 80% of all variants whereas documented sex chromosome fissions are limited to three species. We propose that meiotic drive and drive suppression provide viable explanations for the evolution of many of these variant systems, particularly those involving autosomal fusions. We highlight taxa particularly worthy of further study and provide experimental predictions for testing the role of conflict and its alternatives in generating observed sex chromosome diversity.

SLRfinder: A method to detect candidate sex-linked regions with linkage disequilibrium clustering.

Despite their critical roles in genetic sex determination, sex chromosomes remain unknown in many non-model organisms, especially those having recently evolved sex-linked regions (SLRs). These evolutionarily young and labile sex chromosomes are important for understanding early sex chromosome evolution but are difficult to identify due to the lack of Y/W degeneration and SLRs limited to small genomic regions. Here, we present SLRfinder, a method to identify candidate SLRs using linkage disequilibrium (LD) clustering, heterozygosity and genetic divergence. SLRfinder does not rely on specific sequencing methods or a specific type of reference genome (e.g., from the homomorphic sex). In addition, the input of SLRfinder does not require phenotypic sexes, which may be unknown from population sampling, but sex information can be incorporated and is necessary to validate candidate SLRs. We tested SLRfinder using various published datasets and compared it to the local principal component analysis (PCA) method and the depth-based method Sex Assignment Through Coverage (SATC). As expected, the local PCA method could not be used to identify unknown SLRs. SATC works better on conserved sex chromosomes, whereas SLRfinder outperforms SATC in analysing labile sex chromosomes, especially when SLRs harbour inversions. Power analyses showed that SLRfinder worked better when sampling more populations that share the same SLR. If analysing one population, a relatively larger sample size (around 50) is needed for sufficient statistical power to detect significant SLR candidates, although true SLRs are likely always top-ranked. SLRfinder provides a novel and complementary approach for identifying SLRs and uncovering additional sex chromosome diversity in nature.

The differential virulence of Fusarium strains causing corneal infections and plant diseases is associated with accessory chromosome composition.

Fusarium oxysporum is a cross-kingdom pathogen. While some strains cause disseminated fusariosis and blinding corneal infections in humans, others are responsible for devastating vascular wilt diseases in plants. To better understand the distinct adaptations of F. oxysporum to animal or plant hosts, we conducted a comparative phenotypic and genetic analysis of two strains: MRL8996 (isolated from a keratitis patient) and Fol4287 (isolated from a wilted tomato [Solanum lycopersicum]). In vivo infection of mouse corneas and tomato plants revealed that, while both strains cause symptoms in both hosts, MRL8996 caused more severe corneal ulceration and perforation in mice, whereas Fol4287 induced more pronounced wilting symptoms in tomato. In vitro assays using abiotic stress treatments revealed that the human pathogen MRL8996 was better adapted to elevated temperatures, whereas the plant pathogen Fol4287 was more tolerant of osmotic and cell wall stresses. Both strains displayed broad resistance to antifungal treatment, with MRL8996 exhibiting the paradoxical effect of increased tolerance to higher concentrations of the antifungal caspofungin. We identified a set of accessory chromosomes (ACs) and protein-encoding genes with distinct transposon profiles and functions, respectively, between MRL8996 and Fol4287. Interestingly, ACs from both genomes also encode proteins with shared functions, such as chromatin remodeling and post-translational protein modifications. Our phenotypic assays and comparative genomics analyses lay the foundation for future studies correlating genotype with phenotype and for developing targeted antifungals for agricultural and clinical uses.

SMC-mediated dosage compensation in C. elegans evolved in the presence of an ancestral nematode mechanism.

Mechanisms of X chromosome dosage compensation have been studied extensively in a few model species representing clades of shared sex chromosome ancestry. However, the diversity within each clade as a function of sex chromosome evolution is largely unknown. Here, we anchor ourselves to the nematode Caenorhabditis elegans, for which a well-studied mechanism of dosage compensation occurs through a specialized structural maintenance of chromosomes (SMC) complex, and explore the diversity of dosage compensation in the surrounding phylogeny of nematodes. Through phylogenetic analysis of the C. elegans dosage compensation complex and a survey of its epigenetic signatures, including X-specific topologically associating domains (TADs) and X-enrichment of H4K20me1, we found that the condensin-mediated mechanism evolved recently in the lineage leading to Caenorhabditis through an SMC-4 duplication. Intriguingly, an independent duplication of SMC-4 and the presence of X-specific TADs in Pristionchus pacificus suggest that condensin-mediated dosage compensation arose more than once. mRNA-seq analyses of gene expression in several nematode species indicate that dosage compensation itself is ancestral, as expected from the ancient XO sex determination system. Indicative of the ancestral mechanism, H4K20me1 is enriched on the X chromosomes in Oscheius tipulae, which does not contain X-specific TADs or SMC-4 paralogs. Together, our results indicate that the dosage compensation system in C. elegans is surprisingly new, and condensin may have been co-opted repeatedly in nematodes, suggesting that the process of evolving a chromosome-wide gene regulatory mechanism for dosage compensation is constrained. Significance statement: X chromosome dosage compensation mechanisms evolved in response to Y chromosome degeneration during sex chromosome evolution. However, establishment of dosage compensation is not an endpoint. As sex chromosomes change, dosage compensation strategies may have also changed. In this study, we performed phylogenetic and epigenomic analyses surrounding Caenorhabditis elegans and found that the condensin-mediated dosage compensation mechanism in C. elegans is surprisingly new, and has evolved in the presence of an ancestral mechanism. Intriguingly, condensin-based dosage compensation may have evolved more than once in the nematode lineage, the other time in Pristionchus. Together, our work highlights a previously unappreciated diversity of dosage compensation mechanisms within a clade, and suggests constraints in evolving new mechanisms in the presence of an existing one.

Genetic drift drives faster-Z evolution in the salmon louse Lepeophtheirus salmonis.

Sex chromosome evolution is a complex sub-field of population genetics with unresolved questions about how quickly and adaptively these chromosomes should evolve compared to autosomes. One key limitation to existing knowledge is an intense focus on only a handful of taxa, resulting in uncertainty about whether observed patterns reflect general processes or are idiosyncratic to the more widely-studied clades. In particular, the Z chromosomes of female heterogametic (ZW) systems tend to be quickly but not adaptively evolving in birds, while in butterflies and moths Z chromosomes tend to be evolving adaptively, but not always faster than autosomes. To understand how these two observations fit into broader evolutionary patterns, we explore patterns of Z chromosome evolution outside of these two well-studied clades. We utilize a publicly available high-quality genome, gene expression, population, and outgroup data for the salmon louse Lepeophtheirus salmonis, an important aquacultural pest copepod. We find that the Z chromosome is faster evolving than the autosomes, but that this effect is driven by increased drift rather than adaptive evolution. Due to high rates of female reproductive failure, the Z chromosome exhibits only a slightly lower effective population size than the autosomes which is nonetheless sufficient to decrease efficiency of hemizygous selection acting on the Z. These results highlight the usefulness of organismal life history in calibrating population genetic expectations and demonstrate the value of the ever-expanding wealth of modern publicly available genomic data to help resolve outstanding evolutionary questions.

Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae.

BACKGROUND: The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. RESULTS: We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary. CONCLUSIONS: Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Improved assembly of the Pungitius pungitius reference genome.

The nine-spined stickleback (Pungitius pungitius) has been increasingly used as a model system in studies of local adaptation and sex chromosome evolution but its current reference genome assembly is far from perfect, lacking distinct sex chromosomes. We generated an improved assembly of the nine-spined stickleback reference genome (98.3% BUSCO completeness) with the aid of linked-read mapping. While the new assembly (v8) was of similar size as the earlier version (v7), we were able to assign 4.4 times more contigs to the linkage groups and improve the contiguity of the genome. Moreover, the new assembly contains a ∼22.8 Mb Y-linked scaffold (LG22) consisting mainly of previously assigned X-contigs, putative Y-contigs, putative centromere contigs and highly repetitive elements. The male individual showed an even mapping depth on LG12 (pseudo X chromosome) and LG22 (Y-linked scaffold) in the segregating sites, suggesting near-pure X and Y representation in the v8 assembly. A total of 26,803 genes were annotated, and about 33% of the assembly was found to consist of repetitive elements. The high proportion of repetitive elements in LG22 (53.10%) suggests it can be difficult to assemble the complete sequence of the species' Y chromosome. Nevertheless, the new assembly is a significant improvement over the previous version and should provide a valuable resource for genomic studies of stickleback fishes.