Evaluation of genome sequencing quality in selected plant species using expressed sequence tags

High-quality reference genome sequences of two Cannaceae species provide insights into the evolution of Cannaceae

Canna edulis Ker-Gawl and Canna indica L. are species belonging to the Cannaceae family and both have a very high economic value. Here, we aimed to assemble genomes of C. edulis and C. indica at the chromosome level to generate a reference genome for the Cannaceae family. We also comparatively analyzed the genomes of C. edulis and C. indica and examined the molecular mechanisms responsible for the remarkable differences in plant characteristics in C. edulis varieties. Our results indicated that genome-wide duplication events had recently occurred in C. edulis and C. indica. The comparative analysis of the genomes of C. edulis and C. indica revealed that C. edulis exhibited a remarkable level of replication of genes in the starch and sucrose metabolic pathways, especially during sucrose hydrolysis. This finding is consistent with the fact that the starch content of the C. edulis tuber is higher than that of C. indica. Simplified genome re-sequencing revealed the population structure of 241 C. edulis genes, and a genome-wide association study of leaf traits revealed the location of key genes related to leaf color and morphology. These findings extend our understanding of Cannaceade at the molecular level, and provide an effective theoretical basis for further study and utilization of Cannaceae plants.

A First Genome Survey and Genomic SSR Marker Analysis of Trematomus loennbergii Regan, 1913

Simple Summary

The scaly rockcod (Trematomus loennbergii) is distributed in the Antarctic Ocean and this area is isolated by the Antarctic Circumpolar Current. It is important region to study evolutionary diversity. Trematomus is the main genus, having 11 species, and their habit distribution is well known. However, their genetic and genomic information is not studied. In addition, some species have similar morphology. In this study, a genome survey of T. loennbergii and microsatellite motif analysis were conducted to obtain genomic profile. The fundamental data such as genome size, heterozygosity ratio, duplication ration and microsatellite motifs were obtained. These data will provide a foundation for further whole-genome sequencing and the development of new molecular markers of T. loennbergii.

Abstract

Trematomus loennbergii Regan, 1913, is an evolutionarily important marine fish species distributed in the Antarctic Ocean. However, its genome has not been studied to date. In the present study, whole genome sequencing was performed using next-generation sequencing (NGS) technology to characterize its genome and develop genomic microsatellite markers. The 25-mer frequency distribution was estimated to be the best, and the genome size was predicted to be 815,042,992 bp. The heterozygosity, average rate of read duplication, and sequencing error rates were 0.536%, 0.724%, and 0.292%, respectively. These data were used to analyze microsatellite markers, and a total of 2,264,647 repeat motifs were identified. The most frequent repeat motif was di-nucleotide with 87.00% frequency, followed by tri-nucleotide (10.45%), tetra-nucleotide (1.94%), penta-nucleotide (0.34%), and hexa-nucleotide (0.27%). The AC repeat motif was the most abundant motif among di-nucleotides and among all repeat motifs. Among microsatellite markers, 181 markers were selected and PCR technology was used to validate several markers. A total of 15 markers produced only one band. In summary, these results provide a good basis for further studies, including evolutionary biology studies and population genetics of Antarctic fish species.

The karyotype, genome survey, and assembly of Mud artemisia (Artemisia selengensis)

BACKGROUND

Artemisia selengensis is traditional Chinese medicine and phytochemical analysis indicated that A. selengensis contains essential oils, fatty acids and phenolic acids. The lack of reference genomic information may lead to tardiness in molecular biology research of A. selengensis.

METHOD AND RESULTS

Karyotype analysis, genome survey, and genome assembly was employed to acquire information on the genome structure of A. selengensis. The chromosome number is 2n = 2x = 36, karyotype formula is 28 m + 8Sm, karyotype asymmetry coefficient is 58.8%, and karyotypes were symmetric to Stebbins' type 2A. Besides, the flow cytometry findings reported that the mean peak value of fluorescent intensity is 1,170,677, 2C DNA content is 12 pg and the genome size was estimated to be approximately 5.87 Gb. Furthermore, the genome survey generates 341,478,078 clean reads, unfortunately, after K-mer analysis, no significant peak can be observed, the heterozygosity, repetitive rate and genome size was unable to estimated. It is speculated that this phenomenon might be due to the complexity of genome structure. 37,266 contigs are preliminary assembled with Oxford Nanopore Technology (ONT) sequencing, totaling 804 Mb and GC content was 34.08%. The total length is 804,475,881 bp, N50 is 29,624 bp, and the largest contig length is 239,792 bp.

CONCLUSION

This study reveals the preliminary information of genome size of A. selengensis. These findings may provide supportive information for sequencing and assembly of whole-genome sequencing and encourage the progress of functional gene discovery, genetic improvement, evolutionary study, and structural studies of A. selengensis.

Genome survey of Zanthoxylum bungeanum and development of genomic-SSR markers in congeneric species

Zanthoxylum bungeanum, a spice and medicinal plant, is cultivated in many parts of China and some countries in Southeast Asia; however, data on its genome are lacking. In the present study, we performed a whole-genome survey and developed novel genomic-SSR markers of Z. bungeanum. Clean data (∼197.16 Gb) were obtained and assembled into 11185221 scaffolds with an N50 of 183 bp. K-mer analysis revealed that Z. bungeanum has an estimated genome size of 3971.92 Mb, and the GC content, heterozygous rate, and repeat sequence rate are 37.21%, 1.73%, and 86.04%, respectively. These results indicate that the genome of Z. bungeanum is complex. Furthermore, 27153 simple sequence repeat (SSR) loci were identified from 57288 scaffolds with a minimum length > 1 kb. Mononucleotide repeats (19706) were the most abundant type, followed by dinucleotide repeats (5154). The most common motifs were A/T, followed by AT/AT; these SSRs accounted for 71.42% and 11.84% of all repeats, respectively. A total of 21243 non-repeating primer pairs were designed, and 100 were randomly selected and validated by PCR analysis using DNA from 10 Z. bungeanum individuals and 5 Zanthoxylum armatum individuals. Finally, 36 polymorphic SSR markers were developed with polymorphism information content (PIC) values ranging from 0.16 to 0.75. Cluster analysis revealed that Z. bungeanum and Z. armatum could be divided into two major clusters, suggesting that these newly developed SSR markers are useful for genetic diversity and germplasm resource identification in Z. bungeanum and Z. armatum.

Genome survey sequencing and identification of genomic SSR markers for Rhododendron micranthum

Rhododendron micranthum is an evergreen shrub species widely distributed in China that has high ornamental and medicinal value. However, there is a lack of molecular and genomic data for this plant, which severely restricts the development of its relevant research. The objective of the present study was to conduct a first genomic survey of R. micranthum and determine its whole-genome sequencing scheme. Next-generation sequencing (Illumina Hi-Seq Xten) was used to measure the genome size of R. micranthum, K-mer analysis were employed to investigate its genomic profile. Finally, we conducted bioinformatics methods to performed SSR (simple sequence repeat) prediction based on the genomic data. The genome size of R. micranthum was estimated to be 554.22 Mb. The heterozygosity ratio was 0.93%, and the sequence repeat ratio was calculated to be 49.17%. The clean reads of R. micranthum were assembled into 2281551 scaffolds with a N50 value of 916 bp. A total of 479724 SSR molecular markers were identified in the R. micranthum genome, and 871656 pairs of primers designed for application. Among of them, 100 primer pairs were validated, and 71 primer pairs were successfully amplified. In summary, the R. micranthum genome is complex with high heterozygosity and low repeated sequences. In future whole-genome research in R. micranthum, higher-depth '2+3' (Illumina+PacBio) sequencing may yield better assembly results.

Characterization of DNA methylation variations during fruit development and ripening of Vitis vinifera (cv. 'Fujiminori')

The fruit is the most important economical organ in the grape; accordingly, to investigate the grapevine genomic methylation landscape and examine its functional significance during fruit development, we generated whole genome DNA methylation maps for various developmental stages in the fruit of grapevine. In this study, thirteen DNA methylation-related genes and their expression profiles were identified and analyzed. The methylation levels for mC, mCG, mCHG, and mCHH contexts in 65 days after flowering (65DAF) fruit (véraison stage) were higher than those in 40DAF (green stage) and 90DAF (mature stage) fruits. Relative to methylation in the mC context, methylation levels in the mCHH context were higher than those of mCG and mCHG. The DNA methylation level in the ncRNA regions was significantly higher than that in exon, gene, intron, and mRNA regions. The differentially methylated regions (DMRs) and differentially methylated promoters (DMPs) in 65DAF_vs_40DAF were both higher than those in 90DAF_vs_65DAF and 90DAF_vs_40DAF. Most DMRs (or DMPs) were involved in metabolic processes and cell processes, binding, and catalytic activity. These results indicated that DNA methylation represses gene expression during grape fruit development, and it broadens our understanding of the landscape and function of DNA methylation in grapevine genomes.

Genome survey and SSR analysis of Apocynum venetum

Apocynum venetum is an eco-economic plant that exhibits high stress resistance. In the present paper, we carried out a whole-genome survey of A. venetum in order to provide a foundation for its whole-genome sequencing. High-throughput sequencing technology (Illumina NovaSep) was first used to measure the genome size of A. venetum, and bioinformatics methods were employed for the evaluation of the genome size, heterozygosity ratio, repeated sequences, and GC content in order to provide a foundation for subsequent whole-genome sequencing. The sequencing analysis results indicated that the preliminary estimated genome size of A. venetum was 254.40 Mbp, and its heterozygosity ratio and percentage of repeated sequences were 0.63 and 40.87%, respectively, indicating that it has a complex genome. We used k-mer = 41 to carry out a preliminary assembly and obtained contig N50, which was 3841 bp with a total length of 223949699 bp. We carried out further assembly to obtain scaffold N50, which was 6196 bp with a total length of 227322054 bp. We performed simple sequence repeat (SSR) molecular marker prediction based on the A. venetum genome data and identified a total of 101918 SSRs. The differences between the different types of nucleotide repeats were large, with mononucleotide repeats being most numerous and hexanucleotide repeats being least numerous. We recommend the use of the ‘2+3’ (Illumina+PacBio) sequencing combination to supplement the Hi-C technique and resequencing technique in future whole-genome research in A. venetum.

Genome survey sequencing and genetic background characterization of yellow horn based on next-generation sequencing

Yellowhorn (Xanthoceras sorbifolium Bunge) is an important wood oil tree species, with high ornamental and medicinal value. Nevertheless, genomic information of yellowhorn is currently unavailable. Here, for the first time, we conducted a genome survey of two yellowhorn cultivars, Zhongshi 4 and Zhongshi 9, which had distinct differences on the phenotype and drought resistance, to obtain knowledge on the genomic information by next generation sequencing (NGS). Meanwhile, its genome size was estimated using flow cytometry. As a result, the whole genome survey of Zhongshi 4 and Zhongshi 9 generated 34.40 and 39.55 GB sequence data. The genome size of Zhongshi 4 and Zhongshi 9 estimated were about 536.58 Mb and 569.52 Mb, which were closed to results of flow cytometry. The heterozygosity rates were calculated to be 0.75% and 0.89%, and the repeat rates were 60.08% and 62.00%. These reads were assembled into 1024,373 and 885,404 contigs with a N50 length of 1005 bp and 1219 bp, respectively, which were further assembled into 714,369 and 686,128 scaffolds with scaffold N50 length of ~ 1963 bp and ~ 1938 bp, total length of 386,915 Kb and 391,904 Kb. These results indicated that there was little difference in genome size and complexity among different cultivars. In addition, 63137 and 65271 high-quality genomic simple sequence repeat (SSR) markers in Zhongshi 4 and Zhongshi 9 were generated. We suggest that the technologies combining Illumina and PacBio, assisted by Hi-C and matching assemble software should be used to one of two yellowhorn cultivars genome sequencing. The result will help to design whole genome sequencing strategies for yellowhorn, and provided a large amount of gene resources for further excavation and utilization of yellowhorn.

High-throughput sequencing data and the impact of plant gene annotation quality

The use of draft genomes of different species and re-sequencing of accessions and populations are now common tools for plant biology research. The de novo assembled draft genomes make it possible to identify pivotal divergence points in the plant lineage and provide an opportunity to investigate the genomic basis and timing of biological innovations by inferring orthologs between species. Furthermore, re-sequencing facilitates the mapping and subsequent molecular characterization of causative loci for traits, such as those for plant stress tolerance and development. In both cases high-quality gene annotation-the identification of protein-coding regions, gene promoters, and 5'- and 3'-untranslated regions-is critical for investigation of gene function. Annotations are constantly improving but automated gene annotations still require manual curation and experimental validation. This is particularly important for genes with large introns, genes located in regions rich with transposable elements or repeats, large gene families, and segmentally duplicated genes. In this opinion paper, we highlight the impact of annotation quality on evolutionary analyses, genome-wide association studies, and the identification of orthologous genes in plants. Furthermore, we predict that incorporating accurate information from manual curation into databases will dramatically improve the performance of automated gene predictors.

Experimental Verification and Evolutionary Origin of 5'-UTR Polyadenylation Sites in Arabidopsis thaliana

Messenger RNA (mRNA) polyadenylation is an indispensable step during post-transcriptional pre-mRNA processing for most genes in eukaryotes. The usage of one poly(A) site over another is known as alternative polyadenylation (APA). APA has been implicated in gene expression regulation through its role of selecting the ends of a transcript. Recent studies of polyadenylation profiles in the Arabidopsis database unexpectedly predicted that a portion of the poly(A) sites are located in the 5'-UTR, which remains to be experimentally verified. We selected 16 genes from a dataset of 744, based on criteria designed to minimize problems in interpretation. Here, we experimentally verify 5'-UTR-APA in Arabidopsis for 10 of the 16 selected genes, and show for the first time existence of independent polyadenylated 5'-UTR transcripts, arising due to alternative polyadenylation. We used 3'-RACE and sequencing to validate poly(A) sites and northern blot to show that the observed short upstream transcripts do not arise from the 3'-end of a previously unrecognized convergent gene. Evidence is reported showing that two of the independent upstream open reading frame (uORF) transcripts studied, one containing a complex dual uORF, very likely arose by exon shuffling following duplication of the 5'-end from the downstream major open reading frame (mORF). Finally, results are presented to show that the uORF in this gene may encode two short functional proteins, based on observation of amino acid sequence conservation encoded by the dual uORFs.

Safety, Security, and Policy Considerations for Plant Genome Editing

Genome editing with engineered nucleases (GEEN) is increasingly used as a tool for gene discovery and trait development in crops through generation of targeted changes in endogenous genes. The development of the CRISPR-Cas9 system (clustered regularly interspaced short palindromic repeats with associated Cas9 protein), in particular, has enabled widespread use of genome editing. Research to date has not comprehensively addressed genome-editing specificity and off-target mismatches that may result in unintended changes within plant genomes or the potential for gene drive initiation. Governance and regulatory considerations for bioengineered crops derived from using GEEN will require greater clarity as to target specificity, the potential for mismatched edits, unanticipated downstream effects of off-target mutations, and assurance that genome reagents do not occur in finished products. Since governance and regulatory decision making involves robust standards of evidence extending from the laboratory to the postcommercial marketplace, developers of genome-edited crops must anticipate significant engagement and investment to address questions of regulators and civil society.

Ploidy-dependent changes in the epigenome of symbiotic cells correlate with specific patterns of gene expression

The formation of symbiotic nodule cells in Medicago truncatula is driven by successive endoreduplication cycles and transcriptional reprogramming in different temporal waves including the activation of more than 600 cysteine-rich NCR genes expressed only in nodules. We show here that the transcriptional waves correlate with growing ploidy levels and have investigated how the epigenome changes during endoreduplication cycles. Differential DNA methylation was found in only a small subset of symbiotic nodule-specific genes, including more than half of the NCR genes, whereas in most genes DNA methylation was unaffected by the ploidy levels and was independent of the genes' active or repressed state. On the other hand, expression of nodule-specific genes correlated with ploidy-dependent opening of the chromatin as well as, in a subset of tested genes, with reduced H3K27me3 levels combined with enhanced H3K9ac levels. Our results suggest that endoreduplication-dependent epigenetic changes contribute to transcriptional reprogramming in the differentiation of symbiotic cells.

Achieving Plant CRISPR Targeting that Limits Off-Target Effects

The CRISPR-Cas9 system (clustered regularly interspaced short palindromic repeats with associated Cas9 protein) has been used to generate targeted changes for direct modification of endogenous genes in an increasing number of plant species; but development of plant genome editing has not yet fully considered potential off-target mismatches that may lead to unintended changes within the genome. Assessing the specificity of CRISPR-Cas9 for increasing editing efficiency as well as the potential for unanticipated downstream effects from off-target mutations is an important regulatory consideration for agricultural applications. Increasing genome-editing specificity entails developing improved design methods that better predict the prevalence of off-target mutations as a function of genome composition and design of the engineered ribonucleoprotein (RNP). Early results from CRISPR-Cas9 genome editing in plant systems indicate that the incidence of off-target mutation frequencies is quite low; however, by analyzing CRISPR-edited plant lines and improving both computational tools and reagent design, it may be possible to further decrease unanticipated effects at potential mismatch sites within the genome. This will provide assurance that CRISPR-Cas9 reagents can be designed and targeted with a high degree of specificity. Improved and experimentally validated design tools for discriminating target and potential off-target positions that incorporate consideration of the designed nuclease fidelity and selectivity will help to increase confidence for regulatory decision making for genome-edited plants.

Dissecting the U, M, S and C genomes of wild relatives of bread wheat (Aegilops spp.) into chromosomes and exploring their synteny with wheat

Goat grasses (Aegilops spp.) contributed to the evolution of bread wheat and are important sources of genes and alleles for modern wheat improvement. However, their use in alien introgression breeding is hindered by poor knowledge of their genome structure and a lack of molecular tools. The analysis of large and complex genomes may be simplified by dissecting them into single chromosomes via flow cytometric sorting. In some species this is not possible due to similarities in relative DNA content among chromosomes within a karyotype. This work describes the distribution of GAA and ACG microsatellite repeats on chromosomes of the U, M, S and C genomes of Aegilops, and the use of microsatellite probes to label the chromosomes in suspension by fluorescence in situ hybridization (FISHIS). Bivariate flow cytometric analysis of chromosome DAPI fluorescence and fluorescence of FITC-labelled microsatellites made it possible to discriminate all chromosomes and sort them with negligible contamination by other chromosomes. DNA of purified chromosomes was used as a template for polymerase chain reation (PCR) using Conserved Orthologous Set (COS) markers with known positions on wheat A, B and D genomes. Wheat-Aegilops macrosyntenic comparisons using COS markers revealed significant rearrangements in the U and C genomes, while the M and S genomes exhibited structure similar to wheat. Purified chromosome fractions provided an attractive resource to investigate the structure and evolution of the Aegilops genomes, and the COS markers assigned to Aegilops chromosomes will facilitate alien gene introgression into wheat.

Orphan Crops Browser: a bridge between model and orphan crops

Many important crops have received little attention by the scientific community, either because they are not considered economically important or due to their large and complex genomes. De novo transcriptome assembly, using next-generation sequencing data, is an attractive option for the study of these orphan crops. In spite of the large amount of sequencing data that can be generated, there is currently a lack of tools which can effectively help molecular breeders and biologists to mine this type of information. Our goal was to develop a tool that enables molecular breeders, without extensive bioinformatics knowledge, to efficiently study de novo transcriptome data from any orphan crop (http://www.bioinformatics.nl/denovobrowser/db/species/index). The Orphan Crops Browser has been designed to facilitate the following tasks (1) search and identification of candidate transcripts based on phylogenetic relationships between orthologous sequence data from a set of related species and (2) design specific and degenerate primers for expression studies in the orphan crop of interest. To demonstrate the usability and reliability of the browser, it was used to identify the putative orthologues of 17 known lignin biosynthetic genes from maize and sugarcane in the orphan crop Miscanthus sinensis. Expression studies in miscanthus stem internode tissue differing in maturation were subsequently carried out, to follow the expression of these genes during lignification. Our results showed a negative correlation between lignin content and gene expression. The present data are in agreement with recent findings in maize and other crops, and it is further discussed in this paper.

Expressed Sequence Tags Analysis and Design of Simple Sequence Repeats Markers from a Full-Length cDNA Library in Perilla frutescens (L.)

Perilla frutescens is valuable as a medicinal plant as well as a natural medicine and functional food. However, comparative genomics analyses of P. frutescens are limited due to a lack of gene annotations and characterization. A full-length cDNA library from P. frutescens leaves was constructed to identify functional gene clusters and probable EST-SSR markers via analysis of 1,056 expressed sequence tags. Unigene assembly was performed using basic local alignment search tool (BLAST) homology searches and annotated Gene Ontology (GO). A total of 18 simple sequence repeats (SSRs) were designed as primer pairs. This study is the first to report comparative genomics and EST-SSR markers from P. frutescens will help gene discovery and provide an important source for functional genomics and molecular genetic research in this interesting medicinal plant.

Sugarcane genome sequencing by methylation filtration provides tools for genomic research in the genus Saccharum

Many economically important crops have large and complex genomes that hamper their sequencing by standard methods such as whole genome shotgun (WGS). Large tracts of methylated repeats occur in plant genomes that are interspersed by hypomethylated gene-rich regions. Gene-enrichment strategies based on methylation profiles offer an alternative to sequencing repetitive genomes. Here, we have applied methyl filtration with McrBC endonuclease digestion to enrich for euchromatic regions in the sugarcane genome. To verify the efficiency of methylation filtration and the assembly quality of sequences submitted to gene-enrichment strategy, we have compared assemblies using methyl-filtered (MF) and unfiltered (UF) libraries. The use of methy filtration allowed a better assembly by filtering out 35% of the sugarcane genome and by producing 1.5× more scaffolds and 1.7× more assembled Mb in length compared with unfiltered dataset. The coverage of sorghum coding sequences (CDS) by MF scaffolds was at least 36% higher than by the use of UF scaffolds. Using MF technology, we increased by 134× the coverage of gene regions of the monoploid sugarcane genome. The MF reads assembled into scaffolds that covered all genes of the sugarcane bacterial artificial chromosomes (BACs), 97.2% of sugarcane expressed sequence tags (ESTs), 92.7% of sugarcane RNA-seq reads and 98.4% of sorghum protein sequences. Analysis of MF scaffolds from encoded enzymes of the sucrose/starch pathway discovered 291 single-nucleotide polymorphisms (SNPs) in the wild sugarcane species, S. spontaneum and S. officinarum. A large number of microRNA genes was also identified in the MF scaffolds. The information achieved by the MF dataset provides a valuable tool for genomic research in the genus Saccharum and for improvement of sugarcane as a biofuel crop.

Flow cytometric chromosome sorting from diploid progenitors of bread wheat, T. urartu, Ae. speltoides and Ae. tauschii

Chromosomes 5A (u) , 5S and 5D can be isolated from wild progenitors, providing a chromosome-based approach to develop tools for breeding and to study the genome evolution of wheat. The three subgenomes of hexaploid bread wheat originated from Triticum urartu (A(u)A(u)), from a species similar to Aegilops speltoides (SS) (progenitor of the B genome), and from Ae. tauschii (DD). Earlier studies indicated the potential of chromosome genomics to assist gene transfer from wild relatives of wheat and discover novel genes for wheat improvement. This study evaluates the potential of flow cytometric chromosome sorting in the diploid progenitors of bread wheat. Flow karyotypes obtained by analysing DAPI-stained chromosomes were characterized and the contents of the chromosome peaks were determined. FISH analysis with repetitive DNA probes proved that chromosomes 5A(u), 5S and 5D could be sorted with purities of 78-90 %, while the remaining chromosomes could be sorted in groups of three. Twenty-five conserved orthologous set (COS) markers covering wheat homoeologous chromosome groups 1-7 were used for PCR with DNA amplified from flow-sorted chromosomes and genomic DNA. These assays validated the cytomolecular results as follows: peak I on flow karyotypes contained chromosome groups 1, 4 and 6, peak II represented homoeologous group 5, while peak III consisted of groups 2, 3 and 7. The isolation of individual chromosomes of wild progenitors provides an attractive opportunity to investigate the structure and evolution of the polyploid genome and to deliver tools for wheat improvement.

Natural product biosynthesis in Medicago species

The genus Medicago, a member of the legume (Fabaceae) family, comprises 87 species of flowering plants, including the forage crop M. sativa (alfalfa) and the model legume M. truncatula (barrel medic). Medicago species synthesize a variety of bioactive natural products that are used to engage into symbiotic interactions but also serve to deter pathogens and herbivores. For humans, these bioactive natural products often possess promising pharmaceutical properties. In this review, we focus on the two most interesting and well characterized secondary metabolite classes found in Medicago species, the triterpene saponins and the flavonoids, with a detailed overview of their biosynthesis, regulation, and profiling methods. Furthermore, their biological role within the plant as well as their potential utility for human health or other applications is discussed. Finally, we give an overview of the advances made in metabolic engineering in Medicago species and how the development of novel molecular and omics toolkits can influence a better understanding of this genus in terms of specialized metabolism and chemistry. Throughout, we critically analyze the current bottlenecks and speculate on future directions and opportunities for research and exploitation of Medicago metabolism.