The International Oryza Map Alignment Project (IOMAP): the Americas-past achievements and future directions.

The wild relatives of rice hold unexplored genetic diversity that can be employed to feed an estimated population of 10 billion by 2050. The Oryza Map Alignment Project (OMAP) initiated in 2003 has provided comprehensive genomic resources for comparative, evolutionary, and functional characterization of the wild relatives of rice, facilitating the cloning of >600 rice genes, including those for grain width (GW5) and submergence tolerance (SUB1A). Following in the footsteps of the original project, the goal of 'IOMAP: the Americas' is to investigate the present and historic genetic diversity of wild Oryza species endemic to the Americas through the sequencing of herbaria and in situ specimens. The generation of a large diversity panel describing past and current genetic status and potential erosion of genetic variation in the populations will provide useful knowledge for the conservation of the biodiversity in these species. The wild relatives of rice in the Americas present a wide range of resistance traits useful for crop improvement and neodomestication approaches. In the race against time for a sustainable food future, the neodomestication of the first cereal species recently accomplished in O. alta opens the door to the potential neodomestication of the other wild Oryza species in Americas.

A new genome allows the identification of genes associated with natural variation in aluminium tolerance in Brachiaria grasses.

Toxic concentrations of aluminium cations and low phosphorus availability are the main yield-limiting factors in acidic soils, which represent half of the potentially available arable land. Brachiaria grasses, which are commonly sown as forage in the tropics because of their resilience and low demand for nutrients, show greater tolerance to high concentrations of aluminium cations (Al3+) than most other grass crops. In this work, we explored the natural variation in tolerance to Al3+ between high and low tolerant Brachiaria species and characterized their transcriptional differences during stress. We identified three QTLs (quantitative trait loci) associated with root vigour during Al3+ stress in their hybrid progeny. By integrating these results with a new Brachiaria reference genome, we identified 30 genes putatively responsible for Al3+ tolerance in Brachiaria. We observed differential expression during stress of genes involved in RNA translation, response signalling, cell wall composition, and vesicle location homologous to aluminium-induced proteins involved in limiting uptake or localizing the toxin. However, there was limited regulation of malate transporters in Brachiaria, which suggests that exudation of organic acids and other external tolerance mechanisms, common in other grasses, might not be relevant in Brachiaria. The contrasting regulation of RNA translation and response signalling suggests that response timing is critical in high Al3+-tolerant Brachiaria.

Translocation of a parthenogenesis gene candidate to an alternate carrier chromosome in apomictic Brachiaria humidicola.

BACKGROUND: The apomictic reproductive mode of Brachiaria (syn. Urochloa) forage species allows breeders to faithfully propagate heterozygous genotypes through seed over multiple generations. In Brachiaria, reproductive mode segregates as single dominant locus, the apospory-specific genomic region (ASGR). The AGSR has been mapped to an area of reduced recombination on Brachiaria decumbens chromosome 5. A primer pair designed within ASGR-BABY BOOM-like (BBML), the candidate gene for the parthenogenesis component of apomixis in Pennisetum squamulatum, was diagnostic for reproductive mode in the closely related species B. ruziziensis, B. brizantha, and B. decumbens. In this study, we used a mapping population of the distantly related commercial species B. humidicola to map the ASGR and test for conservation of ASGR-BBML sequences across Brachiaria species. RESULTS: Dense genetic maps were constructed for the maternal and paternal genomes of a hexaploid (2n = 6x = 36) B. humidicola F1 mapping population (n = 102) using genotyping-by-sequencing, simple sequence repeat, amplified fragment length polymorphism, and transcriptome derived single nucleotide polymorphism markers. Comparative genomics with Setaria italica provided confirmation for x = 6 as the base chromosome number of B. humidicola. High resolution molecular karyotyping indicated that the six homologous chromosomes of the sexual female parent paired at random, whereas preferential pairing of subgenomes was observed in the apomictic male parent. Furthermore, evidence for compensated aneuploidy was found in the apomictic parent, with only five homologous linkage groups identified for chromosome 5 and seven homologous linkage groups of chromosome 6. The ASGR mapped to B. humidicola chromosome 1, a region syntenic with chromosomes 1 and 7 of S. italica. The ASGR-BBML specific PCR product cosegregated with the ASGR in the F1 mapping population, despite its location on a different carrier chromosome than B. decumbens. CONCLUSIONS: The first dense molecular maps of B. humidicola provide strong support for cytogenetic evidence indicating a base chromosome number of six in this species. Furthermore, these results show conservation of the ASGR across the Paniceae in different chromosomal backgrounds and support postulation of the ASGR-BBML as candidate genes for the parthenogenesis component of apomixis.

Two Novel DNAs That Enhance Symptoms and Overcome CMD2 Resistance to Cassava Mosaic Disease.

UNLABELLED: Cassava mosaic begomoviruses (CMBs) cause cassava mosaic disease (CMD) across Africa and the Indian subcontinent. Like all members of the geminivirus family, CMBs have small, circular single-stranded DNA genomes. We report here the discovery of two novel DNA sequences, designated SEGS-1 and SEGS-2 (forsequencesenhancinggeminivirussymptoms), that enhance symptoms and break resistance to CMD. The SEGS are characterized by GC-rich regions and the absence of long open reading frames. Both SEGS enhanced CMD symptoms in cassava (Manihot esculentaCrantz) when coinoculated withAfrican cassava mosaic virus(ACMV),East African cassava mosaic Cameroon virus(EACMCV), orEast African cassava mosaic virus-Uganda(EACMV-UG). SEGS-1 also overcame resistance of a cassava landrace carrying the CMD2 resistance locus when coinoculated with EACMV-UG. Episomal forms of both SEGS were detected in CMB-infected cassava but not in healthy cassava. SEGS-2 episomes were also found in virions and whiteflies. SEGS-1 has no homology to geminiviruses or their associated satellites, but the cassava genome contains a sequence that is 99% identical to full-length SEGS-1. The cassava genome also includes three sequences with 84 to 89% identity to SEGS-2 that together encompass all of SEGS-2 except for a 52-bp region, which includes the episomal junction and a 26-bp sequence related to alphasatellite replication origins. These results suggest that SEGS-1 is derived from the cassava genome and facilitates CMB infection as an integrated copy and/or an episome, while SEGS-2 was originally from the cassava genome but now is encapsidated into virions and transmitted as an episome by whiteflies. IMPORTANCE: Cassava is a major crop in the developing world, with its production in Africa being second only to maize. CMD is one of the most important diseases of cassava and a serious constraint to production across Africa. CMD2 is a major CMD resistance locus that has been deployed in many cassava cultivars through large-scale breeding programs. In recent years, severe, atypical CMD symptoms have been observed occasionally on resistant cultivars, some of which carry the CMD2 locus, in African fields. In this report, we identified and characterized two DNA sequences, SEGS-1 and SEGS-2, which produce similar symptoms when coinoculated with cassava mosaic begomoviruses onto a susceptible cultivar or a CMD2-resistant landrace. The ability of SEGS-1 to overcome CMD2 resistance and the transmission of SEGS-2 by whiteflies has major implications for the long-term durability of CMD2 resistance and underscore the need for alternative sources of resistance in cassava.

Molecular analysis of the expression of a crtB transgene and the endogenous psy2-y 1 and psy2-y 2 genes of cassava and their effect on root carotenoid content.

A conventional breeding program was established to transfer the bacterial phytoene synthase transgene-crtB-from a transgenic, white-rooted cassava to yellow-rooted cassava plants carrying the endogenous phytoene synthase alleles named psy2-y 1 and/or psy2-y 2. Combining endogenous phytoene synthase enzymes (PSYs) with CRTB in a single cassava plant would allow the molecular dissection of individual allele contributions to carotenoid synthesis and/or accumulation in cassava roots. The simultaneous expression of the crtB transgene and psy2-y 2 in individuals planted in the field coincided with higher total, HPLC-quantified carotenoid content in roots, although the variability among replications (plants) precluded the detection of statistically significant differences. Nevertheless, the highest total carotenoid content in roots within a family coincided with one individual of the F1 progeny carrying both psy2-y 2 and crtB genes. The results also indicated the presence of at least one more key gene-different from psy or crtB-which too is necessary for the synthesis and/or accumulation of Pro-Vitamin A carotenoids in cassava roots.

An integrated framework for discovery and genotyping of genomic variants from high-throughput sequencing experiments.

Recent advances in high-throughput sequencing (HTS) technologies and computing capacity have produced unprecedented amounts of genomic data that have unraveled the genetics of phenotypic variability in several species. However, operating and integrating current software tools for data analysis still require important investments in highly skilled personnel. Developing accurate, efficient and user-friendly software packages for HTS data analysis will lead to a more rapid discovery of genomic elements relevant to medical, agricultural and industrial applications. We therefore developed Next-Generation Sequencing Eclipse Plug-in (NGSEP), a new software tool for integrated, efficient and user-friendly detection of single nucleotide variants (SNVs), indels and copy number variants (CNVs). NGSEP includes modules for read alignment, sorting, merging, functional annotation of variants, filtering and quality statistics. Analysis of sequencing experiments in yeast, rice and human samples shows that NGSEP has superior accuracy and efficiency, compared with currently available packages for variants detection. We also show that only a comprehensive and accurate identification of repeat regions and CNVs allows researchers to properly separate SNVs from differences between copies of repeat elements. We expect that NGSEP will become a strong support tool to empower the analysis of sequencing data in a wide range of research projects on different species.

Flexible and scalable genotyping-by-sequencing strategies for population studies.

BACKGROUND: Many areas critical to agricultural production and research, such as the breeding and trait mapping in plants and livestock, require robust and scalable genotyping platforms. Genotyping-by-sequencing (GBS) is a one such method highly suited to non-human organisms. In the GBS protocol, genomic DNA is fractionated via restriction digest, then reduced representation is achieved through size selection. Since many restriction sites are conserved across a species, the sequenced portion of the genome is highly consistent within a population. This makes the GBS protocol highly suited for experiments that require surveying large numbers of markers within a population, such as those involving genetic mapping, breeding, and population genomics. We have modified the GBS technology in a number of ways. Custom, enzyme specific adaptors have been replaced with standard Illumina adaptors compatible with blunt-end restriction enzymes. Multiplexing is achieved through a dual barcoding system, and bead-based library preparation protocols allows for in-solution size selection and eliminates the need for columns and gels. RESULTS: A panel of eight restriction enzymes was selected for testing on B73 maize and Nipponbare rice genomic DNA. Quality of the data was demonstrated by identifying that the vast majority of reads from each enzyme aligned to restriction sites predicted in silico. The link between enzyme parameters and experimental outcome was demonstrated by showing that the sequenced portion of the genome was adaptable by selecting enzymes based on motif length, complexity, and methylation sensitivity. The utility of the new GBS protocol was demonstrated by correctly mapping several in a maize F2 population resulting from a B73×Country Gentleman test cross. CONCLUSIONS: This technology is readily adaptable to different genomes, highly amenable to multiplexing and compatible with over forty commercially available restriction enzymes. These advancements represent a major improvement in genotyping technology by providing a highly flexible and scalable GBS that is readily implemented for studies on genome-wide variation.

Provitamin A accumulation in cassava (Manihot esculenta) roots driven by a single nucleotide polymorphism in a phytoene synthase gene.

Cassava (Manihot esculenta) is an important staple crop, especially in the arid tropics. Because roots of commercial cassava cultivars contain a limited amount of provitamin A carotenoids, both conventional breeding and genetic modification are being applied to increase their production and accumulation to fight vitamin A deficiency disorders. We show here that an allelic polymorphism in one of the two expressed phytoene synthase (PSY) genes is capable of enhancing the flux of carbon through carotenogenesis, thus leading to the accumulation of colored provitamin A carotenoids in storage roots. A single nucleotide polymorphism present only in yellow-rooted cultivars cosegregates with colored roots in a breeding pedigree. The resulting amino acid exchange in a highly conserved region of PSY provides increased catalytic activity in vitro and is able to increase carotenoid production in recombinant yeast and Escherichia coli cells. Consequently, cassava plants overexpressing a PSY transgene produce yellow-fleshed, high-carotenoid roots. This newly characterized PSY allele provides means to improve cassava provitamin A content in cassava roots through both breeding and genetic modification.

Genetic bases of resistance to the rice hoja blanca disease deciphered by a quantitative trait locus approach.

Rice hoja blanca (RHB) is one of the most serious diseases in rice-growing areas in tropical Americas. Its causal agent is RHB virus (RHBV), transmitted by the planthopper Tagosodes orizicolus Müir. Genetic resistance is the most effective and environment-friendly way of controlling the disease. So far, only 1 major quantitative trait locus (QTL) of Oryza sativa ssp. japonica origin, qHBV4.1, that alters the incidence of the virus symptoms in 2 Colombian cultivars has been reported. This resistance has already started to be broken, stressing the urgent need for diversifying the resistance sources. In the present study, we performed a search for new QTLs of O. sativa indica origin associated with RHB resistance. We used 4 F2:3-segregating populations derived from indica-resistant varieties crossed with a highly susceptible japonica pivot parent. Besides the standard method for measuring disease incidence, we developed a new method based on computer-assisted image processing to determine the affected leaf area (ALA) as a measure of symptom severity. Based on the disease severity and incidence scores in the F3 families under greenhouse conditions and SNP genotyping of the F2 individuals, we identified 4 new indica QTLs for RHB resistance on rice chromosomes 4, 6, and 11, namely, qHBV4.2WAS208, qHBV6.1PTB25, qHBV11.1, and qHBV11.2, respectively. We also confirmed the wide-range action of qHBV4.1. Among the 5 QTLs, qHBV4.1 and qHBV11.1 had the largest effects on incidence and severity, respectively. These results provide a more complete understanding of the genetic bases of RHBV resistance in the cultivated rice gene pool and can be used to develop marker-aided breeding strategies to improve RHB resistance. The power of joint- and meta-analyses allowed precise mapping and candidate gene identification, providing the basis for positional cloning of the 2 major QTLs qHBV4.1 and qHBV11.1.

In-depth molecular and phenotypic characterization in a rice insertion line library facilitates gene identification through reverse and forward genetics approaches.

We report here the molecular and phenotypic features of a library of 31,562 insertion lines generated in the model japonica cultivar Nipponbare of rice (Oryza sativa L.), called Oryza Tag Line (OTL). Sixteen thousand eight hundred and fourteen T-DNA and 12,410 Tos17 discrete insertion sites have been characterized in these lines. We estimate that 8686 predicted gene intervals--i.e. one-fourth to one-fifth of the estimated rice nontransposable element gene complement--are interrupted by sequence-indexed T-DNA (6563 genes) and/or Tos17 (2755 genes) inserts. Six hundred and forty-three genes are interrupted by both T-DNA and Tos17 inserts. High quality of the sequence indexation of the T2 seed samples was ascertained by several approaches. Field evaluation under agronomic conditions of 27,832 OTL has revealed that 18.2% exhibit at least one morphophysiological alteration in the T1 progeny plants. Screening 10,000 lines for altered response to inoculation by the fungal pathogen Magnaporthe oryzae allowed to observe 71 lines (0.7%) developing spontaneous lesions simulating disease mutants and 43 lines (0.4%) exhibiting an enhanced disease resistance or susceptibility. We show here that at least 3.5% (four of 114) of these alterations are tagged by the mutagens. The presence of allelic series of sequence-indexed mutations in a gene among OTL that exhibit a convergent phenotype clearly increases the chance of establishing a linkage between alterations and inserts. This convergence approach is illustrated by the identification of the rice ortholog of AtPHO2, the disruption of which causes a lesion-mimic phenotype owing to an over-accumulation of phosphate, in nine lines bearing allelic insertions.

Foodomics: A Data-Driven Approach to Revolutionize Nutrition and Sustainable Diets.

Globally, we are failing to meet numerous nutritional, health, and environmental targets linked to food. Defining food composition in its full chemical and quantitative diversity is central to data-driven decision making for supporting nutrition and sustainable diets. "Foodomics"-the application of omics-technology to characterize and quantify biomolecules to improve wellbeing-has the potential to comprehensively elucidate what is in food, how this composition varies across the food system, and how diet composition as an ensemble of foods guides outcomes for nutrition, health, and sustainability. Here, we outline: (i) challenges of evaluating food composition; (ii) state-of-the-art omics technology and innovations for the analysis of food; and (iii) application of foodomics as a complementary data-driven approach to revolutionize nutrition and sustainable diets. Featuring efforts of the Periodic Table of Food Initiative, a participatory effort to create a globally shared foodomics platform, we conclude with recommendations to accelerate foodomics in ways that strengthen the capacity of scientists and benefit all people.

A new method for designing degenerate primers and its use in the identification of sequences in Brachiaria showing similarity to apomixis-associated genes.

MOTIVATION: We developed a technique and a tool for degenerate primer design based on multiple local alignments employing the MEME algorithm supported with electronic PCR. The objective is to find adequate primers starting from sequences with poor global similarity. We show an example of its application in our laboratory to find sequences in Brachiaria with similarity to ESTs related to apomixis.

Identification of a Rice stripe necrosis virus resistance locus and yield component QTLs using Oryza sativa x O. glaberrima introgression lines.

BACKGROUND: Developing new population types based on interspecific introgressions has been suggested by several authors to facilitate the discovery of novel allelic sources for traits of agronomic importance. Chromosome segment substitution lines from interspecific crosses represent a powerful and useful genetic resource for QTL detection and breeding programs. RESULTS: We built a set of 64 chromosome segment substitution lines carrying contiguous chromosomal segments of African rice Oryza glaberrima MG12 (acc. IRGC103544) in the genetic background of Oryza sativa ssp. tropical japonica (cv. Caiapó). Well-distributed simple-sequence repeats markers were used to characterize the introgression events. Average size of the substituted chromosomal segments in the substitution lines was about 10 cM and covered the whole donor genome, except for small regions on chromosome 2 and 4. Proportions of recurrent and donor genome in the substitution lines were 87.59% and 7.64%, respectively. The remaining 4.78% corresponded to heterozygotes and missing data. Strong segregation distortion was found on chromosomes 3 and 6, indicating the presence of interspecific sterility genes. To illustrate the advantages and the power of quantitative trait loci (QTL) detection using substitution lines, a QTL detection was performed for scored traits. Transgressive segregation was observed for several traits measured in the population. Fourteen QTLs for plant height, tiller number per plant, panicle length, sterility percentage, 1000-grain weight and grain yield were located on chromosomes 1, 3, 4, 6 and 9. Furthermore, a highly significant QTL controlling resistance to the Rice stripe necrosis virus was located between SSR markers RM202-RM26406 (44.5-44.8 cM) on chromosome 11. CONCLUSIONS: Development and phenotyping of CSSL libraries with entire genome coverage represents a useful strategy for QTL discovery. Mapping of the RSNV locus represents the first identification of a genetic factor underlying resistance to this virus. This population is a powerful breeding tool. It also helps in overcoming hybrid sterility barriers between species of rice.

In planta gene expression analysis of Xanthomonas oryzae pathovar oryzae, African strain MAI1.

BACKGROUND: Bacterial leaf blight causes significant yield losses in rice crops throughout Asia and Africa. Although both the Asian and African strains of the pathogen, Xanthomonas oryzae pv. oryzae (Xoo), induce similar symptoms, they are nevertheless genetically different, with the African strains being more closely related to the Asian X. oryzae pv. oryzicola (Xoc). RESULTS: Changes in gene expression of the African Xoo strain MAI1 in the susceptible rice cultivar Nipponbare were profiled, using an SSH Xoo DNA microarray. Microarray hybridization was performed comparing bacteria recovered from plant tissues at 1, 3, and 6 days after inoculation (dai) with bacteria grown in vitro. A total of 710 bacterial genes were found to be differentially expressed, with 407 up-regulated and 303 down-regulated. Expression profiling indicated that less than 20% of the 710 bacterial transcripts were induced in the first 24 h after inoculation, whereas 63% were differentially expressed at 6 dai. The 710 differentially expressed genes were one-end sequenced. 535 sequences were obtained from which 147 non-redundant sequences were identified. Differentially expressed genes were related to metabolism, secretion and transport, pathogen adherence to plant tissues, plant cell-wall degradation, IS elements, and virulence. In addition, various other genes encoding proteins with unknown function or showing no similarity to other proteins were also induced. The Xoo MAI1 non-redundant set of sequences was compared against several X. oryzae genomes, revealing a specific group of genes that was present only in MAI1. Numerous IS elements were also found to be differentially expressed. Quantitative real-time PCR confirmed 86% of the identified profile on a set of 14 genes selected according to the microarray analysis. CONCLUSIONS: This is the first report to compare the expression of Xoo genes in planta across different time points during infection. This work shows that as-yet-unidentified and potentially new virulence factors are appearing in an emerging African pathogen. It also confirms that African Xoo strains do differ from their Asian counterparts, even at the transcriptional level.

A universal core genetic map for rice.

To facilitate the creation of easily comparable, low-resolution genetic maps with evenly distributed markers in rice (Oryza sativa L.), we conceived of and developed a Universal Core Genetic Map (UCGM). With this aim, we derived a set of 165 anchors, representing clusters of three microsatellite or simple sequence repeat (SSR) markers arranged into non-recombining groups. Each anchor consists of at least three, closely linked SSRs, located within a distance below the genetic resolution provided by common, segregating populations (<500 individuals). We chose anchors that were evenly distributed across the rice chromosomes, with spacing between 2 and 3.5 Mbp (except in the telomeric regions, where spacing was 1.5 Mbp). Anchor selection was performed using in silico tools and data: the O. sativa cv. Nipponbare rice genome sequence, the CHARM tool, information from the Gramene database and the OrygenesDB database. Sixteen AA-genome accessions of the Oryza genus were used to evaluate polymorphisms for the selected markers, including accessions from O. sativa, O. glaberrima, O. barthii, O. rufipogon, O. glumaepatula and O. meridionalis. High levels of polymorphism were found for the tested O. sativa x O. glaberrima or O. sativa x wild rice combinations. We developed Paddy Map, a simple database that is helpful in selecting optimal sets of polymorphic SSRs for any cross that involves the previously mentioned species. Validation of the UCGM was done by using it to develop three interspecific genetic maps and by comparing genetic SSR locations with their physical positions on the rice pseudomolecules. In this study, we demonstrate that the UCGM is a useful tool for the rice genetics and breeding community, especially in strategies based on interspecific hybridisation.

Efficient Agrobacterium-mediated Transformation of The Elite-Indica Rice Variety Komboka.

Genetic transformation is crucial for both investigating gene functions and for engineering of crops to introduce new traits. Rice (Oryza sativa L.) is an important model in plant research, since it is the staple food for more than half of the world's population. As a result, numerous transformation methods have been developed for both indica and japonica rice. Since breeders continuously develop new rice varieties, transformation protocols have to be adapted for each new variety. Here we provide an optimized transformation protocol with detailed tips and tricks for a new African variety Komboka using immature embryos. In Komboka, we obtained an apparent transformation rate of up to 48% for GUS/GFP reporter gene constructs using this optimized protocol. This protocol is also applicable for use with other elite indica rice varieties.

Multiplying the efficiency and impact of biofortification through metabolic engineering.

Ending all forms of hunger by 2030, as set forward in the UN-Sustainable Development Goal 2 (UN-SDG2), is a daunting but essential task, given the limited timeline ahead and the negative global health and socio-economic impact of hunger. Malnutrition or hidden hunger due to micronutrient deficiencies affects about one third of the world population and severely jeopardizes economic development. Staple crop biofortification through gene stacking, using a rational combination of conventional breeding and metabolic engineering strategies, should enable a leap forward within the coming decade. A number of specific actions and policy interventions are proposed to reach this goal.

Arabidopsis LEC1 and LEC2 Orthologous Genes Are Key Regulators of Somatic Embryogenesis in Cassava.

High genotype-dependent variation in friable embryogenic callus (FEC) induction and subsequent somaclonal variation constitute bottlenecks for the application and scaling of genetic transformation (GT) technology to more farmer- and industry-preferred cassava varieties. The understanding and identification of molecular factors underlying embryogenic development in cassava may help to overcome these constraints. Here, we described the Arabidopsis thaliana LEAFY COTYLEDON (LEC) LEC1 and LEC2 orthologous genes in cassava, designated as MeLEC1 and MeLEC2, respectively. Expression analyses showed that both, MeLEC1 and MeLEC2, are expressed at higher levels in somatic embryogenic (SE) tissues in contrast with differentiated mature tissues. The rapid expression increase of MeLEC genes at early SE induction times strongly suggests that they are involved in the transition from a somatic to an embryonic state, and probably, in the competence acquisition for SE development in cassava. The independent overexpression of the MeLEC genes resulted in different regenerated events with embryogenic characteristics such as MeLEC1OE plants with cotyledon-like leaves and MeLEC2OE plants with somatic-like embryos that emerged over the surface of mature leaves. Transcript increases of other embryo-specific regulating factors were also detected in MeLECOE plants, supporting their mutual interaction in the embryo development coordination. The single overexpression of MeLEC2 was enough to reprogram the vegetative cells and induce direct somatic embryogenesis, which converts this gene into a tool that could improve the recovery of transformed plants of recalcitrant genotypes. The identification of MeLEC genes contributes not only to improve our understanding of SE process in cassava, but also provides viable alternatives to optimize GT and advance in gene editing in this crop, through the development of genotype-independent protocols.