Model For Functional Genomics Research Articles

Construction of ethyl methane sulfonate mutant library in G. arboreum and rapid identification of mutant genes via repeated re-sequencing

Cotton (Gossypium) is the leading economic crop and provides most of renewable spinning fibers for the worldwide textile industry. The diploid Asiatic cotton (G. arboreum), with spinnable fiber and relatively simple genome, is an ideal model for cotton functional genomics and fiber biology. In this study, we constructed an ethyl methane sulfonate (EMS) mutant library of Asiatic cotton. Stable mutants with various phenotypic variations in leaf, flower, fiber and plant architecture were obtained in M1 generation. Genome-wide mutation analysis revealed high-density mutations distributed evenly on whole genome. Using repeated re-sequencing of M1 mutant plants, the candidate genes of 10 out of 15 mutants were identified. The functions of GaYUC4 and GaPDX1, candidate genes for bar blade and yellow leaf vein, respectively, were further confirmed using virus induced gene silencing in Asiatic and/or upland cottons. Our work paved way for expansion of Asiatic cotton mutant library and subsequent applications in cotton functional genomics and cotton breeding.

Chinese cabbage: an emerging model for functional genomics in leafy vegetable crops.

Leafy vegetable crops (LVCs) are consumed worldwide and offer essential nutrients for humans. Unlike model plant species, systematic characterisation of gene function is lacking, although whole-genome sequences (WGSs) are available for various LVCs. Several recent studies in Chinese cabbage have reported high-density mutant populations linking genotype to phenotype, providing blueprints for functional LVC genomics and beyond.

Identification of SNPs and Candidate Genes Associated with Salt Stress in Two Korean Sorghum Cultivars and Understanding Selection Pressures in the Breeding Process

One of the abiotic stresses, salt stress, has an impact on the production and development of crops around the world. Sorghum is a functional genomics model crop of C4 plants due to its small genome size, and it is suitable for providing a clue to the mechanism associated with salt tolerance at the transcriptomic level. However, the mechanism of salt-related genes in sorghum has not been well described. RNA sequencing, using QuantSeq, was performed on two Korean cultivars, ‘Sodamchal’ and ‘Nampungchal’, which are known to have different intensities in response to salt stress, between a control and high-salinity treatment over a different time-course. In addition, physiological responses such as the proline, anthocyanin, chlorophyll, and reducing sugar contents were evaluated under the salt-stress treatment between these two sorghum cultivars. Moreover, differentially expressed genes (DEGs) between the Nampungchal and Sodamchal cultivars were identified in their leaves and roots, respectively. Moreover, the function of DEGs was confirmed through GO classification and KEGG pathway. We also analyzed the correlation between the selection pressure with DEGs by identifying Ka/Ks of DEGs. In the breeding process, the role of positive or negative selected genes was analyzed. Therefore, a new hypothesis on selection pressure was proposed from the breeding perspective of cultivars. A comparative analysis of the two sorghum cultivars provides candidate genes involved in the salt-stress response and may offer a better understanding of the salt-tolerance mechanism in sorghum.

High-Quality Genome Assembly and Comprehensive Transcriptome of the Painted Lady Butterfly Vanessa cardui.

The painted lady butterfly, Vanessa cardui, has the longest migration routes, the widest hostplant diversity, and one of the most complex wing patterns of any insect. Due to minimal culturing requirements, easily characterized wing pattern elements, and technical feasibility of CRISPR/Cas9 genome editing, V. cardui is emerging as a functional genomics model for diverse research programs. Here, we report a high-quality, annotated genome assembly of the V. cardui genome, generated using 84× coverage of PacBio long-read data, which we assembled into 205 contigs with a total length of 425.4 Mb (N50 = 10.3 Mb). The genome was very complete (single-copy complete Benchmarking Universal Single-Copy Orthologs [BUSCO] 97%), with contigs assembled into presumptive chromosomes using synteny analyses. Our annotation used embryonic, larval, and pupal transcriptomes, and 20 transcriptomes across five different wing developmental stages. Gene annotations showed a high level of accuracy and completeness, with 14,437 predicted protein-coding genes. This annotated genome assembly constitutes an important resource for diverse functional genomic studies ranging from the developmental genetic basis of butterfly color pattern, to coevolution with diverse hostplants.

Sorghum breeding in the genomic era: opportunities and challenges.

Key messageThe importance and potential of the multi-purpose crop sorghum in global food security have not yet been fully exploited, and the integration of the state-of-art genomics and high-throughput technologies into breeding practice is required.Sorghum, a historically vital staple food source and currently the fifth most important major cereal, is emerging as a crop with diverse end-uses as food, feed, fuel and forage and a model for functional genetics and genomics of tropical grasses. Rapid development in high-throughput experimental and data processing technologies has significantly speeded up sorghum genomic researches in the past few years. The genomes of three sorghum lines are available, thousands of genetic stocks accessible and various genetic populations, including NAM, MAGIC, and mutagenised populations released. Functional and comparative genomics have elucidated key genetic loci and genes controlling agronomical and adaptive traits. However, the knowledge gained has far away from being translated into real breeding practices. We argue that the way forward is to take a genome-based approach for tailored designing of sorghum as a multi-functional crop combining excellent agricultural traits for various end uses. In this review, we update the new concepts and innovation systems in crop breeding and summarise recent advances in sorghum genomic researches, especially the genome-wide dissection of variations in genes and alleles for agronomically important traits. Future directions and opportunities for sorghum breeding are highlighted to stimulate discussion amongst sorghum academic and industrial communities.

A zebrafish functional genomics model to investigate the role of human A20 variants in vivo

Germline loss-of-function variation in TNFAIP3, encoding A20, has been implicated in a wide variety of autoinflammatory and autoimmune conditions, with acquired somatic missense mutations linked to cancer progression. Furthermore, human sequence data reveals that the A20 locus contains ~ 400 non-synonymous coding variants, which are largely uncharacterised. The growing number of A20 coding variants with unknown function, but potential clinical impact, poses a challenge to traditional mouse-based approaches. Here we report the development of a novel functional genomics approach that utilizes a new A20-deficient zebrafish (Danio rerio) model to investigate the impact of TNFAIP3 genetic variants in vivo. A20-deficient zebrafish are hyper-responsive to microbial immune activation and exhibit spontaneous early lethality. Ectopic addition of human A20 rescued A20-null zebrafish from lethality, while missense mutations at two conserved A20 residues, S381A and C243Y, reversed this protective effect. Ser381 represents a phosphorylation site important for enhancing A20 activity that is abrogated by its mutation to alanine, or by a causal C243Y mutation that triggers human autoimmune disease. These data reveal an evolutionarily conserved role for TNFAIP3 in limiting inflammation in the vertebrate linage and show how this function is controlled by phosphorylation. They also demonstrate how a zebrafish functional genomics pipeline can be utilized to investigate the in vivo significance of medically relevant human TNFAIP3 gene variants.

The evolutionary origin and domestication history of goldfish (Carassius auratus)

Goldfish have been subjected to over 1,000 y of intensive domestication and selective breeding. In this report, we describe a high-quality goldfish genome (2n = 100), anchoring 95.75% of contigs into 50 pseudochromosomes. Comparative genomics enabled us to disentangle the two subgenomes that resulted from an ancient hybridization event. Resequencing 185 representative goldfish variants and 16 wild crucian carp revealed the origin of goldfish and identified genomic regions that have been shaped by selective sweeps linked to its domestication. Our comprehensive collection of goldfish varieties enabled us to associate genetic variations with a number of well-known anatomical features, including features that distinguish traditional goldfish clades. Additionally, we identified a tyrosine-protein kinase receptor as a candidate causal gene for the first well-known case of Mendelian inheritance in goldfish-the transparent mutant. The goldfish genome and diversity data offer unique resources to make goldfish a promising model for functional genomics, as well as domestication.

Multiple molecular defense strategies in Brachypodium distachyon surmount Hessian fly (Mayetiola destructor) larvae-induced susceptibility for plant survival

The Hessian fly is a destructive pest of wheat causing severe economic damage. Numerous genes and associated biological pathways have been implicated in defense against Hessian fly. However, due to limited genetic resources, compounded with genome complexity, functional analysis of the candidate genes are challenging in wheat. Physically, Brachypodium distachyon (Bd) exhibits nonhost resistance to Hessian fly, and with a small genome size, short life cycle, vast genetic resources and amenability to transformation, it offers an alternate functional genomic model for deciphering plant-Hessian fly interactions. RNA-sequencing was used to reveal thousands of Hessian fly-responsive genes in Bd one, three, and five days after egg hatch. Genes encoding defense proteins, stress-regulating transcription factors, signaling kinases, and secondary metabolites were strongly up-regulated within the first 24 hours of larval feeding indicating an early defense, similar to resistant wheat. Defense was mediated by a hypersensitive response that included necrotic lesions, up-regulated ROS-generating and -scavenging enzymes, and H2O2 production. Suppression of cell wall-associated proteins and increased cell permeability in Bd resembled susceptible wheat. Thus, Bd molecular responses shared similarities to both resistant and susceptible wheat, validating its suitability as a model genome for undertaking functional studies of candidate Hessian fly-responsive genes.

Plant Growth and Agrobacterium-mediated Floral-dip Transformation of the Extremophyte Schrenkiella parvula.

Schrenkiella parvula is an extremophyte adapted to various abiotic stresses, including multiple ion toxicity stresses. Despite high-quality genomic resources available to study how plants adapt to environmental stresses, its value as a functional genomics model and tool has been limited by the lack of a feasible transformation system. In this protocol, we report how to generate stable transgenic S. parvula lines using an Agrobacterium-mediated floral-dip method. We modified the transformation protocol used for A. thaliana to account for unique traits of S. parvula, such as an indeterminate flowering habit and a high epicuticular wax content on leaves. Briefly, S. parvula seeds were stratified at 4 °C for five days before planting. Plants were grown at a photoperiod of a 14 h light and 10 h dark and a 130 µmol m-2s-1 light intensity, at 22 °C to 24 °C. Eight to nine week-old plants with multiple inflorescences were selected for transformation. These inflorescences were dipped in an infiltration solution of Agrobacterium tumefaciens GV3101 carrying the pMP90RK plasmid. We performed two rounds of flower dipping with an interval of three to four weeks to increase the transformation efficiency. The T1 seeds were collected and dried for four weeks in a container with desiccants before germination to screen for candidate transformed lines. Resistance to BASTA was used to screen T1 plants. We sprayed the BASTA solution three times with an interval of three days starting at two week-old plants to reduce false positives. A BASTA drop test was performed on surviving individual plants to identify true positive transformants. The transformation efficiency was 0.033%, yielding 3-4 transgenic plants per 10,000 T1 seeds propagated.

Comprehensive functional genomic resource and integrative model for the human brain.

Despite progress in defining genetic risk for psychiatric disorders, their molecular mechanisms remain elusive. Addressing this, the PsychENCODE Consortium has generated a comprehensive online resource for the adult brain across 1866 individuals. The PsychENCODE resource contains ~79,000 brain-active enhancers, sets of Hi-C linkages, and topologically associating domains; single-cell expression profiles for many cell types; expression quantitative-trait loci (QTLs); and further QTLs associated with chromatin, splicing, and cell-type proportions. Integration shows that varying cell-type proportions largely account for the cross-population variation in expression (with >88% reconstruction accuracy). It also allows building of a gene regulatory network, linking genome-wide association study variants to genes (e.g., 321 for schizophrenia). We embed this network into an interpretable deep-learning model, which improves disease prediction by ~6-fold versus polygenic risk scores and identifies key genes and pathways in psychiatric disorders.

Transcriptome-Wide Analysis Reveals Modulation of Human Macrophage Inflammatory Phenotype Through Alternative Splicing.

Human macrophages can shift phenotype across the inflammatory M1 and reparative M2 spectrum in response to environmental challenges, but the mechanisms promoting inflammatory and cardiometabolic disease-associated M1 phenotypes remain incompletely understood. Alternative splicing (AS) is emerging as an important regulator of cellular function, yet its role in macrophage activation is largely unknown. We investigated the extent to which AS occurs in M1 activation within the cardiometabolic disease context and validated a functional genomic cell model for studying human macrophage-related AS events. From deep RNA-sequencing of resting, M1, and M2 primary human monocyte-derived macrophages, we found 3860 differentially expressed genes in M1 activation and detected 233 M1-induced AS events; the majority of AS events were cell- and M1-specific with enrichment for pathways relevant to macrophage inflammation. Using genetic variant data for 10 cardiometabolic traits, we identified 28 trait-associated variants within the genomic loci of 21 alternatively spliced genes and 15 variants within 7 differentially expressed regulatory splicing factors in M1 activation. Knockdown of 1 such splicing factor, CELF1, in primary human macrophages led to increased inflammatory response to M1 stimulation, demonstrating CELF1's potential modulation of the M1 phenotype. Finally, we demonstrated that an induced pluripotent stem cell-derived macrophage system recapitulates M1-associated AS events and provides a high-fidelity macrophage AS model. AS plays a role in defining macrophage phenotype in a cell- and stimulus-specific fashion. Alternatively spliced genes and splicing factors with trait-associated variants may reveal novel pathways and targets in cardiometabolic diseases.

Divergent clonal selection dominates medulloblastoma at recurrence.

The development of targeted anti-cancer therapies through the study of cancer genomes is intended to increase survival rates and decrease treatment-related toxicity. We treated a transposon-driven, functional genomic mouse model of medulloblastoma with 'humanized' in vivo therapy (microneurosurgical tumour resection followed by multi-fractionated, image-guided radiotherapy). Genetic events in recurrent murine medulloblastoma exhibit a very poor overlap with those in matched murine diagnostic samples (<5%). Whole-genome sequencing of 33 pairs of human diagnostic and post-therapy medulloblastomas demonstrated substantial genetic divergence of the dominant clone after therapy (<12% diagnostic events were retained at recurrence). In both mice and humans, the dominant clone at recurrence arose through clonal selection of a pre-existing minor clone present at diagnosis. Targeted therapy is unlikely to be effective in the absence of the target, therefore our results offer a simple, proximal, and remediable explanation for the failure of prior clinical trials of targeted therapy.

A functional genomic model for predicting prognosis in idiopathic pulmonary fibrosis.

BackgroundThe course of disease for patients with idiopathic pulmonary fibrosis (IPF) is highly heterogeneous. Prognostic models rely on demographic and clinical characteristics and are not reproducible. Integrating data from genomic analyses may identify novel prognostic models and provide mechanistic insights into IPF.MethodsTotal RNA of peripheral blood mononuclear cells was subjected to microarray profiling in a training (45 IPF individuals) and two independent validation cohorts (21 IPF/10 controls, and 75 IPF individuals, respectively). To identify a gene set predictive of IPF prognosis, we incorporated genomic, clinical, and outcome data from the training cohort. Predictor genes were selected if all the following criteria were met: 1) Present in a gene co-expression module from Weighted Gene Co-expression Network Analysis (WGCNA) that correlated with pulmonary function (p < 0.05); 2) Differentially expressed between observed “good” vs. “poor” prognosis with fold change (FC) >1.5 and false discovery rate (FDR) < 2 %; and 3) Predictive of mortality (p < 0.05) in univariate Cox regression analysis. “Survival risk group prediction” was adopted to construct a functional genomic model that used the IPF prognostic predictor gene set to derive a prognostic index (PI) for each patient into either high or low risk for survival outcomes. Prediction accuracy was assessed with a repeated 10-fold cross-validation algorithm and independently assessed in two validation cohorts through multivariate Cox regression survival analysis.ResultsA set of 118 IPF prognostic predictor genes was used to derive the functional genomic model and PI. In the training cohort, high-risk IPF patients predicted by PI had significantly shorter survival compared to those labeled as low-risk patients (log rank p < 0.001). The prediction accuracy was further validated in two independent cohorts (log rank p < 0.001 and 0.002). Functional pathway analysis revealed that the canonical pathways enriched with the IPF prognostic predictor gene set were involved in T-cell biology, including iCOS, T-cell receptor, and CD28 signaling.ConclusionsUsing supervised and unsupervised analyses, we identified a set of IPF prognostic predictor genes and derived a functional genomic model that predicted high and low-risk IPF patients with high accuracy. This genomic model may complement current prognostic tools to deliver more personalized care for IPF patients.Electronic supplementary materialThe online version of this article (doi:10.1186/s12890-015-0142-8) contains supplementary material, which is available to authorized users.

Generation of a Transcriptome in a Model Lepidopteran Pest, Heliothis virescens, Using Multiple Sequencing Strategies for Profiling Midgut Gene Expression.

Heliothine pests such as the tobacco budworm, Heliothis virescens (F.), pose a significant threat to production of a variety of crops and ornamental plants and are models for developmental and physiological studies. The efforts to develop new control measures for H. virescens, as well as its use as a relevant biological model, are hampered by a lack of molecular resources. The present work demonstrates the utility of next-generation sequencing technologies for rapid molecular resource generation from this species for which lacks a sequenced genome. In order to amass a de novo transcriptome for this moth, transcript sequences generated from Illumina, Roche 454, and Sanger sequencing platforms were merged into a single de novo transcriptome assembly. This pooling strategy allowed a thorough sampling of transcripts produced under diverse environmental conditions, developmental stages, tissues, and infections with entomopathogens used for biological control, to provide the most complete transcriptome to date for this species. Over 138 million reads from the three platforms were assembled into the final set of 63,648 contigs. Of these, 29,978 had significant BLAST scores indicating orthologous relationships to transcripts of other insect species, with the top-hit species being the monarch butterfly (Danaus plexippus) and silkworm (Bombyx mori). Among identified H. virescens orthologs were immune effectors, signal transduction pathways, olfactory receptors, hormone biosynthetic pathways, peptide hormones and their receptors, digestive enzymes, and insecticide resistance enzymes. As an example, we demonstrate the utility of this transcriptomic resource to study gene expression profiling of larval midguts and detect transcripts of putative Bacillus thuringiensis (Bt) Cry toxin receptors. The substantial molecular resources described in this study will facilitate development of H. virescens as a relevant biological model for functional genomics and for new biological experimentation needed to develop efficient control efforts for this and related Noctuid pest moths.

Neocaridina denticulata: A Decapod Crustacean Model for Functional Genomics.

A decapod crustacean model is needed for understanding the molecular mechanisms underlying physiological processes, such as reproduction, sex determination, molting and growth, immunity, regeneration, and response to stress. Criteria for selection are: life-history traits, adult size, availability and ease of culture, and genomics and genetic manipulation. Three freshwater species are considered: cherry shrimp, Neocaridina denticulata; red swamp crayfish, Procambarus clarkii; and redclaw crayfish, Cherax quadricarinatus. All three are readily available, reproduce year round, and grow rapidly. The crayfish species require more space for culture than does N. denticulata. The transparent cuticle of cherry shrimp provides for direct assessment of reproductive status, stage of molt, and tissue-specific expression of reporter genes, and facilitates screening of mutations affecting phenotype. Moreover, a preliminary genome of N. denticulata is available and efforts toward complete genome sequencing and transcriptome sequencing have been initiated. Neocaridina denticulata possesses the best combination of traits that make it most suitable as a model for functional genomics. The next step is to obtain the complete genome sequence and to develop molecular technologies for the screening of mutants and for manipulating tissue-specific gene expression.

System approaches to study root hairs as a single cell plant model: current status and future perspectives.

Our current understanding of plant functional genomics derives primarily from measurements of gene, protein and/or metabolite levels averaged over the whole plant or multicellular tissues. These approaches risk diluting the response of specific cells that might respond strongly to the treatment but whose signal is diluted by the larger proportion of non-responding cells. For example, if a gene is expressed at a low level, does this mean that it is indeed lowly expressed or is it highly expressed, but only in a few cells? In order to avoid these issues, we adopted the soybean root hair cell, derived from a single, differentiated root epidermal cell, as a single-cell model for functional genomics. Root hair cells are intrinsically interesting since they are major conduits for root water and nutrient uptake and are also the preferred site of infection by nitrogen-fixing rhizobium bacteria. Although a variety of other approaches have been used to study single plant cells or single cell types, the root hair system is perhaps unique in allowing application of the full repertoire of functional genomic and biochemical approaches. In this mini review, we summarize our published work and place this within the broader context of root biology, with a significant focus on understanding the initial events in the soybean-rhizobium interaction.

Inflammatory environment and oxidized LDL convert circulating human proangiogenic cells into functional antigen-presenting cells.

The function of human circulating PACs has been described extensively. However, little focus has been placed on understanding how these cells differ in their functions in the presence of microenvironments mimicking vascular inflammation. We hypothesized that exposure to proinflammatory cytokines or the oxLDL, an autoantigen abundant in advanced atherosclerotic plaques, converts PACs into immune-modulating/proinflammatory cells. Hence, we examined the effect of oxLDL and inflammatory stimuli on their phenotype by use of a functional genomics model based on secretome and whole genome transcriptome profiling. PACs obtained from culturing a PBMC fraction in angiogenic medium were primed with DC differentiation cytokines and then exposed to proinflammatory cytokines or oxLDL. Under these conditions, PACs converted into APCs, expressed maturation markers CD80 and CD83, and showed an increased up-regulation of CD86. APCcy and APCox induced a robust T cell BrdU incorporation. Despite a similar ability to induce lymphocyte proliferation, APCcy and APCox differed for the secretory pathway and mRNA expression. Analysis of the differentially expressed genes identified 4 gene "clusters," showing reciprocal modulation in APCcy vs. APCox, justifying, according to functional genomics analyses, a different putative function of the cells in antigen processing. Together, these data show that treatment with inflammatory cytokines or oxLDL converts human PAC phenotypes and functions into that of APCs with similar lymphocyte-activating ability but distinct maturation degree and paracrine functions.

Establishing a multidisciplinary context for modeling 3D facial shape from DNA.

Establishing a multidisciplinary context for modeling 3D facial shape from DNA.

Novel genomes and genome constitutions identified by GISH and 5S rDNA and knotted1 genomic sequences in the genus Setaria

BackgroundThe Setaria genus is increasingly of interest to researchers, as its two species, S. viridis and S. italica, are being developed as models for understanding C4 photosynthesis and plant functional genomics. The genome constitution of Setaria species has been studied in the diploid species S. viridis, S. adhaerans and S. grisebachii, where three genomes A, B and C were identified respectively. Two allotetraploid species, S. verticillata and S. faberi, were found to have AABB genomes, and one autotetraploid species, S. queenslandica, with an AAAA genome, has also been identified. The genomes and genome constitutions of most other species remain unknown, even though it was thought there are approximately 125 species in the genus distributed world-wide.ResultsGISH was performed to detect the genome constitutions of Eurasia species of S. glauca, S. plicata, and S. arenaria, with the known A, B and C genomes as probes. No or very poor hybridization signal was detected indicating that their genomes are different from those already described. GISH was also performed reciprocally between S. glauca, S. plicata, and S. arenaria genomes, but no hybridization signals between each other were found. The two sets of chromosomes of S. lachnea both hybridized strong signals with only the known C genome of S. grisebachii. Chromosomes of Qing 9, an accession formerly considered as S. viridis, hybridized strong signal only to B genome of S. adherans. Phylogenetic trees constructed with 5S rDNA and knotted1 markers, clearly classify the samples in this study into six clusters, matching the GISH results, and suggesting that the F genome of S. arenaria is basal in the genus.ConclusionsThree novel genomes in the Setaria genus were identified and designated as genome D (S. glauca), E (S. plicata) and F (S. arenaria) respectively. The genome constitution of tetraploid S. lachnea is putatively CCC’C’. Qing 9 is a B genome species indigenous to China and is hypothesized to be a newly identified species. The difference in genome constitution and origin of S. verticillata and S. faberi is also discussed. The new genomes and the genome constitutions of Setaria species identified in this report provide useful information for Setaria germplasm management, foxtail millet breeding, grass evolution and the development of S. viridis and S. italica as a new model for functional genomics.

Genetic transformation of western clover (Trifolium occidentale D. E. Coombe.) as a model for functional genomics and transgene introgression in clonal pasture legume species.

BackgroundWestern clover (Trifolium occidentale) is a perennial herb with characteristics compatible for its development as an attractive model species for genomics studies relating to the forage legume, white clover (Trifolium repens). Its characteristics such as a small diploid genome, self-fertility and ancestral contribution of one of the genomes of T. repens, facilitates its use as a model for genetic analysis of plants transformed with legume or novel genes.ResultsIn this study, a reproducible transformation protocol was established following screening of T. occidentale accessions originating from England, Ireland, France, Spain and Portugal. The protocol is based upon infection of cotyledonary explants dissected from mature seed with the Agrobacterium tumefaciens strain GV3101 carrying vectors which contain the bar selection marker gene. Transformation frequencies of up to 7.5% were achieved in 9 of the 17 accessions tested. Transformed plants were verified by PCR and expression of the gusA reporter gene, while integration of the T-DNA was confirmed by Southern blot hybridisation and segregation of progeny in the T1 generation.ConclusionsDevelopment of this protocol provides a valuable contribution toward establishing T. occidentale as a model species for white clover. This presents opportunities for further improvement in white clover through the application of biotechnology.