A graph-based pan-genome of Brassica oleracea provides new insights into its domestication and morphotype diversification.

The domestication of Brassica oleracea has resulted in diverse morphological types with distinct patterns of organ development. Here we report a graph-based pan-genome of B. oleracea constructed from high-quality genome assemblies of different morphotypes. The pan-genome harbors over 200 structural variant hotspot regions enriched in auxin- and flowering-related genes. Population genomic analyses revealed that early domestication of B. oleracea focused on leaf or stem development. Gene flows resulting from agricultural practices and variety improvement were detected among different morphotypes. Selective-sweep and pan-genome analyses identified an auxin-responsive small auxin up-regulated RNA gene and a CLAVATA3/ESR-RELATED family gene as crucial players in leaf-stem differentiation during the early stage of B. oleracea domestication and the BoKAN1 gene as instrumental in shaping the leafy heads of cabbage and Brussels sprouts. Our pan-genome and functional analyses further revealed that variations in the BoFLC2 gene play key roles in the divergence of vernalization and flowering characteristics among different morphotypes, and variations in the first intron of BoFLC3 are involved in fine-tuning the flowering process in cauliflower. This study provides a comprehensive understanding of the pan-genome of B. oleracea and sheds light on the domestication and differential organ development of this globally important crop species.

TaTPP-7A positively feedback regulates grain filling and wheat grain yield through T6P-SnRK1 signalling pathway and sugar-ABA interaction.

Grain size and filling are two key determinants of grain thousand-kernel weight (TKW) and crop yield, therefore they have undergone strong selection since cereal was domesticated. Genetic dissection of the two traits will improve yield potential in crops. A quantitative trait locus significantly associated with wheat grain TKW was detected on chromosome 7AS flanked by a simple sequence repeat marker of Wmc17 in Chinese wheat 262 mini-core collection by genome-wide association study. Combined with the bulked segregant RNA-sequencing (BSR-seq) analysis of an F2 genetic segregation population with extremely different TKW traits, a candidate trehalose-6-phosphate phosphatase gene located at 135.0 Mb (CS V1.0), designated as TaTPP-7A, was identified. This gene was specifically expressed in developing grains and strongly influenced grain filling and size. Overexpression (OE) of TaTPP-7A in wheat enhanced grain TKW and wheat yield greatly. Detailed analysis revealed that OE of TaTPP-7A significantly increased the expression levels of starch synthesis- and senescence-related genes involved in abscisic acid (ABA) and ethylene pathways. Moreover, most of the sucrose metabolism and starch regulation-related genes were potentially regulated by SnRK1. In addition, TaTPP-7A is a crucial domestication- and breeding-targeted gene and it feedback regulates sucrose lysis, flux, and utilization in the grain endosperm mainly through the T6P-SnRK1 pathway and sugar-ABA interaction. Thus, we confirmed the T6P signalling pathway as the central regulatory system for sucrose allocation and source-sink interactions in wheat grains and propose that the trehalose pathway components have great potential to increase yields in cereal crops.

Resistance of QYm.nau-2D to wheat yellow mosaic virus was derived from an alien introgression into common wheat.

KEY MESSAGE: The QYm.nau-2D locus conferring wheat yellow mosaic virus resistance is an exotic introgression and we developed 11 diagnostic markers tightly linked to QYm.nau-2D. Wheat yellow mosaic virus (WYMV) is a serious disease of winter wheat in China. Breeding resistant varieties is the most effective strategy for WYMV control. A WYMV resistant locus QYm.nau-2D on the chromosome arm 2DL has been repeatedly reported but the mapped region is large. In the present study, we screened recombinants using a biparental population and mapped QYm.nau-2D into an 18.8 Mb physical interval. By genome-wide association studies of 372 wheat varieties for WYMV resistance in four environments, we narrowed down QYm.nau-2D into a 16.4 Mb interval. Haplotype analysis indicated QYm.nau-2D were present as six different states due to recombination during hybridization breeding. QYm.nau-2D was finally mapped into a linkage block of 11.2 Mb. Chromosome painting using 2D specific probes and collinearity analysis among the published sequences corresponding to QYm.nau-2D region indicated the block was an exotic introgression. The Illumina-sequenced reads of four diploid Aegilops species were mapped to the sequence of Fielder, a variety having the introgression. The mapping reads were significantly increased at the putative introgression regions of Fielder. Ae. uniaristata (NN) had the highest mapping reads, suggesting that QYm.nau-2D was possibly an introgression from genome N. We investigated the agronomic performances of different haplotypes and observed no linkage drag of the alien introgression for the 15 tested traits. For marker-assisted selection of QYm.nau-2D, we developed 11 diagnostic markers tightly linked to the locus. This research provided a case study of an exotic introgression, which has been utilized in wheat improvement for WYMV resistance.

A chromosome-scale genome assembly of Dasypyrum villosum provides insights into its application as a broad-spectrum disease resistance resource for wheat improvement.

Dasypyrum villosum is one of the most valuable gene resources in wheat improvement, especially for disease resistance. The mining of favorable genes from D. villosum is frustrated by the lack of a whole genome sequence. In this study, we generated a doubled-haploid line, 91C43DH, using microspore culture and obtained a 4.05-GB high-quality, chromosome-scale genome assembly for D. villosum. The assembly contains39 727 high-confidence genes, and 85.31% of the sequences are repetitive. Two reciprocal translocation events were detected, and 7VS-4VL is a unique translocation in D. villosum. The prolamin seed storage protein-coding genes were found to be duplicated; in particular, the genes encoding low-molecular-weight glutenin at the Glu-V3 locus were significantly expanded. RNA sequencing (RNA-seq) analysis indicated that, after Blumeria graminearum f.sp tritici (Bgt) inoculation, there were more upregulated genes involved in the pattern-triggered immunity and effector-triggered immunity defense pathways in D. villosum than in Triticum urartu. MNase hypersensitive sequencing (MH-seq) identified two Bgt-inducible MH sites (MHSs), one in the promoter and one in the 3' terminal region of the powdery mildew resistance (Pm) gene NLR1-V. Each site had two subpeaks and they were termed MHS1 (MHS1.1/1.2) and MHS2 (MHS2.1/2.2). Bgt-inducible MHS2.2 was uniquely present in D. villosum, and MHS1.1 was more inducible in D. villosum than in wheat, suggesting that MHSs may be critical for regulation of NLR1-V expression and plant defense. In summary, this study provides a valuable genome resource for functional genomics studies and wheat-D. villosum introgression breeding. The identified regulatory mechanisms may also be exploited to develop new strategies for enhancing Pm resistance by optimizing gene expression in wheat.

Fast integration and accumulation of beneficial breeding alleles through an AB-NAMIC strategy in wheat.

Wheat (Triticum aestivum) is among the most important staple crops for safeguarding the food security of the growing world population. To bridge the gap between genebank diversity and breeding programs, we developed an advanced backcross-nested association mapping plus inter-crossed population (AB-NAMIC) by crossing three popular wheat cultivars as recurrent founders to 20 germplasm lines from a mini core collection. Selective backcrossing combined with selection against undesirable traits and extensive crossing within and between sub-populations created new opportunities to detect unknown genes and increase the frequency of beneficial alleles in the AB-NAMIC population. We performed phenotyping of 590 AB-NAMIC lines and a natural panel of 476 cultivars for six consecutive growing seasons and genotyped these 1066 lines with a 660K SNP array. Genome-wide association studies of both panels for plant development and yield traits demonstrated improved power to detect rare alleles and loci with medium genetic effects in AB-NAMIC. Notably, genome-wide association studies in AB-NAMIC detected the candidate gene TaSWEET6-7B (TraesCS7B03G1216700), which has high homology to the rice SWEET6b gene and exerts strong effects on adaptation and yield traits. The commercial release of two derived AB-NAMIC lines attests to its direct applicability in wheat improvement. Valuable information on genome-wide association study mapping, candidate genes, and their haplotypes for breeding traits are available through WheatGAB. Our research provides an excellent framework for fast-tracking exploration and accumulation of beneficial alleles stored in genebanks.

Genomic footprints of Kam Sweet Rice domestication indicate possible migration routes of the Dong people in China and provide resources for future rice breeding.

The Dong people are one of China's 55 recognized ethnic minorities, but there has been a long-standing debate about their origins. In this study, we performed whole-genome resequencing of Kam Sweet Rice (KSR), a valuable, rare, and ancient rice landrace unique to the Dong people. Through comparative genomic analyses of KSR and other rice landraces from south of the Yangtze River Basin in China, we provide evidence that the ancestors of the Dong people likely originated from the southeast coast of China at least 1000 years ago. Alien introgression and admixture in KSR demonstrated multiple migration events in the history of the Dong people. Genomic footprints of domestication demonstrated characteristics of KSR that arose from artificial selection and geographical adaptation by the Dong people. The key genes GS3, Hd1, and DPS1 (related to agronomic traits) and LTG1 and MYBS3 (related to cold tolerance) were identified as domestication targets, reflecting crop improvement and changes in the geographical environment of the Dong people during migration. A genome-wide association study revealed a candidate yield-associated gene, Os01g0923300, a specific haplotype in KSR that is important for regulating grain number per panicle. RNA-sequencing and quantitative reverse transcription-PCR results showed that this gene was more highly expressed in KSR than in ancestral populations, indicating that it may have great value in increasing yield potential in other rice accessions. In summary, our work develops a novel approach for studying human civilization and migration patterns and provides valuable genomic datasets and resources for future breeding of high-yield and climate-resilient rice varieties.

Genome-wide signatures of the geographic expansion and breeding of soybean.

Soybean is a leguminous crop that provides oil and protein. Exploring the genomic signatures of soybean evolution is crucial for breeding varieties with improved adaptability to environmental extremes. We analyzed the genome sequences of 2,214 soybeans and proposed a soybean evolutionary route, i.e., the expansion of annual wild soybean (Glycine soja Sieb. & Zucc.) from southern China and its domestication in central China, followed by the expansion and local breeding selection of its landraces (G. max (L.) Merr.). We observed that the genetic introgression in soybean landraces was mostly derived from sympatric rather than allopatric wild populations during the geographic expansion. Soybean expansion and breeding were accompanied by the positive selection of flowering time genes, including GmSPA3c. Our study sheds light on the evolutionary history of soybean and provides valuable genetic resources for its future breeding.

Pleiotropic Effect of the compactum Gene and Its Combined Effects with Other Loci for Spike and Grain-Related Traits in Wheat.

Club wheat (Triticum aestivum ssp. compactum) with a distinctly compact spike morphology was conditioned by the dominant compactum (C) locus on chromosome 2D and resulted in a redistribution of spike yield components. The disclosure of the genetic basis of club wheat was a prerequisite for the development of widely adapted, agronomically competitive club wheat cultivars. In this study, we used a recombinant inbred line population derived from a cross between club wheat Hiller and modern cultivar Yangmai 158 to construct a genetic linkage map and identify quantitative trait loci associated with 15 morphological traits. The club allele acted in a semi-dominant manner and the C gene was mapped to 370.12-406.29 Mb physical region on the long arm of 2D. Apart from compact spikes, C exhibited a pleiotropic effect on ten other agronomic traits, including plant height, three spike-related traits and six grain-related traits. The compact spike phenotype was correlated with decreased grain size and weight, but with an increase in floret fertility and grain number. These pleiotropic effects make club wheat have compatible spike weight with a normal spike from common wheat. The genetic effects of various gene combinations of C with four yield-related genes, including Ppd-D1, Vrn-D3, Rht-B1b and Rht8, were evaluated. C had no epistatic interaction with any of these genes, indicating that their combinations would have an additive effect on other agronomically important traits. Our research provided a theoretical foundation for the potentially effective deployment of C gene into modern breeding varieties in combination with other favorable alleles.

Genomic insights into historical improvement of heterotic groups during modern hybrid maize breeding.

Single-cross maize hybrids display superior heterosis and are produced from crossing two parental inbred lines belonging to genetically different heterotic groups. Here we assembled 1,604 historically utilized maize inbred lines belonging to various female heterotic groups (FHGs) and male heterotic groups (MHGs), and conducted phenotyping and genomic sequencing analyses. We found that the FHGs and MHGs have undergone both convergent and divergent changes for different sets of agronomic traits. Using genome-wide selection scans and association analyses, we identified a large number of candidate genes that contributed to the improvement of agronomic traits of the FHGs and MHGs. Moreover, we observed increased genetic differentiation between the FHGs and MHGs across the breeding eras, and we found a positive correlation between increasing heterozygosity levels in the differentiated genes and heterosis in hybrids. Furthermore, we validated the function of two selected genes and a differentiated gene. This study provides insights into the genomic basis of modern hybrid maize breeding.

Wheat breeding history reveals synergistic selection of pleiotropic genomic sites for plant architecture and grain yield.

Diversity surveys of crop germplasm are important for gaining insights into the genomic basis for plant architecture and grain yield improvement, which is still poorly understood in wheat. In this study, we exome sequenced 287 wheat accessions that were collected in the past 100 years. Population genetics analysis identified that 6.7% of the wheat genome falls within the selective sweeps between landraces and cultivars, which harbors the genes known for yield improvement. These regions were asymmetrically distributed on the A and B subgenomes with regulatory genes being favorably selected. Genome-wide association study (GWAS) identified genomic loci associated with traits for yield potential, and two underlying genes, TaARF12 encoding an auxin response factor and TaDEP1 encoding the G-protein γ-subunit, were located and characterized to pleiotropically regulate both plant height and grain weight. Elite single-nucleotide haplotypes with increased allele frequency in cultivars relative to the landraces were identified and found to have accumulated over the course of breeding. Interestingly, we found that TaARF12 and TaDEP1 function in epistasis with the classical plant height Rht-1 locus, leading to propose a "Green Revolution"-based working model for historical wheat breeding. Collectively, our study identifies selection signatures that fine-tune the gibberellin pathway during modern wheat breeding and provides a wealth of genomic diversity resources for the wheat research community.

Genome-wide association study-based identification genes influencing agronomic traits in rice (Oryza sativa L.).

Rice is one of the most important cereal crops, providing the daily dietary intake for approximately 50% of the global human population. Here, we re-sequenced 259 rice accessions, generating 1371.65 Gb of raw data. Furthermore, we performed genome-wide association studies (GWAS) on 13 agronomic traits using 2.8 million single nucleotide polymorphisms (SNPs) characterized in 259 rice accessions. Phenotypic data and best linear unbiased prediction (BLUP) values of each of the 13 traits over two years of each trait were used for the GWAS. The results showed that 816 SNP signals were significantly associated with the 13 agronomic traits. Then we detected candidate genes related to target traits within 200 kb upstream and downstream of the associated SNP loci, based on linkage disequilibrium (LD) blocks in the whole rice genome. These candidate genes were further identified through haplotype block constructions. This comprehensive study provides a timely and important genomic resource for breeding high yielding rice cultivars.

Identification of rice (Oryza sativa L.) genes involved in sheath blight resistance via a genome-wide association study.

Rice sheath blight (RSB) is an economically significant disease affecting rice yield worldwide. Genetic resistance to RSB is associated with multiple minor genes, with each providing a minor phenotypic effect, but the underlying dominant resistance genes remain unknown. A genome-wide association study (GWAS) of 259 diverse rice varieties, with genotypes based on a single nucleotide polymorphism (SNP) and haplotype, was conducted to assess their sheath blight reactions at three developmental stages (seedlings, tillering and booting). A total of 653 genes were correlated with sheath blight resistance, of which the disease resistance protein RPM1 (OsRSR1) and protein kinase domain-containing protein (OsRLCK5) were validated by overexpression and knockdown assays. We further found that the coiled-coil (CC) domain of OsRSR1 (OsRSR1-CC) and full-length OsRLCK5 interacted with serine hydroxymethyltransferase 1 (OsSHM1) and glutaredoxin (OsGRX20), respectively. It was found that OsSHM1, which has a role in the reactive oxygen species (ROS) burst, and OsGRX20 enhanced the antioxidation ability of plants. A regulation model of the new RSB resistance though the glutathione (GSH)-ascorbic acid (AsA) antioxidant system was therefore revealed. These results enhance our understanding of RSB resistance mechanisms and provide better gene resources for the breeding of disease resistance in rice.

Triticum population sequencing provides insights into wheat adaptation.

Bread wheat expanded its habitat from a core area of the Fertile Crescent to global environments within ~10,000 years. The genetic mechanisms of this remarkable evolutionary success are not well understood. By whole-genome sequencing of populations from 25 subspecies within the genera Triticum and Aegilops, we identified composite introgression from wild populations contributing to a substantial portion (4-32%) of the bread wheat genome, which increased the genetic diversity of bread wheat and allowed its divergent adaptation. Meanwhile, convergent adaptation to human selection showed 2- to 16-fold enrichment relative to random expectation-a certain set of genes were repeatedly selected in Triticum species despite their drastic differences in ploidy levels and growing zones, indicating the important role of evolutionary constraints in shaping the adaptive landscape of bread wheat. These results showed the genetic necessities of wheat as a global crop and provided new perspectives on transferring adaptive success across species for crop improvement.

Resequencing of 145 Landmark Cultivars Reveals Asymmetric Sub-genome Selection and Strong Founder Genotype Effects on Wheat Breeding in China.

Controlled pedigrees and the multi-decade timescale of national crop plant breeding programs offer a unique experimental context for examining how selection affects plant genomes. More than 3000 wheat cultivars have been registered, released, and documented since 1949 in China. In this study, a set of 145 elite cultivars selected from historical points of wheat breeding in China were re-sequenced. A total of 43.75 Tb of sequence data were generated with an average read depth of 17.94× for each cultivar, and more than 60.92 million SNPs and 2.54 million InDels were captured, based on the Chinese Spring RefSeq genome v1.0. Seventy years of breeder-driven selection led to dramatic changes in grain yield and related phenotypes, with distinct genomic regions and phenotypes targeted by different breeders across the decades. There are very clear instances illustrating how introduced Italian and other foreign germplasm was integrated into Chinese wheat programs and reshaped the genomic landscape of local modern cultivars. Importantly, the resequencing data also highlighted significant asymmetric breeding selection among the three sub-genomes: this was evident in both the collinear blocks for homeologous chromosomes and among sets of three homeologous genes. Accumulation of more newly assembled genes in newer cultivars implied the potential value of these genes in breeding. Conserved and extended sharing of linkage disequilibrium (LD) blocks was highlighted among pedigree-related cultivars, in which fewer haplotype differences were detected. Fixation or replacement of haplotypes from founder genotypes after generations of breeding was related to their breeding value. Based on the haplotype frequency changes in LD blocks of pedigree-related cultivars, we propose a strategy for evaluating the breeding value of any given line on the basis of the accumulation (pyramiding) of beneficial haplotypes. Collectively, our study demonstrates the influence of "founder genotypes" on the output of breeding efforts over many decades and also suggests that founder genotype perspectives are in fact more dynamic when applied in the context of modern genomics-informed breeding.

Ethyl methanesulfonate mutant library construction in Gossypium hirsutum L. for allotetraploid functional genomics and germplasm innovation.

As the gene pool is exposed to both strain on land resources and a lack of diversity in elite allotetraploid cotton, the acquisition and identification of novel alleles has taken on epic importance in facilitating cotton genetic improvement and functional genomics research. Ethyl methanesulfonate (EMS) is an excellent mutagen that induces genome-wide efficient mutations to activate the mutagenic potential of plants with many advantages. The present study established, determined and verified the experimental procedure suitable for EMS-based mutant library construction as the general reference guide in allotetraploid upland cotton. This optimized method and procedure are efficient, and abundant EMS mutant libraries (approximately 12 000) in allotetraploid cotton were successfully obtained. More than 20 mutant phenotypes were observed and screened, including phenotypes of the leaf, flower, fruit, fiber and plant architecture. Through the plants mutant library, high-throughput and high-resolution melting technology-based variation evaluation detected the EMS-induced site mutation. Additionally, based on overall genome-wide mutation analyses by re-sequencing and mutant library assessment, the examination results demonstrated the ideal quality of the cotton EMS-treated mutant library constructed in this study with appropriate high mutation density and saturated genome. What is more, the collection is composed of a broad repertoire of mutants, which is the valuable resource for basic genetic research and functional genomics underlying complex allotetraploid traits, as well as cotton breeding.

Sympatric speciation of wild emmer wheat driven by ecology and chromosomal rearrangements.

In plants, the mechanism for ecological sympatric speciation (SS) is little known. Here, after ruling out the possibility of secondary contact, we show that wild emmer wheat, at the microclimatically divergent microsite of "Evolution Canyon" (EC), Mt. Carmel, Israel, underwent triple SS. Initially, it split following a bottleneck of an ancestral population, and further diversified to three isolated populations driven by disruptive ecological selection. Remarkably, two postzygotically isolated populations (SFS1 and SFS2) sympatrically branched within an area less than 30 m at the tropical hot and dry savannoid south-facing slope (SFS). A series of homozygous chromosomal rearrangements in the SFS1 population caused hybrid sterility with the SFS2 population. We demonstrate that these two populations developed divergent adaptive mechanisms against severe abiotic stresses on the tropical SFS. The SFS2 population evolved very early flowering, while the SFS1 population alternatively evolved a direct tolerance to irradiance by improved ROS scavenging activity that potentially accounts for its evolutionary fate with unstable chromosome status. Moreover, a third prezygotically isolated sympatric population adapted on the abutting temperate, humid, cool, and forested north-facing slope (NFS), separated by 250 m from the SFS wild emmer wheat populations. The NFS population evolved multiple resistant loci to fungal diseases, including powdery mildew and stripe rust. Our study illustrates how plants sympatrically adapt and speciate under disruptive ecological selection of abiotic and biotic stresses.