The mitochondrial genome is a compact, closed circular, double-stranded DNA molecule that sits outside the cell nucleus. In vertebrates it typically spans 16 to 20 kilobases and carries 37 genes: 13 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes, plus a non-coding control region that regulates replication and transcription. Because it is inherited maternally, lacks introns, evolves relatively quickly, and exists in many copies per cell, mitochondrial DNA has become one of the most widely used molecular markers in evolutionary genetics, species delineation, and conservation biology. For fish taxonomists wrestling with groups of closely related and morphologically similar species, a complete mitogenome can be the difference between a resolved phylogeny and a tangle of ambiguous branches.

The sequencing effort began with tissue samples collected from Dal Lake and preserved in 95 percent ethanol. After DNA extraction and quantification with a Qubit 4.0 fluorometer, the team prepared sequencing libraries with an average insert size of 379 base pairs and ran them on an Illumina Novaseq6000 platform using paired-end 150-base-pair chemistry, generating roughly four to five gigabytes of data per sample. The reads were assembled de novo with the GetOrganelle toolkit against an animal mitochondrial database, and the finished sequence was deposited in the NCBI GenBank under accession number PP778510. Annotation of protein-coding, tRNA, and rRNA genes was carried out with MitoFish and MitoAnnotator, while gene arrangement was visualized with the CGView web server.

The resulting mitogenome of S. curvifrons measures 16,592 base pairs and contains the standard vertebrate complement of 37 genes plus one control region. Twelve of the thirteen protein-coding genes sit on the heavy strand, with only ND6 on the light strand, a pattern consistent with other published fish mitogenomes. Fourteen of the 22 tRNA genes are encoded on the heavy strand and eight on the light strand, and both rRNA genes are on the heavy strand. The 12S and 16S ribosomal RNA genes measure 956 and 1,679 base pairs respectively, and the tRNA genes range from 69 to 73 nucleotides in length. Twenty-one of the 22 tRNAs fold into the canonical cloverleaf secondary structure; the exception, trnS1, lacks the D arm, a well-known deviation in animal mitochondria.

Base composition analysis revealed a genome biased toward adenine and thymine, with an overall A plus T content of about 55 percent. Among the protein-coding genes, ATP8 showed the highest A plus T content at 64.2 percent, while ND1 had the lowest at 51.1 percent. Skewness analysis showed negative GC skews and positive AT skews across the genome, a pattern conserved across vertebrates and consistent with other Schizothorax species such as S. prenanti, S. oconnori, and S. waltoni. Notably, the ND6 gene displayed a positive GC skew, unlike the other 12 protein-coding genes. The authors suggest this could reflect an adaptive response to environmental stresses such as high reactive oxygen species levels in oxygen-rich habitats, a condition relevant to high-altitude fish like the schizothoracines.

Codon usage followed familiar trends: the codons for methionine, tryptophan, glutamine, and lysine were the most frequently used, while leucine, arginine, and serine were encoded by the largest number of distinct codons. Most protein-coding genes begin with the standard ATG start codon, but the team found variation at the COXI gene, which begins with GTG, echoing similar observations in other bony fishes. The control region, the most variable segment of the mitogenome and a key player in replication and transcription, spans roughly 936 to 938 base pairs and sits between the tRNA-Pro and tRNA-Phe genes, with an AT content of 66.2 percent. The genome also contains 13 intergenic spacers totaling 69 base pairs and six overlapping regions totaling 22 base pairs, with the longest overlaps occurring among the ATP6, ND4L, and ND5 genes.

To place S. curvifrons in its evolutionary context, the researchers constructed phylogenetic trees from complete mitochondrial genome sequences using maximum likelihood and neighbor-joining methods with the Kimura 2-parameter model in MEGA 11, supported by 1000 bootstrap replicates. The comparison set included other Schizothorax species, Indian major carps such as catla, rohu, and mrigal, salmonids including Atlantic salmon and rainbow trout, zebrafish, and outgroups as distant as human and rat. The analysis showed that S. curvifrons clusters closely with S. niger, sharing a confined evolutionary connection within the same clade, while S. esocinus and S. labiatus form a distinct lineage. S. plagistomus separated into its own clade, distinct from the group containing S. esocinus, S. labiatus, S. niger, and S. curvifrons.

The broader picture that emerges is one of clear evolutionary separation. All Schizothorax species form a clade distinct from other cyprinid genera such as Labeo, Tor, and Cirrhinus, which each form well-supported monophyletic groups of their own. The schizothoracines also split sharply from the Salmonidae, a divergence the authors link in part to contrasting life histories: salmonids are often anadromous, migrating between freshwater and the ocean, while Schizothorax species are strictly freshwater endemics of Himalayan drainages. Outgroup species such as zebrafish, rainbow trout, and Atlantic salmon showed high genetic distances, confirming their distant relationship to the core cyprinid lineage, and the gap extends all the way to higher vertebrates like humans and rats.

Why does this matter beyond the phylogenetic trees? The schizothoracine fishes comprise roughly 100 species across 15 genera, and five species are recognized from Kashmir alone: S. niger, S. esocinus, S. plagiostomus, S. curvifrons, and S. labiatus. Their taxonomic validity has long been debated, and resolving species boundaries is a prerequisite for effective management. Many of these fish are now listed as vulnerable on the IUCN Red List of Threatened Species as overfishing and habitat degradation erode their populations. S. curvifrons itself can reach 56 centimeters and 1.3 kilograms, spawning in tributary streams over rocky and sandy substrates across the Indus River system, including the Chenab, Jhelum, Ravi, and Sutlej rivers and the Dal and Manasbal lakes. A complete reference mitogenome gives conservation geneticists a tool for population structure studies, species identification, and phylogeographic reconstruction.

The study also highlights how much genome length can vary even among close relatives. The S. curvifrons mitogenome differs in size from those of S. esocinus at 16,585 base pairs, S. labiatus at 16,583, S. niger at 16,582, and S. plagiostomus at 16,576, differences the authors attribute to variation in intergenic regions, repetitive sequences, and other structural elements. Such length variation, together with the conserved gene order and base composition shared with other bony fishes, paints a picture of a genome that is architecturally conservative yet subtly divergent in ways that track species boundaries. For a fish that has been commercially and culturally important in Kashmir for generations, the first complete mitochondrial genome is more than a technical milestone. It is a baseline dataset intended to guide future research, resolve lingering ambiguities among the four Schizothorax species of the Kashmir region, and inform conservation strategies aimed at keeping these ancient snow trout swimming in Himalayan waters for generations to come.

Subject of Research: Complete mitochondrial genome sequencing and phylogenetic analysis of the snow trout Schizothorax curvifrons

Article Title: First report of complete mitochondrial genome of Schizothorax curvifrons using next generation sequencing platform (NGS) and phylogentic relationship with other Schizothoracinae

Article References: Agade, A. K., Somasundaram, I., Khan, I. A., Bhat, I. A., Ahmad, I., Qadiri, S. S. N., Bhat, B. A., Bazaz, A. I., Dar, B. N., Aga, A. F., & Rather, M. A. (2025). First report of complete mitochondrial genome of Schizothorax curvifrons using next generation sequencing platform (NGS) and phylogentic relationship with other Schizothoracinae. Blue Biotechnology, 2(1), Article 18. https://doi.org/10.1186/s44315-025-00039-5

Image Credits: AI Generated

DOI: 10.1186/s44315-025-00039-5

Keywords: Schizothorax curvifrons, mitochondrial genome, next-generation sequencing, phylogenetics, Schizothoracinae, snow trout, Kashmir, fish conservation, cyprinidae, mtDNA, molecular taxonomy, Dal Lake

The research, led by Zheng-Hang Liu and Fang-Yuan Xiong, with Zhong-Wei Wang as corresponding author, brought together teams from the Scientific Institute of Pearl River Water Resources Protection in Guangzhou and the State Key Laboratory of Breeding Biotechnology and Sustainable Aquaculture at the Institute of Hydrobiology, Chinese Academy of Sciences in Wuhan, along with conservation centers in Guigang dedicated to protecting rare Hongshui River fish. Their motivation was both scientific and practical. Without knowing how sex is determined in A. longipinnis, breeders cannot develop monosex production strategies, the gold standard for exploiting the growth advantage of females. And without a mapped sex chromosome system, the search for the master switch that sets off male or female development remains a needle-in-haystack problem.

The team’s method of choice was 2b-RAD sequencing, a streamlined genotyping approach that has become a favorite tool for fish geneticists working on species with little or no reference genomic information. The technique exploits a specific type of restriction enzyme to slice the genome at thousands of precisely defined positions, capturing short, uniform tags of DNA adjacent to each cut site. Sequencing these tags in a large panel of individuals generates a dense but manageable snapshot of genetic variation across the genome. Because the fragments are so consistent in length and composition, they can be compared across dozens or hundreds of animals to find sequences that appear in one sex and are completely absent in the other, which is precisely the signature expected of DNA residing on a Y chromosome.

Before any sex-linked screening could begin, however, the researchers needed a genomic scaffold on which to place their findings. They assembled a reference genome from a female individual, using k-mer depth distribution analysis to estimate genome size and assess heterozygosity, smudgeplot analysis to confirm ploidy, and BUSCO assessments to verify that the assembly captured a complete set of conserved single-copy genes. Quality and continuity metrics were visualized for two haplotypes, designated Hap-A and Hap-B, and predicted genes were functionally annotated against multiple databases, with circos plots providing global views of each haplotype’s landscape. The assembled female genome contained 50 chromosomes, consistent with the karyology of many cyprinids, and provided the coordinate system into which any male-specific sequence could ultimately be anchored.

The screen itself proved remarkably productive. Out of the many thousands of 2b-RAD tags compared between XY males and XX females, exactly three male-specific sequences emerged, designated ref65115, ref94872, and ref134458. In genomics, three hits might sound like a meager harvest, but in sex chromosome research it is precisely the kind of sparse, consistent signal that indicates a young or recently differentiated sex chromosome rather than a sprawling region of male-limited DNA. Each of the three sequences was detected only in males, never in females, and their joint presence pointed to a single shared origin on a male-specific segment of the genome.

The decisive step came when the researchers mapped these three sequences against the female reference genome using synteny and coverage depth analyses. All three aligned to a single contiguous window between positions 5,724,913 and 5,869,801 base pairs, a stretch of roughly 145 kilobases. Crucially, the alignment pattern indicated that the sequences did not correspond to an allele present in both sexes but instead represented a Y-specific insertion: a block of DNA that has been added to, or retained exclusively on, the Y chromosome and is entirely missing from the X. In other words, the male sex chromosome of A. longipinnis carries a discrete inserted segment at this genomic address, and that insertion is what separates males from females at the molecular level.

Confirmation followed through rigorous population-level validation. Using PCR primers designed from the flanking regions of the three markers, the team generated four diagnostic amplification assays, named ref65115-IN, ref94872-3, ref134458-1, and ref94872-5. On agarose gels, ref65115-IN produced a male-specific band of 217 base pairs, ref94872-3 a male-specific band of 335 base pairs, ref134458-1 a male-specific band of 265 base pairs, and ref94872-5 a shared band of 336 base pairs in both sexes together with a male-specific band of 211 base pairs. Each assay was then tested across 72 cultured individuals drawn from two geographically distinct populations, 18 males and 18 females each from Guangxi and Guizhou Provinces. The male-specific bands appeared in every male and in no female, providing strong evidence that the Y-specific insertion is shared across the species’ range rather than being a local variant confined to one river system.

Those results settle a question that has quietly constrained both the biology and the farming of this species. An XX/XY system means that sex is determined by which chromosome a sperm carries: females inherit two X chromosomes, males one X and one Y. In evolutionary terms, the finding is also intriguing. Many fish lineages show fluid, rapidly evolving sex determination, with different species even within the same genus employing XY, ZW, or environmental systems. Identifying a Y-specific insertion in A. longipinnis adds another data point to this comparative landscape, suggesting that a discrete, recent insertion event on the Y chromosome may have been the founding mutation of male determination in this lineage. That kind of simple origin story is increasingly common in fish genomics, where sex chromosomes are often young enough that the male-specific region is small and tractable.

For the aquaculture industry, the practical implications are immediate and considerable. Because females of A. longipinnis grow significantly faster than males, a monosex, all-female production system would raise yields, shorten production cycles, and reduce feed costs per kilogram of fish, much as all-female Nile tilapia and all-female grass carp culture have done in other systems. The newly validated PCR markers now give breeders a direct molecular test for genotypic sex: rather than waiting months for fish to reach a size at which gonads can be inspected, hatcheries can sex larvae or juveniles from a tiny tissue sample and identify XX females, XY males, and, in future work, potentially phenotypic sex-reversed individuals. Such marker-assisted selection is the foundation on which sex-control breeding programs are built, allowing the production of neomales, or XX individuals carrying testes, whose sperm carry only X chromosomes and who therefore sire all-female offspring when crossed with ordinary XX females.

The study’s authors are careful to frame their work as a beginning rather than an end. The three male-specific markers pin down the sex-determining region to a Y-specific insertion of about 145 kilobases on the female reference assembly, but the region still requires fine mapping against the male genome, and the identity of the master sex-determining gene within or near that insertion remains unknown. Future work will aim at chromosome-scale assignment of the sex chromosomes, high-resolution mapping of the sex-determining region, and ultimately the cloning of the master sex-determining gene, the single locus whose presence or absence tips embryonic development toward testis or ovary. Because the Y-specific region appears compact, that cloning effort is unusually feasible for a non-model fish, and the newly assembled, well-annotated female genome with its 50 chromosomal scaffolds provides exactly the platform such work requires. All experimental procedures were approved by the Animal Ethics Committee of the Institute of Hydrobiology, Chinese Academy of Sciences, and fish were anesthetized with MS-222 before tissue collection, reflecting the ethical standards applied throughout the study. Funded by the Scientific Institute of Pearl River Water Resources Protection, the research is open access under a Creative Commons license, ensuring that breeders, conservationists, and comparative genomics researchers across the field can build on these markers without restriction. For a species prized on dinner tables across southwestern China, and for the farmers who raise it, a 145-kilobase stretch of Y-chromosomal DNA has just become one of the most valuable sequences in the entire genome.

Subject of Research: Identification of male-specific DNA markers and an XX/XY sex determination system in Acrossocheilus longipinnis

Article Title: Identification of male-specific markers via 2b-RAD sequencing reveals an XX/XY sex determination system in Acrossocheilus longipinnis

Article References: Liu, Z.-H., Xiong, F.-Y., Zhou, Y.-L., E, Z.-C., Li, M., Dong, S.-H., Tian, J., Zhen, P., Gui, J.-F., & Wang, Z.-W. (2026). Identification of male-specific markers via 2b-RAD sequencing reveals an XX/XY sex determination system in Acrossocheilus longipinnis. BMC Genomics. https://doi.org/10.1186/s12864-026-13372-9

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13372-9

Keywords: Acrossocheilus longipinnis, sex determination, XX/XY system, 2b-RAD sequencing, male-specific markers, genome assembly, Y-specific insertion, sex-determining region, monosex breeding, aquaculture, marker-assisted selection, cyprinidae