Genome-wide association analyses using multi-locus models on bananas (Musa spp.) reveal candidate genes related to morphology, fruit quality, and yield.

Bananas (Musa spp.) are an essential fruit worldwide and rank as the fourth most significant food crop for addressing malnutrition due to their rich nutrients and starch content. The potential of their genetic diversity remains untapped due to limited molecular breeding tools. Our study examined a phenotypically diverse group of 124 accessions from the Colombian Musaceae Collection (CMC) conserved in AGROSAVIA. We assessed 12 traits categorized into morphology, fruit quality, and yield, alongside sequence data. Our sequencing efforts provided valuable insights, with an average depth of about 7X per accession, resulting in 187,133 SNPs against Musa acuminata (A genome) and 220,451 against Musa balbisiana (B genome). Population structure analysis grouped samples into four and five clusters based on the reference genome. By using different association models, we identified marker-trait associations (MTAs). The mixed linear model (MLM) revealed four MTAs, while the BLINK and FarmCPU models identified 82 and 70 MTAs, respectively. We identified 38 and 40 candidate genes in linkage proximity to significant MTAs for the A-genome and B-genome, respectively. Our findings provide insights into the genetic underpinnings of morphology, fruit quality, and yield. Once validated, the SNP markers and candidate genes can potentially drive advancements in genomic-guided breeding strategies to enhance banana crop improvement.

Characterizing subgenome recombination and chromosomal imbalances in banana varietal lineages.

BACKGROUND: Bananas and plantains (Musa spp.) are among the most important crops worldwide. The cultivated varieties are vegetatively propagated, so their genetic diversity is essentially fixed over time. Musa acuminata, M. balbisiana and M. schizocarpa have provided the named A, B and S subgenomes that predominantly constitute these varieties. Here we aimed to characterize intergenetic recombination and chromosomal imbalances between these A/B/S subgenomes, which often result in copy-number variants (CNVs) leading to changes in gene dosage and phenotype, in a diverse panel of bananas and plantains. This will allow us to characterize varietal lineages better and identify sources of genetic variation. METHODS: We delimited population structure and clonal lineages in a diverse panel of 188 banana and plantain accessions from the most common cultivars using admixture, principal component and phylogenetic analyses. We used new scalable alignment-based methods, Relative Averaged Alignment (RAA) and Relative Coverage, to infer subgenome composition (AA, AAB, etc.) and interspecific recombination. RESULTS: In our panel, we identified ten varietal lineages composed of somatic clones, plus three groups of tetraploid accessions. We identified chromosomal exchanges resulting in gains/losses in chromosomal segments (CNVs), particularly in AAB and ABB varieties. CONCLUSIONS: We demonstrated alignment-based RAA and Relative Coverage can identify subgenome composition and introgressions with similar results to more complex approaches based on single nucleotide polymorphism (SNP) databases. These ab initio species-agnostic methods can be used without sequencing a panel of wild ancestors to find private SNPs, or in recently diverged pools where private SNPs are uncommon. The extensive A/B/S exchanges and the variation in the length of some introgressions between lineages further support multiple foundational events of hybridization and residual backcrossing. Imbalances between A/B/S may have resulted in CNVs and gene dosage variation. Since most edible banana genomes are fixed on time, these CNVs are stable genetic variations probably associated with phenotypic variation for future genetic studies.

Complex polyploid and hybrid species in an apomictic and sexual tropical forage grass group: genomic composition and evolution in Urochloa (Brachiaria) species.

BACKGROUND AND AIMS: Diploid and polyploid Urochloa (including Brachiaria, Panicum and Megathyrsus species) C4 tropical forage grasses originating from Africa are important for food security and the environment, often being planted in marginal lands worldwide. We aimed to characterize the nature of their genomes, the repetitive DNA and the genome composition of polyploids, leading to a model of the evolutionary pathways within the group including many apomictic species. METHODS: Some 362 forage grass accessions from international germplasm collections were studied, and ploidy was determined using an optimized flow cytometry method. Whole-genome survey sequencing and molecular cytogenetic analysis were used to identify chromosomes and genomes in Urochloa accessions belonging to the 'brizantha' and 'humidicola' agamic complexes and U. maxima. KEY RESULTS: Genome structures are complex and variable, with multiple ploidies and genome compositions within the species, and no clear geographical patterns. Sequence analysis of nine diploid and polyploid accessions enabled identification of abundant genome-specific repetitive DNA motifs. In situ hybridization with a combination of repetitive DNA and genomic DNA probes identified evolutionary divergence and allowed us to discriminate the different genomes present in polyploids. CONCLUSIONS: We suggest a new coherent nomenclature for the genomes present. We develop a model of evolution at the whole-genome level in diploid and polyploid accessions showing processes of grass evolution. We support the retention of narrow species concepts for Urochloa brizantha, U. decumbens and U. ruziziensis, and do not consider diploids and polyploids of single species as cytotypes. The results and model will be valuable in making rational choices of parents for new hybrids, assist in use of the germplasm for breeding and selection of Urochloa with improved sustainability and agronomic potential, and assist in measuring and conserving biodiversity in grasslands.

Genomic regions and candidate genes selected during the breeding of rice in Vietnam.

Vietnam harnesses a rich diversity of rice landraces adapted to a range of conditions, which constitute a largely untapped source of diversity for the continuous improvement of cultivars. We previously identified a strong population structure in Vietnamese rice, which is captured in five Indica and four Japonica subpopulations, including an outlying Indica-5 group. Here, we leveraged that strong differentiation and 672 native rice genomes to identify genomic regions and genes putatively selected during the breeding of rice in Vietnam. We identified significant distorted patterns in allele frequency (XP-CLR) and population differentiation scores (F ST) resulting from differential selective pressures between native subpopulations, and later annotated them with QTLs previously identified by GWAS in the same panel. We particularly focussed on the outlying Indica-5 subpopulation because of its likely novelty and differential evolution, where we annotated 52 selected regions, which represented 8.1% of the rice genome. We annotated the 4576 genes in these regions and selected 65 candidate genes as promising breeding targets, several of which harboured alleles with nonsynonymous substitutions. Our results highlight genomic differences between traditional Vietnamese landraces, which are likely the product of adaption to multiple environmental conditions and regional culinary preferences in a very diverse country. We also verified the applicability of this genome scanning approach to identify potential regions harbouring novel loci and alleles to breed a new generation of sustainable and resilient rice.

Unraveling the Genome of a High Yielding Colombian Sugarcane Hybrid.

Recent developments in High Throughput Sequencing (HTS) technologies and bioinformatics, including improved read lengths and genome assemblers allow the reconstruction of complex genomes with unprecedented quality and contiguity. Sugarcane has one of the most complicated genomes among grassess with a haploid length of 1Gbp and a ploidies between 8 and 12. In this work, we present a genome assembly of the Colombian sugarcane hybrid CC 01-1940. Three types of sequencing technologies were combined for this assembly: PacBio long reads, Illumina paired short reads, and Hi-C reads. We achieved a median contig length of 34.94 Mbp and a total genome assembly of 903.2 Mbp. We annotated a total of 63,724 protein coding genes and performed a reconstruction and comparative analysis of the sucrose metabolism pathway. Nucleotide evolution measurements between orthologs with close species suggest that divergence between Saccharum officinarum and Saccharum spontaneum occurred <2 million years ago. Synteny analysis between CC 01-1940 and the S. spontaneum genome confirms the presence of translocation events between the species and a random contribution throughout the entire genome in current sugarcane hybrids. Analysis of RNA-Seq data from leaf and root tissue of contrasting sugarcane genotypes subjected to water stress treatments revealed 17,490 differentially expressed genes, from which 3,633 correspond to genes expressed exclusively in tolerant genotypes. We expect the resources presented here to serve as a source of information to improve the selection processes of new varieties of the breeding programs of sugarcane.

Global multi-environment resistance QTL for foliar late blight resistance in tetraploid potato with tropical adaptation.

The identification of environmentally stable and globally predictable resistance to potato late blight is challenged by the clonal and polyploid nature of the crop and the rapid evolution of the pathogen. A diversity panel of tetraploid potato germplasm bred for multiple resistance and quality traits was genotyped by genotyping by sequencing (GBS) and evaluated for late blight resistance in three countries where the International Potato Center (CIP) has established breeding work. Health-indexed, in vitro plants of 380 clones and varieties were distributed from CIP headquarters and tuber seed was produced centrally in Peru, China, and Ethiopia. Phenotypes were recorded following field exposure to local isolates of Phytophthora infestans. QTL explaining resistance in four experiments conducted across the three countries were identified in chromosome IX, and environment-specific QTL were found in chromosomes III, V, and X. Different genetic models were evaluated for prediction ability to identify best performing germplasm in each and all environments. The best prediction ability (0.868) was identified with the genomic best linear unbiased predictors (GBLUPs) when using the diploid marker data and QTL-linked markers as fixed effects. Genotypes with high levels of resistance in all environments were identified from the B3, LBHT, and B3-LTVR populations. The results show that many of the advanced clones bred in Peru for high levels of late blight resistance maintain their resistance in Ethiopia and China, suggesting that the centralized selection strategy has been largely successful.

Flow Cytometry-Based Determination of Ploidy from Dried Leaf Specimens in Genomically Complex Collections of the Tropical Forage Grass Urochloa s. l.

Urochloa (including Brachiaria, Megathyrus and some Panicum) tropical grasses are native to Africa and are now, after selection and breeding, planted worldwide, particularly in South America, as important forages with huge potential for further sustainable improvement and conservation of grasslands. We aimed to develop an optimized approach to determine ploidy of germplasm collection of this tropical forage grass group using dried leaf material, including approaches to collect, dry and preserve plant samples for flow cytometry analysis. Our methods enable robust identification of ploidy levels (coefficient of variation of G0/G1 peaks, CV, typically <5%). Ploidy of some 348 forage grass accessions (ploidy range from 2x to 9x), from international genetic resource collections, showing variation in basic chromosome numbers and reproduction modes (apomixis and sexual), were determined using our defined standard protocol. Two major Urochloa agamic complexes are used in the current breeding programs at CIAT and EMBRAPA: the 'brizantha' and 'humidicola' agamic complexes are variable, with multiple ploidy levels. Some U. brizantha accessions have odd level of ploidy (5x), and the relative differences in fluorescence values of the peak positions between adjacent cytotypes is reduced, thus more precise examination of this species is required. Ploidy measurement of U. humidicola revealed aneuploidy.

Allele mining in diverse accessions of tropical grasses to improve forage quality and reduce environmental impact.

BACKGROUND AND AIMS: The C4Urochloa species (syn. Brachiaria) and Megathyrsus maximus (syn. Panicum maximum) are used as pasture for cattle across vast areas in tropical agriculture systems in Africa and South America. A key target for variety improvement is forage quality: enhanced digestibility could decrease the amount of land required per unit production, and enhanced lipid content could decrease methane emissions from cattle. For these traits, loss-of-function (LOF) alleles in known gene targets are predicted to improve them, making a reverse genetics approach of allele mining feasible. We therefore set out to look for such alleles in diverse accessions of Urochloa species and Megathyrsus maximus from the genebank collection held at the CIAT. METHODS: We studied allelic diversity of 20 target genes (11 for digestibility, nine for lipid content) in 104 accessions selected to represent genetic diversity and ploidy levels of U. brizantha, U. decumbens, U. humidicola, U. ruziziensis and M. maximum. We used RNA sequencing and then bait capture DNA sequencing to improve gene models in a U. ruziziensis reference genome to assign polymorphisms with high confidence. KEY RESULTS: We found 953 non-synonymous polymorphisms across all genes and accessions; within these, we identified seven putative LOF alleles with high confidence, including those in the non-redundant SDP1 and BAHD01 genes present in diploid and tetraploid accessions. These LOF alleles could respectively confer increased lipid content and digestibility if incorporated into a breeding programme. CONCLUSIONS: We demonstrated a novel, effective approach to allele discovery in diverse accessions using a draft reference genome from a single species. We used this to find gene variants in a collection of tropical grasses that could help reduce the environmental impact of cattle production.

Resequencing of 672 Native Rice Accessions to Explore Genetic Diversity and Trait Associations in Vietnam.

BACKGROUND: Vietnam possesses a vast diversity of rice landraces due to its geographical situation, latitudinal range, and a variety of ecosystems. This genetic diversity constitutes a highly valuable resource at a time when the highest rice production areas in the low-lying Mekong and Red River Deltas are enduring increasing threats from climate changes, particularly in rainfall and temperature patterns. RESULTS: We analysed 672 Vietnamese rice genomes, 616 newly sequenced, that encompass the range of rice varieties grown in the diverse ecosystems found throughout Vietnam. We described four Japonica and five Indica subpopulations within Vietnam likely adapted to the region of origin. We compared the population structure and genetic diversity of these Vietnamese rice genomes to the 3000 genomes of Asian cultivated rice. The named Indica-5 (I5) subpopulation was expanded in Vietnam and contained lowland Indica accessions, which had very low shared ancestry with accessions from any other subpopulation and were previously overlooked as admixtures. We scored phenotypic measurements for nineteen traits and identified 453 unique genotype-phenotype significant associations comprising twenty-one QTLs (quantitative trait loci). The strongest associations were observed for grain size traits, while weaker associations were observed for a range of characteristics, including panicle length, heading date and leaf width. CONCLUSIONS: We showed how the rice diversity within Vietnam relates to the wider Asian rice diversity by using a number of approaches to provide a clear picture of the novel diversity present within Vietnam, mainly around the Indica-5 subpopulation. Our results highlight differences in genome composition and trait associations among traditional Vietnamese rice accessions, which are likely the product of adaption to multiple environmental conditions and regional preferences in a very diverse country. Our results highlighted traits and their associated genomic regions that are a potential source of novel loci and alleles to breed a new generation of low input sustainable and climate resilient rice.

Differential expression of starch and sucrose metabolic genes linked to varying biomass yield in Miscanthus hybrids.

BACKGROUND: Miscanthus is a commercial lignocellulosic biomass crop owing to its high biomass productivity and low chemical input requirements. Within an interspecific Miscanthus cross, progeny with high biomass yield were shown to have low concentrations of starch and sucrose but high concentrations of fructose. We performed a transcriptional RNA-seq analysis between selected Miscanthus hybrids with contrasting values for these phenotypes to clarify how these phenotypes are genetically controlled. RESULTS: We observed that genes directly involved in the synthesis and degradation of starch and sucrose were down-regulated in high-yielding Miscanthus hybrids. At the same time, glycolysis and export of triose phosphates were up-regulated in high-yielding Miscanthus hybrids. These differentially expressed genes and biological functions were regulated by a well-connected network of less than 25 co-regulated transcription factors. CONCLUSIONS: Our results evidence a direct relationship between high expression of essential enzymatic genes in the starch and sucrose pathways and co-expression with their transcriptional regulators, with high starch concentrations and lower biomass production. The strong interconnectivity between gene expression and regulators, chemotype and agronomic traits opens the door to use the expression of well-characterised genes associated with carbohydrate metabolism, particularly in the starch and sucrose pathway, for the early selection of high biomass-yielding genotypes from large Miscanthus populations.

A Comparison of Differential Gene Expression in Response to the Onset of Water Stress Between Three Hybrid Brachiaria Genotypes.

Brachiaria (Trin.) Griseb. (syn. Urochloa P. Beauv.) is a C4 grass genus belonging to the Panicoideae. Native to Africa, these grasses are now widely grown as forages in tropical areas worldwide and are the subject of intensive breeding, particularly in South America. Tolerance to abiotic stresses such as aluminum and drought are major breeding objectives. In this study, we present the transcriptomic profiling of leaves and roots of three Brachiaria interspecific hybrid genotypes with the onset of water stress, Br12/3659-17 (gt-17), Br12/2360-9 (gt-9), and Br12/3868-18 (gt-18), previously characterized as having good, intermediate and poor tolerance to drought, respectively, in germplasm evaluation programs. RNA was extracted from leaf and root tissue of plants at estimated growing medium water contents (EWC) of 35, 15, and 5%. Differentially expressed genes (DEGs) were compared between different EWCs, 35/15, 15/5, and 35/5 using DESeq2. Overall, the proportions of DEGs enriched in all three genotypes varied in a genotype-dependent manner in relation to EWC comparison, with intermediate and sensitive gt-9 and gt-18 being more similar to each other than to drought tolerant gt-17. More specifically, GO terms relating to carbohydrate and cell wall metabolism in the leaves were enriched by up-regulated DEGs in gt-9 and gt-18, but by down-regulated DEGs in gt-17. Across all genotypes, analysis of DEG enzyme activities indicated an excess of down-regulated putative apoplastic peroxidases in the roots as water stress increased. This suggests that changes in root cell-wall architecture may be an important component of the response to water stress in Brachiaria.

Physiological and transcriptional response to drought stress among bioenergy grass Miscanthus species.

BACKGROUND: Miscanthus is a commercial lignocellulosic biomass crop owing to its high biomass productivity, resilience and photosynthetic capacity at low temperature. These qualities make Miscanthus a particularly good candidate for temperate marginal land, where yields can be limited by insufficient or excessive water supply. Differences in response to water stress have been observed among Miscanthus species, which correlated to origin. In this study, we compared the physiological and molecular responses among Miscanthus species under excessive (flooded) and insufficient (drought) water supply in glasshouse conditions. RESULTS: A significant biomass loss was observed under drought conditions in all genotypes. M. x giganteus showed a lower reduction in biomass yield under drought conditions compared to the control than the other species. Under flooded conditions, biomass yield was as good as or better than control conditions in all species. 4389 of the 67,789 genes (6.4%) in the reference genome were differentially expressed during drought among four Miscanthus genotypes from different species. We observed the same biological processes were regulated across Miscanthus species during drought stress despite the DEGs being not similar. Upregulated differentially expressed genes were significantly involved in sucrose and starch metabolism, redox, and water and glycerol homeostasis and channel activity. Multiple copies of the starch metabolic enzymes BAM and waxy GBSS-I were strongly up-regulated in drought stress in all Miscanthus genotypes, and 12 aquaporins (PIP1, PIP2 and NIP2) were also up-regulated in drought stress across genotypes. CONCLUSIONS: Different phenotypic responses were observed during drought stress among Miscanthus genotypes from different species, supporting differences in genetic adaption. The low number of DEGs and higher biomass yield in flooded conditions supported Miscanthus use in flooded land. The molecular processes regulated during drought were shared among Miscanthus species and consistent with functional categories known to be critical during drought stress in model organisms. However, differences in the regulated genes, likely associated with ploidy and heterosis, highlighted the value of exploring its diversity for breeding.

Draft genome assembly of the biofuel grass crop Miscanthus sacchariflorus.

Miscanthus sacchariflorus (Maxim.) Hack. is a highly productive C4 perennial rhizomatous biofuel grass crop. M. sacchariflorus is among the most widely distributed species in the genus, particularly at cold northern latitudes, and is one of the progenitor species of the commercial M. × giganteus genotypes. We generated a 2.54 Gb whole-genome assembly of the diploid M. sacchariflorus cv. "Robustus 297" genotype, which represented ~59% of the expected total genome size. We later anchored this assembly using the chromosomes from the M. sinensis genome to generate a second assembly with improved contiguity. We annotated 86,767 and 69,049 protein-coding genes in the unanchored and anchored assemblies, respectively. We estimated our assemblies included ~85% of the M. sacchariflorus genes based on homology and core markers. The utility of the new reference for genomic studies was evidenced by a 99% alignment rate of the RNA-seq reads from the same genotype.  The raw data, unanchored and anchored assemblies, and respective gene annotations are publicly available.

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.

Genome biology of the paleotetraploid perennial biomass crop Miscanthus.

Miscanthus is a perennial wild grass that is of global importance for paper production, roofing, horticultural plantings, and an emerging highly productive temperate biomass crop. We report a chromosome-scale assembly of the paleotetraploid M. sinensis genome, providing a resource for Miscanthus that links its chromosomes to the related diploid Sorghum and complex polyploid sugarcanes. The asymmetric distribution of transposons across the two homoeologous subgenomes proves Miscanthus paleo-allotetraploidy and identifies several balanced reciprocal homoeologous exchanges. Analysis of M. sinensis and M. sacchariflorus populations demonstrates extensive interspecific admixture and hybridization, and documents the origin of the highly productive triploid bioenergy crop M. × giganteus. Transcriptional profiling of leaves, stem, and rhizomes over growing seasons provides insight into rhizome development and nutrient recycling, processes critical for sustainable biomass accumulation in a perennial temperate grass. The Miscanthus genome expands the power of comparative genomics to understand traits of importance to Andropogoneae grasses.

Biological Nitrification Inhibition (BNI): Phenotyping of a Core Germplasm Collection of the Tropical Forage Grass Megathyrsus maximus Under Greenhouse Conditions.

Modern intensively managed pastures that receive large external nitrogen (N) inputs account for high N losses in form of nitrate (NO3 -) leaching and emissions of the potent greenhouse gas nitrous oxide (N2O). The natural plant capacity to shape the soil N cycle through exudation of organic compounds can be exploited to favor N retention without affecting productivity. In this study, we estimated the relationship between biological nitrification inhibition (BNI), N2O emissions and plant productivity for 119 germplasm accessions of Guineagrass (Megathyrsus maximus), an important tropical forage crop for livestock production. This relation was tested in a greenhouse experiment measuring BNI as (i) rates of soil nitrification; (ii) abundance of ammonia-oxidizing bacteria (AOB) and archaea (AOA); and (iii) the capacity of root tissue extracts to inhibit nitrification in vitro. We then measured N2O emissions, aboveground biomass and forage nutrition quality parameters. Reductions on nitrification activity ranging between 30 and 70% were found across the germplasm collection of M. maximus. Accessions with low nitrification rates showed a lower abundance of AOB as well as a reduction in N2O emissions compared to accessions of high nitrification rates. The BNI capacity was not correlated to N uptake of plants, suggesting that there may be intraspecific variation in the exploitation of different N sources in this grass species. A group of accessions (cluster) with the most desirable agronomic and environmental traits among the collection was identified for further field validation. These results provide evidence of the ability of M. maximus to suppress soil nitrification and N2O emissions and their relationship with productivity and forage quality, pointing a way to develop N conservative improved forage grasses for tropical livestock production.

Population structure and genetic diversity in red clover (Trifolium pratense L.) germplasm.

Red clover (Trifolium pratense L.) is a highly adaptable forage crop for temperate livestock agriculture. Genetic variation can be identified, via molecular techniques, and used to assess diversity among populations that may otherwise be indistinguishable. Here we have used genotyping by sequencing (GBS) to determine the genetic variation and population structure in red clover natural populations from Europe and Asia, and varieties or synthetic populations. Cluster analysis differentiated the collection into four large regional groups: Asia, Iberia, UK, and Central Europe. The five varieties clustered with the geographical area from which they were derived. Two methods (BayeScan and Samßada) were used to search for outlier loci indicating signatures of selection. A total of 60 loci were identified by both methods, but no specific genomic region was highlighted. The rate of decay in linkage disequilibrium was fast, and no significant evidence of any bottlenecks was found. Phenotypic analysis showed that a more prostrate and spreading growth habit was predominantly found among populations from Iberia and the UK. A genome wide association study identified a single nucleotide polymorphism (SNP) located in a homologue of the VEG2 gene from pea, associated with flowering time. The identification of genetic variation within the natural populations is likely to be useful for enhancing the breeding of red clover in the future.

Comprehensive assembly and analysis of the transcriptome of maritime pine developing embryos.

BACKGROUND: There are clear differences in embryo development between angiosperm and gymnosperm species. Most of the current knowledge on gene expression and regulation during plant embryo development has derived from studies on angiosperms species, in particular from the model plant Arabidopsis thaliana. The few published studies on transcript profiling of conifer embryogenesis show the existence of many putative embryo-specific transcripts without an assigned function. In order to extend the knowledge on the transcriptomic expression during conifer embryogenesis, we sequenced the transcriptome of zygotic embryos for several developmental stages that cover most of Pinus pinaster (maritime pine) embryogenesis. RESULTS: Total RNA samples collected from five zygotic embryo developmental stages were sequenced with Illumina technology. A de novo transcriptome was assembled as no genome sequence is yet published for Pinus pinaster. The transcriptome of reference for the period of zygotic embryogenesis in maritime pine contains 67,429 transcripts, which likely encode 58,527 proteins. The annotation shows a significant percentage, 31%, of predicted proteins exclusively present in pine embryogenesis. Functional categories and enrichment analysis of the differentially expressed transcripts evidenced carbohydrate transport and metabolism over-representation in early embryo stages, as highlighted by the identification of many putative glycoside hydrolases, possibly associated with cell wall modification, and carbohydrate transport transcripts. Moreover, the predominance of chromatin remodelling events was detected in early to middle embryogenesis, associated with an active synthesis of histones and their post-translational modifiers related to increased transcription, as well as silencing of transposons. CONCLUSIONS: Our results extend the understanding of gene expression and regulation during zygotic embryogenesis in conifers and are a valuable resource to support further improvements in somatic embryogenesis for vegetative propagation of conifer species. Specific transcripts associated with carbohydrate metabolism, monosaccharide transport and epigenetic regulation seem to play an important role in pine early embryogenesis and may be a source of reliable molecular markers for early embryogenesis.

A Parthenogenesis Gene Candidate and Evidence for Segmental Allopolyploidy in Apomictic Brachiaria decumbens.

Apomixis, asexual reproduction through seed, enables breeders to identify and faithfully propagate superior heterozygous genotypes by seed without the disadvantages of vegetative propagation or the expense and complexity of hybrid seed production. The availability of new tools such as genotyping by sequencing and bioinformatics pipelines for species lacking reference genomes now makes the construction of dense maps possible in apomictic species, despite complications including polyploidy, multisomic inheritance, self-incompatibility, and high levels of heterozygosity. In this study, we developed saturated linkage maps for the maternal and paternal genomes of an interspecific Brachiaria ruziziensis (R. Germ. and C. M. Evrard) × B. decumbens Stapf. F1 mapping population in order to identify markers linked to apomixis. High-resolution molecular karyotyping and comparative genomics with Setaria italica (L.) P. Beauv provided conclusive evidence for segmental allopolyploidy in B. decumbens, with strong preferential pairing of homologs across the genome and multisomic segregation relatively more common in chromosome 8. The apospory-specific genomic region (ASGR) was mapped to a region of reduced recombination on B. decumbens chromosome 5. The Pennisetum squamulatum (L.) R.Br. PsASGR-BABY BOOM-like (psASGR-BBML)-specific primer pair p779/p780 was in perfect linkage with the ASGR in the F1 mapping population and diagnostic for reproductive mode in a diversity panel of known sexual and apomict Brachiaria (Trin.) Griseb. and P. maximum Jacq. germplasm accessions and cultivars. These findings indicate that ASGR-BBML gene sequences are highly conserved across the Paniceae and add further support for the postulation of the ASGR-BBML as candidate genes for the apomictic function of parthenogenesis.

Red clover (Trifolium pratense L.) draft genome provides a platform for trait improvement.

Red clover (Trifolium pratense L.) is a globally significant forage legume in pastoral livestock farming systems. It is an attractive component of grassland farming, because of its high yield and protein content, nutritional value and ability to fix atmospheric nitrogen. Enhancing its role further in sustainable agriculture requires genetic improvement of persistency, disease resistance, and tolerance to grazing. To help address these challenges, we have assembled a chromosome-scale reference genome for red clover. We observed large blocks of conserved synteny with Medicago truncatula and estimated that the two species diverged ~23 million years ago. Among the 40,868 annotated genes, we identified gene clusters involved in biochemical pathways of importance for forage quality and livestock nutrition. Genotyping by sequencing of a synthetic population of 86 genotypes show that the number of markers required for genomics-based breeding approaches is tractable, making red clover a suitable candidate for association studies and genomic selection.