Evolution of Einkorn wheat centromeres is driven by the mutualistic interplay of two LTR retrotransposons.

BACKGROUND: Centromere function is highly conserved across eukaryotes, but the underlying centromeric DNA sequences vary dramatically between species. Centromeres often contain a high proportion of repetitive DNA, such as tandem repeats and/or transposable elements (TEs). Einkorn wheat centromeres lack tandem repeat arrays and are instead composed mostly of the two long terminal repeat (LTR) retrotransposon families RLG_Cereba and RLG_Quinta which specifically insert in centromeres. However, it is poorly understood how these two TE families relate to each other and if and how they contribute to centromere function and evolution. RESULTS: Based on conservation of diagnostic motifs (LTRs, integrase and primer binding site and polypurine-tract), we propose that RLG_Cereba and RLG_Quinta are a pair of autonomous and non-autonomous partners, in which the autonomous RLG_Cereba contributes all the proteins required for transposition, while the non-autonomous RLG_Quinta contributes GAG protein. Phylogenetic analysis of predicted GAG proteins showed that the RLG_Cereba lineage was present for at least 100 million years in monocotyledon plants. In contrast, RLG_Quinta evolved from RLG_Cereba between 28 and 35 million years ago in the common ancestor of oat and wheat. Interestingly, the integrase of RLG_Cereba is fused to a so-called CR-domain, which is hypothesized to guide the integrase to the functional centromere. Indeed, ChIP-seq data and TE population analysis show only the youngest subfamilies of RLG_Cereba and RLG_Quinta are found in the active centromeres. Importantly, the LTRs of RLG_Quinta and RLG_Cereba are strongly associated with the presence of the centromere-specific CENH3 histone variant. We hypothesize that the LTRs of RLG_Cereba and RLG_Quinta contribute to wheat centromere integrity by phasing and/or placing CENH3 nucleosomes, thus favoring their persistence in the competitive centromere-niche. CONCLUSION: Our data show that RLG_Cereba cross-mobilizes the non-autonomous RLG_Quinta retrotransposons. New copies of both families are specifically integrated into functional centromeres presumably through direct binding of the integrase CR domain to CENH3 histone variants. The LTRs of newly inserted RLG_Cereba and RLG_Quinta elements, in turn, recruit and/or phase new CENH3 deposition. This mutualistic interplay between the two TE families and the plant host dynamically maintains wheat centromeres.

Ruthenium drug BOLD-100 regulates BRAFMT colorectal cancer cell apoptosis through AhR/ROS/ATR signaling axis modulation.

Patients with class I V600EBRAF-mutant (MT) colorectal cancer (CRC) have a poor prognosis and their response to combined anti-BRAF/EGFR inhibition remains limited. There is clearly an unmet need in further understanding the biology of V600EBRAFMT CRC. We have used differential gene expression of BRAFWT and MT CRC cells to identify pathways underpinning BRAFMT CRC. We tested a panel of molecularly/genetically subtyped CRC cells for their sensitivity to the Unfolded Protein Response (UPR) activator BOLD-100. To identify novel combination strategies for BOLD-100, we performed RNA sequencing and high-throughput drug screening. Pathway enrichment analysis identified that the UPR and DNA repair pathways were significantly enriched in BRAFMT CRC. We found that oncogenic BRAF plays a crucial role in mediating response to BOLD-100. Using a systems biology approach, we identified V600EBRAFMT-dependent activation of the replication stress response kinase ATR as a key mediator of resistance to BOLD-100. Further analysis identified acute increases in BRAFMT-dependent-reactive oxygen species (ROS) levels following treatment with BOLD-100 that was demonstrated to promote ATR/CHK1 activation and apoptosis. Furthermore, activation of ROS/ATR/CHK1 following BOLD-100 was found to be mediated through the AHR transcription factor and CYP1A1. Importantly, pharmacological blockade of this resistance pathway with ATR inhibitors synergistically increased BOLD-100-induced apoptosis and growth inhibition in BRAFMT models. These results unveil possible novel therapeutic opportunity for BRAFMT CRC. Implications: BOLD-100 induces BRAFMT-dependent replication stress, and targeted strategies against replication stress (eg. by using ATR inhibitors) in combination with BOLD-100 may serve as a potential novel therapeutic strategy for clinically aggressive BRAFMT CRC.

Plasticity in cell migration modes across development, physiology, and disease.

Cell migration is fundamental to both development and adult physiology, including gastrulation, brain development, angiogenesis, wound healing, bone remodeling, tissue homeostasis, and the immune response. Additionally, misguided cellular migration is implicated in disease pathologies such as cancer metastasis and fibrosis. The microenvironment influences cell migration modes such as mesenchymal, amoeboid, lobopodial, and collective, and these are governed through local signaling by affecting the gene expression and epigenetic alteration of migration-related genes. Plasticity in switching between migration modes is essential for key cellular processes across various contexts. Understanding the mechanisms of cell migration modes and its plasticity is essential for unraveling the complexities of this process and revealing its implications in physiological and pathological contexts. This review focuses on different modes of cell migration, including their aberrant migration in disease pathologies and how they can be therapeutically targeted in disease conditions such as cancer.

Decoding gene regulatory circuitry underlying TNBC chemoresistance reveals biomarkers for therapy response and therapeutic targets.

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype characterised by extensive intratumoral heterogeneity, high rates of metastasis and chemoresistance, leading to poor clinical outcomes. Despite progress, the mechanistic basis of chemotherapy resistance in TNBC patients remains poorly understood. Here, leveraging single-cell transcriptome datasets of matched longitudinal TNBC chemoresponsive and chemoresistant patient cohorts, we unravel distinct cell subpopulations intricately associated with chemoresistance and the signature genes defining these populations. Notably, using genome-wide mapping of the H3K27ac mark, we show that the expression of these chemoresistance genes is driven via a set of TNBC super-enhancers and associated transcription factor networks across TNBC subtypes. Furthermore, genetic screens reveal that a subset of these transcription factors is essential for the survival of TNBC cells, and their loss increases sensitivity to chemotherapeutic agents. Overall, our study has revealed epigenetic and transcription factor networks underlying chemoresistance and suggests novel avenues to stratify and improve the treatment of patients with a high risk of developing resistance.

Basal-epithelial subpopulations underlie and predict chemotherapy resistance in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype, characterized by extensive intratumoral heterogeneity, high metastasis, and chemoresistance, leading to poor clinical outcomes. Despite progress, the mechanistic basis of these aggressive behaviors remains poorly understood. Using single-cell and spatial transcriptome analysis, here we discovered basal epithelial subpopulations located within the stroma that exhibit chemoresistance characteristics. The subpopulations are defined by distinct signature genes that show a frequent gain in copy number and exhibit an activated epithelial-to-mesenchymal transition program. A subset of these genes can accurately predict chemotherapy response and are associated with poor prognosis. Interestingly, among these genes, elevated ITGB1 participates in enhancing intercellular signaling while ACTN1 confers a survival advantage to foster chemoresistance. Furthermore, by subjecting the transcriptional signatures to drug repurposing analysis, we find that chemoresistant tumors may benefit from distinct inhibitors in treatment-naive versus post-NAC patients. These findings shed light on the mechanistic basis of chemoresistance while providing the best-in-class biomarker to predict chemotherapy response and alternate therapeutic avenues for improved management of TNBC patients resistant to chemotherapy.

Ageing impairs the regenerative capacity of regulatory T cells in mouse central nervous system remyelination.

Myelin regeneration (remyelination) is essential to prevent neurodegeneration in demyelinating diseases such as Multiple Sclerosis, however, its efficiency declines with age. Regulatory T cells (Treg) recently emerged as critical players in tissue regeneration, including remyelination. However, the effect of ageing on Treg-mediated regenerative processes is poorly understood. Here, we show that expansion of aged Treg does not rescue age-associated remyelination impairment due to an intrinsically diminished capacity of aged Treg to promote oligodendrocyte differentiation and myelination in male and female mice. This decline in regenerative Treg functions can be rescued by a young environment. We identified Melanoma Cell Adhesion Molecule 1 (MCAM1) and Integrin alpha 2 (ITGA2) as candidates of Treg-mediated oligodendrocyte differentiation that decrease with age. Our findings demonstrate that ageing limits the neuroregenerative capacity of Treg, likely limiting their remyelinating therapeutic potential in aged patients, and describe two mechanisms implicated in Treg-driven remyelination that may be targetable to overcome this limitation.

MicroRNA-335-5p suppresses voltage-gated sodium channel expression and may be a target for seizure control.

There remains an urgent need for new therapies for treatment-resistant epilepsy. Sodium channel blockers are effective for seizure control in common forms of epilepsy, but loss of sodium channel function underlies some genetic forms of epilepsy. Approaches that provide bidirectional control of sodium channel expression are needed. MicroRNAs (miRNA) are small noncoding RNAs which negatively regulate gene expression. Here we show that genome-wide miRNA screening of hippocampal tissue from a rat epilepsy model, mice treated with the antiseizure medicine cannabidiol, and plasma from patients with treatment-resistant epilepsy, converge on a single target-miR-335-5p. Pathway analysis on predicted and validated miR-335-5p targets identified multiple voltage-gated sodium channels (VGSCs). Intracerebroventricular injection of antisense oligonucleotides against miR-335-5p resulted in upregulation of Scn1a, Scn2a, and Scn3a in the mouse brain and an increased action potential rising phase and greater excitability of hippocampal pyramidal neurons in brain slice recordings, consistent with VGSCs as functional targets of miR-335-5p. Blocking miR-335-5p also increased voltage-gated sodium currents and SCN1A, SCN2A, and SCN3A expression in human induced pluripotent stem cell-derived neurons. Inhibition of miR-335-5p increased susceptibility to tonic-clonic seizures in the pentylenetetrazol seizure model, whereas adeno-associated virus 9-mediated overexpression of miR-335-5p reduced seizure severity and improved survival. These studies suggest modulation of miR-335-5p may be a means to regulate VGSCs and affect neuronal excitability and seizures. Changes to miR-335-5p may reflect compensatory mechanisms to control excitability and could provide biomarker or therapeutic strategies for different types of treatment-resistant epilepsy.

A complex epigenome-splicing crosstalk governs epithelial-to-mesenchymal transition in metastasis and brain development.

Epithelial-to-mesenchymal transition (EMT) renders epithelial cells migratory properties. While epigenetic and splicing changes have been implicated in EMT, the mechanisms governing their crosstalk remain poorly understood. Here we discovered that a C2H2 zinc finger protein, ZNF827, is strongly induced during various contexts of EMT, including in brain development and breast cancer metastasis, and is required for the molecular and phenotypic changes underlying EMT in these processes. Mechanistically, ZNF827 mediated these responses by orchestrating a large-scale remodelling of the splicing landscape by recruiting HDAC1 for epigenetic modulation of distinct genomic loci, thereby slowing RNA polymerase II progression and altering the splicing of genes encoding key EMT regulators in cis. Our findings reveal an unprecedented complexity of crosstalk between epigenetic landscape and splicing programme in governing EMT and identify ZNF827 as a master regulator coupling these processes during EMT in brain development and breast cancer metastasis.

Tumour invasion and dissemination.

Activating invasion and metastasis are one of the primary hallmarks of cancer, the latter representing the leading cause of death in cancer patients. Whilst many advances in this area have been made in recent years, the process of cancer dissemination and the underlying mechanisms governing invasion are still poorly understood. Cancer cells exhibit multiple invasion strategies, including switching between modes of invasion and plasticity in response to therapies, surgical interventions and environmental stimuli. The ability of cancer cells to switch migratory modes and their inherent plasticity highlights the critical challenge preventing the successful design of cancer and anti-metastatic therapies. This mini-review presents current knowledge on the critical models of tumour invasion and dissemination. We also discuss the current issues surrounding current treatments and arising therapeutic opportunities. We propose that the establishment of novel approaches to study the key biological mechanisms underlying the metastatic cascade is critical in finding novel targets that could ultimately lead to complete inhibition of cancer cell invasion and dissemination.

Mnt Represses Epithelial Identity To Promote Epithelial-to-Mesenchymal Transition.

The multistep process of epithelial-to-mesenchymal transition (EMT), whereby static epithelial cells become migratory mesenchymal cells, plays a critical role during various developmental contexts, wound healing, and pathological conditions such as cancer metastasis. Despite the established function of basic helix-loop-helix (bHLH) transcription factors (TFs) in cell fate determination, only a few have been examined for their role in EMT. Here, using transcriptome analysis of distinct stages during stepwise progression of transforming growth factor beta (TGFß)-induced EMT in mammary epithelial cells, we revealed distinct categories of bHLH TFs that show differential expression kinetics during EMT. Using a short interfering RNA-mediated functional screen for bHLH TFs during EMT, we found Max network transcription repressor (MNT) to be essential for EMT in mammary epithelial cells. We show that the depletion of MNT blocks TGFß-induced morphological changes during EMT, and this is accompanied by derepression of a large number of epithelial genes. We show that MNT mediates the repression of epithelial identity genes during EMT by recruiting HDAC1 and mediating the loss of H3K27ac and chromatin accessibility. Lastly, we show that MNT is expressed at higher levels in EMT-High breast cancer cells and is required for their migration. Taken together, these findings establish MNT as a critical regulator of cell fate changes during mammary EMT. IMPORTANCE The bHLH TF Mnt promotes epithelial to mesenchymal transition through epigenetic repression of the epithelial gene expression program.

High-resolution mapping of the Mov-1 locus in wheat by combining radiation hybrid (RH) and recombination-based mapping approaches.

KEY MESSAGE: This work reports a quick method that integrates RH mapping and genetic mapping to map the dominant Mov-1 locus to a 1.1-Mb physical interval with a small number of candidate genes. Bread wheat is an important crop for global human population. Identification of genes and alleles controlling agronomic traits is essential toward sustainably increasing crop production. The unique multi-ovary (MOV) trait in wheat holds potential for improving yields and is characterized by the formation of 2-3 grains per spikelet. The genetic basis of the multi-ovary trait is known to be monogenic and dominant in nature. Its precise mapping and functional characterization is critical to utilizing this trait in a feasible manner. Previous mapping efforts of the locus controlling multiple ovary/pistil formation in the hexaploid wheat have failed to produce a consensus for a particular chromosome. We describe a mapping strategy integrating radiation hybrid mapping and high-resolution genetic mapping to locate the chromosomal position of the Mov-1 locus in hexaploid wheat. We used RH mapping approach using a panel of 188 lines to map the Mov-1 locus in the terminal part of long arm of wheat chromosome 2D with a map resolution of 1.67 Mb/cR1500. Then using a genetic population of MOV × Synthetic wheat of F2 lines, we delineated the Mov-1 locus to a 1.1-Mb physical region with a small number of candidate genes. This demonstrates the value of this integrated strategy to mapping dominant genes in wheat.

The Independent Domestication of Timopheev's Wheat: Insights from Haplotype Analysis of the Brittle rachis 1 (BTR1-A) Gene.

Triticum turgidum and T. timopheevii are two tetraploid wheat species sharing T. urartu as a common ancestor, and domesticated accessions from both of these allopolyploids exhibit nonbrittle rachis (i.e., nonshattering spikes). We previously described the loss-of-function mutations in the Brittle Rachis 1 genes BTR1-A and BTR1-B in the A and B subgenomes, respectively, that are responsible for this most visible domestication trait in T. turgidum. Resequencing of a large panel of wild and domesticated T. turgidum accessions subsequently led to the identification of the two progenitor haplotypes of the btr1-A and btr1-B domesticated alleles. Here, we extended the haplotype analysis to other T. turgidum subspecies and to the BTR1 homologues in the related T. timopheevii species. Our results showed that all the domesticated wheat subspecies within T. turgidum share common BTR1-A and BTR1-B haplotypes, confirming their common origin. In T. timopheevii, however, we identified a novel loss-of-function btr1-A allele underlying a partially brittle spike phenotype. This novel recessive allele appeared fixed within the pool of domesticated Timopheev's wheat but was also carried by one wild timopheevii accession exhibiting partial brittleness. The promoter region for BTR1-B could not be amplified in any T. timopheevii accessions with any T. turgidum primer combination, exemplifying the gene-level distance between the two species. Altogether, our results support the concept of independent domestication processes for the two polyploid, wheat-related species.

Unresolved Complexity in the Gene Regulatory Network Underlying EMT.

Epithelial to mesenchymal transition (EMT) is the process whereby a polarized epithelial cell ceases to maintain cell-cell contacts, loses expression of characteristic epithelial cell markers, and acquires mesenchymal cell markers and properties such as motility, contractile ability, and invasiveness. A complex process that occurs during development and many disease states, EMT involves a plethora of transcription factors (TFs) and signaling pathways. Whilst great advances have been made in both our understanding of the progressive cell-fate changes during EMT and the gene regulatory networks that drive this process, there are still gaps in our knowledge. Epigenetic modifications are dynamic, chromatin modifying enzymes are vast and varied, transcription factors are pleiotropic, and signaling pathways are multifaceted and rarely act alone. Therefore, it is of great importance that we decipher and understand each intricate step of the process and how these players at different levels crosstalk with each other to successfully orchestrate EMT. A delicate balance and fine-tuned cooperation of gene regulatory mechanisms is required for EMT to occur successfully, and until we resolve the unknowns in this network, we cannot hope to develop effective therapies against diseases that involve aberrant EMT such as cancer. In this review, we focus on data that challenge these unknown entities underlying EMT, starting with EMT stimuli followed by intracellular signaling through to epigenetic mechanisms and chromatin remodeling.

Assessment and site-specific manipulation of DNA (hydroxy-)methylation during mouse corticogenesis.

Dynamic changes in DNA (hydroxy-)methylation are fundamental for stem cell differentiation. However, the signature of these epigenetic marks in specific cell types during corticogenesis is unknown. Moreover, site-specific manipulation of cytosine modifications is needed to reveal the significance and function of these changes. Here, we report the first assessment of (hydroxy-)methylation in neural stem cells, neurogenic progenitors, and newborn neurons during mammalian corticogenesis. We found that gain in hydroxymethylation and loss in methylation occur sequentially at specific cellular transitions during neurogenic commitment. We also found that these changes predominantly occur within enhancers of neurogenic genes up-regulated during neurogenesis and target of pioneer transcription factors. We further optimized the use of dCas9-Tet1 manipulation of (hydroxy-)methylation, locus-specifically, in vivo, showing the biological relevance of our observations for Dchs1, a regulator of corticogenesis involved in developmental malformations and cognitive impairment. Together, our data reveal the dynamics of cytosine modifications in lineage-related cell types, whereby methylation is reduced and hydroxymethylation gained during the neurogenic lineage concurrently with up-regulation of pioneer transcription factors and activation of enhancers for neurogenic genes.

The centrosome protein AKNA regulates neurogenesis via microtubule organization.

The expansion of brain size is accompanied by a relative enlargement of the subventricular zone during development. Epithelial-like neural stem cells divide in the ventricular zone at the ventricles of the embryonic brain, self-renew and generate basal progenitors1 that delaminate and settle in the subventricular zone in enlarged brain regions2. The length of time that cells stay in the subventricular zone is essential for controlling further amplification and fate determination. Here we show that the interphase centrosome protein AKNA has a key role in this process. AKNA localizes at the subdistal appendages of the mother centriole in specific subtypes of neural stem cells, and in almost all basal progenitors. This protein is necessary and sufficient to organize centrosomal microtubules, and promote their nucleation and growth. These features of AKNA are important for mediating the delamination process in the formation of the subventricular zone. Moreover, AKNA regulates the exit from the subventricular zone, which reveals the pivotal role of centrosomal microtubule organization in enabling cells to both enter and remain in the subventricular zone. The epithelial-to-mesenchymal transition is also regulated by AKNA in other epithelial cells, demonstrating its general importance for the control of cell delamination.

TILL-D: An Aegilops tauschii TILLING Resource for Wheat Improvement.

Aegilops tauschii (2n = 2x = 14, genome DD), also known as Tausch's goatgrass, is the D genome donor of bread or hexaploid wheat Triticum aestivum (2n = 2x = 42, AABBDD genome). It is a rich reservoir of useful genes for biotic and abiotic stress tolerance for wheat improvement. We developed a TILLING (Targeting Induced Local Lesions In Genomes) resource for Ae. tauschii for discovery and validation of useful genes in the D genome of wheat. The population, referred to as TILL-D, was developed with ethyl methanesulfonate (EMS) mutagen. The survival rate in M1 generation was 73%, out of which 22% plants were sterile. In the M2 generation 25% of the planted seeds showed phenotypic mutations such as albinos, chlorinas, no germination, variegated, sterile and partially fertile events, and 2,656 produced fertile M2 plants. The waxy gene was used to calculate the mutation frequency (1/70 kb) of the developed population, which was found to be higher than known mutation frequencies for diploid plants (1/89-1/1000 kb), but lower than that for a polyploid species (1/24-1/51 kb). The TILL-D resource, together with the newly published Ae. tauschii reference genome sequence, will facilitate gene discoveries and validations of agronomically important traits and their eventual fine transfer in bread wheat.

Stage-Specific Transcription Factors Drive Astrogliogenesis by Remodeling Gene Regulatory Landscapes.

A broad molecular framework of how neural stem cells are specified toward astrocyte fate during brain development has proven elusive. Here we perform comprehensive and integrated transcriptomic and epigenomic analyses to delineate gene regulatory programs that drive the developmental trajectory from mouse embryonic stem cells to astrocytes. We report molecularly distinct phases of astrogliogenesis that exhibit stage- and lineage-specific transcriptomic and epigenetic signatures with unique primed and active chromatin regions, thereby revealing regulatory elements and transcriptional programs underlying astrocyte generation and maturation. By searching for transcription factors that function at these elements, we identified NFIA and ATF3 as drivers of astrocyte differentiation from neural precursor cells while RUNX2 promotes astrocyte maturation. These transcription factors facilitate stage-specific gene expression programs by switching the chromatin state of their target regulatory elements from primed to active. Altogether, these findings provide integrated insights into the genetic and epigenetic mechanisms steering the trajectory of astrogliogenesis.

FBXO32 promotes microenvironment underlying epithelial-mesenchymal transition via CtBP1 during tumour metastasis and brain development.

The set of events that convert adherent epithelial cells into migratory cells are collectively known as epithelial-mesenchymal transition (EMT). EMT is involved during development, for example, in triggering neural crest migration, and in pathogenesis such as metastasis. Here we discover FBXO32, an E3 ubiquitin ligase, to be critical for hallmark gene expression and phenotypic changes underlying EMT. Interestingly, FBXO32 directly ubiquitinates CtBP1, which is required for its stability and nuclear retention. This is essential for epigenetic remodeling and transcriptional induction of CtBP1 target genes, which create a suitable microenvironment for EMT progression. FBXO32 is also amplified in metastatic cancers and its depletion in a NSG mouse xenograft model inhibits tumor growth and metastasis. In addition, FBXO32 is essential for neuronal EMT during brain development. Together, these findings establish that FBXO32 acts as an upstream regulator of EMT by governing the gene expression program underlying this process during development and disease.

ERK signalling modulates epigenome to drive epithelial to mesenchymal transition.

The series of events that allow the conversion from adherent epithelial cells into migratory cells is collectively known as epithelial-mesenchymal transition (EMT). EMT is employed during embryonic development such as for gastrulation and neural crest migration and is misused in diseases, such as cancer metastasis. ERK signalling is known to be essential for EMT, however its influence on the epigenetic and transcriptional programme underlying EMT is poorly understood. Here, using a comprehensive genome-wide analysis of H3K27ac mark and gene expression in mammary epithelial cells undergoing EMT, we found that ERK signalling is essential for the epigenetic reprogramming underlying hallmark gene expression and phenotypic changes of EMT. We show that the chemical inhibition of Erk signalling during EMT prevents the loss and gain of the H3K27ac mark at regulatory regions of epithelial and mesenchymal genes, respectively, and results in a transcriptome and epigenome closer to those of epithelial cells. Further computational analyses identified a distinct set of transcription factor motifs enriched at distal regulatory regions that are epigenetically remodelled by ERK signalling. Altogether, our findings reveal an ERK-dependent epigenetic remodelling of regulatory elements that results in a gene expression programme essential for driving EMT.

GSP: a web-based platform for designing genome-specific primers in polyploids.

MOTIVATION: The sequences among subgenomes in a polyploid species have high similarity, making it difficult to design genome-specific primers for sequence analysis. RESULTS: We present GSP, a web-based platform to design genome-specific primers that distinguish subgenome sequences in a polyploid genome. GSP uses BLAST to extract homeologous sequences of the subgenomes in existing databases, performs a multiple sequence alignment, and design primers based on sequence variants in the alignment. An interactive primers diagram, a sequence alignment viewer and a virtual electrophoresis are displayed as parts of the primer design result. GSP also designs specific primers from multiple sequences uploaded by users. AVAILABILITY AND IMPLEMENTATION: GSP is a user-friendly and efficient web platform freely accessible at http://probes.pw.usda.gov/GSP Source code and command-line application are available at https://github.com/bioinfogenome/GSP CONTACTS: yong.gu@ars.usda.gov or devin.coleman-derr@ars.usda.gov SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.