Genome, Transcriptome, and Germplasm Sequencing Uncovers Functional Variation in the Warm-Season Grain Legume Horsegram Macrotyloma uniflorum (Lam.) Verdc.

Horsegram is a grain legume with excellent nutritional and remedial properties and good climate resilience, able to adapt to harsh environmental conditions. Here, we used a combination of short- and long-read sequencing technologies to generate a genome sequence of 279.12Mb, covering 83.53% of the estimated total size of the horsegram genome, and we annotated 24,521 genes. De novo prediction of DNA repeats showed that approximately 25.04% of the horsegram genome was made up of repetitive sequences, the lowest among the legume genomes sequenced so far. The major transcription factors identified in the horsegram genome were bHLH, ERF, C2H2, WRKY, NAC, MYB, and bZIP, suggesting that horsegram is resistant to drought. Interestingly, the genome is abundant in Bowman-Birk protease inhibitors (BBIs), which can be used as a functional food ingredient. The results of maximum likelihood phylogenetic and estimated synonymous substitution analyses suggested that horsegram is closely related to the common bean and diverged approximately 10.17 million years ago. The double-digested restriction associated DNA (ddRAD) sequencing of 40 germplasms allowed us to identify 3,942 high-quality SNPs in the horsegram genome. A genome-wide association study with powdery mildew identified 10 significant associations similar to the MLO and RPW8.2 genes. The reference genome and other genomic information presented in this study will be of great value to horsegram breeding programs. In addition, keeping the increasing demand for food with nutraceutical values in view, these genomic data provide opportunities to explore the possibility of horsegram for use as a source of food and nutraceuticals.

Narrow genetic base shapes population structure and linkage disequilibrium in an industrial oilseed crop, Brassica carinata A. Braun.

Ethiopian mustard (Brassica carinata A. Braun) is an emerging sustainable source of vegetable oil, in particular for the biofuel industry. The present study exploited genome assemblies of the Brassica diploids, Brassica nigra and Brassica oleracea, to discover over 10,000 genome-wide SNPs using genotype by sequencing of 620 B. carinata lines. The analyses revealed a SNP frequency of one every 91.7 kb, a heterozygosity level of 0.30, nucleotide diversity levels of 1.31 × 10-05, and the first five principal components captured only 13% molecular variation, indicating low levels of genetic diversity among the B. carinata collection. Genome bias was observed, with greater SNP density found on the B subgenome. The 620 lines clustered into two distinct sub-populations (SP1 and SP2) with the majority of accessions (88%) clustered in SP1 with those from Ethiopia, the presumed centre of origin. SP2 was distinguished by a collection of breeding lines, implicating targeted selection in creating population structure. Two selective sweep regions on B3 and B8 were detected, which harbour genes involved in fatty acid and aliphatic glucosinolate biosynthesis, respectively. The assessment of genetic diversity, population structure, and LD in the global B. carinata collection provides critical information to assist future crop improvement.

Sequencing Analysis of Genetic Loci for Resistance for Late Leaf Spot and Rust in Peanut (Arachis hypogaea L.).

The aim of this study was to identify candidate resistance genes for late leaf spot (LLS) and rust diseases in peanut (Arachis hypogaea L.). We used a double-digest restriction-site associated DNA sequencing (ddRAD-Seq) technique based on next-generation sequencing (NGS) for genotyping analysis across the recombinant inbred lines (RILs) derived from a cross between a susceptible line, TAG 24, and a resistant line, GPBD 4. A total of 171 SNPs from the ddRAD-Seq together with 282 markers published in the previous studies were mapped on a genetic map covering 1510.1 cM. Subsequent quantitative trait locus (QTL) analysis revealed major genetic loci for LLS and rust resistance on chromosomes A02 and A03, respectively. Heterogeneous inbred family-derived near isogenic lines and the pedigree of the resistant gene donor, A. cardenasii Krapov. & W.C. Greg., including the resistant derivatives of ICGV 86855 and VG 9514 as well as GPBD 4, were employed for whole-genome resequencing analysis. The results indicated the QTL candidates for LLS and rust resistance were located in 1.4- and 2.7-Mb genome regions on A02 and A03, respectively. In these regions, four and six resistance-related genes with deleterious mutations were selected as candidates for LLS and rust resistance, respectively. These delimited genomic regions may be beneficial in breeding programs aimed at improving disease resistance and enhancing peanut productivity.

Differential expression of two galactinol synthase isoforms LcGolS1 and LcGolS2 in developing lentil (Lens culinaris Medik. cv CDC Redberry) seeds.

Galactinol synthase (GS, EC 2.4.1.123) catalyzes the transfer of a galactosyl residue from UDP-galactose to myo-inositol to synthesize galactinol, a precursor for raffinose family oligosaccharides (RFO) biosynthesis. Screening, a cDNA library constructed with RNA isolated from developing lentil seeds, with partial GS genes resulted in identification of cDNA clones for two isoforms of GS, LcGolS1 (1336 bp, ORF-1002 bp, 334 amino acids) and LcGolS2 (1324bp, ORF-975bp, 325 amino acids) with predicted molecular weights of 38.7 kDa and 37.6 kDa, respectively. During lentil seed development, LcGolS1 transcripts showed higher accumulation during 26-32 days after flowering (DAF) corresponding to seed desiccation, while LcGolS2 showed maximum accumulation at 24 DAF, prior to increase in LcGolS1 transcripts. GS enzyme activity was maximum at 26 and 28 DAF and corresponded to galactinol accumulation, which also increased rapidly at 22 DAF with maximum accumulation at 26 DAF. Substrates for GS activity, myo-inositol and glucose/galactose were present in high concentrations during early stages of seed development but gradually decreased from 20 DAF to 32 DAF when galactinol concentration increased coinciding with increased GS enzyme activity.

The developmental transcriptome atlas of the biofuel crop Camelina sativa.

Camelina sativa is currently being embraced as a viable industrial bio-platform crop due to a number of desirable agronomic attributes and the unique fatty acid profile of the seed oil that has applications for food, feed and biofuel. The recent completion of the reference genome sequence of C. sativa identified a young hexaploid genome. To complement this work, we have generated a genome-wide developmental transcriptome map by RNA sequencing of 12 different tissues covering major developmental stages during the life cycle of C. sativa. We have generated a digital atlas of this comprehensive transcriptome resource that enables interactive visualization of expression data through a searchable database of electronic fluorescent pictographs (eFP browser). An analysis of this dataset supported expression of 88% of the annotated genes in C. sativa and provided a global overview of the complex architecture of temporal and spatial gene expression patterns active during development. Conventional differential gene expression analysis combined with weighted gene expression network analysis uncovered similarities as well as differences in gene expression patterns between different tissues and identified tissue-specific genes and network modules. A high-quality census of transcription factors, analysis of alternative splicing and tissue-specific genome dominance provided insight into the transcriptional dynamics and sub-genome interplay among the well-preserved triplicated repertoire of homeologous loci. The comprehensive transcriptome atlas in combination with the reference genome sequence provides a powerful resource for genomics research which can be leveraged to identify functional associations between genes and understand the regulatory networks underlying developmental processes.

A high-density SNP genotyping array for Brassica napus and its ancestral diploid species based on optimised selection of single-locus markers in the allotetraploid genome.

KEY MESSAGE: The Brassica napus Illumina array provides genome-wide markers linked to the available genome sequence, a significant tool for genetic analyses of the allotetraploid B. napus and its progenitor diploid genomes. A high-density single nucleotide polymorphism (SNP) Illumina Infinium array, containing 52,157 markers, was developed for the allotetraploid Brassica napus. A stringent selection process employing the short probe sequence for each SNP assay was used to limit the majority of the selected markers to those represented a minimum number of times across the highly replicated genome. As a result approximately 60 % of the SNP assays display genome-specificity, resolving as three clearly separated clusters (AA, AB, and BB) when tested with a diverse range of B. napus material. This genome specificity was supported by the analysis of the diploid ancestors of B. napus, whereby 26,504 and 29,720 markers were scorable in B. oleracea and B. rapa, respectively. Forty-four percent of the assayed loci on the array were genetically mapped in a single doubled-haploid B. napus population allowing alignment of their physical and genetic coordinates. Although strong conservation of the two positions was shown, at least 3 % of the loci were genetically mapped to a homoeologous position compared to their presumed physical position in the respective genome, underlying the importance of genetic corroboration of locus identity. In addition, the alignments identified multiple rearrangements between the diploid and tetraploid Brassica genomes. Although mostly attributed to genome assembly errors, some are likely evidence of rearrangements that occurred since the hybridisation of the progenitor genomes in the B. napus nucleus. Based on estimates for linkage disequilibrium decay, the array is a valuable tool for genetic fine mapping and genome-wide association studies in B. napus and its progenitor genomes.

A reliable and rapid method for soluble sugars and RFO analysis in chickpea using HPAEC-PAD and its comparison with HPLC-RI.

A high performance anion exchange chromatography (HPAEC) coupled with pulsed amperometric detection (PAD) was optimised to separate with precision, accuracy and high reproducibility soluble sugars including oligosaccharides present in pulse meal samples. The optimised method within 20min separated myo-inositol, galactinol, glucose, fructose, sucrose, raffinose, stachyose and verbascose in chickpea seed meal extracts. Gradient method of eluting solvent (sodium hydroxide) resulted in higher sensitivity and rapid detection compared to similar analytical methods. Peaks asymmetry equivalent to one and resolution value ⩾1.5 support column's precision and accuracy for quantitative determinations of soluble sugars in complex mixtures. Intermediate precision determined as relative standard deviation (1.8-3.5%) for different soluble sugars confirms reproducibility of the optimised method. The developed method has superior sensitivity to detect even scarcely present verbascose in chickpea. It also quantifies myo-inositol and galactinol making it suitable both for RFO related genotype screening and biosynthetic studies.

Genotype and growing environment interaction shows a positive correlation between substrates of raffinose family oligosaccharides (RFO) biosynthesis and their accumulation in chickpea ( Cicer arietinum L.) seeds.

To develop genetic improvement strategies to modulate raffinose family oligosaccharides (RFO) concentration in chickpea ( Cicer arietinum L.) seeds, RFO and their precursor concentrations were analyzed in 171 chickpea genotypes from diverse geographical origins. The genotypes were grown in replicated trials over two years in the field (Patancheru, India) and in the greenhouse (Saskatoon, Canada). Analysis of variance revealed a significant impact of genotype, environment, and their interaction on RFO concentration in chickpea seeds. Total RFO concentration ranged from 1.58 to 5.31 mmol/100 g and from 2.11 to 5.83 mmol/100 g in desi and kabuli genotypes, respectively. Sucrose (0.60-3.59 g/100 g) and stachyose (0.18-2.38 g/100 g) were distinguished as the major soluble sugar and RFO, respectively. Correlation analysis revealed a significant positive correlation between substrate and product concentration in RFO biosynthesis. In chickpea seeds, raffinose, stachyose, and verbascose showed a moderate broad sense heritability (0.25-0.56), suggesting the use of a multilocation trials based approach in chickpea seed quality improvement programs.