tGBS® genotyping-by-sequencing enables reliable genotyping of heterozygous loci

Introduction to the Special Issue of Plants on "Advances in Plant Reproductive Ecology and Conservation Biology"

Plant reproductive ecology explores aspects of the biology and ecology of plants ranging from breeding systems, plant-pollinator interactions, seed germination, floral traits, and much more [...].

Genetic Diversity within a Collection of Italian Maize Inbred Lines: A Resource for Maize Genomics and Breeding

Genetic diversity is fundamental for studying the complex architecture of the traits of agronomic importance, controlled by major and minor loci. Moreover, well-characterized germplasm collections are essential tools for dissecting and analyzing genetic and phenotypic diversity in crops. A panel of 360 entries, a subset of a larger collection maintained within the GenBank at CREA Bergamo, which includes the inbreds derived from traditional Italian maize open-pollinated (OP) varieties and advanced breeding ones (Elite Inbreds), was analyzed to identify SNP markers using the tGBS® genotyping-by-sequencing technology. A total of 797,368 SNPs were found during the initial analysis. Imputation and filtering processes were carried out based on the percentage of missing data, redundant markers, and rarest allele frequencies, resulting in a final dataset of 15,872 SNP markers for which a physical map position was identified. Using this dataset, the inbred panel was characterized for linkage disequilibrium (LD), genetic diversity, population structure, and genetic relationships. LD decay at a genome-wide level indicates that the collection is a suitable resource for association mapping. Population structure analyses, which were carried out with different clustering methods, showed stable grouping statistics for four groups, broadly corresponding to 'Insubria', 'Microsperma', and 'Scagliolino' genotypes, with a fourth group composed prevalently of elite accessions derived from Italian and US breeding programs. Based on these results, the CREA Italian maize collection, genetically characterized in this study, can be considered an important tool for the mapping and characterization of useful traits and associated loci/alleles, to be used in maize breeding programs.

Novel Alleles from Cicer reticulatum L. for Genetic Improvement of Cultivated Chickpeas Identified through Genome Wide Association Analysis

The availability of wild chickpea (Cicer reticulatum L.) accessions has the potential to be used for the improvement of important traits in cultivated chickpeas. The main objectives of this study were to evaluate the phenotypic and genetic variations of chickpea progeny derived from interspecific crosses between C. arietinum and C. reticulatum, and to establish the association between single nucleotide polymorphism (SNP) markers and a series of important agronomic traits in chickpea. A total of 486 lines derived from interspecific crosses between C. arietinum (CDC Leader) and 20 accessions of C. reticulatum were evaluated at different locations in Saskatchewan, Canada in 2017 and 2018. Significant variations were observed for seed weight per plant, number of seeds per plant, thousand seed weight, and plant biomass. Path coefficient analysis showed significant positive direct effects of the number of seeds per plant, thousand seed weight, and biomass on the total seed weight. Cluster analysis based on the agronomic traits generated six groups that allowed the identification of potential heterotic groups within the interspecific lines for yield improvement and resistance to ascochyta blight disease. Genotyping of the 381 interspecific lines using a modified genotyping by sequencing (tGBS) generated a total of 14,591 SNPs. Neighbour-joining cluster analysis using the SNP data grouped the lines into 20 clusters. The genome wide association analysis identified 51 SNPs that had significant associations with different traits. Several candidate genes associated with early flowering and yield components were identified. The candidate genes and the significant SNP markers associated with different traits have a potential to aid the trait introgression in the breeding program.

Quantitative Trait Locus Mapping for Fire Blight Resistance in an F2 Population of Malus fusca MAL0045 Uncovers Novel Resistance Loci

Several fire blight resistance loci in Malus genotypes map on different linkage groups (LGs) representing chromosomes of the domesticated apple. Prior genetics studies primarily focused on F1 populations. A strong resistance quantitative trait locus (QTL) explained up to 66% of phenotypic variance in an F1 progeny derived from crossing the highly resistant wild apple genotype Malus fusca MAL0045 and the highly susceptible apple cultivar 'Idared', which was previously mapped on LG10 (Mfu10) of MAL0045. Strains of the causative bacterial pathogen Erwinia amylovora, notably those that show a single nucleotide polymorphism in the avrRpt2EA effector protein sequence at position 156 (e.g., Ea3049), are more virulent and overcome some known fire blight resistance donors and their QTLs. However, MAL0045 is resistant to Ea3049 and Mfu10 is not overcome, but most of the F1 progeny were highly susceptible to this strain. This phenomenon led to the assumption that other putative resistance factors not segregating in the F1 progeny might be present in the genome of MAL0045. Here, we crossed F1 progeny together to obtain 135 F2 individuals. Facilitated by genotyping-by-sequencing and phenotypic assessments, we identified and mapped two novel resistance QTLs in these F2 individuals on LGs 4 and 15, which were not identified in the F1. To our knowledge, these are the first resistance QTLs mapped in F2 progeny in Malus. In addition, we report that neither MAL0045 nor Mfu10 is broken down by a highly aggressive U.S. strain, LA635, after analyses in the original F1 individuals. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Genomic data of two chickpea populations sharing a potential Ascochyta blight resistance region

Chickpea (Cicer arietinum L.) is one of the most important crops worldwide and a valuable nutritional source. The availability of data from different genotypes and populations is important for the comprehension of the biology and trait control of chickpea. Tissue from young leaves was collected from adult plants and sequenced using an Illumina HiSeq X platform, which provided sequencing data for a total of 169 individuals from two different populations. Furthermore, functional annotation was performed with BLAST2GO software in a candidate region for resistance to Ascochyta blight, a devastating disease that produces huge yield reductions if the growth conditions are optimal for the fungus. A total of 273 different genes in a region spanning ∼4.67 Mb in chromosome 4 were fully annotated. The raw DNA sequences and functional annotation data can be reused by the scientific community for the analysis of different agronomic traits of interest in chickpea.

Genetic determinants of lettuce resistance to drop caused by Sclerotinia minor identified through genome-wide association mapping frequently co-locate with loci regulating anthocyanin content

KEY MESSAGE

GWAS identified 19 QTLs for resistance to Sclerotinia minor, 11 of them co-locating with red leaf color. Lower disease incidence was observed in red and dark red accessions. Lettuce (Lactuca sativa L.), one of the most economically important vegetables grown primarily in moderate climates around the world, is susceptible to many diseases including lettuce drop caused by the soilborne fungus Sclerotinia minor. Complete resistance to S. minor has not been identified in cultivated lettuce or its wild relatives. We conducted five experiments over 4 years with the diversity panel of almost 500 lettuce accessions to evaluate their response to the pathogen in an artificially infested field. The lowest disease incidence (DI) was observed in cultivars Eruption, Infantry, and Annapolis (median DI of 12.1-17.5%), while the highest DI was recorded for cultivars Reine des Glaces, Wayahead, and line FL. 43007 (median DI of 81.0-95.2%). Overall, significantly lower DI was observed in red and dark red accessions compared to those with a lower anthocyanin content. Genome-wide association mapping identified 19 QTLs for resistance to S. minor, 21 for the presence of red leaf color or its variations caused by the anthocyanin content, and one for the green color intensity. Eleven of the QTLs for disease resistance were located within 10 Mb of the loci associated with red color or anthocyanin content identified in this diversity panel. The frequent, non-random co-location of QTLs, together with the lower DI observed in red and dark red accessions suggests that lettuce interaction with S. minor may be partly influenced by anthocyanins. We have identified RLL2 and ANS, the genes of the anthocyanin biosynthesis pathway that co-locate with resistance QTLs, as candidates for functional studies to ascertain the involvement of anthocyanins in lettuce resistance against S. minor. Resistance QTLs closely linked with QTLs for anthocyanin content could be used to develop lettuce with a relatively high partial resistance and red color, while those not associated with anthocyanins could be used to develop partially resistant cultivars of green color.

Exome-based new allele-specific PCR markers and transferability for sodicity tolerance in bread wheat (Triticum aestivum L.)

Targeted exome-based genotype by sequencing (t-GBS), a sequencing technology that tags SNPs and haplotypes in gene-rich regions was used in previous genome-wide association studies (GWAS) for sodicity tolerance in bread wheat. Thirty-nine novel SNPs including 18 haplotypes for yield and yield-components were identified. The present study aimed at developing SNP-derived markers by precisely locating new SNPs on ~180 bp allelic sequence of t-GBS, marker validation, and SNP functional characterization based on its exonic location. We identified unknown locations of significant SNPs/haplotypes by aligning allelic sequences on to IWGSC RefSeqv1.0 on respective chromosomes. Eighteen out of the target 39 SNP locations fulfilled the criteria for producing PCR markers, among which only eight produced polymorphic signals. These eight markers associated with yield, plants m-2, heads m-2, and harvest index, including a pleiotropic marker for yield, harvest index, and grains/head were validated for its amplification efficiency and phenotypic effects in focused identification germplasm strategy (FIGS) wheat set and a doubled haploid (DH) population (Scepter/IG107116). The phenotypic variation explained by these markers are in the range of 4.1-37.6 in the FIGS population. High throughput PCR-based genotyping using new markers and association with phenotypes in FIGS wheat set and DH population validated the effect of functional SNP on closely associated genes-calcineurin B-like- and dirigent protein, basic helix-loop-helix (BHLH-), plant homeodomain (PHD-) and helix-turn-helix myeloblastosis (HTH myb) type -transcription factor. Further, genome-wide SNP annotation using SnpEff tool confirmed that these SNPs are in gene regulatory regions (upstream, 3'-UTR, and intron) modifying gene expression and protein-coding. This integrated approach of marker design for t-GBS alleles, SNP functional annotation, and high-throughput genotyping of functional SNP offers translation solutions across crops and complex traits in crop improvement programs.

Dataset on the single nucleotide variation in diversity panel of 500 lettuce accessions genotyped with tunable genotyping-by-sequencing (tGBS) method

Lettuce (Lactuca sativa L.) is an important leafy vegetable cultivated in moderate climates around the world. Based on phenotypic characteristics, such as formation and size of 'heads' (leaves arranged in a dense rosette), size, and texture of leaves, length of stems, and size and composition of seeds, lettuce cultivars can be classified into several distinct horticultural types (Batavia, butterhead, iceberg, Latin, leaf, oilseed, romaine, and stem). L. serriola, a wild progenitor of cultivated lettuce, is frequently utilized in breeding programs to introgress desirable genes and alleles (e.g., resistance to diseases) into lettuce gene pool. A diversity panel of ∼500 lettuce accessions was genotyped with tunable genotyping-by-sequencing (tGBS) method to identify single nucleotide polymorphism (SNP) sites. 115,261 SNPs were positioned on lettuce genome using the reference genome of cultivar Salinas. The described diversity panel together with the set of SNP markers can be used for mapping quantitative trait loci (QTL) and to develop marker assays for marker assisted selection (MAS). Identified SNP sites can also be used to identify F1 hybrids, genotype gene bank collections, and in other areas of lettuce genetics and breeding.

The Conservation Genetics of Iris lacustris (Dwarf Lake Iris), a Great Lakes Endemic

Iris lacustris, a northern Great Lakes endemic, is a rare species known from 165 occurrences across Lakes Michigan and Huron in the United States and Canada. Due to multiple factors, including habitat loss, lack of seed dispersal, patterns of reproduction, and forest succession, the species is threatened. Early population genetic studies using isozymes and allozymes recovered no to limited genetic variation within the species. To better explore genetic variation across the geographic range of I. lacustris and to identify units for conservation, we used tunable Genotyping-by-Sequencing (tGBS) with 171 individuals across 24 populations from Michigan and Wisconsin, and because the species is polyploid, we filtered the single nucleotide polymorphism (SNP) matrices using polyRAD to recognize diploid and tetraploid loci. Based on multiple population genetic approaches, we resolved three to four population clusters that are geographically structured across the range of the species. The species migrated from west to east across its geographic range, and minimal genetic exchange has occurred among populations. Four units for conservation are recognized, but nine adaptive units were identified, providing evidence for local adaptation across the geographic range of the species. Population genetic analyses with all, diploid, and tetraploid loci recovered similar results, which suggests that methods may be robust to variation in ploidy level.

The potentialities of omics resources for millet improvement

Millets are nutrient-rich (nutri-rich) cereals with climate resilience attributes. However, its full productive potential is not realized due to the lack of a focused yield improvement approach, as evidenced by the available literature. Also, the lack of well-characterized genomic resources significantly limits millet improvement. But the recent availability of genomic data and advancement in omics tools has shown its enormous potential to enhance the efficiency and precision faced by conventional breeding in millet improvement. The development of high throughput genotyping platforms based on next-generation sequencing (NGS) has provided a low-cost method for genomic information, specifically for neglected nutri-rich cereals with the availability of a limited number of reference genome sequences. NGS has created new avenues for millet biotechnological interventions such as mutation-based study, GWAS, GS, and other omics technologies. The simultaneous discovery of high-throughput markers and multiplexed genotyping platform has aggressively aided marker-assisted breeding for millet improvement. Therefore, omics technology offers excellent opportunities to explore and combine useful variations for targeted traits that could impart high nutritional value to high-yielding cultivars under changing climatic conditions. In millet improvement, an in-depth account of NGS, integrating genomics data with different biotechnology tools, is reviewed in this context.

Genome diversity and highland-adaptative variation in Tibet barley landrace population of China

Barley landraces accumulated variation in adapting to extreme highland environments during long-term domestication in Tibet, but little is known about their population structure and genomic selection traces. In this study, tGBS (tunable genotyping by sequencing) sequencing, molecular marker and phenotypic analyses were conducted on 1,308 highland and 58 inland barley landraces in China. The accessions were divided into six sub-populations and clearly distinguished most six-rowed, naked barley accessions (Qingke in Tibet) from inland barley. Genome-wide differentiation was observed in all five sub-populations of Qingke and inland barley accessions. High genetic differentiation in the pericentric regions of chromosomes 2H and 3H contributed to formation of five types of Qingke. Ten haplotypes of the pericentric regions of 2H, 3H, 6H and 7H were further identified as associated with ecological diversification of these sub-populations. There was genetic exchange between eastern and western Qingke but they shared the same progenitor. The identification of 20 inland barley types indicated multiple origins of Qingke in Tibet. The distribution of the five types of Qingke corresponded to specific environments. Two predominant highland-adaptative variations were identified for low temperature tolerance and grain color. Our results provide new insights into the origin, genome differentiation, population structure and highland adaptation in highland barley which will benefit both germplasm enhancement and breeding of naked barley.

HvGST plays a key role in anthocyanin accumulation in colored barley

Blue aleurone of barley is caused by the accumulation of delphinidin-based derivatives. Although these compounds are ideal nutrients for human health, they are undesirable contaminants in malt brewing. Therefore, the ability to add and remove this trait easily would facilitate breeding barley for different purposes. Here we identified a glutathione S-transferase gene (HvGST) that was responsible for the blue aleurone trait in Tibetan qingke barley by performing a genome-wide association study and RNA-sequencing analysis. Gene variation and expression analysis indicated that HvGST also participates in the transport and accumulation of anthocyanin in purple barley. Haplotype and the geographic distribution analyses of HvGST alleles revealed two independent natural variants responsible for the emergence of white aleurone: a 203-bp deletion causing premature termination of translation in qingke barley and two key single nucleotide polymorphisms in the promoter resulting in low transcription in Western barley. This study contributes to a better understanding of mechanisms of colored barley formation, and provides a comprehensive reference for marker-assisted barley breeding.

Identification of Quantitative Trait Loci Associated with Bacterial Leaf Spot Resistance in Baby Leaf Lettuce

Spring mix is a popular packaged salad that contains lettuce (Lactuca sativa L.) as one of its main ingredients. Plants for baby leaf lettuce (BLL) production are grown at very high densities, which enhances the occurrence of bacterial leaf spot (BLS) caused by Xanthomonas hortorum pv. vitians (Xhv), a disease that can make the crop unmarketable. The market demands disease-free, high-quality BLL all year round. Growing highly BLS-resistant cultivars will reduce loss of yield and quality, thus minimizing economic detriment to lettuce and spring mix growers. The research objectives were to identify lettuce accessions resistant to BLS and associated quantitative trait loci (QTL). A total of 495 lettuce accessions were screened with six isolates (BS0347, BS2861, BS3127, L7, L44, and Sc8B) of Xhv. Accessions showing overall high-level resistance to all tested Xhv isolates were 'Bunte Forellen', PI 226514, 'La Brillante', ARM09-161-10-1-4, 'Grenadier', 'Bella', PI 491210, 'Delight', and 'Romana Verde del Mercado'. Genome-wide association studies of BLS resistance by mixed linear model analyses identified significant QTLs on four lettuce chromosomes (2, 4, 6, and 8). The most significant QTL was on Chromosome 8 (P = 1.42 × 10-7), which explained 6.7% of total phenotypic variation for the disease severity. Accessions with a high level of resistance detected in this study are valuable resources for lettuce germplasm improvement. Molecular markers closely linked to QTLs can be considered for marker-assisted selection to develop new BLL lettuce cultivars with resistance to multiple races of Xhv.

Genotyping by Sequencing Advancements in Barley

Barley is considered an ideal crop to study cereal genetics due to its close relationship with wheat and diploid ancestral genome. It plays a crucial role in reducing risks to global food security posed by climate change. Genetic variations in the traits of interest in crops are vital for their improvement. DNA markers have been widely used to estimate these variations in populations. With the advancements in next-generation sequencing, breeders could access different types of genetic variations within different lines, with single-nucleotide polymorphisms (SNPs) being the most common type. However, genotyping barley with whole genome sequencing (WGS) is challenged by the higher cost and computational demand caused by the large genome size (5.5GB) and a high proportion of repetitive sequences (80%). Genotyping-by-sequencing (GBS) protocols based on restriction enzymes and target enrichment allow a cost-effective SNP discovery by reducing the genome complexity. In general, GBS has opened up new horizons for plant breeding and genetics. Though considered a reliable alternative to WGS, GBS also presents various computational difficulties, but GBS-specific pipelines are designed to overcome these challenges. Moreover, a robust design for GBS can facilitate the imputation to the WGS level of crops with high linkage disequilibrium. The complete exploitation of GBS advancements will pave the way to a better understanding of crop genetics and offer opportunities for the successful improvement of barley and its close relatives.

Genome-wide association mapping reveals genomic regions frequently associated with lettuce field resistance to downy mildew

KEY MESSAGE

GWAS identified 63 QTLs for resistance to downy mildew. Though QTLs were distributed across all chromosomes, the genomic regions frequently associated with resistance were located on chromosomes 4 and 5. Lettuce downy mildew is one of the most economically important diseases of cultivated lettuce worldwide. We have applied the genome-wide association mapping (GWAS) approach to detect QTLs for field resistance to downy mildew in the panel of 496 accessions tested in 21 field experiments. The analysis identified 131 significant marker-trait associations that could be grouped into 63 QTLs. At least 51 QTLs were novel, while remaining 12 QTLs overlapped with previously described QTLs for lettuce field resistance to downy mildew. Unlike race-specific, dominant Dm genes that mostly cluster on three out of nine lettuce chromosomes, QTLs (qDMR loci) for polygenic resistance are randomly distributed across all nine chromosomes. The genomic regions frequently associated with lettuce field resistance to downy mildew are located on chromosomes 4 and 5 and could be used for detailed study of the mechanism of polygenic resistance. The most resistant accessions identified in the current study (cvs. Auburn, Grand Rapids, Romabella, PI 226514, and PI 249536) are being incorporated into our breeding program. Markers closely linked to the resistance QTLs could be potentially used for marker-assisted selection, or in combination with other markers in the genome, for a combined genomic and marker-assisted selection. Up to date this is the most comprehensive study of QTLs for field resistance to downy mildew and the first study that uses GWAS for mapping disease resistance loci in lettuce.

Taxonomic re-appraisal for toothfish (Dissostichus: Notothenioidea) across the Antarctic Polar Front using genomic and morphological studies

The Patagonian toothfish, Dissostichus eleginoides, is one of the largest predatory fishes inhabiting Southern Ocean waters spanning the Antarctic Polar Front (APF), a prominent biogeographic boundary restricting gene flow and driving species divergence between Antarctic and sub-Antarctic waters. In the light of emerging threats to toothfish conservation and sustainability, this study investigated genetic [mtDNA sequences and genome wide nuclear single nucleotide polymorphisms (SNPs)] and morphological data to critically evaluate the taxonomic status of toothfish north (Chile and Patagonian shelf) and south (South Georgia and South Sandwich Islands) of the APF. mtDNA revealed reciprocally monophyletic lineages on either side of the APF with coalescent analysis indicating these diverged during the Pleistocene. Integration with data from other sources suggests the Chilean/Patagonian lineage is endemic. SNP analysis confirmed restricted nuclear gene flow between both groups and revealed a consensus suite of positive outlier SNPs compatible with adaptive divergence between these groups. Finally, several morphological features permit unequivocal assignment of individuals to either of the clades. Based on the genetic, phenotypic and ecological divergence, the authors propose that toothfish on either side of the APF be recognised as distinct species, with the name D. eleginoides used for toothfish occurring in South American waters north of the APF and toothfish south of the APF being classified using the new name D. australis reflecting their southern distribution.

SNP-Based Genome-Wide Association Mapping of Pollen Viability Under Heat Stress in Tropical Zea mays L. Inbred Lines

Global environmental changes with more extreme episodes of heat waves are major threats to agricultural productivity. Heat stress in spring affects the reproductive stage of maize, resulting in tassel blast, pollen abortion, poor pollination, reduced seed set, barren ears and ultimately yield loss. As an aneamophelous crop, maize has a propensity for pollen abortion under heat stress conditions. To overcome the existing challenges of heat stress and pollen abortion, this study utilized a broad genetic base of maize germplasm to identify superior alleles to be utilized in breeding programs. A panel of 375 inbred lines was morpho-physiologically screened under normal and heat stress conditions in two locations across two consecutive planting seasons, 2017 and 2018. The exposure of pollen to high temperature showed drastic decline in pollen germination percentage. The average pollen germination percentage (PGP) at 35 and 45°C was 40.3% and 9.7%, respectively, an average decline of 30.6%. A subset of 275 inbred lines were sequenced using tunable genotyping by sequencing, resulting in 170,098 single nucleotide polymorphisms (SNPs) after filtration. Genome wide association of PGP in a subset of 122 inbred lines resulted in ten SNPs associated with PGP35°C (p ≤ 10−5), nine with PGP45°C (p ≤ 10−6–10−8) and ten SNPs associated with PGP ratio (p ≤ 10−5). No SNPs were found to be in common across PGP traits. The number of favorable alleles possessed by each inbred line for PGP35°C, PGP45°C, and the PGP ratio ranged between 4 and 10, 3–13 and 5–13, respectively. In contrast, the number of negative alleles for these traits ranged between 2 and 8, 3–13 and 3–13, respectively. Genetic mapping of yield (adjusted weight per plant, AWP−1) and flowering time (anthesis-silking interval, ASI) in 275 lines revealed five common SNPs: three shared for AWP−1 between normal and heat stress conditions, one for ASI between conditions, and one SNP, CM007648.1-86615409, was associated with both ASI and AWP−1. Variety selection can be performed based on these favorable alleles for various traits. These marker trait associations identified in the diversity panel can be utilized in breeding programs to improve heat stress tolerance in maize.

OUP accepted manuscript

Substantial morphological variation in land plants remains inaccessible to genetic analysis because current models lack variation in important ecological and agronomic traits. The genus Gilia was historically a model for biosystematics studies and includes variation in morphological traits that are poorly understood at the genetic level. We assembled a chromosome-scale reference genome of G. yorkii and used it to investigate genome evolution in the Polemoniaceae. We performed QTL (quantitative trait loci) mapping in a G. yorkii×G. capitata interspecific population for traits related to inflorescence architecture and flower color. The genome assembly spans 2.75 Gb of the estimated 2.80-Gb genome, with 96.7% of the sequence contained in the nine largest chromosome-scale scaffolds matching the haploid chromosome number. Gilia yorkii experienced at least one round of whole-genome duplication shared with other Polemoniaceae after the eudicot paleohexaploidization event. We identified QTL linked to variation in inflorescence architecture and petal color, including a candidate for the major flower color QTL—a tandem duplication of flavanol 3′,5′-hydroxylase. Our results demonstrate the utility of Gilia as a forward genetic model for dissecting the evolution of development in plants including the causal loci underlying inflorescence architecture transitions.

Addressing Iron and Zinc Micronutrient Malnutrition Through Nutrigenomics in Pearl Millet: Advances and Prospects

Iron (Fe) and zinc (Zn) micronutrient deficiencies are significant health concerns, particularly among the underprivileged and resource-poor people in the semi-arid tropics globally. Pearl millet is regarded as a climate-smart crop with low water and energy footprints. It thrives well under adverse agro-ecologies such as high temperatures and limited rainfall. Pearl millet is regarded as a nutri-cereal owing to health-promoting traits such as high grain Fe and Zn content, metabolizable energy, high antioxidant and polyphenols, high proportion of slowly digestible starches, dietary fibers, and favorable essential amino acid profile compared to many cereals. Higher genetic variability for grain Fe and Zn content has facilitated considerable progress in mapping and mining QTLs, alleles and genes underlying micronutrient metabolism. This has been made possible by developing efficient genetic and genomic resources in pearl millet over the last decade. These include genetic stocks such as bi-parental RIL mapping populations, association mapping panels, chromosome segment substitution lines (CSSLs) and TILLING populations. On the genomics side, considerable progress has been made in generating genomic markers, such as SSR marker repository development. This was followed by the development of a next-generation sequencing-based genome-wide SNP repository. The circa 1,000 genomes re-sequencing project played a significant role. A high-quality reference genome was made available by re-sequencing of world diversity panel, mapping population parents and hybrid parental lines. This mini-review attempts to provide information on the current developments on mapping Fe and Zn content in pearl millet and future outlook.

Genetic diversity assessment of Georgia peanut cultivars developed during ninety years of breeding

The cultivated peanut(Arachis hypogaea L.) has experienced severe genetic bottlenecks over the course of its evolution and domestication. Most genetic diversity studies in peanut have focused on global genetic stocks, wild accessions, and related species, but few have focused on elite cultivars. The objective of this project was to assess the genetic diversity of 32 peanut cultivars developed by the University of Georgia breeding program since its inception in 1931. Quantifying genetic similarity (GSIM) among these cultivars will provide a better understanding of their relationships and aid in breeding for cultivar development. Genotyping-by-sequencing (GBS), in concert with the recently published A. hypogaea genome sequence, was used to identify a total of 27,142 single nucleotide polymorphisms (SNPs) among these cultivars. Coefficients of parentage (CoP) were calculated based on publicly available pedigree information and compared with SNP-based GSIM estimates; the resulting correlations were low, ranging from R² = 0.212-0.279. Although genetic diversity is generally low in cultivated peanut, our data indicate that the genetic diversity of Georgia cultivars has actually increased since the early days of the breeding program, likely a result of the incorporation of diverse germplasm and breeding lines into the program. The results reported here provide a valuable understanding of genetic variation among elite Georgia peanut cultivars that have had a significant impact on the peanut industry within the United States.

Population genetics and genome-wide association studies provide insights into the influence of selective breeding on genetic variation in lettuce

Genetic diversity is an important resource in crop breeding to improve cultivars with desirable traits. Selective breeding can lead to a reduction of genetic diversity. However, our understanding on this subject remains limited in lettuce (Lactuca sativa L.). Genotyping-by-sequencing (GBS) can provide a reduced version of the genome as a cost-effective method to identify genetic variants across the genome. We genotyped a diverse set of 441 lettuce accessions using the GBS method. Phylogenetic and population genetic analyses indicated substantial genetic divergence among four horticultural types of lettuce: butterhead, crisphead, leaf, and romaine. Genetic-diversity estimates between and within the four types indicated that the crisphead type was the most differentiated from other types, whereas its population was the most homogenous with the slowest linkage disequilibrium (LD) decay among the four types. These results suggested that crisphead lettuces had relatively less genetic variation across the genome as well as low gene flow from other types. We identified putative selective sweep regions that showed low genetic variation in the crisphead type. Genome-wide association study (GWAS) and quantitative trait loci (QTL) analyses provided evidence that these genomic regions were, in part, associated with delayed bolting, implicating the positive selection of delayed bolting in reducing variation. Our findings enhance the current understanding of genetic diversity and the impacts of selective breeding on patterning genetic variation in lettuce.

Genomic Prediction of Yield Traits in Single-Cross Hybrid Rice (Oryza sativa L.)

Hybrid rice varieties can outyield the best inbred varieties by 15 – 30% with appropriate management. However, hybrid rice requires more inputs and management than inbred rice to realize a yield advantage in high-yielding environments. The development of stress-tolerant hybrid rice with lowered input requirements could increase hybrid rice yield relative to production costs. We used genomic prediction to evaluate the combining abilities of 564 stress-tolerant lines used to develop Green Super Rice with 13 male sterile lines of the International Rice Research Institute for yield-related traits. We also evaluated the performance of their F₁ hybrids. We identified male sterile lines with good combining ability as well as F₁ hybrids with potential further use in product development. For yield per plant, accuracies of genomic predictions of hybrid genetic values ranged from 0.490 to 0.822 in cross-validation if neither parent or up to both parents were included in the training set, and both general and specific combining abilities were modeled. The accuracy of phenotypic selection for hybrid yield per plant was 0.682. The accuracy of genomic predictions of male GCA for yield per plant was 0.241, while the accuracy of phenotypic selection was 0.562. At the observed accuracies, genomic prediction of hybrid genetic value could allow improved identification of high-performing single crosses. In a reciprocal recurrent genomic selection program with an accelerated breeding cycle, observed male GCA genomic prediction accuracies would lead to similar rates of genetic gain as phenotypic selection. It is likely that prediction accuracies of male GCA could be improved further by targeted expansion of the training set. Additionally, we tested the correlation of parental genetic distance with mid-parent heterosis in the phenotyped hybrids. We found the average mid-parent heterosis for yield per plant to be consistent with existing literature values at 32.0%. In the overall population of study, parental genetic distance was significantly negatively correlated with mid-parent heterosis for yield per plant (r = −0.131) and potential yield (r = −0.092), but within female families the correlations were non-significant and near zero. As such, positive parental genetic distance was not reliably associated with positive mid-parent heterosis.

Chromosome-level genome assembly of a regenerable maize inbred line A188

BACKGROUND

The maize inbred line A188 is an attractive model for elucidation of gene function and improvement due to its high embryogenic capacity and many contrasting traits to the first maize reference genome, B73, and other elite lines. The lack of a genome assembly of A188 limits its use as a model for functional studies.

RESULTS

Here, we present a chromosome-level genome assembly of A188 using long reads and optical maps. Comparison of A188 with B73 using both whole-genome alignments and read depths from sequencing reads identify approximately 1.1 Gb of syntenic sequences as well as extensive structural variation, including a 1.8-Mb duplication containing the Gametophyte factor1 locus for unilateral cross-incompatibility, and six inversions of 0.7 Mb or greater. Increased copy number of carotenoid cleavage dioxygenase 1 (ccd1) in A188 is associated with elevated expression during seed development. High ccd1 expression in seeds together with low expression of yellow endosperm 1 (y1) reduces carotenoid accumulation, accounting for the white seed phenotype of A188. Furthermore, transcriptome and epigenome analyses reveal enhanced expression of defense pathways and altered DNA methylation patterns of the embryonic callus.

CONCLUSIONS

The A188 genome assembly provides a high-resolution sequence for a complex genome species and a foundational resource for analyses of genome variation and gene function in maize. The genome, in comparison to B73, contains extensive intra-species structural variations and other genetic differences. Expression and network analyses identify discrete profiles for embryonic callus and other tissues.

Meta-Analysis Identifies Pleiotropic Loci Controlling Phenotypic Trade-offs in Sorghum

Community association populations are composed of phenotypically and genetically diverse accessions. Once these populations are genotyped, the resulting marker data can be reused by different groups investigating the genetic basis of different traits. Because the same genotypes are observed and scored for a wide range of traits in different environments, these populations represent a unique resource to investigate pleiotropy. Here we assembled a set of 234 separate trait datasets for the Sorghum Association Panel, a group of 406 sorghum genotypes widely employed by the sorghum genetics community. Comparison of genome wide association studies conducted with two independently generated marker sets for this population demonstrate that existing genetic marker sets do not saturate the genome and likely capture only 35-43% of potentially detectable loci controlling variation for traits scored in this population. While limited evidence for pleiotropy was apparent in cross-GWAS comparisons, a multivariate adaptive shrinkage approach recovered both known pleiotropic effects of existing loci and new pleiotropic effects, particularly significant impacts of known dwarfing genes on root architecture. In addition, we identified new loci with pleiotropic effects consistent with known trade-offs in sorghum development. These results demonstrate the potential for mining existing trait datasets from widely used community association populations to enable new discoveries from existing trait datasets as new, denser genetic marker datasets are generated for existing community association populations.

Recent and Historical Gene Flow in Cultivars, Landraces, and a Wild Taxon of Cucurbita pepo in Mexico

Gene flow among crops and their wild relatives is an active study area in evolutionary biology and horticulture, because genetic exchange between them may impact their evolutionary trajectories and increase the genetic variation of the cultivated lineages. Mexico is a center of diversity for the genus Cucurbita that includes pumpkins, squash and gourds. Gene flow between domesticated and wild species has been reported as common in Cucurbita; but gene flow among populations of C. pepo ssp. pepo from Mexico and its wild relative has not been studied. We used 2,061 SNPs, derived from tunable genotyping by sequencing (tGBS) to estimate gene flow among 14 Mexican traditional landraces of C. pepo ssp. pepo, also including individuals from five improved cultivars of C. pepo ssp. pepo and C. pepo ssp. ovifera var. ovifera, and individuals of their wild relative C. pepo ssp. fraterna. We found moderate to high levels of genetic diversity, and low to moderate genetic differentiation. In the test of introgression between lineages, we found that all possible arrangements for ancestral and derived sites between the lineages showed similar frequencies; thus, incomplete lineage sorting, but also gene flow, might be taking place in C. pepo. Overall, our results suggest that gene flow between these subspecies and cultigens, incomplete lineage sorting and the retention of ancestral characters shaped the evolutionary trajectory of C. pepo in its area of origin and diversification. In addition, we found evidence of the use of Mexican landraces as genetic material for the improvement of commercial cultivars. The landraces of Mexico are an important source of genetic diversity for C. pepo, which has been preserved both by management practices of small farmers and by the natural gene flow that exists between the different crop fields of the region.

The domestication of Cucurbita argyrosperma as revealed by the genome of its wild relative

Despite their economic importance and well-characterized domestication syndrome, the genomic impact of domestication and the identification of variants underlying the domestication traits in Cucurbita species (pumpkins and squashes) is currently lacking. Cucurbita argyrosperma, also known as cushaw pumpkin or silver-seed gourd, is a Mexican crop consumed primarily for its seeds rather than fruit flesh. This makes it a good model to study Cucurbita domestication, as seeds were an essential component of early Mesoamerican diet and likely the first targets of human-guided selection in pumpkins and squashes. We obtained population-level data using tunable Genotype by Sequencing libraries for 192 individuals of the wild and domesticated subspecies of C. argyrosperma across Mexico. We also assembled the first high-quality wild Cucurbita genome. Comparative genomic analyses revealed several structural variants and presence/absence of genes related to domestication. Our results indicate a monophyletic origin of this domesticated crop in the lowlands of Jalisco. We found evidence of gene flow between the domesticated and wild subspecies, which likely alleviated the effects of the domestication bottleneck. We uncovered candidate domestication genes that are involved in the regulation of growth hormones, plant defense mechanisms, seed development, and germination. The presence of shared selected alleles with the closely related species Cucurbita moschata suggests domestication-related introgression between both taxa.

Genotyping by Sequencing Reveals Genetic Relatedness of Southwestern U.S. Blue Maize Landraces

Maize has played a key role in the sustenance and cultural traditions of the inhabitants of the southwestern USA for many centuries. Blue maize is an important component of the diverse landraces still cultivated in the region but the degree to which they are related is unknown. This research was designed to ascertain the genotypic, morphological, and phenotypic diversity of six representative southwestern blue maize landraces. Their genotypic diversity was examined using tunable genotyping-by-sequencing (tGBS™). A total of 81,038 high quality SNPs were identified and obtained through tGBS. A total of 45 morphological and biochemical traits were evaluated at two locations in New Mexico. The varieties Los Lunas High and Flor del Rio were genetically less related with other southwestern landraces whereas diffusion between Navajo Blue, Hopi Blue, Yoeme Blue, and Taos Blue demonstrated that these landraces were genetically related. Phenotypic variability was highest for kernel traits and least for plant traits. Plant, ear, and kernel traits were fairly consistent within and across locations. Principal component analysis and tGBS showed that Corn Belt variety ‘Ohio Blue’ was distinctly different from southwestern landraces. Genotypic analysis displayed that southwestern landraces are genetically closely related, but selection has resulted in differing phenotypes. This study has provided additional insight into the genetic relatedness of southwestern blue maize landraces.

Genome-wide analyses reveal footprints of divergent selection and popping-related traits in CIMMYT's maize inbred lines

Popcorn (Zea mays L. var. Everta) is the most ancient type of cultivated maize. However, there is little known about the genetics of popping-related traits based on genotyping-by-sequencing (GBS) technology. Here, we characterized the phenotypic variation for seven popping-related traits in maize kernels among 526 CIMMYT inbred lines (CMLs). In total, 155 083 high-quality single nucleotide polymorphism (SNP) markers were identified by a GBS approach. Several trait-associated loci were detected by genome-wide association study for color, popping expansion volume, shape, pericarp, flotation index, floury/vitreous, and protein content, explaining a majority of the observed phenotypic variance, and these were validated by a diverse panel comprising 764 tropical landrace accessions. Sixty two of the identified loci were recognized to have undergone selection. On average, there was a 55.27% frequency for alleles that promote popping in CMLs. Our work not only pinpoints previously unknown loci for popping-related traits, but also reveals that many of these loci have undergone selection. Beyond establishing a new benchmark for the genetics of popcorn, our study provides a foundation for gene discovery and breeding. It also presents evidence to investigate the role of a gradual loss of popping ability as a by-product of diversification of culinary uses throughout the evolution of teosinte-to-modern maize.

Characterization, Genetic Analyses, and Identification of QTLs Conferring Metabolic Resistance to a 4-Hydroxyphenylpyruvate Dioxygenase Inhibitor in Sorghum (Sorghum bicolor)

Postemergence grass weed control continues to be a major challenge in grain sorghum [Sorghum bicolor (L.) Moench], primarily due to lack of herbicide options registered for use in this crop. The development of herbicide-resistant sorghum technology to facilitate broad-spectrum postemergence weed control can be an economical and viable solution. The 4-hydroxyphenylpyruvate dioxygenase-inhibitor herbicides (e.g., mesotrione or tembotrione) can control a broad spectrum of weeds including grasses, which, however, are not registered for postemergence application in sorghum due to crop injury. In this study, we identified two tembotrione-resistant sorghum genotypes (G-200, G-350) and one susceptible genotype (S-1) by screening 317 sorghum lines from a sorghum association panel (SAP). These tembotrione-resistant and tembotrione-susceptible genotypes were evaluated in a tembotrione dose-response [0, 5.75, 11.5, 23, 46, 92 (label recommended dose), 184, 368, and 736 g ai ha^-1] assay. Compared with S-1, the genotypes G-200 and G-350 exhibited 10- and seven fold more resistance to tembotrione, respectively. To understand the inheritance of tembotrione-resistant trait, crosses were performed using S-1 and G-200 or G-350 to generate F₁ and F₂ progeny. The F₁ and F₂ progeny were assessed for their response to tembotrione treatment. Genetic analyses of the F₁ and F₂ progeny demonstrated that the tembotrione resistance in G-200 and G-350 is a partially dominant polygenic trait. Furthermore, cytochrome P450 (CYP)-inhibitor assay using malathion and piperonyl butoxide suggested possible CYP-mediated metabolism of tembotrione in G-200 and G-350. Genotype-by-sequencing based quantitative trait loci (QTL) mapping revealed QTLs associated with tembotrione resistance in G-200 and G-350 genotypes. Overall, the genotypes G-200 and G-350 confer a high level of metabolic resistance to tembotrione and controlled by a polygenic trait. There is an enormous potential to introgress the tembotrione resistance into breeding lines to develop agronomically desirable sorghum hybrids.

Identification of fertility-related genes for maize CMS-S via Bulked Segregant RNA-Seq

Cytoplasmic male sterility (CMS) is extensively used in maize hybrid production, and identification of genes related to fertility restoration for CMS is important for hybrid breeding. The fertility restoration of S type CMS is governed by several loci with major and minor effects, while the mechanism of fertility restoration for CMS-S is still unknown. In this study, BSR-Seq was conducted with two backcrossing populations with the fertility restoration genes, Rf3 and Rf10, respectively. Genetic mapping via BSR-Seq verified the positions of the two loci. A total of 353 and 176 differentially expressed genes (DEGs) between the male fertility and male sterile pools were identified in the populations with Rf3 and Rf10, respectively. In total, 265 DEGs were co-expressed in the two populations, which were up-regulated in the fertile plants, and they might be related to male fertility involving in anther or pollen development. Moreover, 35 and seven DEGs were specifically up-regulated in the fertile plants of the population with Rf3 and Rf10, respectively. Function analysis of these DEGs revealed that jasmonic acid (JA) signal pathway might be involved in the Rf3 mediated fertility restoration for CMS-S, while the small ubiquitin-related modifier system could play a role in the fertility restoration of Rf10.

Understanding Host-Pathogen Interactions in Brassica napus in the Omics Era

Brassica napus (canola/oilseed rape/rapeseed) is an economically important crop, mostly found in temperate and sub-tropical regions, that is cultivated widely for its edible oil. Major diseases of Brassica crops such as Blackleg, Clubroot, Sclerotinia Stem Rot, Downy Mildew, Alternaria Leaf Spot and White Rust have caused significant yield and economic losses in rapeseed-producing countries worldwide, exacerbated by global climate change, and, if not remedied effectively, will threaten global food security. To gain further insights into the host-pathogen interactions in relation to Brassica diseases, it is critical that we review current knowledge in this area and discuss how omics technologies can offer promising results and help to push boundaries in our understanding of the resistance mechanisms. Omics technologies, such as genomics, proteomics, transcriptomics and metabolomics approaches, allow us to understand the host and pathogen, as well as the interaction between the two species at a deeper level. With these integrated data in multi-omics and systems biology, we are able to breed high-quality disease-resistant Brassica crops in a more holistic, targeted and accurate way.

Construction of a dense genetic map of the Malus fusca fire blight resistant accession MAL0045 using tunable genotyping-by-sequencing SNPs and microsatellites

Although, the Pacific crabapple, Malus fusca, is a hardy and disease resistant species, studies relating to the genetics of its unique traits are very limited partly due to the lack of a genetic map of this interesting wild apple. An accession of M. fusca (MAL0045) of Julius Kühn-Institut collection in Germany is highly resistant to fire blight disease, incited by different strains of the causative pathogen—Erwinia amylovora. This is the most destructive bacterial disease of Malus of which most of the domesticated apples (Malus domestica) are susceptible. Using a scarcely dense genetic map derived from a population of 134 individuals of MAL0045 × ‘Idared’, the locus (Mfu10) controlling fire blight resistance mapped on linkage group 10 (LG10) and explained up to 66% of the phenotypic variance with different strains. Although the development of robust and tightly linked molecular markers on LG10 through chromosome walking approach led to the identification of a major candidate gene, any minor effect locus remained elusive possibly due to the lack of marker density of the entire genetic map. Therefore, we have developed a dense genetic map of M. fusca using tunable genotyping-by-sequencing (tGBS) approach. Of thousands of de novo SNPs identified, 2677 were informative in M. fusca and 90.5% of these successfully mapped. In addition, integration of SNP data and microsatellite (SSR) data resulted in a final map comprising 17 LGs with 613 loci spanning 1081.35 centi Morgan (cM). This map will serve as a template for mapping using different strains of the pathogen.

Applications of genotyping-by-sequencing (GBS) in maize genetics and breeding

Genotyping-by-Sequencing (GBS) is a low-cost, high-throughput genotyping method that relies on restriction enzymes to reduce genome complexity. GBS is being widely used for various genetic and breeding applications. In the present study, 2240 individuals from eight maize populations, including two association populations (AM), backcross first generation (BC1), BC1F2, F2, double haploid (DH), intermated B73 × Mo17 (IBM), and a recombinant inbred line (RIL) population, were genotyped using GBS. A total of 955,120 of raw data for SNPs was obtained for each individual, with an average genotyping error of 0.70%. The rate of missing genotypic data for these SNPs was related to the level of multiplex sequencing: ~ 25% missing data for 96-plex and ~ 55% for 384-plex. Imputation can greatly reduce the rate of missing genotypes to 12.65% and 3.72% for AM populations and bi-parental populations, respectively, although it increases total genotyping error. For analysis of genetic diversity and linkage mapping, unimputed data with a low rate of genotyping error is beneficial, whereas, for association mapping, imputed data would result in higher marker density and would improve map resolution. Because imputation does not influence the prediction accuracy, both unimputed and imputed data can be used for genomic prediction. In summary, GBS is a versatile and efficient SNP discovery approach for homozygous materials and can be effectively applied for various purposes in maize genetics and breeding.

Identifying barley pan-genome sequence anchors using genetic mapping and machine learning

We identified 1.844 million barley pan-genome sequence anchors from 12,306 genotypes using genetic mapping and machine learning. There is increasing evidence that genes from a given crop genotype are far to cover all genes in that species; thus, building more comprehensive pan-genomes is of great importance in genetic research and breeding. Obtaining a thousand-genotype scale pan-genome using deep-sequencing data is currently impractical for species like barley which has a huge and highly repetitive genome. To this end, we attempted to identify barley pan-genome sequence anchors from a large quantity of genotype-by-sequencing (GBS) datasets by combining genetic mapping and machine learning algorithms. Based on the GBS sequences from 11,166 domesticated and 1140 wild barley genotypes, we identified 1.844 million pan-genome sequence anchors. Of them, 532,253 were identified as presence/absence variation (PAV) tags. Through aligning these PAV tags to the genome of hulless barley genotype Zangqing320, our analysis resulted in a validation of 83.6% of them from the domesticated genotypes and 88.6% from the wild barley genotypes. Association analyses against flowering time, plant height and kernel size showed that the relative importance of the PAV and non-PAV tags varied for different traits. The pan-genome sequence anchors based on GBS tags can facilitate the construction of a comprehensive pan-genome and greatly assist various genetic studies including identification of structural variation, genetic mapping and breeding in barley.

Characterization of Genetic Diversity and Genome-Wide Association Mapping of Three Agronomic Traits in Qingke Barley (Hordeum Vulgare L.) in the Qinghai-Tibet Plateau

Barley (Hordeum vulgare L.) is one of the most important cereal crops worldwide. In the Qinghai-Tibet Plateau, six-rowed hulless (or naked) barley, called “qingke” in Chinese or “nas” in Tibetan, is produced mainly in Tibet. The complexity of the environment in the Qinghai-Tibet Plateau has provided unique opportunities for research on the breeding and adaptability of qingke barley. However, the genetic architecture of many important agronomic traits for qingke barley remains elusive. Heading date (HD), plant height (PH), and spike length (SL) are three prominent agronomic traits in barley. Here, we used genome-wide association (GWAS) mapping and GWAS with eigenvector decomposition (EigenGWAS) to detect quantitative trait loci (QTL) and selective signatures for HD, PH, and SL in a collection of 308 qingke barley accessions. The accessions were genotyped using a newly-developed, proprietary genotyping-by-sequencing (tGBS) technology, that yielded 14,970 high quality single nucleotide polymorphisms (SNPs). We found that the number of SNPs was higher in the varieties than in the landraces, which suggested that Tibetan varieties and varieties in the Tibetan area may have originated from different landraces in different areas. We have identified 62 QTLs associated with three important traits, and the observed phenotypic variation is well-explained by the identified QTLs. We mapped 114 known genes that include, but are not limited to, vernalization, and photoperiod genes. We found that 83.87% of the identified QTLs are located in the non-coding regulatory regions of annotated barley genes. Forty-eight of the QTLs are first reported here, 28 QTLs have pleotropic effects, and three QTL are located in the regions of the well-characterized genes HvVRN1, HvVRN3, and PpD-H2. EigenGWAS analysis revealed that multiple heading-date-related loci bear signatures of selection. Our results confirm that the barley panel used in this study is highly diverse, and showed a great promise for identifying the genetic basis of adaptive traits. This study should increase our understanding of complex traits in qingke barley, and should facilitate genome-assisted breeding for qingke barley improvement.

Genetic Analysis of Domestication Parallels in Annual and Perennial Sunflowers (Helianthus spp.): Routes to Crop Development

Parallels exist between the domestication of new species and the improvement of various crops through selection on traits which favor the sowing, harvest and retention of yield potential and the directed efforts to improve their agronomics, disease resistance and quality characteristics. Common selection pressures may result in the parallel selection of orthologs underlying these traits and homologies between crop species can be exploited by plant breeders to improve germplasm. Perennial grains and oilseeds are a class of proposed crops for improving the diversity and sustainability of agricultural systems. Maximilian sunflower (Helianthus maximiliani Schrad.) is a perennial crop wild relative of sunflower (Helianthus annuus L.) and a candidate perennial oilseed species. Understanding parallels between cultivated H. annuus and H. maximiliani may provide new tools for the development of Maximilian sunflower and other wild relatives of sunflower as crops to enhance functional diversity in cropping systems. F₂ populations of Maximilian sunflower segregating for traits associated with the domestication ideotype of cultivated sunflower including branching architecture, capitulum morphology and flowering time were developed to investigate parallels between H. maximiliani and H. annuus. Genotype-by-sequencing (GBS) was employed to genotype novel Maximilian sunflower populations and perform quantitative-trait-loci (QTL) analysis. A total of 11 QTL in five regions were identified across 21 linkage groups using 4142 GBS derived single nucleotide polymorphism markers called using the sunflower reference genome as a guide. A major QTL on linkage group 17b, associated with aspects of floral development and apical dominance, was discovered and corresponds with a known domestication QTL hotspot in H. annuus and candidate genes were identified. This suggests the potential to exploit orthologs for neo-domestication of H. maximiliani for traits such as branching architecture, timing of anthesis, and capitulum size and morphology for the development of a perennial oilseed crop from wild relatives of cultivated sunflower.

Genome-wide association mapping reveals loci for shelf life and developmental rate of lettuce

KEY MESSAGE

Two major QTL, one for shelf life that corresponds to qSL4 and one, qDEV7, for developmental rate, were identified. Associated markers will be useful in breeding for improved fresh-cut lettuce. Fresh-cut lettuce in packaged salad can have short shelf life, and visible deterioration may start within a week after processing. Yield and developmental rate are an important aspect of lettuce production. Genetic diversity and genome-wide association studies (GWAS) were performed on 493 accessions with the genotypic data of 4615 high-quality single nucleotide polymorphism markers. Population structure (Q), principal component (PC), and phylogenetic analyses displayed genetic relationships associated with lettuce types and geographic distribution. Data for shelf life, yield, developmental rate, and their stability indices were used for statistical analysis, and GWAS was performed by general and mixed linear models. The genetic relationship among the individuals was incorporated into the models using kinship matrix, PC, and Q. Broad-sense heritability (H2) across environments was 0.43 for shelf life, 0.36 for yield, and 0.60 for developmental rate. There was a negative correlation between yield and developmental rate. Significant marker-trait association (SMTA) was detected for shelf life on chromosome 4. The most significant quantitative trait locus (QTL,  qSL4, P = 2.23E-17) explained 24% of the total phenotypic variation (R2). The major QTL for developmental rate was detected on chromosome 7 (qDEV7, P = 2.43E-16, R2 = 17%), while additional QTLs with smaller effect were found in all chromosomes. No SMTA was detected for yield. The study identified lettuce accessions with extended and stable shelf life, stable yield, and desirable developmental rate. Molecular markers closely linked to traits can be applied for selection of preferable genotypes and for identification of genes associated with these traits.

A Domestication-Associated Gene GmPRR3b Regulates the Circadian Clock and Flowering Time in Soybean

Improved soybean cultivars have been adapted to grow at a wide range of latitudes, enabling expansion of cultivation worldwide. However, the genetic basis of this broad adaptation is still not clear. Here, we report the identification of GmPRR3b as a major flowering time regulatory gene that has been selected during domestication and genetic improvement for geographic expansion. Through a genome-wide association study of a diverse soybean landrace panel consisting of 279 accessions, we identified 16 candidate quantitative loci associated with flowering time and maturity time. The strongest signal resides in the known flowering gene E2, verifying the effectiveness of our approach. We detected strong signals associated with both flowering and maturity time in a genomic region containing GmPRR3b. Haplotype analysis revealed that GmPRR3bH6 is the major form of GmPRR3b that has been utilized during recent breeding of modern cultivars. mRNA profiling analysis showed that GmPRR3bH6 displays rhythmic and photoperiod-dependent expression and is preferentially induced under long-day conditions. Overexpression of GmPRR3bH6 increased main stem node number and yield, while knockout of GmPRR3bH6 using CRISPR/Cas9 technology delayed growth and the floral transition. GmPRR3bH6 appears to act as a transcriptional repressor of multiple predicted circadian clock genes, including GmCCA1a, which directly upregulates J/GmELF3a to modulate flowering time. The causal SNP (Chr12:5520945) likely endows GmPRR3bH6 a moderate but appropriate level of activity, leading to early flowering and vigorous growth traits preferentially selected during broad adaptation of landraces and improvement of cultivars.

Genetic Dissection of Phomopsis Stem Canker Resistance in Cultivated Sunflower Using High Density SNP Linkage Map

Phomopsis stem canker (PSC) caused by Diaporthe helianthi is increasingly becoming a global threat for sunflower production. In this study, the genetic basis of PSC resistance was investigated in a recombinant inbred line (RIL) population developed from a cross between HA 89 (susceptible) and HA-R3 (resistant). The RIL population was evaluated for PSC disease incidence (DI) in seven screening trials at multiple locations during 2016–2018. The distribution of PSC DI in the RIL population was continuous, confirming a polygenic inheritance of the trait. A moderately high broad-sense heritability (H², 0.76) was estimated for the trait across environments. In the combined analysis, both the genotype and the genotype × environment interactions were highly significant. A linkage map spanning 1505.33 cM was constructed using genotyping-by-sequencing derived markers. Marker–trait association analysis identified a total of 15 quantitative trait loci (QTL) associated with PSC resistance on 11 sunflower chromosomes, each explaining between 5.24 and 17.39% of the phenotypic variation. PSC resistance QTL were detected in two genomic regions each on chromosomes 3, 5, 13, and 17, while one QTL each was detected in the remaining seven chromosomes. Tightly linked single nucleotide polymorphism (SNP) markers flanking the PSC resistance QTL will facilitate marker-assisted selection in PSC resistance sunflower breeding.

Evaluation of Brassica oleracea accessions for resistance to Plasmodiophora brassicae and identification of genomic regions associated with resistance

Clubroot disease caused by Plasmodiophora brassicae is a challenge to Brassica crop production. Breakdown of resistance controlled by major genes of the Brassica A genome has been reported. Therefore, identification of resistance in the Brassica C genome is needed to broaden the genetic base of resistance in Brassica napus canola. In this study, we evaluated 135 Brassica oleracea accessions, belonging to eight variants of this species to identify resistant accessions as well as to identify the genomic regions associated with resistance to two recently evolved P. brassicae pathotypes, F3-14 (3A) and F-359-13 (5X L-G2). Resistance to these pathotypes was observed more frequently in var. acephala (kale) followed by var. capitata (cabbage); few accessions also carried resistance to both pathotypes. Association mapping using single nucleotide polymorphism (SNP) markers developed through genotyping by sequencing technique identified 10 quantitative trait loci (QTL) from six C-genome chromosomes to be associated with resistance to these pathotypes; among these, two QTL associated with resistance to 3A and one QTL associated with resistance to 5X L-G2 carried ≥3 SNP markers. The 10 QTL identified in this study individually accounted for 8%–18% of the total phenotypic variance. Thus, the results from this study can be used in molecular breeding of Brassica crops for resistance to this disease.

Construction of High-Density Genetic Map and Mapping of Sex-Related Loci in the Yellow Catfish (Pelteobagrus fulvidraco)

The yellow catfish (Pelteobagrus fulvidraco) is a very important aquaculture species distributed in freshwater area of China. All-male yellow catfish is very popular in aquaculture because of their significant sex dimorphism phenomena. The males grow much faster than females in full-sibling family. However, the sex dimorphism mechanism is still unclear in yellow catfish. In order to better understand the genetic basis of yellow catfish sexual dimorphism, it is vital to map the sex-related traits and localize the candidate genes across yellow catfish whole genome. Here, we constructed a high-density linkage map of yellow catfish using genotyping-by-sequencing (GBS) strategy. A total of 5705 single-nucleotide polymorphism (SNP) markers were mapped to 26 different linkage groups (LGs) using 184 F1 offspring. The total genetic map length was 3071.59 cM, with an average interlocus distance of 0.54 cM. Eleven significant sex-related QTLs in yellow catfish were identified. Six sex-related genes were identified from the region of reference genome near these QTLs including amh, gnrhr, vasa, lnnr1, foxl2, and bmp15. The high-density genetic linkage map provides valuable resources for yellow catfish molecular assistant breeding and elucidating sex differentiation process. Moreover, the comparative genomic study was analyzed among yellow catfish, channel catfish, and zebrafish. It revealed highly conserved chromosomal distribution between yellow catfish and channel catfish.

Identification of loci controlling adaptation in Chinese soya bean landraces via a combination of conventional and bioclimatic GWAS

Landraces often contain genetic diversity that has been lost in modern cultivars, including alleles that confer enhanced local adaptation. To comprehensively identify loci associated with adaptive traits in soya bean landraces, for example flowering time, a population of 1938 diverse landraces and 97 accessions of the wild progenitor of cultivated soya bean, Glycine soja was genotyped using tGBS® . Based on 99 085 high-quality SNPs, landraces were classified into three sub-populations which exhibit geographical genetic differentiation. Clustering was inferred from STRUCTURE, principal component analyses and neighbour-joining tree analyses. Using phenotypic data collected at two locations separated by 10 degrees of latitude, 17 trait-associated SNPs (TASs) for flowering time were identified, including a stable locus Chr12:5914898 and previously undetected candidate QTL/genes for flowering time in the vicinity of the previously cloned flowering genes, E1 and E2. Using passport data associated with the collection sites of the landraces, 27 SNPs associated with adaptation to three bioclimatic variables (temperature, daylength, and precipitation) were identified. A series of candidate flowering genes were detected within linkage disequilibrium (LD) blocks surrounding 12 bioclimatic TASs. Nine of these TASs exhibit significant differences in flowering time between alleles within one or more of the three individual sub-populations. Signals of selection during domestication and/or subsequent landrace diversification and adaptation were detected at 38 of the 44 flowering and bioclimatic TASs. Hence, this study lays the groundwork to begin breeding for novel environments predicted to arise following global climate change.

Enhancing Hybrid Prediction in Pearl Millet Using Genomic and/or Multi-Environment Phenotypic Information of Inbreds

Genomic selection (GS) is an emerging methodology that helps select superior lines among experimental cultivars in plant breeding programs. It offers the opportunity to increase the productivity of cultivars by delivering increased genetic gains and reducing the breeding cycles. This methodology requires inexpensive and sufficiently dense marker information to be successful, and with whole genome sequencing, it has become an important tool in many crops. The recent assembly of the pearl millet genome has made it possible to employ GS models to improve the selection procedure in pearl millet breeding programs. Here, three GS models were implemented and compared using grain yield and dense molecular marker information of pearl millet obtained from two different genotyping platforms (C [conventional GBS RAD-seq] and T [tunable GBS tGBS]). The models were evaluated using three different cross-validation (CV) schemes mimicking real situations that breeders face in breeding programs: CV2 resembles an incomplete field trial, CV1 predicts the performance of untested hybrids, and CV0 predicts the performance of hybrids in unobserved environments. We found that (i) adding phenotypic information of parental inbreds to the calibration sets improved predictive ability, (ii) accounting for genotype-by-environment interaction also increased the performance of the models, and (iii) superior strategies should consider the use of the molecular markers derived from the T platform (tGBS).

A Genome-Wide Association Study Identifies Quantitative Trait Loci Affecting Hematological Traits in Camelus bactrianus

Simple Summary

Bactrian camels can adapt to harsh natural environments. This unique tolerance of camels is tightly linked to their hematological traits, which are related to their immune, metabolic, and disease status. Therefore, mapping genomic regions that affect blood cell traits can help identify genomic characteristics that can be used as biomarkers of immune, metabolic, and disease states. This knowledge will further our understanding of the camel’s tolerance mechanisms.

Abstract

Bactrian camels (Camelus bactrianus) are one of the few large livestock species that can survive in the Gobi Desert. Animal immunity and disease resistance are related to hematological traits, which are also associated with tolerance observed in Bactrian camels. However, no genome-wide association studies have examined the genetic mechanism of the immune capability of Bactrian camels. In the present study, we used genotyping-by-sequencing data generated from 366 Bactrian camel accessions to perform a genome-wide association study for 17 hematological traits. Of the 256,616 single-nucleotide polymorphisms (SNPs) obtained, 1,635 trait–SNP associations were among the top quantitative trait locus candidates. Lastly, 664 candidate genes associated with 13 blood traits were identified. The most significant were ZNF772, MTX2, ESRRG, MEI4, IL11, FRMPD4, GABPA, NTF4, CRYBG3, ENPP5, COL16A1, and CD207. The results of our genome-wide association study provide a list of significant SNPs and candidate genes, which offer valuable information for further dissection of the molecular mechanisms that regulate the camel’s hematological traits to ultimately reveal their tolerance mechanisms.

Amand P, Le HQ, Su Z, Bai G. Multiplex restriction amplicon sequencing: a novel next-generation sequencing-based marker platform for high-throughput genotyping

To enable rapid selection of traits in marker-assisted breeding, markers must be technically simple, low-cost, high-throughput and randomly distributed in a genome. We developed such a technology, designated as Multiplex Restriction Amplicon Sequencing (MRASeq), which reduces genome complexity by polymerase chain reaction (PCR) amplification of amplicons flanked by restriction sites. The first PCR primers contain restriction site sequences at 3'-ends, preceded by 6-10 bases of specific or degenerate nucleotide sequences and then by a unique M13-tail sequence which serves as a binding site for a second PCR that adds sequencing primers and barcodes to allow sample multiplexing for sequencing. The sequences of restriction sites and adjacent nucleotides can be altered to suit different species. Physical mapping of MRASeq SNPs from a biparental population of allohexaploid wheat (Triticum aestivum L.) showed a random distribution of SNPs across the genome. MRASeq generated thousands of SNPs from a wheat biparental population and natural populations of wheat and barley (Hordeum vulgare L.). This novel, next-generation sequencing-based genotyping platform can be used for linkage mapping to screen quantitative trait loci (QTL), background selection in breeding and many other genetics and breeding applications of various species.

Genome Wide Association Study and Genomic Selection of Amino Acid Concentrations in Soybean Seeds

Soybean is a major source of protein for human consumption and animal feed. Releasing new cultivars with high nutritional value is one of the major goals in soybean breeding. To achieve this goal, genome-wide association studies of seed amino acid contents were conducted based on 249 soybean accessions from China, US, Japan, and South Korea. The accessions were evaluated for 15 amino acids and genotyped by sequencing. Significant genetic variation was observed for amino acids among the accessions. Among the 231 single nucleotide polymorphisms (SNPs) significantly associated with variations in amino acid contents, fifteen SNPs localized near 14 candidate genes involving in amino acid metabolism. The amino acids were classified into two groups with five in one group and seven amino acids in the other. Correlation coefficients among the amino acids within each group were high and positive, but the correlation coefficients of amino acids between the two groups were negative. Twenty-five SNP markers associated with multiple amino acids can be used to simultaneously improve multi-amino acid concentration in soybean. Genomic selection analysis of amino acid concentration showed that selection efficiency of amino acids based on the markers significantly associated with all 15 amino acids was higher than that based on random markers or markers only associated with individual amino acid. The identified markers could facilitate selection of soybean varieties with improved seed quality.

Single primer enrichment technology as a tool for massive genotyping: a benchmark on black poplar and maize

BACKGROUND AND AIMS

The advent of molecular breeding is advocated to improve the productivity and sustainability of second-generation bioenergy crops. Advanced molecular breeding in bioenergy crops relies on the ability to massively sample the genetic diversity. Genotyping-by-sequencing has become a widely adopted method for cost-effective genotyping. It basically requires no initial investment for design as compared with array-based platforms which have been shown to offer very robust assays. The latter, however, has the drawback of being limited to analyse only the genetic diversity accounted during selection of a set of polymorphisms and design of the assay. In contrast, genotyping-by-sequencing with random sampling of genomic loci via restriction enzymes or random priming has been shown to be fast and convenient but lacks the ability to target specific regions of the genome and to maintain high reproducibility across laboratories.

METHODS

Here we present a first adoption of single-primer enrichment technology (SPET) which provides a highly efficient and scalable system to obtain targeted sequence-based large genotyping data sets, bridging the gaps between array-based systems and traditional sequencing-based protocols. To fully explore SPET performance, we conducted a benchmark study in ten Zea mays lines and a large-scale study of a natural black poplar population of 540 individuals with the aim of discovering polymorphisms associated with biomass-related traits.

KEY RESULTS

Our results showed the ability of this technology to provide dense genotype information on a customized panel of selected polymorphisms, while yielding hundreds of thousands of untargeted variable sites. This provided an ideal resource for association analysis of natural populations harbouring unexplored allelic diversities and structure such as in black poplar.

CONCLUSION

The improvement of sequencing throughput and the development of efficient library preparation methods has made it feasible to carry out targeted genotyping-by-sequencing experiments cost-competitively with either random complexity reduction systems or traditional array-based platforms, while maintaining the key advantages of both technologies.

Single primer enrichment technology as a tool for massive genotyping: a benchmark on black poplar and maize

BACKGROUND AND AIMS

The advent of molecular breeding is advocated to improve the productivity and sustainability of second-generation bioenergy crops. Advanced molecular breeding in bioenergy crops relies on the ability to massively sample the genetic diversity. Genotyping-by-sequencing has become a widely adopted method for cost-effective genotyping. It basically requires no initial investment for design as compared with array-based platforms which have been shown to offer very robust assays. The latter, however, has the drawback of being limited to analyse only the genetic diversity accounted during selection of a set of polymorphisms and design of the assay. In contrast, genotyping-by-sequencing with random sampling of genomic loci via restriction enzymes or random priming has been shown to be fast and convenient but lacks the ability to target specific regions of the genome and to maintain high reproducibility across laboratories.

METHODS

Here we present a first adoption of single-primer enrichment technology (SPET) which provides a highly efficient and scalable system to obtain targeted sequence-based large genotyping data sets, bridging the gaps between array-based systems and traditional sequencing-based protocols. To fully explore SPET performance, we conducted a benchmark study in ten Zea mays lines and a large-scale study of a natural black poplar population of 540 individuals with the aim of discovering polymorphisms associated with biomass-related traits.

KEY RESULTS

Our results showed the ability of this technology to provide dense genotype information on a customized panel of selected polymorphisms, while yielding hundreds of thousands of untargeted variable sites. This provided an ideal resource for association analysis of natural populations harbouring unexplored allelic diversities and structure such as in black poplar.

CONCLUSION

The improvement of sequencing throughput and the development of efficient library preparation methods has made it feasible to carry out targeted genotyping-by-sequencing experiments cost-competitively with either random complexity reduction systems or traditional array-based platforms, while maintaining the key advantages of both technologies.

A Novel QTL Conferring Fusarium Crown Rot Resistance Located on Chromosome Arm 6HL in Barley

Fusarium crown rot (FCR), caused primarily by Fusarium pseudograminearum, is a devastating disease for cereal production in semi-arid regions worldwide. To identify and characterize loci conferring FCR resistance, we assessed a landrace AWCS799 which is among the best lines identified from a systematic screening of more than 1,000 genotypes. Genetic control of its resistance was investigated by generating and analyzing two populations of recombinant inbred lines with AWCS799 as the common parent. One of the populations was used for QTL detection and the other for validation. A novel QTL, located on the long arm of chromosome 6H (designated as Qcrs.caf-6H), was consistently detected in each of the four FCR severity tests conducted against the mapping population. The QTL explained up to 28.3% of the phenotypic variance, and its effect was confirmed in the validation population. Significant interaction between this resistance locus and either plant height or heading date was not detected, further facilitating its manipulation in breeding programs.