Construction of a SSR-based genetic map and identification of QTL for domestication traits using recombinant inbred lines from a cross between wild and cultivated cowpea (<Emphasis Type="Italic">V. unguiculata</Emphasis> (L.) Walp.)

Cowpea (Vigna unguiculata (L.) Walp.) is a grain legume commonly grown and consumed in many parts of the tropics and subtropics. A genetic linkage map was constructed using simple sequence repeat (SSR) markers and a recombinant inbred (RI) population of159 individuals derived from a cross between the breeding line 524B, a California Blackeye, and 219-01, a perennial wild cowpea from Kenya. Out of 912 primer combinations predicted to amplify SSRs in cowpea, 639 reliably produced amplification products in PCR assays and 202 (31.6%) were polymorphic between the two parents. These polymorphic SSRs were used to construct a genetic map consisting of 11 linkage groups (LGs) spanning 677 cM, with an average distance between markers of 3 cM. Agronomic traits related to domestication (seed weight, pod shattering) were analyzed together with the genotypic data. Six quantitative trait loci (QTL) for seed size were revealed with the phenotypic variation ranging from 8.9 to 19.1%. Four QTL for pod shattering were identified with the phenotypic variation ranging from 6.4 to 17.2%. The QTL for seed size and pod shattering mainly cluster in two areas of LGs 1 and 10, facilitating the use of marker-assisted selection to eliminate undesirable wild phenotypes in breeding activities involving introgression of traits from wild germplasm. The generation of an SSR-based molecular map and additional trait-linked markers also contributes to the expanding tool kit available to cowpea breeders, especially in Africa.     

Pod shattering in grain legumes: emerging genetic and environment-related patterns

A reduction in pod shattering is one of the main components of grain legume domestication. Despite this, many domesticated legumes suffer serious yield losses due to shattering, particularly under arid conditions. Mutations related to pod shattering modify the twisting force of pod walls or the structural strength of the dehiscence zone in pod sutures. At a molecular level, a growing body of evidence indicates that these changes are controlled by a relatively small number of key genes that have been selected in parallel across grain legume species, supporting partial molecular convergence. Legume homologs of Arabidopsis thaliana silique shattering genes play only minor roles in legume pod shattering. Most domesticated grain legume species contain multiple shattering-resistance genes, with mutants of each gene typically showing only partial shattering resistance. Hence, crosses between varieties with different genes lead to transgressive segregation of shattering alleles, producing plants with either enhanced shattering resistance or atavistic susceptibility to the trait. The frequency of these resistance pod-shattering alleles is often positively correlated with environmental aridity. The continued development of pod-shattering-related functional information will be vital for breeding crops that are suited to the increasingly arid conditions expected in the coming decades.

Recent genetic, genomic, and phenotypic studies of pod shattering in grain legumes lay the foundation for breeding crops suited for increasingly arid conditions.     

Sequence polymorphisms in wild,weedy, and cultivated rice suggest seed‐shattering locus sh4 played a minor role in Asian rice domestication

The predominant view regarding Asian rice domestication is that the initial origin of nonshattering involved a single gene of large effect, specifically, the sh4 locus via the evolutionary replacement of a dominant allele for shattering with a recessive allele for reduced shattering. Data have accumulated to challenge this hypothesis. Specifically, a few studies have reported occasional seed‐shattering plants from populations of the wild progenitor of cultivated rice (Oryza rufipogon complex) being homozygous for the putative “nonshattering” sh4 alleles. We tested the sh4 hypothesis for the domestication of cultivated rice by obtaining genotypes and phenotypes for a diverse set of samples of wild, weedy, and cultivated rice accessions. The cultivars were fixed for the putative “nonshattering” allele and nonshattering phenotype, but wild rice accessions are highly polymorphic for the putative “nonshattering” allele (frequency ~26%) with shattering phenotype. All weedy rice accessions are the “nonshattering” genotype at the sh4 locus but with shattering phenotype. These data challenge the widely accepted hypothesis that a single nucleotide mutation (“G”/“T”) of the sh4 locus is the major driving force for rice domestication. Instead, we hypothesize that unidentified shattering loci are responsible for the initial domestication of cultivated rice through reduced seed shattering.     

More than one way to evolve a weed: parallel evolution of US weedy rice through independent genetic mechanisms

Many different crop species were selected for a common suite of ‘domestication traits’, which facilitates their use for studies of parallel evolution. Within domesticated rice (Oryza sativa), there has also been independent evolution of weedy strains from different cultivated varieties. This makes it possible to examine the genetic basis of parallel weed evolution and the extent to which this process occurs through shared genetic mechanisms. We performed comparative QTL mapping of weediness traits using two recombinant inbred line populations derived from crosses between an indica crop variety and representatives of each of the two independently evolved weed strains found in US rice fields, strawhull (S) and blackhull awned (B). Genotyping‐by‐sequencing provided dense marker coverage for linkage map construction (average marker interval <0.25 cM), with 6016 and 13 730 SNPs mapped in F₅ lines of the S and B populations, respectively. For some weediness traits (awn length, hull pigmentation and pericarp pigmentation), QTL mapping and sequencing of underlying candidate genes confirmed that trait variation was largely attributable to individual loci. However, for more complex quantitative traits (including heading date, panicle length and seed shattering), we found multiple QTL, with little evidence of shared genetic bases between the S and B populations or across previous studies of weedy rice. Candidate gene sequencing revealed causal genetic bases for 8 of 27 total mapped QTL. Together these findings suggest that despite the genetic bottleneck that occurred during rice domestication, there is ample genetic variation in this crop to allow agricultural weed evolution through multiple genetic mechanisms.     

Whole-genome assembly and resequencing reveal genomic imprint and key genes of rapid domestication in narrow-leafed lupin

Shifting from a livestock-based protein diet to a plant-based protein diet has been proposed as an essential requirement to maintain global food sustainability, which requires the increased production of protein-rich crops for direct human consumption. Meanwhile, the lack of sufficient genetic diversity in crop varieties is an increasing concern for sustainable food supplies. Countering this concern requires a clear understanding of the domestication process and dynamics. Narrow-leafed lupin (Lupinus angustifolius L.) has experienced rapid domestication and has become a new legume crop over the past century, with the potential to provide protein-rich seeds. Here, using long-read whole-genome sequencing, we assembled the third-generation reference genome for the narrow-leafed lupin cultivar Tanjil, comprising 20 chromosomes with a total genome size of 615.8 Mb and contig N50 = 5.65 Mb. We characterized the original mutation and putative biological pathway resulting in low seed alkaloid level that initiated the recent domestication of narrow-leafed lupin. We identified a 1133-bp insertion in the cis-regulatory region of a putative gene that may be associated with reduced pod shattering (lentus). A comparative analysis of genomic diversity in cultivars and wild types identified an apparent domestication bottleneck, as precisely predicted by the original model of the bottleneck effect on genetic variability in populations. Our results identify the key domestication genetic loci and provide direct genomic evidence for a domestication bottleneck, and open up the possibility of knowledge-driven de novo domestication of wild plants as an avenue to broaden crop plant diversity to enhance food security and sustainable low-carbon emission agriculture.     

Construction and comparison of three reference‐quality genome assemblies for soybean

We report reference‐quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single‐nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan‐gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40–42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.     

Comprehensive genomic analyses of Vigna unguiculata provide insights into population differentiation and the genetic basis of key agricultural traits

Vigna unguiculata is an important legume crop worldwide. The subsp. sesquipedalis and unguiculata are the two major types grown; the former is mainly grown in Asia to produce fresh pods, while the latter is mainly grown in Africa to produce seeds. Here, a chromosome-scale genome for subsp. sesquipedalis was generated by combining high-fidelity (HiFi) long-read sequencing with high-throughput chromosome conformation capture (Hi-C) technology. The genome size for all contigs and N50 were 594 and 18.5 Mb, respectively. The Hi-C interaction map helped cluster 91% of the contigs into 11 chromosomes. Genome comparisons between subsp. sesquipedalis and unguiculata revealed extensive genomic variations, and some variations resulted in gene loss. A germplasm panel with 315 accessions of V. unguiculata was resequenced, and a genomic variation map was constructed. Population structure and phylogenetic analyses suggested that subsp. sesquipedalis originated from subsp. unguiculata. Highly differentiated genomic regions were also identified, and a number of genes functionally enriched in adaptations were located in these regions. Two traits, pod length (PL) and pod width (PW), were observed for this germplasm, and genome-wide association analysis of these traits was performed. The quantitative trait loci (QTLs) for these two traits were identified, and their candidate genes were uncovered. Interestingly, genomic regions of PL QTLs also showed strong signals of artificial selection. Taken together, the results of this study provide novel insights into the population differentiation and genetic basis of key agricultural traits in V. unguiculata.     

Fine mapping of a major quantitative trait locus that regulates pod shattering in soybean

Legumes represent the second most important family of crop plants, accounting for ~27 % of the world’s crop production. While some legumes are grown as forages or vegetables, most crop legumes are grown for harvesting their nutritious seeds. The legume seeds are contained in the pod, which is composed of a single seed-bearing carpel that, when matures, splits open along two seams, a process called pod dehiscence or pod shattering. Pod shattering before or during harvest causes yield losses of grain legumes. Moreover, the dominant shattering trait of the wild progenitors is a limiting factor for efficient introgression of value-added traits into elite breeding lines. Knowledge of the genetic mechanisms underlying pod shattering will facilitate breeding of shattering-resistant varieties, expedite introgression of agronomically favorable traits from wild species to elite breeding lines, and enrich our understanding of the evolution of seed dispersal and crop domestication in diverse crop species. Here we report fine mapping of a major quantitative trait locus (designated as qPDH1) that regulates pod shattering in soybean (Glycine max). A combination of linkage and association mapping allowed us to delimit the qPDH1 locus within a 47-kb region on chromosome 16. The data reported here will facilitate positional cloning of the underlying gene and the development of breeder-friendly genetic markers for marker-assisted selection in soybean.     

Multiple genetic pathways for seed shattering in the grasses

Shattering is an essential seed dispersal mechanism in wild species. It is believed that independent mutations at orthologous loci led to convergent domestication of cereal crops. To investigate genetic relationships of Triticeae shattering genes with those of other grasses, we mapped spike-, barrel- (B-type), and wedge-type (W-type) spikelet disarticulation genes in wheat and its wild relatives. The Br1 gene for W-type disarticulation was mapped to a region delimited by Xpsr598 and Xpsr1196 on the short arm of chromosomes 3A in Triticum timopheevii and 3S in Aegilops speltoides. The spike- and W-type disarticulation genes are allelic at Br1 in Ae. speltoides. The B-type disarticulation gene, designated as Br2, was mapped to an interval of 4.4 cM between Xmwg2013 and Xpsr170 on the long arm of chromosome 3D in Aegilops tauschii, the D-genome donor of common wheat. Therefore, B- and W-type disarticulations are governed by two different orthologous loci on group-3 chromosomes. Based on map position, orthologs of Br1 and Br2 were not detected in barley, maize, rice, and sorghum, indicating multiple genetic pathways for shattering in grasses. The implications of the mapping results are discussed with regard to the evolution of polyploid wheat and domestication of cereals.Supplementary material is available in the online version of this article at     

Largely unlinked gene sets targeted by selection for domestication syndrome phenotypes in maize and sorghum

The domestication of diverse grain crops from wild grasses was a result of artificial selection for a suite of overlapping traits producing changes referred to in aggregate as ‘domestication syndrome’. Parallel phenotypic change can be accomplished by either selection on orthologous genes or selection on non‐orthologous genes with parallel phenotypic effects. To determine how often artificial selection for domestication traits in the grasses targeted orthologous genes, we employed resequencing data from wild and domesticated accessions of Zea (maize) and Sorghum (sorghum). Many ‘classic’ domestication genes identified through quantitative trait locus mapping in populations resulting from wild/domesticated crosses indeed show signatures of parallel selection in both maize and sorghum. However, the overall number of genes showing signatures of parallel selection in both species is not significantly different from that expected by chance. This suggests that while a small number of genes will extremely large phenotypic effects have been targeted repeatedly by artificial selection during domestication, the optimization part of domestication targeted small and largely non‐overlapping subsets of all possible genes which could produce equivalent phenotypic alterations.     

Transcriptional profile of Pseudomonas syringae pv. phaseolicola NPS3121 in response to tissue extracts from a susceptible Phaseolus vulgaris L. cultivar

Background  

Pseudomonas syringae pv. phaseolicola is a Gram-negative plant-pathogenic bacterium that causes "halo blight" disease of beans (Phaseolus vulgaris L.). This disease affects both foliage and pods, and is a major problem in temperate areas of the world. Although several bacterial genes have been determined as participants in pathogenesis, the overall process still remains poorly understood, mainly because the identity and function of many of the genes are largely unknown. In this work, a genomic library of P. syringae pv. phaseolicola NPS3121 was constructed and PCR amplification of individual fragments was carried out in order to print a DNA microarray. This microarray was used to identify genes that are differentially expressed when bean leaf extracts, pod extracts or apoplastic fluid were added to the growth medium.     

Development of sequence-specific PCR markers linked to the Tardus gene that reduces pod shattering in narrow-leafed lupin (Lupinus angustifolius L.)

Seed pods of wild-type narrow-leafed lupins (Lupinus angustifolius L.) shatter upon maturity, dispersing their seeds. Recessive alleles of the genes Tardus and Lentus that confer reduced pod shattering have been incorporated into domesticated cultivars to facilitate harvesting. Tardus was mapped in an F₈ recombinant inbred population of a cross between domesticated and wild lupins. A microsatellite–anchored fragment length polymorphism marker (TaM1), which mapped 2.1 cM from Tardus, was converted to a locus-specific PCR assay. Marker TaM2, a restriction fragment length polymorphism marker was converted to a PCR assay and mapped to 3.9 cM on the other side of Tardus. Marker TaM3, a cleaved amplified polymorphic sequence marker, was positioned along-side marker TaM1 at 3.9 cM from Tardus. One or more markers was polymorphic in 70% of possible pairwise crosses between Australian domesticated lines and wild accessions tested, indicating wide applicability of the markers in crosses between wild and domesticated germplasm.     

The genome of cowpea (Vigna unguiculata [L.] Walp.)

Cowpea (Vigna unguiculata [L.] Walp.) is a major crop for worldwide food and nutritional security, especially in sub‐Saharan Africa, that is resilient to hot and drought‐prone environments. An assembly of the single‐haplotype inbred genome of cowpea IT97K‐499‐35 was developed by exploiting the synergies between single‐molecule real‐time sequencing, optical and genetic mapping, and an assembly reconciliation algorithm. A total of 519 Mb is included in the assembled sequences. Nearly half of the assembled sequence is composed of repetitive elements, which are enriched within recombination‐poor pericentromeric regions. A comparative analysis of these elements suggests that genome size differences between Vigna species are mainly attributable to changes in the amount of Gypsy retrotransposons. Conversely, genes are more abundant in more distal, high‐recombination regions of the chromosomes; there appears to be more duplication of genes within the NBS‐LRR and the SAUR‐like auxin superfamilies compared with other warm‐season legumes that have been sequenced. A surprising outcome is the identification of an inversion of 4.2 Mb among landraces and cultivars, which includes a gene that has been associated in other plants with interactions with the parasitic weed Striga gesnerioides. The genome sequence facilitated the identification of a putative syntelog for multiple organ gigantism in legumes. A revised numbering system has been adopted for cowpea chromosomes based on synteny with common bean (Phaseolus vulgaris). An estimate of nuclear genome size of 640.6 Mbp based on cytometry is presented.     

Molecular mechanism underlying pseudopeloria in Habenaria radiata (Orchidaceae)

Habenaria radiata (Orchidaceae) has two whorls of perianth, comprising three greenish sepals, two white petals and one lip (labellum). By contrast, the pseudopeloric (with a decreased degree of zygomorphy) mutant cultivar of H. radiata , ‘Hishou’, has changes in the identities of the dorsal sepal to a petaloid organ and the two ventral sepals to lip‐like organs. Here, we isolated four DEFICIENS ‐ like and two AGL 6 ‐like genes from H. radiata , and characterized their expression. Most of these genes revealed similar expression patterns in the wild type and in the ‘Hishou’ cultivar, except Hr DEF ‐C3. The Hr DEF ‐C3 gene was expressed in petals and lip in the wild type but was ectopically expressed in sepal, petals, lip, leaf, root and bulb in ‘Hishou’. Sequence analysis of the Hr DEF ‐C3 loci revealed that the ‘Hishou’ genome harbored two types of Hr DEF ‐C3 genes: one identical to wild‐type Hr DEF ‐C3 and the other carrying a retrotransposon insertion in its promoter. Genetic linkage analysis of the progeny derived from an intraspecific cross between ‘Hishou’ and the wild type demonstrated that the mutant pseudopeloric trait was dominantly inherited and was linked to the Hr DEF ‐C3 gene carrying the retrotransposon. These results indicate that the pseudopeloric phenotype is caused by retrotransposon insertion in the Hr DEF ‐C3 promoter, resulting in the ectopic expression of Hr DEF ‐C3 . As the expression of Hr AGL 6‐C2 was limited to lateral sepals and lip, the overlapping expression of Hr DEF ‐C3 and Hr AGL 6‐C2 is likely to be responsible for the sepal to lip‐like identity in the lateral sepals of the ‘Hishou’ cultivar.     

Signatures of de‐domestication in autochthonous pig breeds and of domestication in wild boar populations from MC1R and NR6A1 allele distribution

Autochthonous pig breeds are usually reared in extensive or semi‐extensive production systems that might facilitate contact with wild boars and, thus, reciprocal genetic exchanges. In this study, we analysed variants in the melanocortin 1 receptor (MC1R) gene (which cause different coat colour phenotypes) and in the nuclear receptor subfamily 6 group A member 1 (NR6A1) gene (associated with increased vertebral number) in 712 pigs of 12 local pig breeds raised in Italy (Apulo‐Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano and Sarda) and south‐eastern European countries (Kr?kopolje from Slovenia, Black Slavonian and Turopolje from Croatia, Mangalitsa and Moravka from Serbia and East Balkan Swine from Bulgaria) and compared the data with the genetic variability at these loci investigated in 229 wild boars from populations spread in the same macro‐geographic areas. None of the autochthonous pig breeds or wild boar populations were fixed for one allele at both loci. Domestic and wild‐type alleles at these two genes were present in both domestic and wild populations. Findings of the distribution of MC1R alleles might be useful for tracing back the complex genetic history of autochthonous breeds. Altogether, these results indirectly demonstrate that bidirectional introgression of wild and domestic alleles is derived and affected by the human and naturally driven evolutionary forces that are shaping the Sus scrofa genome: autochthonous breeds are experiencing a sort of ‘de‐domestication’ process, and wild resources are challenged by a ‘domestication’ drift. Both need to be further investigated and managed.