Assessment of Genetic Diversity and Structure of Large Garlic (Allium sativum) Germplasm Bank, by Diversity Arrays Technology "Genotyping-by-Sequencing" Platform (DArTseq)

DArTseq-Based, High-Throughput Identification of Novel Molecular Markers for the Detection of Blackleg (Leptosphaeria Spp.) Resistance in Rapeseed

Blackleg disease, caused by Leptosphaeria spp. fungi, is one of the most important diseases of Brassica napus, responsible for severe yield losses worldwide. Blackleg resistance is controlled by major R genes and minor quantitative trait loci (QTL). Due to the high adaptation ability of the pathogen, R-mediated resistance can be easily broken, while the resistance mediated via QTL is believed to be more durable. Thus, the identification of novel molecular markers linked to blackleg resistance for B. napus breeding programs is essential. In this study, 183 doubled haploid (DH) rapeseed lines were assessed in field conditions for resistance to Leptosphaeria spp. Subsequently, DArTseq-based Genome-Wide Association Study (GWAS) was performed to identify molecular markers linked to blackleg resistance. A total of 133,764 markers (96,121 SilicoDArT and 37,643 SNP) were obtained. Finally, nine SilicoDArT and six SNP molecular markers were associated with plant resistance to Leptosphaeria spp. at the highest significance level, p < 0.001. Importantly, eleven of these fifteen markers were found within ten genes located on chromosomes A06, A07, A08, C02, C03, C06 and C08. Given the immune-related functions of the orthologues of these genes in Arabidopsis thaliana, the identified markers hold great promise for application in rapeseed breeding programs.

Ecological genetics of isolated loach populations indicate compromised adaptive potential

Many endangered species live in fragmented and isolated populations with low genetic variability, signs of inbreeding, and small effective population sizes - all features elevating their extinction risk. The flat-headed loach (Oreonectes platycephalus), a small noemacheilid fish, is widely across southern China, but only in the headwaters of hillstreams; as a result, they are spatially isolated from conspecific populations. We surveyed single nucleotide polymorphisms in 16 Hong Kong populations of O. platycephalus to determine whether loach populations from different streams were genetically isolated from each other, showed low levels of genetic diversity, signs of inbreeding, and had small contemporary effective population sizes. Estimates of average observed heterozygosity (HO = 0.0473), average weighted nucleotide diversity (πw = 0.0546) and contemporary effective population sizes (Ne = 10.2 ~ 129.8) were very low, and several populations showed clear signs of inbreeding as judged from relatedness estimates. The degree of genetic differentiation among populations was very high (average FST = 0.668), even over short geographic distances (<1.5 km), with clear patterns of isolation by distance. These results suggest that Hong Kong populations of O. platycephalus have experienced strong genetic drift and loss of genetic variability because sea-level rise after the last glaciation reduced connectedness among paleodrainages, isolating populations in headwaters. All this, together with the fact that the levels of genetic diversity and contemporary effective population sizes within O. platycephalus populations are lower than most other freshwater fishes, suggests that they face high local extinction risk and have limited capacity for future adaptation.

Genetic Heterogeneity in Cowpea Genotypes (Vigna unguiculata L. Walp) Using DArTseq (GBS)-Derived Single Nucleotide Polymorphisms

Cowpeas (Vigna unguiculata L. Walp) have been credible constituents of nutritious food and forage in human and animal diets since the Neolithic era. The modern technique of Diversity Array Technology (DArTseq) is both cost-effective and rapid in producing thousands of high-throughputs, genotyped, single nucleotide polymorphisms (SNPs) in wide-genomic analyses of genetic diversity. The aim of this study was to assess the heterogeneity in cowpea genotypes using DArTseq-derived SNPs. A total of 92 cowpea genotypes were selected, and their fourteen-day-old leaves were freeze-dried for five days. DNA was extracted using the CTAB protocol, genotyped using DArTseq, and analysed using DArTsoft14. A total of 33,920 DArTseq-derived SNPs were recalled for filtering analysis, with a final total of 16,960 SNPs. The analyses were computed using vcfR, poppr, and ape in R Studio v1.2.5001-3 software. The heatmap revealed that the TVU 9596 (SB26), Orelu (SB72), 90K-284-2 (SB55), RV 403 (SB17), and RV 498 (SB16) genotypes were heterogenous. The mean values for polymorphic information content, observed heterozygosity, expected heterozygosity, major allele frequency, and the inbreeding coefficient were 0.345, 0.386, 0.345, 0.729, and 0.113, respectively. Moreover, they validated the diversity of the evaluated cowpea genotypes, which could be used for potential breeding programmes and management of cowpea germplasm.

Phenotypic characterization and analysis of genetic diversity between commercial crossbred and indigenous chickens from three different agro-ecological zones using DArT-Seq technology

Indigenous and were used to study genetic diversity and population structure analyses. Polymorphism information content (PIC) values ranged from 0.0 to 0.5, with 21,285 SNP markers (35%) being in the lowest PIC value range (0 to 0.15) while 13,511 (commercial chickens have developed unique adaptations to their environments, which may include nutrition, pathogens, and thermal stress. Besides, environmental pressures and artificial selection have generated significant genome-wide divergence in chickens, as those selection pressures contribute a considerable evolutionary force to phenotypic and genotypic differentiation. Herein, we determined genomic diversity of indigenous chickens from semi-deciduous rainforest (SDR), coastal savannah (CS) and Guinea savannah (GS) agro-ecological zones (AEZs) in Ghana and commercial crossbreds (CC) reared at the Kwame Nkrumah University of Science and Technology (KNUST). We generated SNP markers from 82 chickens (62 indigenous chicken ecotypes and 26 commercial crossbred ecotype) using DArT-Seq technology. A total of 85,396 SNP markers were generated and after filtering the data, 58,353 markers 21%) were in the highest PIC value range (0.45 to 0.50). The CC were more genetically diverse than the indigenous birds, with the highest expected heterozygosity value of 0.220. Between the commercial crossbreds population and the indigenous ecotypes, pairwise FST values were estimated to be 0.105 between CS, 0.096 between SDF, and 0.133 between GS. Furthermore, PCA analysis showed that the CC, SDF and GS chickens clustered together and are genetically distant from the commercial crossbred. We herein show that chickens from the AEZs studied can be considered as one population. However, due the abundance of agro-byproducts in the SDR compared to the CS and GS, chickens from the SDR AEZ had better growth compared to their counterparts. It is suggested that the genetic diversity within the local ecotypes could form the basis for genetic improvement.

NMR Metabolite Profiling for the Characterization of Vessalico Garlic Ecotype and Bioactivity against Xanthomonas campestris pv. campestris

The Italian garlic ecotype "Vessalico" possesses distinct characteristics compared to its French parent cultivars Messidor and Messidrôme, used for sowing, as well as other ecotypes in neighboring regions. However, due to the lack of a standardized seed supply method and cultivation protocol among farmers in the Vessalico area, a need to identify garlic products that align with the Vessalico ecotype arises. In this study, an NMR-based approach followed by multivariate analysis to analyze the chemical composition of Vessalico garlic sourced from 17 different farms, along with its two French parent cultivars, was employed. Self-organizing maps allowed to identify a homogeneous subset of representative samples of the Vessalico ecotype. Through the OPLS-DA model, the most discriminant metabolites based on values of VIP (Variable Influence on Projections) were selected. Among them, S-allylcysteine emerged as a potential marker for distinguishing the Vessalico garlic from the French parent cultivars by NMR screening. Additionally, to promote sustainable agricultural practices, the potential of Vessalico garlic extracts and its main components as agrochemicals against Xanthomonas campestris pv. campestris, responsible for black rot disease, was explored. The crude extract exhibited a MIC of 125 μg/mL, and allicin demonstrated the highest activity among the tested compounds (MIC value of 31.25 μg/mL).

Effectiveness of DArTseq markers application in genetic diversity and population structure of indigenous chickens in Eastern Province of Rwanda

BACKGROUND

The application of biotechnologies which make use of genetic markers in chicken breeding is developing rapidly. Diversity Array Technology (DArT) is one of the current Genotyping-By-Sequencing techniques allowing the discovery of whole genome sequencing. In livestock, DArT has been applied in cattle, sheep, and horses. Currently, there is no study on the application of DArT markers in chickens. The aim was to study the effectiveness of DArTSeq markers in the genetic diversity and population structure of indigenous chickens (IC) and SASSO in the Eastern Province of Rwanda.

METHODS

In total 87 blood samples were randomly collected from 37 males and 40 females of indigenous chickens and 10 females of SASSO chickens purposively selected from 5 sites located in two districts of the Eastern Province of Rwanda. Genotyping by Sequencing (GBS) using DArTseq technology was employed. This involved the complexity reduction method through digestion of genomic DNA and ligation of barcoded adapters followed by PCR amplification of adapter-ligated fragments.

RESULTS

From 45,677 DArTseq SNPs and 25,444 SilicoDArTs generated, only 8,715 and 6,817 respectively remained for further analysis after quality control. The average call rates observed, 0.99 and 0.98 for DArTseq SNPs and SilicoDArTs respectively were quite similar. The polymorphic information content (PIC) from SilicoDArTs (0.33) was higher than that from DArTseq SNPs (0.22). DArTseq SNPs and SilicoDArTs had 34.4% and 34% of the loci respectively mapped on chromosome 1. DArTseq SNPs revealed distance averages of 0.17 and 0.15 within IC and SASSO chickens respectively while the respective averages observed with SilicoDArTs were 0.42 and 0.36. The average genetic distance between IC and SASSO chickens was moderate for SilicoDArTs (0.120) compared to that of DArTseq SNPs (0.048). The PCoA and population structure clustered the chicken samples into two subpopulations (1 and 2); 1 is composed of IC and 2 by SASSO chickens. An admixture was observed in subpopulation 2 with 12 chickens from subpopulation 1.

CONCLUSIONS

The application of DArTseq markers have been proven to be effective and efficient for genetic relationship between IC and separated IC from exotic breed used which indicate their suitability in genomic studies. However, further studies using all chicken genetic resources available and large big sample sizes are required.

Geographical differentiation of garlic based on HS-GC-IMS combined with multivariate statistical analysis

Garlic is famous for its unique flavor and health benefits. An effective means of authenticating garlic's origin is through the implementation of the Protected Geographical Indication (PGI) scheme. However, the prevalence of fraudulent behavior raises concerns regarding the reliability of this system. In this study, garlic samples from six distinct production areas (G1: Cangshan garlic, G2: Qixian garlic, G3: Dali single clove garlic, G4: Jinxiang garlic, G5: Yongnian garlic, and G6: Badong garlic) underwent analysis using HS-GC-IMS. A total of 26 VOCs were detected in the samples. The differences in VOCs among the different garlic samples were visually presented in a two-dimensional topographic map and fingerprint map. Principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA) were employed to demonstrate the capacity of the HS-GC-IMS method for effectively distinguishing garlic samples from different geographical sources. Further screening based on the p-value and VIP score threshold identified 12 different aroma substances, which can be utilized for the identification of garlic from different producing areas. The fusion of HS-GC-IMS with multivariate statistical analysis proved to be a rapid, intuitive, and efficient approach for identifying and categorizing garlic VOCs, offering a novel strategy for ascertaining garlic origin and ensuring quality control.

Characterization of the Moroccan Barley Germplasm Preserved in the Polish Genebank as a First Step towards Selecting Forms with Increased Drought Tolerance

In marginal, arid, and semi-arid areas of Morocco, crops are often exposed to multiple abiotic and biotic stresses that have a major impact on yield. Farmer-maintained Moroccan landraces have been shaped by the impact of very strong selection pressures, gradually adapting to the local ecosystem and obsolete low-input agricultural practices without improvement towards high yield and quality. Considering the increasing threat of drought in Poland, it is necessary to introduce germplasm with tolerance to water deficit into barley breeding programs. The aim of this research was a DArTseq-based genetic characterization of a collection of germplasm of Moroccan origin, conserved in the Polish genebank. The results showed that all conserved landraces have a high level of heterogeneity and their gene pool is different from the material developed by Polish breeders. Based on the analysis of eco-geographical data, locations with extremely different intensities of drought stress were selected. A total of 129 SNPs unique to accessions from these locations were identified. In the neighborhood of the clusters of unique SNPs on chromosomes 5H and 6H, genes that may be associated with plant response to drought stress were identified. The results obtained may provide a roadmap for further research to support Polish barley breeding for increased drought tolerance.

DREB1 and DREB2 Genes in Garlic (Allium sativum L.): Genome-Wide Identification, Characterization, and Stress Response

Dehydration-responsive element-binding (DREB) transcription factors (TFs) of the A1 and A2 subfamilies involved in plant stress responses have not yet been reported in Allium species. In this study, we used bioinformatics and comparative transcriptomics to identify and characterize DREB A1 and A2 genes redundant in garlic (Allium sativum L.) and analyze their expression in A. sativum cultivars differing in the sensitivity to cold and Fusarium infection. Eight A1 (AsaDREB1.1-1.8) and eight A2 (AsaDREB2.1-2.8) genes were identified. AsaDREB1.1-1.8 genes located in tandem on chromosome 1 had similar expression patterns, suggesting functional redundancy. AsaDREB2.1-2.8 were scattered on different chromosomes and had organ- and genotype-specific expressions. AsaDREB1 and AsaDREB2 promoters contained 7 and 9 hormone- and stress-responsive cis-regulatory elements, respectively, and 13 sites associated with TF binding and plant development. In both Fusarium-resistant and -sensitive cultivars, fungal infection upregulated the AsaDREB1.1-1.5, 1.8, 2.2, 2.6, and 2.8 genes and downregulated AsaDREB2.5, but the magnitude of response depended on the infection susceptibility of the cultivar. Cold exposure strongly upregulated the AsaDREB1 genes, but downregulated most AsaDREB2 genes. Our results provide the foundation for further functional analysis of the DREB TFs in Allium crops and could contribute to the breeding of stress-tolerant varieties.

Utilizing Genomics to Characterize the Common Oat Gene Pool—The Story of More than a Century of Polish Breeding

This study was undertaken to investigate the diversity and population structure of 487 oat accessions, including breeding lines from the ongoing programs of the three largest Polish breeding companies, along with modern and historical Polish and foreign cultivars. The analysis was based on 7411 DArTseq-derived SNPs distributed among three sub-genomes (A, C, and D). The heterogeneity of the studied material was very low, as only cultivars and advanced breeding lines were examined. Principal component analysis (PCA), principal coordinate analysis (PCoA), and cluster and STRUCTURE analyses found congruent results, which show that most of the examined cultivars and materials from Polish breeding programs formed major gene pools, that only some accessions derived from Strzelce Plant Breeding, and that foreign cultivars were outside of the main group. During the 120 year oat breeding process, only 67 alleles from the old gene pool were lost and replaced by 67 new alleles. The obtained results indicate that no erosion of genetic diversity was observed within the Polish native oat gene pool. Moreover, current oat breeding programs have introduced 673 new alleles into the gene pool relative to historical cultivars. The analysis also showed that most of the changes in relation to historical cultivars occurred within the A sub-genome with emphasis on chromosome 6A. The targeted changes were the rarest in the C sub-genome. This study showed that Polish oat breeding based mainly on traditional breeding methods—although focused on improving traits typical to this crop, i.e., enhancing the grain yield and quality and improving adaptability—did not significantly narrow the oat gene pool and in fact produced cultivars that are not only competitive in the European market but are also reservoirs of new alleles that were not found in the analyzed foreign materials.

Physiological specialization of Puccinia triticina and genome-wide association mapping provide insights into the genetics of wheat leaf rust resistance in Iran

Leaf rust caused by Puccinia triticina Erikss. (Pt) is the most widely distributed and important wheat disease worldwide. The objective of the present study was to determine the frequency of Iranian Pt races, their virulence to key resistance genes and map quantitative trait loci (QTL) for resistance to different Pt races from 185 globally diverse wheat genotypes using a genome-wide association study (GWAS) approach. The virulence pattern of the 33 Pt isolates from various wheat-growing areas of Iran on 55 wheat differentials showed that the FKTPS and FKTTS were relatively frequent pathotypes among the 18 identified races. The weighted average frequency of virulence on the resistance genes Lrb, Lr3bg, Lr14b, Lr16, Lr24, Lr3ka, Lr11 and Lr20 were high (> 90%). However, low virulence on the resistant genes Lr2a, Lr9, Lr19, Lr25, Lr28 and Lr29 indicates that these genes are still effective against the pathogen population in Iran at present. GWAS on a panel of 185 wheat genotypes against 10 Pt races resulted into 62 significant marker-trait associations (MTAs) belonged to 34 quantitative trait loci (QTL) across 16 chromosomes. Among them, 10 QTLs on chromosomes 1A, 1B, 3B, 3D, 4A, 6D, 7A and 7D were identified as potential novel QTLs, of which four QTLs (QLr.iau-3B-2, QLr.iau-7A-2, QLr.iau-7A-3 and QLr.iau-7D-2) are more interesting, as they are associated with resistance to two or more Pt races. The known and novel QTLs associated with different Pt races found here, can be used in future wheat breeding programs to recombine different loci for durable resistance against leaf rust races.

Turning Garlic into a Modern Crop: State of the Art and Perspectives

Garlic is cultivated worldwide for the value of its bulbs, but its cultivation is challenged by the infertility of commercial cultivars and the accumulation of pathogens over time, which occurs as a consequence of vegetative (clonal) propagation. In this review, we summarize the state of the art of garlic genetics and genomics, highlighting recent developments that will lead to its development as a modern crop, including the restoration of sexual reproduction in some garlic strains. The set of tools available to the breeder currently includes a chromosome-scale assembly of the garlic genome and multiple transcriptome assemblies that are furthering our understanding of the molecular processes underlying important traits like the infertility, the induction of flowering and bulbing, the organoleptic properties and resistance to various pathogens.

Phenotypic diversity and distinctiveness of the Belltall garlic landrace

Among the Mediterranean horticultural landraces, garlic is one of the crops most threatened by genetic erosion. Due to its sexual sterility and to the incidence of seed-borne diseases, historical varieties have been widely replaced by commercial cultivars. In Catalonia, despite the historical relevance of the crop, solely the Belltall garlic landrace is cultivated for commercial purposes. To assess the genotypic and phenotypic diversity within the Belltall garlic, we evaluated sixteen local accessions and five recognized traditional and modern varieties as controls. Genetic analysis with SSR and InDel markers showed low genetic diversity within the Belltall population, grouping modern and traditional varieties separately. Farmers and consumers were involved in the definition of the landrace ideotype and classified the materials by means of projective mapping. Scant phenotypic diversity was found within the Belltall landrace, which is characterized by its color profile and the small size of bulb and cloves. The Belltall landrace grown outside its area of origin lost the distinctive quality signals that differentiate the landrace from the commercial cultivars (clove appearance), indicating that the high quality of the landrace is under genotype-by-environment effects (i.e. local adaptation). Moreover, the size of the Belltall sowing clove had a strong effect on the harvested bulb size. Our research represents a case study for the description of the variability within garlic landraces and an approach to quantify the phenomenon of local adaptation that currently drives their conservation.

Genome-Wide Association Study Revealed SNP Alleles Associated with Seed Size Traits in African Yam Bean (Sphenostylis stenocarpa (Hochst ex. A. Rich.) Harms)

Seed size is an important yield and quality-determining trait in higher plants and is also crucial to their evolutionary fitness. In African yam bean (AYB), seed size varies widely among different accessions. However, the genetic basis of such variation has not been adequately documented. A genome-wide marker-trait association study was conducted to identify genomic regions associated with four seed size traits (seed length, seed width, seed thickness, and 100-seed weight) in a panel of 195 AYB accessions. A total of 5416 SNP markers were generated from the diversity array technology sequence (DArTseq) genotype-by-sequencing (GBS)- approach, in which 2491 SNPs were retained after SNP quality control and used for marker-trait association analysis. Significant phenotypic variation was observed for the traits. Broad-sense heritability ranged from 50.0% (seed width) to 66.4% (seed length). The relationships among the traits were positive and significant. Genome-wide association study (GWAS) using the general linear model (GLM) and the mixed linear model (MLM) approaches identified 12 SNP markers significantly associated with seed size traits across the six test environments. The 12 makers explained 6.5-10.8% of the phenotypic variation. Two markers (29420334|F|0-52:C>G-52:C>G and 29420736|F|0-57:G>T-57:G>T) with pleiotropic effects associated with seed width and seed thickness were found. A candidate gene search identified five significant markers (100026424|F|0-37:C>T-37:C>T, 100041049|F|0-42:G>C-42:G>C, 100034480|F|0-31:C>A-31:C>A, 29420365|F|0-55:C>G-55:C>G, and 29420736|F|0-57:G>T-57:G>T) located close to 43 putative genes whose encoding protein products are known to regulate seed size traits. This study revealed significant makers not previously reported for seed size in AYB and could provide useful information for genomic-assisted breeding in AYB.

SSR markers development and their application in genetic diversity evaluation of garlic (Allium sativum) germplasm

Garlic (Allium sativum), an asexually propagated vegetable and medicinal crop, has abundant genetic variation. Genetic diversity evaluation based on molecular markers has apparent advantages since their genomic abundance, environment insensitivity, and non-tissue specific features. However, the limited number of available DNA markers, especially SSR markers, are insufficient to conduct related genetic diversity assessment studies in garlic. In this study, 4372 EST-SSR markers were newly developed, and 12 polymorphic markers together with other 17 garlic SSR markers were used to assess the genetic diversity and population structure of 127 garlic accessions. The averaged polymorphism information content (PIC) of these 29 SSR markers was 0.36, ranging from 0.22 to 0.49. Seventy-nine polymorphic loci were detected among these accessions, with an average of 3.48 polymorphic loci per SSR. Both the clustering analyses based on either the genotype data of SSR markers or the phenotypic data of morphological traits obtained genetic distance divided the 127 garlic accessions into three clusters. Moreover, the Mantel test showed that genetic distance had no significant correlations with geographic distance, and weak correlations were found between genetic distance and the phenotypic traits. AMOVA analysis showed that the main genetic variation of this garlic germplasm collection existed in the within-population or cluster. Results of this study will be of great value for the genetic/breeding studies in garlic and enhance the utilization of these garlic germplasms.

Genetic analysis of novel resistance sources and genome-wide association mapping identified novel QTLs for resistance to Zymoseptoria tritici, the causal agent of septoria tritici blotch in wheat

Septoria tritici blotch (STB) caused by Zymoseptoria tritici is one of the most important foliar diseases of wheat causing significant yield losses worldwide. In this study, a panel of bread wheat genotypes comprised 185 globally diverse genotypes were tested against 10 Z. tritici isolates at the seedling stage. Genome-wide association study (GWAS) using high-throughput DArTseq markers was performed and further gene expression analysis of significant markers trait association (MTAs) associated with resistance to STB was analyzed. Disease severity level showed significant differences among wheat genotypes for resistance to different Z. tritici isolates. We found novel landrace genotypes that showed highly resistance spectra to all tested isolates. GWAS analysis resulted in 19 quantitative trait loci (QTLs) for resistance to STB that were located on 14 chromosomes. Overall, 14 QTLs were overlapped with previously known QTLs or resistance genes, as well as five potentially novel QTLs on chromosomes 1A, 4A, 5B, 5D, and 6D. Identified novel resistance sources and also novel QTLs for resistance to different Z. tritici isolates can be used for gene pyramiding and development of durable resistance cultivars in future wheat breeding programs.

Genetic Diversity and Population Structure of Cowpea [Vigna unguiculata (L.) Walp.] Germplasm Collected from Togo Based on DArT Markers

Crop genetic diversity is a sine qua non for continuous progress in the development of improved varieties, hence the need for germplasm collection, conservation and characterization. Over the years, cowpea has contributed immensely to the nutrition and economic life of the people in Togo. However, the bulk of varieties grown by farmers are landraces due to the absence of any serious genetic improvement activity on cowpea in the country. In this study, the genetic diversity and population structure of 255 cowpea accessions collected from five administrative regions and the agricultural research institute of Togo were assessed using 4600 informative diversity array technology (DArT) markers. Among the regions, the polymorphic information content (PIC) ranged from 0.19 to 0.27 with a mean value of 0.25. The expected heterozygosity (He) varied from 0.22 to 0.34 with a mean value of 0.31, while the observed heterozygosity (Ho) varied from 0.03 to 0.07 with an average of 0.05. The average inbreeding coefficient (FIS) varied from 0.78 to 0.89 with a mean value of 0.83, suggesting that most of the accessions are inbred. Cluster analysis and population structure identified four groups with each comprising accessions from the six different sources. Weak to moderate differentiation was observed among the populations with a genetic differentiation index varying from 0.014 to 0.117. Variation was highest (78%) among accessions within populations and lowest between populations (7%). These results revealed a moderate level of diversity among the Togo cowpea germplasm. The findings of this study constitute a foundation for genetic improvement of cowpea in Togo.

High-density DArT-based SilicoDArT and SNP markers for genetic diversity and population structure studies in cassava (Manihot esculenta Crantz)

Cassava (Manihot esculenta Crantz) is an important industrial and staple crop due to its high starch content, low input requirement, and resilience which makes it an ideal crop for sustainable agricultural systems and marginal lands in the tropics. However, the lack of genomic information on local genetic resources has impeded efficient conservation and improvement of the crop and the exploration of its full agronomic and breeding potential. This work was carried out to obtain information on population structure and extent of genetic variability among some local landraces conserved at the Plant Genetic Resources Research Institute, Ghana and exotic cassava accessions with Diversity Array Technology based SilicoDArT and SNP markers to infer how the relatedness in the genetic materials can be used to enhance germplasm curation and future breeding efforts. A total of 10521 SilicoDArT and 10808 SNP markers were used with varying polymorphic information content (PIC) values. The average PIC was 0.36 and 0.28 for the SilicoDArT and SNPs respectively. Population structure and average linkage hierarchical clustering based on SNPs revealed two distinct subpopulations and a large number of admixtures. Both DArT platforms identified 22 landraces as potential duplicates based on Gower's genetic dissimilarity. The expected heterozygosity which defines the genetic variation within each subpopulation was 0.008 for subpop1 which were mainly landraces and 0.391 for subpop2 indicating the homogeneous and admixture nature of the two subpopulations. Further analysis upon removal of the duplicates increased the expected heterozygosity of subpop1 from 0.008 to 0.357. A mantel test indicated strong interdependence (r = 0.970; P < 0.001) between SilicoDArT and DArTSeq SNP genotypic data suggesting both marker platforms as a robust system for genomic studies in cassava. These findings provide important information for efficient ex-situ conservation of cassava, future heterosis breeding, and marker-assisted selection (MAS) to enhance cassava improvement.

Physical Mapping of QTL in Four Spring Wheat Populations under Conventional and Organic Management Systems. I. Earliness

In previous studies, we reported quantitative trait loci (QTL) associated with the heading, flowering, and maturity time in four hard red spring wheat recombinant inbred line (RIL) populations but the results are scattered in population-specific genetic maps, which is challenging to exploit efficiently in breeding. Here, we mapped and characterized QTL associated with these three earliness traits using the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.0 physical map. Our data consisted of (i) 6526 single nucleotide polymorphisms (SNPs) and two traits evaluated at five conventionally managed environments in the 'Cutler' × 'AC Barrie' population; (ii) 3158 SNPs and two traits evaluated across three organic and seven conventional managements in the 'Attila' × 'CDC Go' population; (iii) 5731 SilicoDArT and SNP markers and the three traits evaluated at four conventional and organic management systems in the 'Peace' × 'Carberry' population; and (iv) 1058 SNPs and two traits evaluated across two conventionally and organically managed environments in the 'Peace' × 'CDC Stanley' population. Using composite interval mapping, the phenotypic data across all environments, and the IWGSC RefSeq v2.0 physical maps, we identified a total of 44 QTL associated with days to heading (11), flowering (10), and maturity (23). Fifteen of the 44 QTL were common to both conventional and organic management systems, and the remaining QTL were specific to either the conventional (21) or organic (8) management systems. Some QTL harbor known genes, including the Vrn-A1, Vrn-B1, Rht-A1, and Rht-B1 that regulate photoperiodism, flowering time, and plant height in wheat, which lays a solid basis for cloning and further characterization.

Comprehensive Metabolite Profiling in Genetic Resources of Garlic (Allium sativum L.) Collected from Different Geographical Regions

Garlic (Allium sativum) is the second most important Allium crop that has been used as a vegetable and condiment from ancient times due to its characteristic flavor and taste. Although garlic is a sterile plant that reproduces vegetatively through cloves, garlic shows high biodiversity, as well as phenotypic plasticity and environmental adaptation capacity. To determine the possible mechanism underlying this phenomenon and to provide new genetic materials for the development of a novel garlic cultivar with useful agronomic traits, the metabolic profiles in the leaf tissue of 30 garlic accessions collected from different geographical regions, with a special focus on the Asian region, were investigated using LC/MS. In addition, the total saponin and fructan contents in the roots and cloves of the investigated garlic accessions were also evaluated. Total saponin and fructan contents did not separate the garlic accessions based on their geographical origin, implying that saponin and fructan contents were clone-specific and agroclimatic changes have affected the quantitative and qualitative levels of saponins in garlic over a long history of cultivation. Principal component analysis (PCA) and dendrogram clustering of the LC/MS-based metabolite profiling showed two major clusters. Specifically, many Japanese and Central Asia accessions were grouped in cluster I and showed high accumulations of flavonol glucosides, alliin, and methiin. On the other hand, garlic accessions grouped in cluster II exhibited a high accumulation of anthocyanin glucosides and amino acids. Although most of the accessions were not separated based on country of origin, the Central Asia accessions were clustered in one group, implying that these accessions exhibited distinct metabolic profiles. The present study provides useful information that can be used for germplasm selection and the development of new garlic varieties with beneficial biotic and abiotic stress-adaptive traits.

Recent Large-Scale Genotyping and Phenotyping of Plant Genetic Resources of Vegetatively Propagated Crops

Several recent national and international projects have focused on large-scale genotyping of plant genetic resources in vegetatively propagated crops like fruit and berries, potatoes and woody ornamentals. The primary goal is usually to identify true-to-type plant material, detect possible synonyms, and investigate genetic diversity and relatedness among accessions. A secondary goal may be to create sustainable databases that can be utilized in research and breeding for several years ahead. Commonly applied DNA markers (like microsatellite DNA and SNPs) and next-generation sequencing each have their pros and cons for these purposes. Methods for large-scale phenotyping have lagged behind, which is unfortunate since many commercially important traits (yield, growth habit, storability, and disease resistance) are difficult to score. Nevertheless, the analysis of gene action and development of robust DNA markers depends on environmentally controlled screening of very large sets of plant material. Although more time-consuming, co-operative projects with broad-scale data collection are likely to produce more reliable results. In this review, we will describe some of the approaches taken in genotyping and/or phenotyping projects concerning a wide variety of vegetatively propagated crops.

Omics approaches in Allium research: Progress and way ahead

BACKGROUND

The genus Allium (Family: Amaryllidaceae) is an economically important group of crops cultivated worldwide for their use as a vegetable and spices. Alliums are also well known for their nutraceutical properties. Among alliums, onion, garlic, leek, and chives cultivated worldwide. Despite their substantial economic and medicinal importance, the genome sequence of any of the Allium is not available, probably due to their large genome sizes. Recently evolved omics technologies are highly efficient and robust in elucidating molecular mechanisms of several complex life processes in plants. Omics technologies, such as genomics, transcriptomics, proteomics, metabolomics, metagenomics, etc. have the potential to open new avenues in research and improvement of allium crops where genome sequence information is limited. A significant amount of data has been generated using these technologies for various Allium species; it will help in understanding the key traits in Allium crops such as flowering, bulb development, flavonoid biosynthesis, male sterility and stress tolerance at molecular and metabolite level. This information will ultimately assist us in speeding up the breeding in Allium crops.

METHOD

In the present review, major omics approaches, and their progress, as well as potential applications in Allium crops, could be discussed in detail.

RESULTS

Here, we have discussed the recent progress made in Allium research using omics technologies such as genomics, transcriptomics, micro RNAs, proteomics, metabolomics, and metagenomics. These omics interventions have been used in alliums for marker discovery, the study of the biotic and abiotic stress response, male sterility, organ development, flavonoid and bulb color, micro RNA discovery, and microbiome associated with Allium crops. Further, we also emphasized the integrated use of these omics platforms for a better understanding of the complex molecular mechanisms to speed up the breeding programs for better cultivars.

CONCLUSION

All the information and literature provided in the present review throws light on the progress and potential of omics platforms in the research of Allium crops. We also mentioned a few research areas in Allium crops that need to be explored using omics technologies to get more insight. Overall, alliums are an under-studied group of plants, and thus, there is tremendous scope and need for research in Allium species.

Insights into phylogeny, age and evolution of Allium (Amaryllidaceae) based on the whole plastome sequences

BACKGROUND AND AIMS

The genus Allium L., one of the largest monocotyledonous genera and one that includes many economically important crops with nutritional and medicinal value, has been the focus of classification or phylogeny studies for centuries. Recent studies suggested that the genus can be divided into 15 subgenera and 72 sections, which were further classified into three evolutionary lineages. However, the phylogenetic relationships reconstructed by one or two loci showed weaker support, especially for the third evolutionary lineage, which might not show the species relationships very clearly and could hinder further adaptive and evolutionary study.

METHODS

In this study, a total of 39 complete chloroplast genomes of Allium (covering 12 Allium subgenera) were collected, and combining these with 125 species of plastomes from 19 other families of monocots, we reconstructed the phylogeny of the genus Allium, estimated the origin and divergence time of the three evolutionary lineages and investigated the adaptive evolution in this genus and related families.

RESULTS

Our phylogenetic analysis confirmed the monophyly and three evolutionary lineages of Allium, while new species relationships were detected within the third evolutionary lineage. The divergence time of the three evolutionary lineages was estimated to be in the early Eocene to the middle Miocene, and numerous positive selected genes (PSGs) and PSGs with high average Ka/Ks values were found in Allium species.

CONCLUSIONS

Our results detected a well-supported phylogenetic relationship of Allium. The PSGs and PSGs with high Ka/Ks values, as well as diversified morphologies, complicated chromosome characteristics and unique reproductive modes may play important roles in the adaptation and evolution of Allium species. This is the first study that conducted phylogenetic and evolutionary analyses on the genus Allium combined with the plastome and morphological and cytological data. We hope that this study can contribute to further analysis of Allium for other researchers.

Discovery of thermophilic Bacillales using reduced-representation genotyping for identification

Background

This study demonstrates the use of reduced-representation genotyping to provide preliminary identifications for thermophilic bacterial isolates. The approach combines restriction enzyme digestion and PCR with next-generation sequencing to provide thousands of short-read sequences from across the bacterial genomes. Isolates were obtained from compost, hot water systems, and artesian bores of the Great Artesian Basin. Genomic DNA was double-digested with two combinations of restriction enzymes followed by PCR amplification, using a commercial provider of DArTseq™, Diversity Arrays Technology Pty Ltd. (Canberra, Australia). The resulting fragments which formed a reduced-representation of approximately 2.3% of the genome were sequenced. The sequence tags obtained were aligned against all available RefSeq bacterial genome assemblies by BLASTn to identify the nearest reference genome.

Results

Based on the preliminary identifications, a total of 99 bacterial isolates were identified to species level, from which 8 isolates were selected for whole-genome sequencing to assess the identification results. Novel species and strains were discovered within this set of isolates. The preliminary identifications obtained by reduced-representation genotyping, as well as identifications obtained by BLASTn alignment of the 16S rRNA gene sequence, were compared with those derived from the whole-genome sequence data, using the same RefSeq sequence database for the three methods. Identifications obtained with reduced-representation sequencing agreed with the identifications provided by whole-genome sequencing in 100% of cases. The identifications produced by BLASTn alignment of 16S rRNA gene sequence to the same database differed from those provided by whole-genome sequencing in 37.5% of cases, and produced ambiguous identifications in 50% of cases.

Conclusions

Previously, this method has been successfully demonstrated for use in bacterial identification for medical microbiology. This study demonstrates the first successful use of DArTseq™ for preliminary identification of thermophilic bacterial isolates, providing results in complete agreement with those obtained from whole-genome sequencing of the same isolates. The growing database of bacterial genome sequences provides an excellent resource for alignment of reduced-representation sequence data for identification purposes, and as the available sequenced genomes continue to grow, the technique will become more effective.

Multiple Small-Effect Alleles of Indica Origin Enhance High Iron-Associated Stress Tolerance in Rice Under Field Conditions in West Africa

Understanding the genetics of field-based tolerance to high iron-associated (HIA) stress in rice can accelerate the development of new varieties with enhanced yield performance in West African lowland ecosystems. To date, few field-based studies have been undertaken to rigorously evaluate rice yield performance under HIA stress conditions. In this study, two NERICA × O. sativa bi-parental rice populations and one O.sativa diversity panel consisting of 296 rice accessions were evaluated for grain yield and leaf bronzing symptoms over multiple years in four West African HIA stress and control sites. Mapping of these traits identified a large number of QTLs and single nucleotide polymorphisms (SNPs) associated with stress tolerance in the field. Favorable alleles associated with tolerance to high levels of iron in anaerobic rice soils were rare and almost exclusively derived from the indica subpopulation, including the most favorable alleles identified in NERICA varieties. These findings highlight the complex genetic architecture underlying rice response to HIA stress and suggest that a recurrent selection program focusing on an expanded indica genepool could be productively used in combination with genomic selection to increase the efficiency of selection in breeding programs designed to enhance tolerance to this prevalent abiotic stress in West Africa.

Single Nucleotide Polymorphisms (SNPs) Reveal Sibship Among Founders of a Bangladeshi Rohu (Labeo rohita) Breeding Population

Rohu (Labeo rohita) is a significant freshwater aquaculture species with approximately 1.8 Mt produced annually. Fin clips obtained from the founders of a newly established Bangladesh-based breeding population (∼140 fish from each of the Halda, Jamuna, and Padma rivers) were used to identify 9157 SNPs and 14 411 silicoDArT markers using the Diversity Arrays Technology (DArT) genotyping-by-sequencing platform known as DArTseq. After quality control, 1985 SNPs were retained and used to examine population structure within and among river systems. Examination of genomic relationships revealed evidence of full- and half-sibling relationships among founders. Accordingly, sibship and dummy parents were assigned within each river population using a maximum likelihood approach with COLONY software. Founders that had no dummy parents in common were then identified for population genetic analyses. Only 40 unique dummy parents and 17 founders with no common dummy parents were identified from the Halda river, compared with 206 (96) from the Jamuna and 184 (83) from the Padma. Overall pairwise FST estimates among rivers were low (<0.005) and the optimum number of clusters using unsupervised K-means clustering was one, indicating little genetic divergence among the river populations in our SNPs. These results suggest that observed sibship among founders should be accounted for in future pedigree-based analyses and it cannot be assumed that fertilized spawn collections are representative samples of river populations.

Identification of bacterial isolates from a public hospital in Australia using complexity-reduced genotyping

Bacterial identification methods used in routine identification of pathogens in medical microbiology include a combination approach of biochemical tests, mass spectrometry or molecular biology techniques. Extensive publicly-available databases of DNA sequence data from pathogenic bacteria have been amassed in recent years; this provides an opportunity for using bacterial genome sequencing for identification purposes. Whole genome sequencing is increasing in popularity, although at present it remains a relatively expensive approach to bacterial identification and typing. Complexity-reduced bacterial genome sequencing provides an alternative. We evaluate genomic complexity-reduction using restriction enzymes and sequencing to identify bacterial isolates. A total of 165 bacterial isolates from hospital patients in the Australian Capital Territory, between 2013 and 2015 were used in this study. They were identified and typed by the Microbiology Department of Canberra Public Hospital, and represented 14 bacterial species. DNA extractions from these samples were processed using a combination of the restriction enzymes PstI with MseI, PstI with HpaII and MseI with HpaII. The resulting sequences (length 30-69 bp) were aligned against publicly available bacterial genome and plasmid sequences. Results of the alignment were processed using a bioinformatics pipeline developed for this project, Currito3.1 DNA Fragment Analysis Software. All 165 samples were correctly identified to genus and species by each of the three combinations of restriction enzymes. A further 35 samples typed to the level of strain identified and compared for consistency with MLST typing data and in silico MLST data derived from the nearest sequenced candidate reference. The high level of agreement between bacterial identification using complexity-reduced genome sequencing and standard hospital identifications indicating that this new approach is a viable alternative for identification of bacterial isolates derived from pathology specimens. The effectiveness of species identification and in particular, strain typing, depends on access to a comprehensive and taxonomically accurate bacterial genome sequence database containing relevant bacterial species and strains.

Genome sequences of horticultural plants: past, present, and future

Horticultural plants play various and critical roles for humans by providing fruits, vegetables, materials for beverages, and herbal medicines and by acting as ornamentals. They have also shaped human art, culture, and environments and thereby have influenced the lifestyles of humans. With the advent of sequencing technologies, there has been a dramatic increase in the number of sequenced genomes of horticultural plant species in the past decade. The genomes of horticultural plants are highly diverse and complex, often with a high degree of heterozygosity and a high ploidy due to their long and complex history of evolution and domestication. Here we summarize the advances in the genome sequencing of horticultural plants, the reconstruction of pan-genomes, and the development of horticultural genome databases. We also discuss past, present, and future studies related to genome sequencing, data storage, data quality, data sharing, and data visualization to provide practical guidance for genomic studies of horticultural plants. Finally, we propose a horticultural plant genome project as well as the roadmap and technical details toward three goals of the project.

SNP discovery of Korean short day onion inbred lines using double digest restriction site-associated DNA sequencing

Onion (Allium cepa L.) is an economically important vegetable crop around the world. Genetic and genomic research into various onion accessions will provide insights into the onion genome to enhance breeding strategies and improve crops. However, the onion’s large genome size means that studies of molecular markers are limited in onion. This study aimed to discover high quality single nucleotide polymorphisms (SNPs) from 192 onion inbred lines relating to short-day cultivation in Korea. Paired-end (PE) double digested restriction site-associated DNA sequencing (ddRAD-seq) was used to discover SNPs in onion. A total of 538,973,706 reads (25.9 GB), with an average of 2,658,491 high-quality reads, were generated using ddRAD-seq. With stringent filtering, 1904 SNPs were discovered based on onion reference scaffolds. Further, population structure and genetic relationship studies suggested that two well-differentiated sub-populations exist in onion lines. SNP-associated flanking sequences were also compared with a public non-redundant database for gene ontology and pathway analysis. To our knowledge, this is the first report to identify high-quality SNPs in onion based on reference sequences using the ddRAD-seq platform. The SNP markers identified will be useful for breeders and the research community to deepen their understanding, enhance breeding programs, and support the management of onion genomic resources.

Genetic Variation and Population Structure of Oryza glaberrima and Development of a Mini-Core Collection Using DArTseq

The sequence variation present in accessions conserved in genebanks can best be used in plant improvement when it is properly characterized and published. Using low cost and high density single nucleotide polymorphism (SNP) assays, the genetic diversity, population structure, and relatedness between pairs of accessions can be quickly assessed. This information is relevant for different purposes, including creating core and mini-core sets that represent the maximum possible genetic variation contained in the whole collection. Here, we studied the genetic variation and population structure of 2,179 Oryza glaberrima Steud. accessions conserved at the AfricaRice genebank using 27,560 DArTseq-based SNPs. Only 14% (3,834 of 27,560) of the SNPs were polymorphic across the 2,179 accessions, which is much lower than diversity reported in other Oryza species. Genetic distance between pairs of accessions varied from 0.005 to 0.306, with 1.5% of the pairs nearly identical, 8.0% of the pairs similar, 78.1% of the pairs moderately distant, and 12.4% of the pairs very distant. The number of redundant accessions that contribute little or no new genetic variation to the O. glaberrima collection was very low. Using the maximum length sub-tree method, we propose a subset of 1,330 and 350 accessions to represent a core and mini-core collection, respectively. The core and mini-core sets accounted for ~61 and 16%, respectively, of the whole collection, and captured 97-99% of the SNP polymorphism and nearly all allele and genotype frequencies observed in the whole O. glaberrima collection available at the AfricaRice genebank. Cluster, principal component and model-based population structure analyses all divided the 2,179 accessions into five groups, based roughly on country of origin but less so on ecology. The first, third and fourth groups consisted of accessions primarily from Liberia, Nigeria, and Mali, respectively; the second group consisted primarily of accessions from Togo and Nigeria; and the fifth and smallest group was a mixture of accessions from multiple countries. Analysis of molecular variance showed between 10.8 and 28.9% of the variation among groups with the remaining 71.1-89.2% attributable to differences within groups.