Chenyang Hao

Genetic Diversity and Linkage Disequilibrium in Chinese Bread Wheat (Triticum aestivum L.) Revealed by SSR Markers

Two hundred and fifty bread wheat lines, mainly Chinese mini core accessions, were assayed for polymorphism and linkage disequilibrium (LD) based on 512 whole-genome microsatellite loci representing a mean marker density of 5.1 cM. A total of 6,724 alleles ranging from 1 to 49 per locus were identified in all collections. The mean PIC value was 0.650, ranging from 0 to 0.965. Population structure and principal coordinate analysis revealed that landraces and modern varieties were two relatively independent genetic sub-groups. Landraces had a higher allelic diversity than modern varieties with respect to both genomes and chromosomes in terms of total number of alleles and allelic richness. 3,833 (57.0%) and 2,788 (41.5%) rare alleles with frequencies of <5% were found in the landrace and modern variety gene pools, respectively, indicating greater numbers of rare variants, or likely new alleles, in landraces. Analysis of molecular variance (AMOVA) showed that A genome had the largest genetic differentiation and D genome the lowest. In contrast to genetic diversity, modern varieties displayed a wider average LD decay across the whole genome for locus pairs with r2>0.05 (P<0.001) than the landraces. Mean LD decay distance for the landraces at the whole genome level was <5 cM, while a higher LD decay distance of 5–10 cM in modern varieties. LD decay distances were also somewhat different for each of the 21 chromosomes, being higher for most of the chromosomes in modern varieties (<5∼25 cM) compared to landraces (<5∼15 cM), presumably indicating the influences of domestication and breeding. This study facilitates predicting the marker density required to effectively associate genotypes with traits in Chinese wheat genetic resources.

Identifying Loci Influencing 1,000-Kernel Weight in Wheat by Microsatellite Screening for Evidence of Selection during Breeding

Chinese wheat mini core collection (262 accessions) was genotyped at 531 microsatellite loci representing a mean marker density of 5.1 cM. One-thousand-kernel weights (TKW) of lines were measured in five trials (three environments in four growing seasons). Structure analysis based on 42 unlinked SSR loci indicated that the materials formed two sub-populations, viz., landraces and modern varieties. A large difference in TKW (7.08 g, P<0.001) was found between the two sub-groups. Therefore, TKW is a major yield component that was improved in the past 6 decades; it increased from a mean 31.5 g in the 1940s to 44.64 g in the 2000s, representing a 2.19 g increase in each decade. Analyses based on a mixed linear model (MLM), population structure (Q) and relative kinship (K) revealed 22 SSR loci that were significantly associated with mean TKW (MTKW) of the five trials estimated by the best linear unbiased predictor (BLUP) method. They were mainly distributed on chromosomes of homoeologous groups 1, 2, 3, 5 and 7. Six loci, cfa2234-3A, gwm156-3B, barc56-5A, gwm234-5B, wmc17-7A and cfa2257-7A individually explained more than 11.84% of the total phenotypic variation. Favored alleles for breeding at the 22 loci were inferred according to their estimated effects on MTKW based on mean difference of varieties grouped by genotypes. Statistical simulation showed that these favored alleles have additive genetic effects. Frequency changes of alleles at loci associated with TKW are much more dramatic than those at neutral loci between the sub-groups. The numbers of favored alleles in modern varieties indicate there is still considerable genetic potential for their use as markers for genome selection of TKW in wheat breeding. Alleles that can be used globally to increase TKW were inferred according to their distribution by latitude and frequency of changes between landraces and the modern varieties.

Global Selection on Sucrose Synthase Haplotypes during a Century of Wheat Breeding

Endosperm starch synthesis pathway is a major target of indirect selection in global wheat breeding for higher yield. Spike number per unit area, number of grains per spike, and thousand kernel weight (TKW) are important yield components. In China, increases in wheat (Triticum aestivum) yields are mainly due to increases in grain number per spike and TKW. TKW mainly depends on starch content, as starch accounts for about 70% of the grain endosperm. Sucrose synthase catalysis is the first step in the conversion of sucrose to starch, that is, the conversion of sucrose to fructose and UDP-glucose by the wheat sucrose synthase genes (TaSus1 and TaSus2) that are located on chromosomes 7A/7B/7D and 2A/2B/2D, respectively. A total of 1,520 wheat accessions were genotyped at the six loci. Two, two, five, and two haplotypes were identified at the TaSus2-2A, TaSus2-2B, TaSus1-7A, and TaSus1-7B loci, respectively. Their main variations were detected within the introns. Significant differences between the haplotypes correlated with TKW differences among 348 modern Chinese cultivars from the core collection. Frequency changes for favored haplotypes showed gradual increases in cultivars released since beginning of the last century in China, Europe, and North America. Geographic distributions and time changes of favored haplotypes were characterized in six major wheat production regions worldwide. Strong selection bottlenecks to haplotype variations occurred at polyploidization and domestication and during breeding of wheat. Genetic-effect differences between haplotypes at the same locus influence the selection time and intensity. This work shows that the endosperm starch synthesis pathway is a major target of indirect selection in global wheat breeding for higher yield.

TaGS5‐3A, a grain size gene selected during wheat improvement for larger kernel and yield

Grain size is a dominant component of grain weight in cereals. Earlier studies have shown that OsGS5 plays a major role in regulating both grain size and weight in rice via promotion of cell division. In this study, we isolated TaGS5 homoeologues in wheat and mapped them on chromosomes 3A, 3B and 3D. Temporal and spatial expression analysis showed that TaGS5 homoeologues were preferentially expressed in young spikes and developing grains. Two alleles of TaGS5-3A, TaGS5-3A-T and TaGS5-3A-G were identified in wheat accessions, and a functional marker was developed to discriminate them. Association analysis revealed that TaGS5-3A-T was significantly correlated with larger grain size and higher thousand kernel weight. Biochemical assays showed that TaGS5-3A-T possesses a higher enzymatic activity than TaGS5-3A-G. Transgenic rice lines overexpressing TaGS5-3A-T also exhibited larger grain size and higher thousand kernel weight than TaGS5-3A-G lines, and the transcript levels of cell cycle-related genes in TaGS5-3A-T lines were higher than those in TaGS5-3A-G lines. Furthermore, systematic evolution analysis in diploid, tetraploid and hexaploid wheat showed that TaGS5-3A underwent strong artificial selection during wheat polyploidization events and the frequency changes of two alleles demonstrated that TaGS5-3A-T was favoured in global modern wheat cultivars. These results suggest that TaGS5-3A is a positive regulator of grain size and its favoured allele TaGS5-3A-T exhibits a larger potential application in wheat high-yield breeding.

A yield-associated gene TaCWI, in wheat: its function, selection and evolution in global breeding revealed by haplotype analysis

Key messageWheat anther-specific invertase genes were haplotyped in wheat. Strong allelic selection occurred during wheat polyploidization, domestication and breeding because of their association with yield traits.AbstractPlant invertase hydrolyzes sucrose into glucose and fructose. Cell wall invertase (CWI), one of the three types of invertase, is essential for plant development. Based on isolated TaCWI genes from chromosomes 4A, 5B and 5D, two SNPs were detected in the promoter region of TaCWI-4A, and four SNPs and two Indels were present in the TaCWI-5D gene. No polymorphism was detected in TaCWI-5B coding or promoter regions. CAPS markers caps4A and caps5D were developed to discriminate haplotypes of TaCWI-4A and TaCWI-5D. Marker/trait association analysis indicated that Hap-5D-C at TaCWI-5D was significantly associated with higher thousand kernel weight (TKW) in 348 Chinese modern cultivars grown in multiple environments. Geographic distributions and changes over time of favored haplotypes showed that Hap-5D-C was the most frequent haplotype in modern cultivars and was strongly positively selected in six major wheat production regions worldwide. However, selection for haplotypes at TaCWI-4A was not so evident, possibly due to balancing effects of the two haplotypes on TKW and grain number per spike (GN). In rainfed production regions, Hap-4A-C was favored because it brought more seeds, but in well irrigated conditions, Hap-4A-T was favored in modern breeding because of higher TKW. Evolutionary analysis among wheat and its relatives showed that genetic diversity of TaCWI genes on chromosomes 4A and 5D declined dramatically in progression from the diploid level to modern polyploid cultivars. There was strong allelic selection during polyploidization, domestication and breeding.

Identifying loci influencing grain number by microsatellite screening in bread wheat (Triticum aestivum L.)

Grain number (GN) is one of three major yield-related components in wheat. We used the Chinese wheat mini core collection to undertake a genome-wide association analysis of grain number using 531 SSR markers randomly located on all 21 chromosomes. Grain numbers of all accessions were measured in four trials, i.e. two environments in four growing seasons. Association analysis based on a mixed linear model (MLM) revealed that 27 SSR loci were significantly associated with mean GN (MGN) estimated by the best linear unbiased predictor (BLUP) method. These included numerous breeder favorable alleles with strong positive effects at 23 loci. Significant or extremely significant differences were detected on MGN between varieties conveying favored allele and varieties with other alleles. Moreover, statistical simulation showed that the favored alleles have additive genetic effects. Although modern varieties combined larger numbers of favored alleles, the numbers of favored alleles were not significantly different from those in landraces, especially those alleles contributing mostly to the phenotypic variation. These results indicate that there is still considerable genetic potential for use of markers for genome selection of GN for high yield in wheat.

Hitchhiking Effect Mapping: A New Approach for Discovering Agronomic Important Genes

Besides the natural selection, the crops cultivated today have experienced two episodes of strong artificial selection, domestic and modern breeding. Domestication led to giant genetic structure differentiation between cultivars and their wild species, while modern breeding made further genetic structure differentiation between the modern varieties and the landraces. In a population, diversity of the loci under strong selection is significantly lower than that of other loci. At the same time, diversity in the genomic regions flanking these selected loci also declines in the process of selection. This phenomenon is called hitchhiking effects or selection sweep in genetics. Genomic regions with selection sweep (haplotype block) could be detected after draft genome scanning (genome typing) with molecular markers in a number of released varieties or natural populations. Marker/trait association analysis in these regions would detect the loci (or QTLs) even the favored alleles (genes) in breeding or natural adaptation. Fine scanning of these genomic regions would help to determine the sizes of haplotype blocks and to discover the key genes, thereby providing very valuable information for isolation of the key genes and molecular design of new varieties. Establishment of high density genetic linkage maps in the major crops and availability of high throughput genotyping platform make it possible to discover agronomic important genes through marker/trait association analysis. On the basis of available publications, we give a brief introduction of the hitchhiking effect mapping approach in this paper using plant height, 1000-grain weight, and phosphorus-deficiency tolerance as examples in wheat.

Association Mapping and Haplotype Analysis of a 3.1-Mb Genomic Region Involved in Fusarium Head Blight Resistance on Wheat Chromosome 3BS

A previous study provided an in-depth understanding of molecular population genetics of European and Asian wheat gene pools using a sequenced 3.1-Mb contig (ctg954) on chromosome 3BS. This region is believed to carry the Fhb1 gene for response to Fusarium head blight. In this study, 266 wheat accessions were evaluated in three environments for Type II FHB response based on the single floret inoculation method. Hierarchical clustering (UPGMA) based on a Manhattan dissimilarity matrix divided the accessions into eight groups according to five FHB-related traits which have a high correlation between them; Group VIII comprised six accessions with FHB response levels similar to variety Sumai 3. Based on the compressed mixed linear model (MLM), association analysis between five FHB-related traits and 42 molecular markers along the 3.1-Mb region revealed 12 significant association signals at a threshold of P<0.05. The highest proportion of phenotypic variation (6.2%) in number of diseased spikelets (NDS) occurred at locus cfb6059, and the physical distance was about 2.9 Kb between umn10 and this marker. Haplotype block (HapB) analysis using a sliding window LD of 5 markers, detected six HapBs in the 3.1-Mb region at r2>0.1 and P<0.001 between random closely linked markers. F-tests among Haps with frequencies >0.05 within each HapB at r2>0.1 and P<0.001 showed significant differences between the Hap carried by FHB resistant resources, such as Sumai 3 and Wangshuibai, and susceptible genotypes in HapB3 and HapB6. These results suggest that Fhb1 is located within HapB6, with the possibility that another gene is located at or near HapB3. SSR markers and Haps detected in this study will be helpful in further understanding the genetic basis of FHB resistance, and provide useful information for marker-assisted selection of Fhb1 in wheat breeding.

Small RNA and Degradome Sequencing Reveal Complex Roles of miRNAs and Their Targets in Developing Wheat Grains

Plant microRNAs (miRNAs) have been shown to play critical roles in plant development. In this study, we employed small RNA combined with degradome sequencing to survey development-related miRNAs and their validated targets during wheat grain development. A total of 186 known miRNAs and 37 novel miRNAs were identified in four small RNA libraries. Moreover, a miRNA-like long hairpin locus was first identified to produce 21~22-nt phased siRNAs that act in trans to cleave target mRNAs. A comparison of the miRNAomes revealed that 55 miRNA families were differentially expressed during the grain development. Predicted and validated targets of these development-related miRNAs are involved in different cellular responses and metabolic processes including cell proliferation, auxin signaling, nutrient metabolism and gene expression. This study provides insight into the complex roles of miRNAs and their targets in regulating wheat grain development.

TaGW2, a Good Reflection of Wheat Polyploidization and Evolution

Hexaploid wheat consists of three subgenomes, namely, A, B, and D. These well-characterized ancestral genomes also exist at the diploid and tetraploid levels, thereby rendering wheat as a good model species for studying polyploidization. Here, we performed intra- and inter-species comparative analyses of wheat and its relatives to dissect polymorphism and differentiation of the TaGW2 genes. Our results showed that genetic diversity of TaGW2 decreased with progression from the diploids to tetraploids and hexaploids. The strongest selection occurred in the promoter regions of TaGW2-6A and TaGW2-6B. Phylogenetic trees clearly indicated that Triticum urartu and Ae. speltoides were the donors of the A and B genomes in tetraploid and hexaploid wheats. Haplotypes detected among hexaploid genotypes traced back to the tetraploid level. Fst and π values revealed that the strongest selection on TaGW2 occurred at the tetraploid level rather than in hexaploid wheat. This infers that grain size enlargement, especially increased kernel width, mainly occurred in tetraploid genotypes. In addition, relative expression levels of TaGW2s significantly declined from the diploid level to tetraploids and hexaploids, further indicating that these genes negatively regulate kernel size. Our results also revealed that the polyploidization events possibly caused much stronger differentiation than domestication and breeding.