The research, led by Chen Chen and Geng Zhou of the Hunan Vegetable Research Institute at the Hunan Academy of Agricultural Sciences, together with colleagues at Yuelushan Laboratory and Hunan Agricultural University, focused on the gene CsMYB60. This gene encodes an R2R3-MYB transcription factor, a type of DNA-binding protein that acts as a master switch for plant pigment production. Earlier work had already established that CsMYB60 drives the biosynthesis of flavonoids and proanthocyanidins, two classes of plant secondary metabolites, in cucumber fruit spines. More recent studies extended that role to the mature fruit skin, showing that the gene helps control flavonoid accumulation and, with it, the color the peel takes on as the fruit ripens.

To dissect the genetic basis of peel color with precision, the team turned to a powerful experimental resource: a pair of near-isogenic lines. These are two cucumber inbred lines, L14 and L19, that are genetically almost identical except at the region governing peel color. L14 develops a yellow mature peel, while L19, a mutant derived from it, remains white. Because the two lines share nearly all of their genetic background, any difference in peel color can be attributed with high confidence to the small stretch of DNA where they differ, eliminating the noise that usually complicates genetic analysis in diverse varieties.

The researchers combined bulked segregant analysis with high-resolution fine mapping. In bulked segregant analysis, scientists cross the two contrasting lines, then pool the DNA of many offspring that share the same trait, in this case yellow peel, and sequence those pools alongside pools from offspring with the opposite trait. Regions of the genome where the pools diverge sharply are candidate homes for the causal gene. This approach, originally developed for rapid mapping in rice, allowed the team to narrow the peel color locus to a 121-kilobase interval on chromosome 4. Within that interval sat exactly one compelling candidate: CsMYB60.

What the team found when they compared the DNA sequences of the two alleles was the study’s central surprise. The yellow-peel line L14 carries not one but two Mutator-like element insertions, known as MULE1 and MULE2, within the second intron of CsMYB60. The white-peel mutant L19, by contrast, retains only MULE1. Mutator-like elements are a family of transposable elements, mobile DNA sequences that can copy or cut themselves out of the genome and reinsert elsewhere. Far from being genomic junk, transposable elements are now recognized as major drivers of plant evolution, frequently rewiring gene regulation by inserting into promoters, introns, or nearby regulatory regions and altering when, where, and how strongly a gene is expressed.

Here, the double insertion appeared to matter. During the peel color-transition stage, when the fruit’s skin is actively changing, the MULE1 plus MULE2 haplotype showed higher CsMYB60 transcript abundance than the MULE1-only haplotype in the near-isogenic materials examined. In other words, the line carrying both insertions produced more messenger RNA from the pigment-regulating gene, and that line was the one that turned yellow. The authors are careful to note that the causal regulatory contribution of MULE2 itself remains to be determined; the association between the structural haplotype and both expression level and peel color is clear, but precisely how the second insertion boosts transcription is a question for future work.

The team did not stop at correlation. Using CRISPR/Cas9 genome editing, they knocked out CsMYB60 directly in the L14 background. The result was unambiguous: the edited plants lost their yellow pigmentation entirely, confirming that CsMYB60 is required for yellow mature peel formation in this genetic background. This kind of functional validation is the gold standard in plant genetics, because it demonstrates necessity rather than mere association. Delete the gene, and the trait disappears; restore or retain it with the right structural haplotype, and the trait appears.

The findings fit into a broader and rapidly evolving picture of how fruit color is controlled in cucurbits. Cucumber peel and flesh color are governed by a constellation of loci affecting chloroplast development, carotenoid accumulation, and flavonoid biosynthesis. Previous studies have implicated genes such as CsMYB36 in yellow-green immature peel, APRR2 in green immature fruit color, and various chloroplast-related genes in light green phenotypes. The R2R3-MYB family, meanwhile, is famous across the plant kingdom for its role in the MYB-bHLH-WDR transcriptional complexes that regulate flavonoid pathways. What makes the new study distinctive is that it identifies a structural haplotype, a specific physical arrangement of DNA, rather than a simple single-nucleotide change, as the molecular signature associated with the trait.

That distinction has practical consequences for breeding. Structural variation like transposon insertions is often invisible to standard SNP-based genotyping, which means breeders can inadvertently lose or overlook valuable haplotypes when crossing varieties. The authors point out that the molecular marker associated with the MULE1 plus MULE2 haplotype may facilitate tracking of this structural variant in relevant cucumber breeding materials. In practice, that means seed companies and public breeding programs could screen seedlings for the double-insertion haplotype and confidently predict which plants will produce the yellow mature peel that certain markets demand, without waiting for the fruit to ripen.

The study also adds to a growing appreciation of transposable elements as engines of agricultural diversity. In crops ranging from maize to grapevine to tomato, transposon insertions have been shown to create or modulate traits as varied as seed coat color, fruit shape, and ripening behavior. The cucumber peel story is a textbook example of how a mobile element can act as a regulatory rheostat: by landing inside an intron, the MULE insertions appear to tune the output of a transcription factor that sits atop an entire pigment biosynthesis pathway. Whether MULE2 acts by altering chromatin structure, changing intron splicing, or providing novel regulatory sequences remains open, and answering it will likely require additional experiments such as reporter assays and splicing analyses.

For now, the work delivers a complete arc from field observation to molecular mechanism. Two nearly identical cucumber lines, one yellow and one white, led researchers to a 121-kilobase window on chromosome 4, then to a single MYB gene, then to two nested transposon insertions, and finally to a CRISPR experiment that sealed the case. The identification of a previously unrecognized MULE1 plus MULE2 structural haplotype associated with distinct CsMYB60 transcript abundance and mature peel coloration gives cucumber geneticists a new landmark on their maps and gives breeders a practical tool. It also serves as a reminder that some of the most important instructions in a crop’s genome are not written in its genes at all, but in the restless, mobile DNA that has been reshaping those genomes for millions of years.

Subject of Research: Transposon-mediated structural haplotypes at the CsMYB60 gene and their association with mature peel coloration in cucumber

Article Title: Transposon-mediated structural haplotypes at CsMYB60 are associated with peel coloration in cucumber

Article References: Chen, C., Chen, J., Chen, H., Liu, X., Lu, X., Tian, Y., & Zhou, G. (2026). Transposon-mediated structural haplotypes at CsMYB60 are associated with peel coloration in cucumber. Theoretical and Applied Genetics, 139(10), Article 285. https://doi.org/10.1007/s00122-026-05388-2

Image Credits: AI Generated

DOI: 10.1007/s00122-026-05388-2

Keywords: cucumber, peel color, CsMYB60, transposable elements, MULE, structural haplotype, flavonoid biosynthesis, R2R3-MYB transcription factor, CRISPR/Cas9, fine mapping, bulked segregant analysis, plant genetics

The software, called BSArray, was developed by M. Habib Widyawan and Zenglu Li of the University of Georgia’s Institute of Plant Breeding, Genetics, and Genomics, together with Qijian Song of the USDA Agricultural Research Service’s Beltsville Agricultural Research Center. It packages a classical genetic technique known as bulked segregant analysis, or BSA, into an interactive R Shiny application, meaning that researchers can run the full analysis through a graphical interface rather than a command line. The developers designed it specifically for data generated on Illumina Infinium SNP genotyping arrays, the fixed-content genotyping platforms that many crop genetics labs already use routinely.

Bulked segregant analysis is an elegant shortcut in quantitative genetics. Instead of genotyping every individual in a breeding population and running a full linkage analysis, researchers create two bulks: one pooling the offspring that show the trait of interest most strongly, and one pooling those that show it least. Any genomic region that truly influences the trait will show a systematic shift in allele frequency between the two bulks, because the individuals in the high-trait bulk share, by selection, the parental chromosome segment carrying the favorable allele. Scanning the genome for such frequency differences points directly to the quantitative trait loci, or QTL, responsible. The approach trades some statistical power for enormous savings in cost and time, which is why it remains popular for early-stage gene discovery in crops.

Where sequencing-based BSA has attracted a rich ecosystem of automated pipelines, the array-based version has lagged behind. Illumina Infinium arrays produce fluorescence intensity measurements for each SNP, and converting those signals into reliable allele frequency estimates for bulks involves several judgment calls: how to normalize intensities, how to weight genotypes by their confidence, and how to filter out markers with poor performance. In practice, researchers have had to compare parental and bulk genotypes by hand, write custom filtering scripts, and assemble results in spreadsheets. Each manual step is an opportunity for error, and none of them scale well when a breeding program wants to analyze dozens of populations across multiple traits and seasons.

BSArray addresses this by implementing the full pipeline end to end. The application starts with the normalized intensity metrics from the array and estimates allele frequencies for each bulk directly from those signals. Its central innovation is a quality-weighted statistic, denoted W, which computes the allele frequency difference between bulks while integrating two pieces of information: the magnitude of the frequency difference itself and the confidence of the underlying genotype calls. Markers that show a large difference but rest on shaky genotype calls are discounted relative to markers where the signal is both large and reliable. This weighting scheme helps suppress false positives arising from noisy or poorly performing array probes, a persistent nuisance in array-based genotyping.

From the per-marker W statistic, the pipeline moves to a sliding window framework that scans the genome in overlapping segments. Rather than relying on a single noisy marker, the window approach asks whether a stretch of consecutive markers collectively shows an elevated frequency difference, which is the spatial signature of a genuine QTL. Significance is assessed against empirical genome-wide thresholds, so the user gets a statistically grounded call on which regions are putative QTL rather than an arbitrary visual cutoff. Once significant regions are identified, the pipeline automatically prioritizes the most informative candidate markers within them and annotates the nearby gene models, giving breeders a direct bridge from a statistical signal to the genes they might pursue for marker-assisted selection.

The developers demonstrated the tool using the SoySNP50K array, a 50,000-SNP platform developed for soybean, in three implementation scenarios that map neatly onto the realistic situations breeders encounter. In the first, BSArray detected a previously reported major-effect locus, confirming that the pipeline recovers known biology. In the second, it mapped a locus that had not previously been placed on the soybean genome, showing its value for genuinely new discovery. In the third, it detected a novel minor-effect locus within a QTL region that was already known, illustrating the tool’s ability to refine and dissect complex trait architecture rather than only finding the largest, most obvious signals. Together, the three cases cover the spectrum from validation to refinement to novel discovery.

The practical significance for crop improvement lies in speed and reproducibility. Marker-assisted selection, the practice of using DNA markers to guide breeding decisions, depends on having well-mapped QTL and reliable flanking markers. Every step that slows QTL discovery—manual genotype comparison, ad hoc scripting, irreproducible filtering—delays the delivery of markers that breeders can deploy in the field. By reducing manual intervention and computational requirements, BSArray offers what its developers describe as a rapid, reproducible, and user-friendly framework that serves as a practical and scalable alternative to sequencing-based approaches, particularly for laboratories whose primary genotyping investment is already in SNP arrays rather than whole-genome sequencing.

The cost logic is an important part of the appeal. Whole-genome resequencing of bulks delivers dense marker coverage but carries substantial per-sample expense and generates data volumes that demand serious computational infrastructure. SNP arrays, by contrast, offer a fixed, curated marker set at a fraction of the cost, and for many crop species the array density is more than sufficient for bulked segregant analysis, which needs enough markers to detect regional allele frequency shifts rather than the causal variant itself. A tool that makes array-based BSA as turnkey as its sequencing-based counterparts could therefore broaden access to QTL discovery for smaller breeding programs, national agricultural research systems, and university labs in developing countries—settings where the soybean work’s funding context, including the United Soybean Board and support for an Indonesian doctoral researcher through the Fulbright FIRST program, reflects the international reach of modern crop genomics.

BSArray is published as open-access software in BMC Bioinformatics, with earlier versions tested by collaborators including Francismar C. Marcelino-Guimarães of the Brazilian Agricultural Research Corporation. The design choices—an interactive web-style interface, automated quality weighting, empirical significance thresholds, and built-in gene annotation—reflect a broader trend in bioinformatics toward tools that embed expert judgment into software so that domain scientists can perform analyses that once required a dedicated computational specialist. For a field racing to develop climate-resilient, disease-resistant crop varieties faster than pathogens and weather can evolve, shaving weeks off the gene-discovery pipeline is not a trivial convenience. It is the difference between a breeding cycle that waits on analysis and one that moves at the speed of the growing season.

Subject of Research: Automated bulked segregant analysis pipeline for QTL discovery using SNP genotyping arrays

Article Title: BSArray: An automated bulked segregant analysis pipeline for QTL discovery using Illumina Infinium SNP genotyping array

Article References: Widyawan, M. H., Song, Q., & Li, Z. (2026). BSArray: An automated bulked segregant analysis pipeline for QTL discovery using Illumina Infinium SNP genotyping array. BMC Bioinformatics. https://doi.org/10.1186/s12859-026-06670-2

Image Credits: AI Generated

DOI: 10.1186/s12859-026-06670-2

Keywords: bulked segregant analysis, QTL discovery, Illumina Infinium, SNP array, soybean, SoySNP50K, marker-assisted selection, R Shiny, bioinformatics pipeline, quantitative trait loci, crop breeding, genetic mapping