Chromosome Evolution Research Articles

ZW sex chromosome structure in Amborella trichopoda.

Sex chromosomes have evolved hundreds of times across the flowering plant tree of life; their recent origins in some members of this clade can shed light on the early consequences of suppressed recombination, a crucial step in sex chromosome evolution. Amborella trichopoda, the sole species of a lineage that is sister to all other extant flowering plants, is dioecious with a young ZW sex determination system. Here we present a haplotype-resolved genome assembly, including highly contiguous assemblies of the Z and W chromosomes. We identify a ~3-megabase sex-determination region (SDR) captured in two strata that includes a ~300-kilobase inversion that is enriched with repetitive sequences and contains a homologue of the Arabidopsis METHYLTHIOADENOSINE NUCLEOSIDASE (MTN1-2) genes, which are known to be involved in fertility. However, the remainder of the SDR does not show patterns typically found in non-recombining SDRs, such as repeat accumulation and gene loss. These findings are consistent with the hypothesis that dioecy is derived in Amborella and the sex chromosome pair has not significantly degenerated.

Phylogenomics resolves key relationships in Rumex and uncovers a dynamic history of independently evolving sex chromosomes

Abstract Sex chromosomes have evolved independently many times across eukaryotes. Despite a considerable body of literature on sex chromosome evolution, the causes and consequences of variation in their formation, degeneration, and turnover remain poorly understood. Chromosomal rearrangements are thought to play an important role in these processes by promoting or extending the suppression of recombination on sex chromosomes. Sex chromosome variation may also contribute to barriers to gene flow, limiting introgression among species. Comparative approaches in groups with sexual system variation can be valuable for understanding these questions. Rumex is a diverse genus of flowering plants harboring significant sexual system and karyotypic variation, including hermaphroditic and dioecious clades with XY (and XYY) sex chromosomes. Previous disagreement in the phylogenetic relationships among key species has rendered the history of sex chromosome evolution uncertain. Resolving this history is important for investigating the interplay of chromosomal rearrangements, introgression, and sex chromosome evolution in the genus. Here, we use new transcriptome assemblies from 11 species representing major clades in the genus, along with a whole-genome assembly generated for a key hermaphroditic species. Using phylogenomic approaches, we find evidence for the independent evolution of sex chromosomes across two major clades, and introgression from unsampled lineages likely predating the formation of sex chromosomes in the genus. Comparative genomic approaches revealed high rates of chromosomal rearrangement, especially in dioecious species, with evidence for a complex origin of the sex chromosomes through multiple chromosomal fusions. However, we found no evidence of elevated rates of fusion on the sex chromosomes in comparison with autosomes, providing no support for an adaptive hypothesis of sex chromosome expansion due to sexually antagonistic selection. Overall, our results highlight a complex history of karyotypic evolution in Rumex, raising questions about the role that chromosomal rearrangements might play in the evolution of large heteromorphic sex chromosomes.

Stochastic Epigenetic Modification and Evolution of Sex Determination in Vertebrates.

In this report, we propose a novel mathematical model of the origin and evolution of sex determination in vertebrates that is based on the stochastic epigenetic modification (SEM) mechanism. We have previously shown that SEM, with rates consistent with experimental observation, can both increase the rate of gene fixation and decrease pseudogenization, thus dramatically improving the efficacy of evolution. Here, we present a conjectural model of the origin and evolution of sex determination wherein the SEM mechanism alone is sufficient to parsimoniously trigger and guide the evolution of heteromorphic sex chromosomes from the initial homomorphic chromosome configuration, without presupposing any allele frequency differences. Under this theoretical model, the SEM mechanism (i) predated vertebrate sex determination origins and evolution, (ii) has been conveniently and parsimoniously co-opted by the vertebrate sex determination systems during the evolutionary transitioning to the extant vertebrate sex determination, likely acting "on top" of these systems, and (iii) continues existing, alongside all known vertebrate sex determination systems, as a universal pan-vertebrate sex determination modulation mechanism.

Evolution of Sex-linked Genes and the Role of Pericentromeric Regions in Sex Chromosomes: Insights from Diploid Willows.

The evolution of sex chromosomes can involve recombination suppression sometimes involving structural changes, such as inversions, allowing subsequent rearrangements, including inversions and gene transpositions. In the two major genus Salix clades, Salix and Vetrix, almost all species are dioecious, and sex-linked regions have evolved on chromosome 7 and 15, with either male or female heterogamety. We used chromosome conformation capture (Hi-C) and PacBio HiFi (high-fidelity) reads to assemble chromosome-level, gap-free X and Y chromosomes from both clades, S. triandra (15XY system), a basal species in the Vetrix clade, and the Salix clade species S. mesnyi (7XY system). Combining these with other available genome assemblies, we found inversions within the sex-linked regions, which are likely to be pericentromeric and probably recombined rarely in the ancestral species, before sex-linkage evolved. The Y-linked regions in all 15XY and 7XY species include partial duplicates containing exon 1 of an ARR17-like gene similar to male-determining factors in other Salicaceae species. We also found duplicates of a Y-specific gene, which we named MSF. The derived Salix clade 7XY chromosome systems appear to have evolved when these two genes transposed from the 15Y to the 7Y. Additionally, the 7Y chromosomes in S. dunnii and S. chaenomeloides probably evolved from the ancestral 7X of the Salix clade, involving a similar transposition, and loss of the ancestral 7Y. We suggest that pericentromeric regions that recombine infrequently may facilitate the evolution of sex linkage.

One-step leaping evolution from an autosomal pair to the heteromorphic sex chromosomes.

Sex chromosomes evolve from an autosomal pair after the acquisition of a sex-determining gene. The primary sex chromosomes are homomorphic in both sexes and often undergo heteromorphism in either sex (XY in males or ZW in females) in association with chromosome rearrangements such as inversion, which creates a non-recombining region, called a stratum. Then, multiple strata may form by sequential inversions and extend the non-recombining region, where gene divergence accelerates, and degeneration of the Y or W chromosome progressively occurs. Contrast to the conventional theory, we propose a shortcut in heteromorphic sex chromosome evolution, where an autosomal pair directly evolves into a heteromorphic sex chromosome pair. We illustrate this with two frog cases where Y chromosome or autosome, which is morphologically inverted, was introgressed from another species through interspecific hybridization, instantly forming a new heteromorphic sex chromosome pair. This event resulted in a distinct non-recombining region immediately after hybridization. The introduction of an inverted chromosome from a different species may be associated with benefits in morphology, breeding behavior, hybrid viability, sex determination, and recovery of the sex ratio of the hybrids. We discuss the molecular mechanisms driving preferential mutations in the introduced, inverted chromosome through interspecific hybridization.

Rapid chromosome evolution and acquisition of thermosensitive stochastic sex determination in nematode androdioecious hermaphrodites

The factors contributing to evolution of androdioecy, the coexistence of hermaphrodites and males such as in Caenorhabditis elegans, remains poorly known. However, nematodes exhibit androdioecy in at last 13 genera with the predatory genus Pristionchus having seven independent transitions towards androdioecy. Nonetheless, associated genomic architecture and sex determination mechanisms are largely known from Caenorhabditis. Here, studying 47 Pristionchus species, we observed repeated chromosome evolution which abolished the ancestral XX/XO sex chromosome system. Two phylogenetically unrelated androdioecious Pristionchus species have no genomic differences between sexes and mating hermaphrodites with males resulted in hermaphroditic offspring only. We demonstrate that stochastic sex determination is influenced by temperature in P. mayeri and P. entomophagus, and CRISPR engineering indicated a conserved role of the transcription factor TRA-1 in P. mayeri. Chromosome-level genome assemblies and subsequent genomic analysis of related Pristionchus species revealed stochastic sex determination to be derived from XY sex chromosome systems through sex chromosome-autosome fusions. Thus, rapid karyotype evolution, sex chromosome evolution and evolvable sex determination mechanisms are general features of this genus, and represent a dynamic background against which androdioecy has evolved recurrently. Future studies might indicate that stochastic sex determination is more common than currently appreciated.

Fully phased genome assemblies and graph-based genetic variants of the olive flounder, Paralichthys olivaceus

The olive flounder, Paralichthys olivaceus, also known as the Korean halibut, is an economically important flatfish in East Asian countries. Here, we provided four fully phased genome assemblies of two different olive flounder individuals using high-fidelity long-read sequencing and their parental short-read sequencing data. We obtained 42–44 Gb of ~15-kb and ~Q30 high-fidelity long reads, and their assembly quality values were ~53. We annotated ~30 K genes, ~170-Mb repetitive sequences, and ~3 M 5-methylcytosine positions for each genome assembly, and established a graph-based draft pan-genome of the olive flounder. We identified 5 M single-nucleotide variants and 100 K structural variants with their genotype information, where ~13% of the variants were possibly fixed in the two Korean individuals. Based on our chromosome-level genome assembly, we also explored chromosome evolution in the Pleuronectiformes family, as reported earlier. Our high-quality genomic resources will contribute to future genomic selection for accelerating the breeding process of the olive flounder.

Improved chromosome-level donkey (Equus asinus) genome provides insights into genome and chromosome evolution

Improved chromosome-level donkey (Equus asinus) genome provides insights into genome and chromosome evolution

Development and evolution of Drosophila chromatin landscape in a 3D genome context

Little is known about how the epigenomic states change during development and evolution in a 3D genome context. Here we use Drosophila pseudoobscura with complex turnover of sex chromosomes as a model to address this, by collecting massive epigenomic and Hi-C data from five developmental stages and three adult tissues. We reveal that over 60% of the genes and transposable elements (TE) exhibit at least one developmental transition of chromatin state. Transitions on specific but not housekeeping enhancers are associated with specific chromatin loops and topologically associated domain borders (TABs). While evolutionarily young TEs are generally silenced, old TEs more often have been domesticated as interacting TABs or specific enhancers. But on the recently evolved X chromosome, young TEs are instead often active and recruited as TABs, due to acquisition of dosage compensation. Overall we characterize how Drosophila epigenomic landscapes change during development and in response to chromosome evolution, and highlight the important roles of TEs in genome organization and regulation.

Shifts in Chromosome Evolution Rates Shape the Karyotype Patterns of Leafcutting Ants.

Trait evolution has become a central focus in evolutionary biology, with phylogenetic comparative methods offering a framework to study how and why traits vary among species. Identifying variations in trait evolution rates within phylogenies is important for uncovering the mechanisms behind these differences. Karyotype variation, which is substantial across eukaryotic organisms, plays an essential role in species diversification. This study investigates karyotype variation within the leafcutting ant clade, focusing on chromosome number and morphology. We aim to determine whether karyotypic traits are phylogenetically dependent and how different evolutionary models explain karyotype diversity. Previous models have been insufficient in explaining these variations. To address these gaps, we employ modern phylogenetic methods to assess the impact of chromosomal fissions and fusions on karyotype evolution. By evaluating various evolutionary models-particularly the Brownian motion model, which suggests neutral chromosomal changes-we pursue for the further understanding the mode and tempo of karyotype evolution in ants. Our research examines how shifts in chromosomal change rates contribute to divergence among leafcutting ant species and assesses the role of chromosomal changes in the clade's evolutionary trajectory. Comparative analysis of leafcutting ant ideograms suggests that shared karyotype traits are strongly related to species relationships. This implies that karyotype diversification in leafcutting ants follows a phylogenetic trajectory at varying rates, with differences in karyotype traits reflecting the evolutionary distance between lineages. Particularly, the increase in the chromosome number of Acromyrmex is likely due to fission rearrangements rather than demi or polyploidization. We discuss and provide insights into the mechanisms driving karyotype variation and its implications for leafcutting ant diversification.

Delineating the W sex chromosome in the clam shrimp, Eulimnadia texana

Introduction: Sex chromosomes have evolved independently across various lineages, often showing convergent degradation of the sex-limited chromosome. While extensively studied in model organisms with ancient sex chromosomal systems, the evolution of early-stage sex chromosomes remains poorly understood. Eulimnadia texana, a freshwater crustacean with a unique androdioecious breeding system (ZZ, ZW, and viable WW genotypes), provides a rare opportunity to study early sex chromosome evolution. This study examines E. texana's W chromosome for evidence of a small localized non-recombining region, characterized by a transposable element (TE) "hotspot", low gene density, and low GC content. Methods: Sex-linked markers were mapped onto the W chromosome (Scaffold 1). TEs in the WW genome were identified using RepeatModeler and RepeatMasker. Statistical analyses compared TE distribution between the genome and Scaffold 1, which was then divided into 20 equal-sized "bins" for finer-scale statistical analyses. Gene density and GC content were analyzed across these bins. Results: While no significant TE accumulation was found across the entire W chromosome compared to the remaining genome, a specific region (6.6-8.8 Mb, fourth bin) showed significantly higher TE accumulation. This region also exhibited low gene density and low GC content, indicative of reduced recombination. Conclusion: Our findings suggest that E. texana's W chromosome contains a smaller region of crossover suppression, supporting the hypothesis that it is in a proto-sex chromosome in early evolutionary development. This study provides valuable insights into early sex chromosome evolution and establishes E. texana as an ideal model for further investigation of evolutionary processes driving proto-sex chromosome differentiation.

The burst of satellite DNA in Leptidea wood white butterflies and their putative role in karyotype evolution.

Satellite DNAs (satDNAs) are abundant components of eukaryotic genomes, playing pivotal roles in chromosomal organization, genome stability, and evolution. Here, we combined cytogenetic and genomic methods to characterize the satDNAs in the genomes of Leptidea butterflies. Leptidea is characterized by the presence of a high heterochromatin content, large genomes, and extensive chromosomal reshuffling as well as the occurrence of cryptic species. We show that, in contrast to other Lepidoptera, satDNAs constitute a considerable proportion of Leptidea genomes, ranging between 4.11% and 11.05%. This amplification of satDNAs, together with the hyperactivity of transposable elements, contributes to the substantial genome expansion in Leptidea. Using chromosomal mapping, we show that, particularly LepSat01-100 and LepSat03-167 satDNAs, are preferentially localized in heterochromatin exhibiting variable distribution that may have contributed to the highly diverse karyotypes within the genus. The satDNAs also exhibit W-chromosome accumulation, suggesting their involvement in sex chromosome evolution. Our results provide insights into the dynamics of satDNAs in Lepidoptera genomes and highlight their role in genome expansion and chromosomal organization, which could influence the speciation process. The high proportion of repetitive DNAs in the genomes of Leptidea underscores the complex evolutionary dynamics revealing the interplay between repetitive DNAs and genomic architecture in the genus.

Whole-chromosome oligo-painting in licorice unveils interspecific chromosomal evolutionary relationships and possible origin of triploid genome species.

Licorice is one of the most extensively studied medicinal plants in the world, whose roots and rhizomes have long been used as both a sweetener and an essential component in numerous herbal preparations. However, the genus Glycyrrhiza has a complex composition, and the interspecies chromosomal relationships, origin, and evolution are still largely unclear. Here, we develop a set of whole-chromosome painting probes that allowed identification of all eight chromosomes of licorice on same metaphase chromosomes. Comparative chromosome painting analyses in seven different Glycyrrhiza species revealed that the genus Glycyrrhiza maintained extraordinarily conserved chromosomal synteny after about 3-12 million years of divergence. No cytologically visible inter-chromosomal rearrangements were identified in any species. By comparative chromosomal karyotype analyses, we revealed interspecific chromosome evolutionary relationships and dramatic variable chromosomal karyotype after independent divergence and demonstrated that G. prostrate was the most closely related to the ancestral type among the seven Glycyrrhiza species. Furthermore, we also discovered a G. glandulosa seed with distinct triploid-genome for the first time in China, suggesting the existence of a polyploid evolutionary pathway in the genus Glycyrrhiza, which challenges the previous notion that only diploids of licorice existed in nature. This study expands our knowledge of the chromosome evolution of licorice and will lay an important foundation for the genome origin and evolution studies in the genus Glycyrrhiza.

Expanding the Triangle of U: Comparative analysis of the Hirschfeldia incana genome provides insights into chromosomal evolution, phylogenomics and high photosynthesis-related traits.

The Brassiceae tribe encompasses many economically important crops and exhibits high intraspecific and interspecific phenotypic variation. After a shared whole-genome triplication (WGT) event (Br-α, ~15.9 million years ago), differential lineage diversification and genomic changes contributed to an array of divergence in morphology, biochemistry, and physiology underlying photosynthesis-related traits. Here, the C3 species Hirschfeldia incana is studied as it displays high photosynthetic rates under high-light conditions. Our aim was to elucidate the evolution that gave rise to the genome of H. incana and its high-photosynthesis traits. We reconstructed a chromosome-level genome assembly for H. incana (Nijmegen, v2.0) using nanopore and chromosome conformation capture (Hi-C) technologies, with 409Mb in size and an N50 of 52Mb (a 10× improvement over the previously published scaffold-level v1.0 assembly). The updated assembly and annotation was subsequently employed to investigate the WGT history of H. incana in a comparative phylogenomic framework from the Brassiceae ancestral genomic blocks and related diploidized crops. Hirschfeldia incana (x=7) shares extensive genome collinearity with Raphanus sativus (x=9). These two species share some commonalities with Brassica rapa and B. oleracea (A genome, x=10 and C genome, x=9, respectively) and other similarities with B. nigra (B genome, x=8). Phylogenetic analysis revealed that H. incana and R. sativus form a monophyletic clade in between the Brassica A/C and B genomes. We postulate that H. incana and R. sativus genomes are results of hybridization or introgression of the Brassica A/C and B genome types. Our results might explain the discrepancy observed in published studies regarding phylogenetic placement of H. incana and R. sativus in relation to the "Triangle of U" species. Expression analysis of WGT retained gene copies revealed sub-genome expression divergence, likely due to neo- or sub-functionalization. Finally, we highlighted genes associated with physio-biochemical-anatomical adaptive changes observed in H. incana which likely facilitate its high-photosynthesis traits under high light. The improved H. incana genome assembly, annotation and results presented in this work will be a valuable resource for future research to unravel the genetic basis of its ability to maintain a high photosynthetic efficiency in high-light conditions and thereby improve photosynthesis for enhanced agricultural production.

High prevalence of chromosome 17 in breast cancer micronuclei: a means to get rid of tumor suppressors?

Micronuclei (MN), defined as small extra-nuclear chromatin bodies enclosed by a nuclear envelope, serve as noticeable markers of chromosomal instability (CIN). The MN have been used for breast cancer (BC) screening, diagnosis, and prognosis. However, more recently they have gained attention as seats for active chromosomal rearrangements. BC subtypes exhibit differential CIN levels and aggressiveness. This study aimed to investigate MN chromosomal contents across BC subtypes, exploring its potential role in aggressiveness and pathogenesis. Immunostaining of BC cells was performed with anti-centromeric antibody followed by confocal microscopy. Further, fluorescence in situ hybridization (FISH) was done to check the presence of specific chromosomes in the MN. The real time PCR was also done from the RNA isolated from MN to check the expression of TP53 gene. BC cell lines (CLs) showed the presence of both centromere-positive ( +) and -negative ( -) MN, with significant variation in frequency among hormone and human epidermal growth factor receptor positive and triple-negative (TN) BC cells. FISH targeting chromosomes 1, 3, 8, 11, and 17 detected centromeric signals for all the above chromosomes in MN with a relatively higher prevalence of chromosome 17 in all the CLs. Out of all the CLs, TNBC cells demonstrated the highest frequency of centromere + and chromosome 17 + MN. TP53 expression could also be demonstrated inside the MN by FISH and real time PCR. Patient sample imprints also confirmed the presence of chromosome 17 in MN with polysomy of the same in corresponding nuclei. The high prevalence of chromosome 17 in BC MN may connote the importance of its rearrangements in the pathogenesis of BC. Further, the higher prevalence of chromosome 17 and 1 signals in TNBC MN point towards the significance of pathogenetic events involving the genes located in these chromosomes in evolution of this more aggressive phenotype.

Sex Chromosome Evolution: Hallmarks and Question Marks.

Sex chromosomes are widespread in species with separate sexes. They have evolved many times independently and display a truly remarkable diversity. New sequencing technologies and methodological developments have allowed the field of molecular evolution to explore this diversity in a large number of model and nonmodel organisms, broadening our vision on the mechanisms involved in their evolution. Diverse studies have allowed us to better capture the common evolutionary routes that shape sex chromosomes; however, we still mostly fail to explain why sex chromosomes are so diverse. We review over half a century of theoretical and empirical work on sex chromosome evolution and highlight pending questions on their origins, turnovers, rearrangements, degeneration, dosage compensation, gene content, and rates of evolution. We also report recent theoretical progress on our understanding of the ultimate reasons for sex chromosomes' existence.

Comparative analysis of chromosome-level genomes provides insights into chromosomal evolution in Chiroptera.

Chiroptera (bats) presents a fascinating model due to its remarkable variation in chromosome numbers, which range from 14 to 62. This astonishing diversity makes bats an excellent subject for studying chromosome evolution. The black-bearded tomb bat (Taphozous melanopogon) occupies a pivotal phylogenetic position within Chiroptera, emphasizing its crucial role in the systematic examination of bat chromosome evolution. In this study, we present the first chromosome-level genome of T. melanopogon within the family Emballonuridae. Together with previously published genomes, we construct a strongly supported phylogenetic tree of bats, which supports that Emballonuridae forms a basal group within Yangochiroptera. Furthermore, we reconstruct ancestral karyotypes at key nodes along the bat phylogeny and conduct a synteny analysis among the genomes of 12 bat species. Our findings identified evolutionary breakpoint regions (EBRs) that are of particular interest. Notably, some bat genomes exhibit an enrichment of genes related to host defense against microbial pathogens within EBRs. Remarkably, one species possesses multiple copies of some β-defensin genes, while six other species have experienced the loss of some β-defensin genes due to EBRs. Furthermore, some olfactory receptor genes are located in EBRs of 12 species, 4 of which have a significant enrichment in sensory perception of smell. Together, our comparative genomic analysis underscores the potential link between chromosome rearrangements and the adaptation of bats to defend against microbial pathogens.

Frequencies of house fly proto-Y chromosomes across populations are predicted by temperature heterogeneity within populations.

Sex chromosomes often differ between closely related species and can even be polymorphic within populations. Species with multifactorial sex determination segregate for multiple different sex determining loci within populations, making them uniquely informative of the selection pressures that drive the evolution of sex chromosomes. The house fly (Musca domestica) is a model species for studying multifactorial sex determination because male determining genes have been identified on all six of the chromosomes, which means that any chromosome can be a "proto-Y". Natural populations of house fly also segregate for a recently derived female-determining locus, meaning house flies also have a proto-W chromosome. The different proto-Y chromosomes are distributed along latitudinal clines on multiple continents, their distributions can be explained by seasonality in temperature, and they have temperature-dependent effects on physiological and behavioral traits. It is not clear, however, how the clinal distributions interact with the effect of seasonality on the frequencies of house fly proto-Y and proto-W chromosomes across populations. To address this question, we measured the frequencies of house fly proto-Y and proto-W chromosomes across nine populations in the United States of America. We confirmed the clinal distribution along the eastern coast of North America, but it is limited to the eastern coast. In contrast, annual mean daily temperature range predicts proto-Y chromosome frequencies across the entire continent. Our results therefore suggest that temperature heterogeneity can explain the distributions of house fly proto-Y chromosomes in a way that does not depend on the cline.

Long-term hematopoietic dysfunction in patients with large-scale mitochondrial DNA deletion syndromes.

Pearson syndrome (PS) and Kearns-Sayre syndrome (KSS) are single large-scale mitochondrial DNA deletion (SLSMD) syndromes. PS is characterized by severe, transient childhood cytopenia, whereas KSS typically manifests later in life without hematologic abnormalities. Despite distinct clinical presentations, both share a common mitochondrial DNA deletion. Recent observations suggest a potential link between PS progression and myeloid malignancy development, indicating that bone marrow failure (BMF) may be a key aspect of PS pathology and potentially universal across SLSMDs. This study explores longitudinal hematological manifestations of SLSMD syndromes, focusing on bone marrow (BM) dysfunction. Sixteen patients with SLSMDs (13 PS and 3 KSS) were followed, of whom 75% experienced cytopenia, necessitating blood transfusions in 56%. Despite achieving transfusion independence at a median age of 24 months, persistent hematological abnormalities were noted. Comprehensive longitudinal BM studies were conducted in 62% of subjects and consistently revealed signs of marrow dysfunction, even without concurrent cytopenia. Median BM cellularity at a median age of four years and eight months was 50%, with histological signs of dyserythropoiesis, abnormal megakaryocytes, and signs suggesting myelodysplasia. Reduced CD34+ counts and BM colony-forming unit capacity were noted, alongside chromosome 7 aberrations in 16% of patients on cytogenetic studies. Our findings establish BM dysfunction as a persistent hallmark of SLSMD syndromes, posing a risk of clonal evolution and acquisition of chromosome 7 aberrations. This aligns with recent literature, emphasizing enduring BMF in SLSMD syndromes and advocating for tailored hematological monitoring guidelines for this unique patient cohort.

Genome-Wide Association Studies (GWAS) and Transcriptome Analysis Reveal Male Heterogametic Sex-Determining Regions and Candidate Genes in Northern Snakeheads (Channa argus).

The Northern snakehead (Channa argus) is a significant economic aquaculture species in China. Exhibiting sexual dimorphism in the growth rate between females and males, mono-sex breeding holds substantial value for aquaculture. This study employed GWAS and transcriptome analysis were applied to identify sex determination genomic regions and develop sex-specific markers. A total of 270 single-nucleotide polymorphisms (SNPs) and 31 insertion-deletions (InDels) were identified as being sexually dimorphic through GWAS and fixation index (Fst) scanning. Based on GWAS results, two sex-specific InDel markers were developed, effectively distinguishing genetic sex for XX females, XY males, and YY super-males via (polymerase chain reaction) PCR amplification. A major genomic segment of approximately 115 kb on chromosome 3 (Chr 03) was identified as the sex-determination region. A comparative transcriptome analysis of gonads for three sexes identified 158 overlapping differentially expressed genes (DEGs). Additionally, three sex-related candidate genes were identified near the sex determination region, including id2, sox11, and rnf144a. Further studies are required to elucidate the functions of these genes. Overall, two sex-specific InDel markers support a male heterogametic XX/XY sex-determination system in Northern snakeheads and three candidate genes offer new insights into sex determination and the evolution of sex chromosomes in teleost fish.