Genomic and genetic insights into speciation and pigment pattern diversification in Danio fishes

Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

Roles for asip2b have been found in cichlid stripe patterning and fin leucophores of Danio rerio.

“Roles for asip2b have been found in cichlid stripe patterning52 as well as fin leucophores of Danio rerio53.”

• cites [52] sourceDisputed

  Citation 52 attaches to the cichlid stripe-patterning claim, not the separate claim about Danio rerio fin leucophores. The paper identifies agouti-related peptide 2 (agrp2), not asip2b, as the gene controlling cichlid stripes; its main text also distinguishes agrp2 from its paralogs (section S; page not provided). (LLM, conf=0.88)

  Checked by a reasoning model — may be wrong.

  Experimental — accuracy not yet measured
• cites [53] 2025, Agouti and BMP signaling drive a naturally occurring fate conversion of melanophores to leucophores in zebrafish, Proceedings of the National Academy of Sciences sourceSupported

Probably not a substantial claim: may not pass S3 (load-bearing).

Igsf11 is an immunoglobulin superfamily adhesion protein required for migration and survival of melanophores.

“Igsf11 is an immunoglobulin superfamily adhesion protein required for migration and survival of melanophores14.”

• cites [14] sourceSupported

Probably not a substantial claim: may not pass S3 (load-bearing).

Gja4 and Gja5b are gap junction proteins that are required autonomously by melanophores and xanthophores for homotypic and heterotypic interactions.

“Gja4 and Gja5b are gap junction proteins that are required autonomously by melanophores and xanthophores for homotypic and heterotypic interactions, i.e., interactions between the same or different pigment cell types, respectively15,16.”

• cites [15] sourcePartly supported

  The Results section on Cx41.8 and Cx39.4 reports that leo and luc function is required in xanthophores and melanophores; transplantation experiments specifically support a cell-autonomous requirement in xanthophores. The Discussion supports homotypic melanophore interactions but says melanophore–xanthophore interaction via gap junctions might occur and that direct interaction in vivo is not known. The claim presents that heterotypic interaction as established, and our copy does not identify the connexins by the names Gja4 and Gja5b. (LLM, conf=0.85)

  Experimental — accuracy not yet measured
• cites [16] sourcePartly supported

  The Results (page unavailable) show that Cx41.8 and Cx39.4 are expressed in both melanophores and xanthophores, while the Discussion says they function together in interactions between those cell types. But the Results explicitly say GJA4 connexins were not analyzed, and the paper does not establish an autonomous requirement for either cell type or a role in interactions between pigment cells of the same type. Its ‘homotypic’ experiments concern connexin channels in oocytes, not same-type pigment-cell interactions. (LLM, conf=0.88)

  Experimental — accuracy not yet measured

May combine two claims.

kcnj13 encodes a potassium channel required in melanophores for interactions with xanthophores and iridophores, regulating the shapes of all three pigment cell types.

“The kcnj13 gene encodes a potassium channel required in melanophores for interactions with xanthophores and iridophores, regulating the shapes of all three pigment cell types10–13.”

• cites [10] sourceNot checked — no full text of the source
• cites [11] sourcePartly supported

  The Results identify Kir7.1 as a potassium channel expressed in melanophores and xanthophores. The Discussion says it is thought to contribute to interactions between those cells, but does not establish that it is required for those interactions. The cited passages do not address interactions with iridophores or regulation of the shapes of all three pigment cell types. Page numbers are unavailable in our copy. (LLM, conf=0.88)

  Experimental — accuracy not yet measured
• cites [12] sourceNot checked — no full text of the source
• cites [13] sourceSupported

May combine two claims.

The igsf11, gja4 and gja5b genes are essential for pigment patterning in D. aesculapii as in zebrafish but do not contribute to the patterning differences between these species.

“The igsf11, gja4 and gja5b genes are essential for pigment patterning in D. aesculapii as in zebrafish but do not contribute to the patterning differences between these species9.”

• cites [9] sourceSupported

May combine two claims.

D-statistics suggested possible gene flow between D. rerio and D. aesculapii, D. rerio and D. nigrofasciatus, and also D. kyathit and D. nigrofasciatus.

“D-statistics (four-taxon)38 suggested possible gene flow between D. rerio and D. aesculapii, D. rerio and D. nigrofasciatus, and also D. kyathit and D. nigrofasciatus (Fig. 3D).”

• cites [38] sourceSupported

Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

The intracellular loop domain of Gja5b may regulate channel gating and trafficking, and other interactions.

“the intracellular loop domain of Gja5b may regulate channel gating and trafficking, and other interactions58.”

• cites [58] sourceCould not verify

  The Introduction describes intracellular loops in G-protein-coupled receptors, and the Life Cycle of FSH Receptor section discusses receptor trafficking. Neither passage addresses Gja5b or whether its intracellular loop regulates channel gating, trafficking, or other interactions (searched: Gja5b, connexin, channel gating, trafficking). (LLM, conf=0.76)

  Experimental — accuracy not yet measured

May combine two claims.

The ancestral pigmentation pattern of the Danio genus is likely horizontal stripes, while vertical bars and spotted patterns recurrently evolved in different Danio lineages over the course of evolution.

“Our results, together with previous research8, provide evidence that the ancestral pigmentation pattern of the Danio genus is likely horizontal stripes, while vertical bars and spotted patterns recurrently evolved in different Danio lineages over the course of evolution.”

• cites [8] sourceSupported

May combine two claims.

Wild-type hybrids between D. rerio and other Danio species generally display a pattern of horizontal stripes, similar to the D. rerio pattern, with stripe aberrations sometimes evident likely depending on genetic backgrounds of different species isolates.

“Wild-type hybrids between D. rerio and other Danio species generally display a pattern of horizontal stripes6,9,46,47, similar to the D. rerio pattern, with stripe aberrations sometimes evident likely depending on !genetic backgrounds of different species isolates, as depicted here for D. tinwini and D. albolineatus hybrids (Fig. 5A).”

• cites [6] sourceNot checked — no full text of the source
• cites [9] sourcePartly supported

  Split into 2 atoms — 1 confirmed, 1 partially supported, 0 disputed, 0 misattributed, 0 unverifiable. Expand the claim to see per-atom verdicts.

  Experimental — accuracy not yet measured
• cites [46] sourcePartly supported

  Split into 2 atoms — 1 confirmed, 1 partially supported, 0 disputed, 0 misattributed, 0 unverifiable. Expand the claim to see per-atom verdicts.

  Experimental — accuracy not yet measured
• cites [47] sourcePartly supported

  Split into 2 atoms — 1 confirmed, 1 partially supported, 0 disputed, 0 misattributed, 0 unverifiable. Expand the claim to see per-atom verdicts.

  Experimental — accuracy not yet measured

Classified as probably not a substantial claim. May combine two claims.

The C-terminal domain of Gja5b mediates anchoring and localization by stabilizing the gap junction plaque and facilitating interactions with regulatory proteins like ZO-1, and also modulates channel gating.

“The C-terminal domain of Gja5b mediates anchoring and localization by stabilizing the gap junction plaque and facilitating interactions with regulatory proteins like ZO-1, and also modulates channel gating55–57”

• cites [55] sourcePartly supported

  The Discussion (page unavailable) says the C-terminal domain of Cx43 is thought to regulate gap-junction assembly and gating and to interact with ZO-1; it also describes a proposed role for ZO-1 in recruiting connexons to plaques. The paper concerns Cx43, not Gja5b, and does not establish that Gja5b’s C-terminal domain anchors or stabilizes gap-junction plaques. (LLM, conf=0.86)

  Experimental — accuracy not yet measured
• cites [56] sourceNot checked — no full text of the source
• cites [57] sourcePartly supported

  The Abstract and the section comparing Cx40CT and Cx43CT interactions with ZO-1 (page not provided) support a role for the Cx40 C-terminal domain in channel regulation and show that it binds ZO-1. They do not establish that Gja5b’s C-terminal domain anchors or localizes channels by stabilizing plaques; the plaque discussion concerns Cx43 and describes a possible effect on plaque formation and size, not stabilization. (LLM, conf=0.86)

  Experimental — accuracy not yet measured

May combine two claims.

The difficulty finding a single, unambiguous phylogeny for this group is consistent with a complex evolutionary history of the Rerio clade involving a hybrid origin of D. rerio, high rates of incomplete lineage sorting, and speciation with gene flow.

“The difficulty finding a single, unambiguous phylogeny for this group is consistent with a complex evolutionary history of the Rerio clade involving a hybrid origin of D. rerio, high rates of incomplete lineage sorting, and speciation with gene flow as suggested previously2,31.”

• cites [2] sourcePartly supported

  The Conclusions and Multilocus Phylogeny sections describe multiple gene-tree topologies and a complex history involving rapid speciation and introgression. The introgression tests report gene flow involving the D. rerio lineage after it diverged from D. aesculapii, but do not establish a hybrid origin for D. rerio. The Multilocus Phylogeny section says it remains unclear how often incomplete lineage sorting explains the differing topologies, so it does not support the claim of high rates. (LLM, conf=0.88)

  Experimental — accuracy not yet measured
• cites [31] sourcePartly supported

  Split into 2 atoms — 0 confirmed, 2 partially supported, 0 disputed, 0 misattributed, 0 unverifiable. Expand the claim to see per-atom verdicts.

  Experimental — accuracy not yet measured

May combine two claims.

Gap junctions containing Gja5b are critical for stripe patterning in D. rerio, with mutants developing spots similar to the wild-type pattern of D. tinwini.

“Gap junctions containing Gja5b are critical for stripe patterning in D. rerio with mutants developing spots similar to the wild-type pattern of D. tinwini15,16,54.”

• cites [15] sourcePartly supported

  The Introduction (page unavailable) identifies connexin 41.8/GJA5 as a gap-junction subunit and says leopard mutants in Danio rerio have stripes broken into dark spots. The Abstract (page unavailable) likewise reports spotted patterns from leo mutations. Our copy does not identify the spots as similar to the wild-type pattern of D. tinwini, so that comparison is unsupported. (LLM, conf=0.86)

  Experimental — accuracy not yet measured
• cites [16] sourcePartly supported

  The Results and the ‘Connexins in Zebrafish Skin Patterning’ section report that Cx41.8 mutations change zebrafish stripes into spots, supporting the patterning claim. Our copy does not compare those spots with wild-type D. tinwini or explicitly identify Cx41.8 as Gja5b. (LLM, conf=0.84)

  Experimental — accuracy not yet measured
• cites [54] sourcePartly supported

  The Introduction reports that Cx41.8 is involved in zebrafish stripe formation and that leopard mutants with cx41.8 mutations have spotted skin. The paper calls the connexin Cx41.8 rather than Gja5b, and the passages found do not compare the mutants’ spots with wild-type D. tinwini; that comparison is unsupported by this source. (LLM, conf=0.85)

  Experimental — accuracy not yet measured

Probably not a substantial claim: may not pass S3 (load-bearing).

Recessive loss-of-function alleles in heterozygous F1 fish (c.16_25delCTGCTGGGGA, p.Lys6ThrfsX2) were selected to establish the homozygous mutant line gja5bt21mp developing the characteristic mutant phenotype similar to gja5bt1.

“Recessive loss-of-function alleles in heterozygous F1 fish (c.16_25delCTGCTGGGGA, p.Lys6ThrfsX2) were selected to establish the homozygous mutant line gja5bt21mp developing the characteristic mutant phenotype similar to gja5bt1 109,110.”

• cites [109] sourceNot checked — no full text of the source
• cites [110] sourcePartly supported

  The 'Cx41.8 is the leopard gene' section identifies leo t1 as a recessive null allele, and 'Zebrafish leopard phenotype' reports a spotted pattern in its homozygotes. The paper does not report the specified deletion, the gja5bt21mp line, or a comparison of that line with gja5bt1; it also does not use the gja5bt1 designation. (LLM, conf=0.86)

  Experimental — accuracy not yet measured

May combine two claims.

Of the remaining combinations, six (16.7 % of total) showed variant stripe pattern, but with enough variability that phenotypes could also overlap with those of control hybrids, suggesting genetic polymorphisms in modifier loci within species.

“Of the remaining combinations, six (16.7 % of total) showed variant stripe pattern (circles with dashed white outlines in Fig. 5A), but with enough variability that phenotypes could also overlap with those of control hybrids, suggesting genetic polymorphisms in modifier loci within species48.”

• cites [48] sourcePartly supported

  Split into 2 atoms — 1 confirmed, 0 partially supported, 0 disputed, 0 misattributed, 1 unverifiable. Expand the claim to see per-atom verdicts.

  Experimental — accuracy not yet measured

May combine two claims.

The Danio species can be grouped into two clades: six species belong to the Rerio clade and four belong to the Choprae clade.

“The Danio species can be grouped into two clades: six species belong to the Rerio clade (D. rerio, D. aesculapii, D. kyathit, D. tinwini, D. nigrofasciatus, D. albolineatus), and four belong to the Choprae clade (D. jaintianensis, D. choprae, D. margaritatus, D. erythromicron) (Fig. 1).”

May combine two claims.

The assembled genomes have lengths ranging from 1.4 - 1.7 Gbp for species in the Rerio clade, 1.0 - 1.1 Gbp for those in the Choprae clade, and around 0.7 Gbp for the two Danionella species, all with relatively high continuity and completeness.

“The assembled genomes have lengths ranging from 1.4 – 1.7 Gbp for species in the Rerio clade, 1.0 – 1.1 Gbp for those in the Choprae clade, and around 0.7 Gbp for the two Danionella species, all with relatively high continuity and completeness.”

May combine two claims.

Leveraging newly available genome sequences and transcriptomes of other Danio species identified signs of selection in these genes as well as changes to inferred protein sequences.

“Leveraging the newly available genome sequences and transcriptomes of other Danio species identified signs of selection in these genes as well as changes to inferred protein sequences.”

May combine two claims.

These regions experienced the highest number of non-synonymous substitutions, with the greatest extent of change observed within the Choprae clade.

“These regions experienced the highest number of non-synonymous substitutions, with the greatest extent of change observed within the Choprae clade.”

May combine two claims.

Globally, 37.1 % of regional trees supported D. kyathit as sister to D. rerio, while 18.6 % of them supported D. aesculapii as sister to D. rerio.

“Globally, 37.1 % of regional trees supported D. kyathit as sister to D. rerio, while 18.6 % of them supported D. aesculapii as sister to D. rerio.”

May combine two claims.

The SNP profiles of Danio rerio strains exhibited similar genomic distributions, whereas those of D. aesculapii and D. kyathit showed pronounced divergence in both SNP density and localization.

“In most cases, the SNP profiles of Danio rerio strains exhibited similar genomic distributions, whereas those of D. aesculapii and D. kyathit showed pronounced divergence in both SNP density and localization.”

May combine two claims.

The quartet probability of reconstructing the three species as sister species were very close to one third (random), with less than 5 % difference in relative frequency.

“The quartet probability of reconstructing the three species as sister species were very close to one third (random), with less than 5 % difference in relative frequency (Fig. 3B).”

May combine two claims.

The Ne of D. rerio reached its peak in the Late Pleistocene interglacial stage (about 126-115 thousand years ago (kya)) and began to decrease sharply from the last glacial period (about 115 - 11.7 kya).

“The Ne of D. rerio reached its peak in the Late Pleistocene interglacial stage (about 126–115 thousand years ago (kya)) and began to decrease sharply from the last glacial period (about 115 – 11.7 kya).”

Probably not a substantial claim: may not pass S3 (load-bearing).

There was a positive correlation between gene number and genome size.

“There was a positive correlation between gene number and genome size.”

May combine two claims.

The dominant topologies transitioned from centromere to telomere: near the centromere region, the dominant tree was (((D. rerio, D. kyathit), D. aesculapii), (D. tinwini, D. nigrofasciatus)); for the intermediate region: (((D. rerio, D. aesculapii), D. kyathit), (D. tinwini, D. nigrofasciatus)); and for the telomere region: ((D. rerio, D. aesculapii), (D. kyathit, (D. tinwini, D. nigrofasciatus))

“The dominant topologies transitioned from centromere to telomere: 1) near the centromere region, the dominant tree was (((D. rerio, D. kyathit), D. aesculapii), (D. tinwini, D. nigrofasciatus)); 2) for the intermediate region: (((D. rerio, D. aesculapii), D. kyathit), (D. tinwini, D. nigrofasciatus)); and 3) for the telomere region: ((D. rerio, D. aesculapii), (D. kyathit, (D. tinwini, D. nigrofas”

May combine two claims.

The assemblies of D. tinwini, D. jaintianensis, D. albolineatus, D. choprae, and D. cerebrum show limited contiguity (Scaffold N50 = 4,632 - 18,166 kb) but achieve reasonably complete gene space coverage, with BUSCO scores ranging from 78 to 91.9 %.

“The assemblies of D. tinwini, D. jaintianensis, D. albolineatus, D. choprae, and D. cerebrum were generated using 10X Genomics data. Although these assemblies show limited contiguity (Scaffold N50 = 4,632 – 18,166 kb), they achieve reasonably complete gene space coverage, with BUSCO scores ranging from 78 to 91.9 %.”

May combine two claims.

The dN/dS values of pigment-related PSGs were notably higher than those of other PSGs at node 1 and also the ancestor node of Choprae clade (node 2) but were lower at the ancestor node of Rerio clade (node 6).

“The dN/dS values of pigment-related PSGs were notably higher than those of other PSGs at node 1 and also the ancestor node of Choprae clade (node 2) but were lower at the ancestor node of Rerio clade (node 6) (Fig. 4C).”

Probably not a substantial claim: may not pass S3 (load-bearing), S4 (based).

Although assemblies relying on short reads or linked reads exhibit lower contiguity, they still achieve satisfactory gene-space completeness.

“Although assemblies relying on short reads or linked reads exhibit lower contiguity, they still achieve satisfactory gene-space completeness.”

May combine two claims.

The transmembrane helix domains of Gja5b are highly conserved, whereas the intracellular loop and the C-terminus, both of which are disordered regions, were more variable among species.

“We found that the transmembrane helix domains are highly conserved, whereas the intracellular loop and the C-terminus, both of which are disordered regions, were more variable among species (Fig. 6A and B).”

May combine two claims.

D. albolineatus experienced population expansion during the Upper Pleistocene (126 - 11.7 kya), during which the population was relatively stable in size from 50 - 20 kya.

“The third clade (D. albolineatus): D. albolineatus experienced population expansion during the Upper Pleistocene (126 – 11.7 kya), during which the population was relatively stable in size from 50 – 20 kya.”

May combine two claims.

In most cases (25 out of 36, 69.4 %), we observed no differences between wild-type and hemizygous test hybrids, suggesting that functions of these genes are largely conserved between D. rerio and these other species.

“In most cases (25 out of 36, 69.4 %), we observed no differences between wild-type and hemizygous test hybrids, suggesting that functions of these genes are largely conserved between D. rerio and these other species.”

May combine two claims.

The highest average dN/dS value for pigment-related PSGs was observed in D. margaritatus, with ectodysplasin A receptor (edar, dN/dS = 19.9) and engrailed homeobox 1b (en1b, dN/dS = 4.0) showing relatively high dN/dS values.

“The highest average dN/dS value for pigment-related PSGs was observed in D. margaritatus, with ectodysplasin A receptor (edar, dN/dS = 19.9) and engrailed homeobox 1b (en1b, dN/dS = 4.0) showing relatively high dN/dS values.”

May combine two claims.

Structural modeling based on protein sequences suggests that the C-terminal domain of Gja5b may adopt a similar folded conformation in D. margaritatus and D. tinwini, whereas it appears structurally conserved across other species.

“Structural modeling based on protein sequences suggests that this domain may adopt a similar folded conformation in both species, whereas it appears to be structurally conserved across other species (Fig. 6D and E).”

May combine two claims.

The p-value of kcnj13 and gja4 were all less than 0.05 at nodes 2/4/5/6, but kcnj13 underwent positive selection at nodes 4/5/6, while gja4 at nodes 2/4/5.

“The p-value of kcnj13 and gja4 were all less than 0.05 at nodes 2/4/5/6, but kcnj13 underwent positive selection at nodes 4/5/6, while gja4 at nodes 2/4/5 (Table S8).”

May combine two claims.

Of 253 pigment-related SGs, the average dN/dS value was lowest in the horizontally striped species while spotted species had higher values.

“Of 253 pigment-related SGs, the average dN/dS value was lowest in the horizontally striped species while spotted species had higher values (Fig. 4D).”

Probably not a substantial claim: may not pass S3 (load-bearing).

Reads from mitochondrial genes mapped almost exclusively to the maternal D. rerio genome.

“reads from mitochondrial genes mapped almost exclusively to the maternal D. rerio genome as anticipated (Fig. S26C)”

May combine two claims.

D. margaritatus and D. tinwini each carry one unique amino acid residue in the C-terminus of Gja5b.

“D. margaritatus and D. tinwini each carry one unique amino acid residue in the C-terminus.”

Probably not a substantial claim: may not pass S3 (load-bearing).

Across all libraries, many genes displayed a slight bias for the D. rerio allele.

“Across all libraries, many genes displayed a slight bias for the D. rerio allele”

May combine two claims.

Complementation tests suggest functional divergence in igsf11 between D. rerio and D. margaritatus, but less clearly between D. rerio and D. erythromicron or D. choprae.

“In particular, these complementation tests suggest functional divergence in igsf11 between D. rerio and D. margaritatus, but less clearly between D. rerio and D. erythromicron or D. choprae, the closest relatives of D. margaritatus.”

May combine two claims.

The D. rerio strains' assemblies resulted in high scaffold continuity (Scaffold N50 = 49,100 - 51,269 kb) but intermediate BUSCO completeness scores (72.9 - 79.4 %), indicating potential residual errors or gaps despite successful chromosomal assignment.

“This process resulted in high scaffold continuity (Scaffold N50 = 49,100 – 51,269 kb) but intermediate BUSCO completeness scores (72.9 – 79.4 %), indicating potential residual errors or gaps despite successful chromosomal assignment.”

May combine two claims.

Using the Danio rerio (GRCz11) reference genome, the CB strain exhibited 37 % more SNPs compared to the AB strain, and 30 % more SNPs compared to the NA strain.

“Using the Danio rerio (GRCz11) reference genome, the CB strain exhibited 37 % more SNPs compared to the AB strain, and 30 % more SNPs compared to the NA strain”

May combine two claims.

D. aesculapii and D. kyathit populations exhibited expansion-contraction dynamics, but their contraction phases began later than D. rerio's, at 40 - 30 kya and 20 kya, respectively.

“Similarly, D. aesculapii and D. kyathit populations exhibited expansion-contraction dynamics, but their contraction phases began later than D. rerio’s, at 40 – 30 kya and 20 kya, respectively.”

May combine two claims.

The right arm of Dre4 exhibited extensive repetitive matches to the other two genome sequences, yet harbored few orthologs or paralogs on any chromosome of D. aesculapii or D. kyathit.

“The right arm of Dre4 exhibited extensive repetitive matches to the other two genome sequences, yet harbored few orthologs or paralogs on any chromosome of D. aesculapii or D. kyathit.”

May combine two claims.

The D. rerio PAVs were generally shorter in length, while the number of PAVs was almost equal to or greater than that of the other Danio species.

“More specifically, the D. rerio PAVs were generally shorter in length, while the number of PAVs was almost equal to or greater than that of the other Danio species (Table S1).”

May combine two claims.

The genome of D. dracula exhibited excellent contiguity, with a Contig N50 of 2,300 kb and a Scaffold N50 of 10,288 kb, along with high BUSCO completeness (90.3 %).

“The genome of D. dracula, also assembled from PacBio long reads, exhibited excellent contiguity, with a Contig N50 of 2,300 kb and a Scaffold N50 of 10,288 kb, along with high BUSCO completeness (90.3 %).”

May combine two claims.

These trees revealed strong positional effects and that topology frequencies were significantly different at different locations along chromosomes.

“These trees revealed strong positional effects and that topology frequencies were significantly different at different locations along chromosomes (Fig. 3C).”

May combine two claims.

D. erythromicron population expanded during 200-50 kya before contracting at 50 - 40 kya, whereas D. margaritatus population expanded later (100 - 20 kya) with contraction initiating at 20 kya.

“The fourth clade (D. erythromicron, D. margaritatus): D. erythromicron population expanded during 200–50 kya (spanning the Late Pleistocene interglacial to early Last Glacial periods) before contracting at 50 – 40 kya, whereas D. margaritatus population expanded later (100 – 20 kya) with contraction initiating at 20 kya.”

May combine two claims.

Bidirectional gene flow signals were detected between (D. tinwini, D. nigrofasciatus) and D. rerio, also between D. kyathit and D. aesculapii.

“Bidirectional gene flow signals were detected between (D. tinwini, D. nigrofasciatus) and D. rerio, also between D. kyathit and D. aesculapii (Fig. 3D, Fig. S12).”

May combine two claims.

In the WGA, first codon nucleotide, and first+second codon nucleotide trees, D. kyathit was inferred as the sister clade to (D. rerio, D. aesculapii), but distinguished from (D. tinwini, D. nigrofasciatus).

“In the WGA, first codon nucleotide, and first+second codon nucleotide trees, D. kyathit was inferred as the sister clade to (D. rerio, D. aesculapii), but distinguished from (D. tinwini, D. nigrofasciatus).”

May combine two claims.

Expression patterns of asip2b and mc5r were suggestive of xanthophores, and mdkb of iridophores.

“Expression patterns of other genes in this set were suggestive of xanthophores (asip2b, mc5r) and iridophores (mdkb).”

May combine two claims.

Populations of D. tinwini and D. nigrofasciatus began to expand 20 kya, slightly in D. tinwini, while sharply in D. nigrofasciatus.

“Populations of these two species began to expand 20 kya, slightly in D. tinwini, while sharply in D. nigrofasciatus.”

May combine two claims.

The tree topologies ((D. rerio, D. aesculapii), (D. kyathit, (D. tinwini, D. nigrofasciatus)) and ((D. rerio, D. aesculapii), ((D. kyathit, D. tinwini), D. nigrofasciatus)) are both consistent with the existing distribution of the Danio species and are mainly recovered in regions of high recombination near the ends of Danio chromosomes.

“The tree topologies ((D. rerio, D. aesculapii), (D. kyathit, (D. tinwini, D. nigrofasciatus)) and ((D. rerio, D. aesculapii), ((D. kyathit, D. tinwini), D. nigrofasciatus)) are both consistent with the existing distribution of the Danio species and are mainly recovered in regions of high recombination near the ends of Danio chromosomes (Fig. 3C and E).”

May combine two claims.

Assemblies for D. nigrofasciatus, D. margaritatus, and D. erythromicron are structurally fragmented (Scaffold N50 = 152 - 659 kb; scaffold count >1.38 million), yet exhibit moderate BUSCO completeness (79.7 - 84.9 %).

“Assemblies for D. nigrofasciatus, D. margaritatus, and D. erythromicron were constructed using Platanus, primarily based on short-read data. These assemblies are structurally fragmented (Scaffold N50 = 152 – 659 kb; scaffold count >1.38 million), yet exhibit moderate BUSCO completeness (79.7 – 84.9 %), indicating adequate coverage of conserved gene regions.”

May combine two claims.

D. rerio inherited ancestral horizontal stripes, whereas relatives repeatedly evolved spots and vertical bars.

“Whereas D. rerio inherited ancestral horizontal stripes, relatives repeatedly evolved spots and vertical bars.”

• AssumptionRests on an assumption: Pigment patterns play key roles in camouflage, kin recognition, and mate choice, and are targets for natural and sexual selection.

May combine two claims.

The uplift time of the Indo-Burman Ranges overlaps with the divergence time of Rerio clade species, and the divergence of D. choprae and D. jaintianensis.

“This uplift time overlaps with the divergence time of Rerio clade species, and the divergence of D. choprae and D. jaintianensis (Fig. S13-S15).”

• cites [41] sourceNot checked — no full text of the source

Classified as probably not a substantial claim.

To avoid any issues with differing annotation quality, we used a two-step mapping strategy.

“To avoid any issues with differing annotation quality, we used a two-step mapping strategy.”

Classified as probably not a substantial claim.

To minimize the effect of annotation, sequencing, and assembly errors, pseudogenes, non-orthologous alignments, and non-conserved gene structures on subsequent evolutionary rate analyses, a series of rigorous filtering criteria was adopted.

“To minimize the effect of annotation, sequencing, and assembly errors, pseudogenes, non-orthologous alignments, and non-conserved gene structures on subsequent evolutionary rate analyses, a series of rigorous filtering criteria was adopted:”

Classified as probably not a substantial claim.

Phenotypic variability observed between hybrids of different genotypes could be caused by the influence of novel genetic backgrounds from specific species pairs rather than by a functional divergence of the tested genes.

“Phenotypic variability observed between hybrids of different genotypes could be caused by the influence of novel genetic backgrounds from specific species pairs rather than by a functional divergence of the tested genes.”

May combine two claims.

Pigment-related genes, including those regulating direct cell-cell interactions between pigment cells, exhibited signs of positive selection, and these genes evolved more rapidly in some species than others.

“These results indicated that pigment-related genes, i.e., including those regulating direct cell-cell interactions between pigment cells, exhibited signs of positive selection, and these genes evolved more rapidly in some species than others.”

May combine two claims.

A cis-acting regulatory change has increased the expression of aox6 in D. rerio and affected xanthophore pigment intensity.

“This finding suggests that a cis-acting regulatory change has increased the expression of aox6 in D. rerio and affected xanthophore pigment intensity.”

May combine two claims.

Amino acid changes in the intracellular loop of Gja5b may have led to reduced function and discontinuous stripes, with divergence in kcnj13 of D. choprae contributing to reduced contrast between dark bars and light intervening regions, similar to D. aesculapii.

“Here, amino acid changes in the intracellular loop of Gja5b may have led to reduced function and discontinuous stripes, with divergence in kcnj13 of D. choprae contributing to reduced contrast between dark bars and light intervening regions, similar to D. aesculapii.”

• AssumptionRests on an assumption: Pigment pattern formation is a self-organizing process based on interactions between pigment cell types.

May combine two claims.

D. rerio, D. kyathit, and D. nigrofasciatus retained ancestral striped patterns, whereas structural variants in the C-terminus of Gja5b impaired normal intercellular interactions between melanophores and xanthophores, disrupting horizontal stripes and leading to spots in D. tinwini.

“Our findings suggest a model in which D. rerio, D. kyathit, and D. nigrofasciatus retained ancestral striped patterns, whereas structural variants in the C-terminus of Gja5b impaired normal intercellular interactions between melanophores and xanthophores, disrupting horizontal stripes and leading to spots in D. tinwini.”

• AssumptionRests on an assumption: Pigment pattern formation is a self-organizing process based on interactions between pigment cell types.

May combine two claims.

The reporter genes were upregulated in the expected cell types, verifying enrichment of hybrid cells.

“The reporter genes were upregulated in the expected cell types, verifying enrichment of hybrid cells (Fig. S26A).”

May combine two claims.

Gja5b may have functionally converged in D. margaritatus and D. tinwini, which are situated in different clades within Danio, and this convergence might have mediated changes in Gja5b proteins resulting in the convergent evolution of their spotted patterns.

“Our results also suggest that gja5b may have functionally converged in D. margaritatus and D. tinwini, which are situated in different clades within Danio, and that this convergence might have mediated changes in Gja5b proteins resulting in the convergent evolution of their spotted patterns.”

May combine two claims.

The detected gene flow signals between D. rerio and D. aesculapii, and among D. kyathit, D. tinwini and D. nigrofasciatus are consistent with the current geographic distribution of these Danio species.

“The detected gene flow signals between D. rerio and D. aesculapii, and among D. kyathit, D. tinwini and D. nigrofasciatus are consistent with the current geographic distribution of these Danio species (Fig. 3E).”

Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

Fish with different pigment patterns vary considerably in contrast, thus requiring different settings for aperture and exposure time, which can result in slightly different color representations in the pictures.

“Fish with different pigment patterns vary considerably in contrast, thus requiring different settings for aperture and exposure time, which can result in slightly different color representations in the pictures.”

Classified as probably not a substantial claim.

RNA-seq in hybrids requires discerning species of origin as well as gene of origin for each read.

“RNA-seq in hybrids requires discerning species of origin as well as gene of origin for each read.”

May combine two claims.

Wild-type alleles from these other Danio species cannot fully complement the loss of function, and otherwise recessive, D. rerio alleles, and are thus likely to be hypomorphic relative to wild-type D. rerio alleles.

“These results suggest that wild-type alleles from these other Danio species cannot fully complement the loss of function, and otherwise recessive, D. rerio alleles, and are thus likely to be hypomorphic relative to wild-type D. rerio alleles.”

May combine two claims.

Geographical isolation may have accelerated the divergence of D. choprae and D. jaintianensis, after which gene flow occurred among species in nearby basins, respectively.

“Geographical isolation may have accelerated the divergence of D. choprae and D. jaintianensis, after which gene flow occurred among species in nearby basins, respectively (Table S6).”

Classified as probably not a substantial claim.

This complex speciation history resulted in the genetic material of D. rerio and D. kyathit, originally more closely related, being now more closely related to their sympatric Danio species.

“This complex speciation history resulted in the genetic material of D. rerio and D. kyathit, originally more closely related, being now more closely related to their sympatric Danio species.”

Classified as probably not a substantial claim.

To mitigate the possibility of mapping bias toward the D. rerio reference in regions with poor sequence conservation, we excluded reads mapping to 3' and 5' untranslated regions.

“To mitigate the possibility of mapping bias toward the D. rerio reference in regions with poor sequence conservation, we excluded reads mapping to 3’ and 5’ untranslated regions.”

May combine two claims.

The phylogenomic analyses point to a complex history of incomplete lineage sorting, introgression, and possible genomic rearrangements contributing to topological discordances between Danio species.

“The phylogenomic analyses point to a complex history of incomplete lineage sorting, introgression, and possible genomic rearrangements contributing to topological discordances between Danio species.”

Classified as probably not a substantial claim.

Differential expression was assessed using default parameters: Wald test with adjusted p-values using Benjamini-Hochberg correction and FDR < 0.05.

“Differential expression was assessed using default parameters: Wald test with adjusted p-values using Benjamini-Hochberg correction and FDR < 0.05.”

May combine two claims.

There was introgression of alleles from D. tinwini into D. kyathit, and from D. kyathit into D. nigrofasciatus.

“This evidence suggested introgression of alleles from D. tinwini into D. kyathit, and from D. kyathit into D. nigrofasciatus (Table S5, Fig. 3D, Fig. S12).”

Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

The remaining 368 orthogroups likely contain multiple gene copies or may even be completely absent in one or more species, potentially resulting from lineage-specific duplications or losses such as those mediated by the teleost-specific whole-genome duplication (TGD).

“Beyond the 253 pigment-related SGs analyzed above, the remaining 368 orthogroups likely contain multiple gene copies or may even be completely absent in one or more species, potentially resulting from lineage-specific duplications or losses such as those mediated by the teleost-specific whole-genome duplication (TGD).”

Probably not a substantial claim: may not pass S3 (load-bearing).

The very low minimum identity accounts for sequence divergence between species when remapping to GRCz11.

“The very low minimum identity accounts for sequence divergence between species when remapping to GRCz11.”

May combine two claims.

The observed non-complementation phenotypes of hemizygous hybrids were not as severe as those of homozygous D. rerio mutants, implying that these loci have some activity in pigment cells of the other species and that functions of these genes in pigment pattern predate the origin of the Danio genus.

“Our finding that the observed non-complementation phenotypes of hemizygous hybrids were not as severe as those of homozygous D. rerio mutants (Fig. 5A, bottom) further implies that these loci have some activity in pigment cells of the other species and that functions of these genes in pigment pattern predate the origin of the Danio genus.”

Probably not a substantial claim: may not pass S3 (load-bearing), S4 (based).

The integration of long-read technologies with multi-platform scaffolding is essential for producing chromosome-level assemblies.

“In summary, the integration of long-read technologies with multi-platform scaffolding is essential for producing chromosome-level assemblies.”