Two Novel AP2/EREBP Transcription Factor Genes TaPARG Have Pleiotropic Functions on Plant Architecture and Yield-Related Traits in Common Wheat

Genome-wide identification and evolutionary analysis of the AP2/EREBP, COX and LTP genes in Zea mays L. under drought stress

AP2 (APETALA2)/EREBP (ethylene-responsive element-binding protein), cytochrome c oxidase (COX) and nonspecific lipid transfer proteins (LTP) play important roles in the response to drought stress. This is the first study to identify the COX gene in Zea mays L. via genome-wide analysis. The qRT‒PCR results indicated that AP2/EREBP, COX and LTP were downregulated, with fold changes of 0.84, 0.53 and 0.31, respectively, after 12 h of drought stress. Genome-wide analysis identified 78 AP2/EREBP, 6 COX and 10 LTP genes in Z. mays L. Domain analysis confirmed the presence of the AP2 domain, Cyt_c_Oxidase_Vb domain and nsLTP1 in the AP2/EREBP, COX and LTP proteins, respectively. The AP2/EREBP protein family (AP2) includes five different domain types: the AP2/ERF domain, the EREBP-like factor (EREBP), the ethylene responsive factor (ERF), the dehydration responsive element binding protein (DREB) and the SHN SHINE. Synteny analysis of the AP2/EREBP, COX and LTP genes revealed collinearity orthologous relationships in O. sativa, H. vulgare and A. thaliana. AP2/EREBP genes were found on the 10 chromosomes of Z. mays L. COX genes were found on chromosomes 1, 3, 4, 5, 7 and 8. LTP genes were found on chromosomes 1, 3, 6, 8, 9 and 10. In the present study, the Ka/Ks ratios of the AP2/EREBP paralogous pairs indicated that the AP2/EREBP genes were influenced primarily by purifying selection, which indicated that the AP2/EREBP genes received strong environmental pressure during evolution. The Ka/Ks ratios of the COX-3/COX-4 paralogous pairs indicate that the COX-3/COX-4 genes were influenced primarily by Darwinian selection (driving change). For the LTP genes, the Ka/Ks ratios of the LTP-1/LTP-10, LTP-5/LTP-3 and LTP-4/LTP-8 paralogous pairs indicate that these genes were influenced primarily by purifying selection, while the Ka/Ks ratios of the LTP-2/LTP-6 paralogous pairs indicate that these genes were influenced primarily by Darwinian selection. The duplication time of the AP2/EREBP paralogous gene pairs in Z. mays L. ranged from approximately 9.364 to 100.935 Mya. The duplication time of the COX-3/COX-4 paralogous gene pair was approximately 5.217 Mya. The duplication time of the LTP paralogous gene pairs ranged from approximately 19.064 to 96.477 Mya. The major focus of research is to identify the genes that are responsible for drought stress tolerance to improve maize for drought stress tolerance. The results of the present study will improve the understanding of the functions of the AP2/EREBP, COX and LTP genes in response to drought stress.

TaNAC020 homoeologous genes are associated with higher thousand kernel weight and kernel length in Chinese wheat

NAC proteins constitute one of the largest plant-specific transcription factor (TF) families and play significant roles in plant growth and development. In the present study, three TaNAC020 homoeologous genes located on chromosomes 7A, 7B, and 7D were isolated from wheat (Triticum aestivum L.). TaNAC020s were predominantly expressed in developing grains. The developed transgenic rice lines for TaNAC020-B showed higher starch density and lower amylose contents than those of the wild type (WT). Sequence polymorphism studies showed seven and eight SNPs in TaNAC020-A/B, making three and two haplotypes, respectively. No sequence polymorphism was identified in TaNAC020-D. Association analysis revealed that HAP-2 of TaNAC020-A and TaNAC020-B was the favored haplotype for higher thousand kernel weight and length. Geographic distribution and allelic frequency showed that our favored haplotype experienced strong selection in China, and likewise, diversity increased in TaNAC020s during wheat polyploidization. The results obtained in this study demonstrate that TaNAC020s positively influence starch synthesis and accumulation and are one of the key regulators of the kernel (seed) size and kernel number and have the potential for utilization in wheat breeding to improve grain yield. Molecular markers developed in this study stand instrumental in marker-assisted selection for genetic improvement and germplasm enhancement in wheat.

QTL Analysis for Bread Wheat Seed Size, Shape and Color Characteristics Estimated by Digital Image Processing

The size, shape, and color of wheat seeds are important traits that are associated with yield and flour quality (size, shape), nutritional value, and pre-harvest sprouting (coat color). These traits are under multigenic control, and to dissect their molecular and genetic basis, quantitative trait loci (QTL) analysis is used. We evaluated 114 recombinant inbred lines (RILs) in a bi-parental RIL mapping population (the International Triticeae Mapping Initiative, ITMI/MP) grown in 2014 season. We used digital image analysis for seed phenotyping and obtained data for seven traits describing seed size and shape and 48 traits of seed coat color. We identified 212 additive and 34 pairs of epistatic QTLs on all the chromosomes of wheat genome except chromosomes 1A and 5D. Many QTLs were overlapping. We demonstrated that the overlap between QTL regions was low for seed size/shape traits and high for coat color traits. Using the literature and KEGG data, we identified sets of genes in Arabidopsis and rice from the networks controlling seed size and color. Further, we identified 29 and 14 candidate genes for seed size-related loci and for loci associated with seed coat color, respectively.

Identification of TaBADH-A1 allele for improving drought resistance and salt tolerance in wheat (Triticum aestivum L.)

Drought and salt stress can strongly affect the growth and development of wheat. Wheat adapts to drought and salt stress through osmotic regulation. Betaine aldehyde dehydrogenase (BADH) is a key enzyme in the synthesis of betaine, an osmotic regulator. We cloned a region of the TaBADH-A1 promoter and genomic DNA that included the introns and exons, from four Chinese wheat cultivars. Following the analysis of TaBADH-A1 genomic DNA and promoter sequence polymorphisms of 4 cloned and 15 cultivars from the database, 7 haplotypes of TaBADH-A1 gene were identified. We divided the 7 haplotypes with a 254 bp insertion or deletion (indel) into two main alleles, BADH-A1a and BADH-A1b. Meanwhile, a molecular marker was developed based on the 254 bp indel of the third intron of TaBADH-A1 gene. Expression levels of BADH-A1b were found to be significantly higher than those of BADH-A1a under drought and salt stress conditions. Betaine accumulation was significantly higher in wheat containing BADH-A1b compared to BADH-A1a under drought and salt stress. We also identified that the average relative germination and survival rates of wheat with the BADH-A1b allele were significantly higher than wheat with the BADH-A1a allele. The results reveal that wheat containing BADH-A1b has stronger drought and salt tolerance than wheat with BADH-A1a. Meanwhile, the geographic distribution and frequency of the TaBADH-A1 locus alleles indicate that BADH-A1a has been preferred in Chinese wheat breeding programs, while BADH-A1b, associated with favorable stress tolerance, has been neglected. The results of this study provide evidence for an excellent candidate allele for marker-assisted selection of new wheat cultivars with increased salt tolerance and drought resistance.

Comparative Genetic Analysis of Durum Wheat Landraces and Cultivars Widespread in Tunisia

The durum wheat (Triticum turgidum L. ssp. durum Desf.) landraces constitute a useful natural germplasm to increase the genetic diversity in the modern durum cultivars. The Tunisian durum germplasm constitutes 28 accessions conserved in Genebank of Tunisia, which are still unexplored. In this study, a comparative genetic analysis was performed to investigate the relationships between the Tunisian durum lines and the modern cultivars and detect divergent loci involved in breeding history. The genetic diversity analyses carried out using nine morphological descriptors and the 25K single-nucleotide polymorphism (SNP) array allowed us to distinguish two groups of Tunisian landraces and one of durum cultivars. The analysis of molecular variance and diversity indices confirmed the genetic variability among the groups. A total of 529 SNP loci were divergent between Tunisian durum landraces and modern cultivars. Candidate genes related to plant and spike architecture, including FLOWERING LOCUS T (FT-B1), zinc finger CONSTANS, and AP2/EREBPs transcription factors, were identified. In addition, divergent genes involved in grain composition and biotic stress nucleotide-binding site and leucine-reach repeats proteins and disease resistance proteins (NBS-LRR and RPM) were found, suggesting that the Tunisian durum germplasm may represent an important source of favorable alleles to be used in future durum breeding programs for developing well-adapted and resilient cultivars.

Grain Transcriptome Dynamics Induced by Heat in Commercial and Traditional Bread Wheat Genotypes

High temperature (HT) events have negative impact on wheat grains yield and quality. Transcriptome profiles of wheat developing grains of commercial genotypes (Antequera and Bancal) and landraces (Ardito and Magueija) submitted to heatwave-like treatments during grain filling were evaluated. Landraces showed significantly more differentially expressed genes (DEGs) and presented more similar responses than commercial genotypes. DEGs were more associated with transcription and RNA and protein synthesis in Antequera and with metabolism alterations in Bancal and landraces. Landraces upregulated genes encoding proteins already described as HT responsive, like heat shock proteins and cupins. Apart from the genes encoding HSP, two other genes were upregulated in all genotypes, one encoding for Adenylate kinase, essential for the cellular homeostasis, and the other for ferritin, recently related with increased tolerance to several abiotic stress in Arabidopsis. Moreover, a NAC transcription factor involved in plant development, known to be a negative regulator of starch synthesis and grain yield, was found to be upregulated in both commercial varieties and downregulated in Magueija landrace. The detected diversity of molecular processes involved in heat response of commercial and traditional genotypes contribute to understand the importance of genetic diversity and relevant pathways to cope with these extreme events.

Identification of the yield traits related haplotype combinations of transcription factor genes TaHDZ34 in common wheat

A predominant objective in wheat breeding is improving yield-related traits. The homeodomain-leucine zipper (HD-Zip) transcription factor plays a significant role in plant growth and development. In this study, we cloned all homeologs of TaHDZ34, which is a member of the HD-Zip class IV transcription factor family in wheat (Triticum aestivum L.). Sequence polymorphism analysis showed that TaHDZ-A34, TaHDZ-B34, and TaHDZ-D34 formed five, six, and six haplotypes, respectively, and the genes were divided into two main haplotype groups. We also developed functional molecular markers. The TaHDZ34 genes were divided into eight main haplotype combinations. Association analysis and distinct population validation preliminarily indicated that TaHDZ34 genes modulate grain number per spike, effective spikelet number per spike, thousand kernel weight, and flag leaf area per plant in wheat. Hap-ABD was the most effective haplotype combination of TaHDZ34. Subcellular localization showed that TaHDZ-A34 was localized to the nucleus. The interacting proteins of TaHDZ-A34 were involved in protein synthesis/degradation, energy production and transportation, and photosynthesis. Geographic distribution and frequencies of TaHDZ34 haplotype combinations suggested that Hap-Abd and Hap-AbD were preferentially selected in Chinese wheat breeding programs. The high-yield-related haplotype combination Hap-ABD provided beneficial genetic resources for the marker-assisted selection of new wheat cultivars.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01298-5.

Identification of Candidate Chromosome Region Related to Melon (Cucumis melo L.) Fruit Surface Groove Trait Through Biparental Genetic Mapping and Genome-Wide Association Study

The melon fruit surface groove (fsg) not only affects peel structure and causes stress-induced fruit cracking but also fits consumers' requirements in different regions. In this study, genetic inheritance analysis of three F₂ populations derived from six parental lines revealed that the fsg trait is controlled by a simple recessive inherited gene. Through bulked segregant analysis sequencing (BSA-seq), the Cmfsg locus was detected in an 8.96 Mb interval on chromosome 11 and then initially mapped to a region of approximately 1.15 Mb. Further fine mapping with a large F₂ population including 1,200 plants narrowed this region to 207 kb containing 11 genes. A genome-wide association study (GWAS) with 187 melon accessions also produced the same chromosome region for the Cmfsg locus. Due to the rare molecular markers and lack of mutations in the coding and promoter regions of the 11 candidate genes in the fine-mapped interval, we conducted in silico BSA to explore the natural melon panel to predict candidate genes for the Cmfsg locus. A 1.07 kb segment upstream of MELO3C019694.2 (annotated as the AGAMOUS MADS-box transcription factor) exhibited a correlation with the grooved and non-grooved accessions among the F₂ individuals, and a natural panel consisted of 17 melon accessions. The expression level of MELO3C019694.2 in the pericarp was higher in grooved lines than in non-grooved lines and was specifically expressed in fruit compared with other tissues (female flower, male flower, root, and leaf). This work provides fundamental information for further research on melon fsg trait formation and molecular markers for melon breeding.

Transcriptomic Identification of Wheat AP2/ERF Transcription Factors and Functional Characterization of TaERF-6-3A in Response to Drought and Salinity Stresses

AP2/ERF (APETALA2/ethylene responsive factor) is a family of plant-specific transcription factors whose members are widely involved in many biological processes, such as growth, development, and biotic and abiotic stress responses. Here, 20 AP2/ERF genes were identified based on wheat RNA-seq data before and after drought stress, and classified as AP2, ERF, DREB, and RAV. The analysis of gene structure revealed that about 85% of AP2/ERF family members had lost introns, which are presumed to have been lost during the formation and evolution of the wheat genome. The expression of 20 AP2/ERF family genes could be verified by qRT-PCR, which further supported the validity of the RNA-seq data. Subsequently, subcellular localization and transcriptional activity experiments showed that the ERF proteins were mainly located in the nucleus and were self-activating, which further supports their functions as transcription factors. Furthermore, we isolated a novel ERF gene induced by drought, salt, and cold stresses and named it TaERF-6-3A. TaERF-6-3A overexpression increased sensitivity to drought and salt stresses in Arabidopsis, which was supported by physiological and biochemical indices. Moreover, the expression of stress- and antioxidant-related genes was downregulated in TaERF-6-3A-overexpressing plants. Overall, these results contribute to the further understanding of the TaERF-6-3A gene function in wheat.

Identification of Single Nucleotide Polymorphism in TaSBEIII and Development of KASP Marker Associated With Grain Weight in Wheat

Manipulation of genes involved in starch synthesis could significantly affect wheat grain weight and yield. The starch-branching enzyme (SBE) catalyzes the formation of branch points by cleaving the α-1,4 linkage in polyglucans and reattaching the chain via an α-1,6 linkage. Three types of SBE isoforms (SBEI, SBEII, and SBEIII) exist in higher plants, with the number of SBE isoforms being species-specific. In this study, the coding sequence of the wheat TaSBEIII gene was amplified. After the multiple sequence alignment of TaSBEIII genome from 20 accessions in a wheat diversity panel, one SNP was observed in TaSBEIII-A, which formed the allelic marker allele-T. Based on this SNP at 294 bp (C/T), a KASP molecular marker was developed to distinguish allelic variation among the wheat genotypes for thousand grain weight (TGW). The results were validated using 262 accessions of mini core collection (MCC) from China, 153 from Pakistan, 53 from CIMMYT, and 17 diploid and 18 tetraploid genotypes. Association analysis between TaSBEIII-A allelic variation and agronomic traits found that TaSBEIII-A was associated with TGW in mini core collection of China (MCC). The accessions possessing Allele-T had higher TGW than those possessing Allele-C; thus, Allele-T was a favorable allelic variation. By analyzing the frequency of the favorable allelic variation Allele-T in MCC, it increased from pre-1950 (25%) to the 1960s (45%) and increased continuously from 1960 to 1990 (80%). The results suggested that the KASP markers can be utilized in grain weight improvement, which ultimately improves wheat yield by marker-assisted selection in wheat breeding. The favorable allelic variation allele-T should be valuable in enhancing grain yield by improving the source and sink simultaneously. Furthermore, the newly developed KASP marker validated in different genetic backgrounds could be integrated into a breeding kit for screening high TGW wheat.

Genome-Wide Association Mapping and Gene Expression Analysis Reveal the Negative Role of OsMYB21 in Regulating Bacterial Blight Resistance in Rice

BACKGROUND

Bacterial blight (BB), caused by Xanthomonas oryzae pv. oryzae (Xoo), is one of the most devastating diseases in rice all over the world. Due to the diversity and rapid evolution of Xoo, identification and use of the non-race specific quantitative resistance QTLs has been considered the preferred strategy for effective control of this disease. Although numerous QTLs for BB resistance have been identified, they haven't been effectively used for improvement of BB resistance in rice due to their small effects and lack of knowledge on the function of genes underlying the QTLs.

RESULTS

In the present study, a genome-wide association study of BB resistance was performed in a rice core collection from South China. A total of 17 QTLs were identified to be associated with BB resistance. Among them, 13 QTLs were newly identified in the present study and the other 4 QTLs were co-localized with the previously reported QTLs or Xa genes that confer qualitative resistance to Xoo strains. Particularly, the qBBR11-4 on chromosome 11 explained the largest phenotypic variation in this study and was co-localized with the previously identified QTLs for BB and bacterial leaf streak (BLS) resistance against diverse strains in three studies, suggesting its broad-spectrum resistance and potential value in rice breeding. Through combined analysis of differential expression and annotations of the predicted genes within qBBR11-4 between two sets of rice accessions selected based on haplotypes and disease phenotypes, we identified the transcription factor OsMYB21 as the candidate gene for qBBR11-4. The OsMYB21 overexpressing plants exhibited decreased resistance to bacterial blight, accompanied with down-regulation of several defense-related genes compared with the wild-type plants.

CONCLUSION

The results suggest that OsMYB21 negatively regulates bacterial blight resistance in rice, and this gene can be a promising target in rice breeding by using the gene editing method. In addition, the potential candidate genes for the 13 novel QTLs for BB resistance were also analyzed in this study, providing a new source for cloning of genes associated with BB resistance and molecular breeding in rice.

Development and Exploitation of KASP Assays for Genes Underpinning Drought Tolerance Among Wheat Cultivars From Pakistan

High-throughput genotyping for functional markers offers an excellent opportunity to effectively practice marker-assisted selection (MAS) while breeding cultivars. We developed kompetitive allele-specific PCR (KASP) assays for genes conferring drought tolerance in common wheat (Triticum aestivum L.). In total, 11 KASP assays developed in this study and five already reported assays were used for their application in wheat breeding. We investigated alleles at 16 loci associated with drought tolerance among 153 Pakistani hexaploid wheat cultivars released during 1953-2016; 28 diploid wheat accessions (16 for AA and 12 for BB) and 19 tetraploid wheat (AABB) were used to study the evolutionary history of the studied genes. Superior allelic variations of the studied genes were significantly associated with higher grain yield. Favored haplotypes of TaSnRK2.3-1A, TaSnRK2.3-1B, TaSnRK2.9-5A, TaSAP-7B, and TaLTPs-1A predominated in Pakistani wheat germplasm indicating unconscious pyramiding and selection pressure on favorable haplotypes during selection breeding. TaSnRK2.8-5A, TaDreb-B1, 1-feh w3, TaPPH-7A, TaMOC-7A, and TaPARG-2A had moderate to low frequencies of favorable haplotype among Pakistani wheat germplasm pointing toward introgression of favorable haplotypes by deploying functional markers in marker-assisted breeding. The KASP assays were compared with gel-based markers for reliability and phenotypically validated among 62 Pakistani wheat cultivars. Association analyses showed that the favorable allelic variations were significantly associated with grain yield-contributing traits. The developed molecular marker toolkit of the genes can be instrumental for the wheat breeding in Pakistan.

Identification of a novel ERF gene, TaERF8, associated with plant height and yield in wheat

BACKGROUND

Ethylene Responsive Factor (ERF) is involved in various processes of plant development and stress responses. In wheat, several ERFs have been identified and their roles in mediating biotic or abiotic stresses have been elucidated. However, their effects on wheat plant architecture and yield-related traits remain poorly studied.

RESULTS

In this study, TaERF8, a new member of the ERF family, was isolated in wheat (Triticum aestivum L.). Three homoeologous TaERF8 genes, TaERF8-2A, TaERF8-2B and TaERF8-2D (named according to sub-genomic origin), were cloned from the common wheat cultivar Chinese Spring. The three homoeologs showed highly similar protein sequences, with identical AP2 domain. Whereas homoeologs sequence polymorphism analysis allowed the establishment of ten, two and three haplotypes, respectively. Expression analysis revealed that TaERF8s were constitutively expressed through entire wheat developmental stages. Analysis of related agronomic traits of TaERF8-2B overexpressing transgenic lines showed that TaERF8-2B plays a role in regulating plant architecture and yield-related traits. Association analysis between TaERF8-2B haplotypes (Hap-2B-1 and Hap-2B-2) and agronomic traits showed that TaERF8-2B was associated with plant height, heading date and 1000 kernel weight (TKW). The TaERF8-2B haplotypes distribution analysis revealed that Hap-2B-2 frequency increased in domesticated emmer wheat and modern varieties, being predominant in five major China wheat producing zones.

CONCLUSION

These results indicated that TaERF8s are differentially involved in the regulation of wheat growth and development. Haplotype Hap-2B-2 was favored during domestication and in Chinese wheat breeding. Unveiling that the here described molecular marker TaERF8-2B-InDel could be used for marker-assisted selection, plant architecture and TKW improvement in wheat breeding.

Association mapping for yield traits in Paeonia rockii based on SSR markers within transcription factors of comparative transcriptome

BACKGROUND

Allelic variation underlying the quantitative traits in plants is caused by the extremely complex regulation process. Tree peony originated in China is a peculiar ornamental, medicinal and oil woody plant. Paeonia rockii, one of tree peony species, is a precious emerging woody oil crop. However, in this valuable plant, the study of functional loci associated with yield traits has rarely been identified. Therefore, to explore the genetic architecture of 24 yield quantitative traits, the association mapping was first reported in 420 unrelated cultivated P. rockii individuals based on the next-generation sequencing (NGS) and single-molecule long-read sequencing (SMLRS).

RESULTS

The developed 58 pairs of polymorphic expressed sequence tag-simple sequence repeat (EST-SSR) markers from 959 candidate transcription factors (TFs) associated with yield were used for genotyping the 420 P. rockii accessions. We observed a high level of genetic diversity (polymorphic information content, PIC = 0.514) and low linkage disequilibrium (LD) between EST-SSRs. Moreover, four subpopulations in the association population were revealed by STRUCTURE analyses. Further, single-marker association analysis identified 141 significant associations, involving 17 quantitative traits and 41 EST-SSRs. These loci were mainly from AP2, TCP, MYB, HSF, bHLH, GATA, and B3 gene families and showed a small proportion of the phenotypic variance (3.79 to 37.45%).

CONCLUSIONS

Our results summarize a valuable collection of functional loci associated with yield traits in P. rockii, and provide a precious resource that reveals allelic variation underlying quantitative traits in Paeonia and other woody oil crops.

Genome-wide association study of Striga resistance in early maturing white tropical maize inbred lines

BACKGROUND

Striga hermonthica (Benth.) parasitism militates against increased maize production and productivity in savannas of sub-Saharan Africa (SSA). Identification of Striga resistance genes is important in developing genotypes with durable resistance. So far, there is only one report on the existence of QTL for Striga resistance on chromosome 6 of maize. The objective of this study was to identify genomic regions significantly associated with grain yield and other agronomic traits under artificial Striga field infestation. A panel of 132 early-maturing maize inbreds were phenotyped for key agronomic traits under Striga-infested and Striga-free conditions. The inbred lines were also genotyped using 47,440 DArTseq markers from which 7224 markers were retained for population structure analysis and genome-wide association study (GWAS).

RESULTS

The inbred lines were grouped into two major clusters based on structure analysis as well as the neighbor-joining hierarchical clustering. A total of 24 SNPs significantly associated with grain yield, Striga damage at 8 and 10 weeks after planting (WAP), ears per plant and ear aspect under Striga infestation were detected. Under Striga-free conditions, 11 SNPs significantly associated with grain yield, number of ears per plant and ear aspect were identified. Three markers physically located close to the putative genes GRMZM2G164743 (bin 10.05), GRMZM2G060216 (bin 3.06) and GRMZM2G103085 (bin 5.07) were detected, linked to grain yield, Striga damage at 8 and 10 WAP and number of ears per plant under Striga infestation, explaining 9 to 42% of the phenotypic variance. Furthermore, the S9_154,978,426 locus on chromosome 9 was found at 2.61 Mb close to the ZmCCD1 gene known to be associated with the reduction of strigolactone production in the maize roots.

CONCLUSIONS

Presented in this study is the first report of the identification of significant loci on chromosomes 9 and 10 of maize that are closely linked to ZmCCD1 and amt5 genes, respectively and may be related to plant defense mechanisms against Striga parasitism. After validation, the identified loci could be targets for breeders for marker-assisted selection (MAS) to accelerate genetic enhancement of maize for Striga resistance in the tropics, particularly in SSA, where the parasitic weed is endemic.

Genome-Wide Identification and Analysis of the AP2 Transcription Factor Gene Family in Wheat (Triticum aestivum L.)

The AP2 transcription factors play important roles in regulating plant growth and development. However, limited data are available on the contributions of AP2 transcription factors in wheat (Triticum aestivum L.). In the present study, a total of 62 AP2 genes were identified in wheat from a genome-wide search against the latest wheat genome data. Phylogenetic and sequence alignment analyses divided the wheat AP2 genes into 3 clusters, euAP2, euANT, and basalANT. Chromosomal distribution, gene structure and duplication, and motif composition were subsequently investigated. The 62 TaAP2 genes were unevenly distributed on 21 chromosomes. Twenty-four homologous gene sets among A, B, and D sub-genomes were detected, which contributed to the expansion of the wheat AP2 gene family. The expression levels of TaAP2 genes were examined using the WheatExp database; most detected genes exhibited tissue-specific expression patterns. The transcript levels of 9 randomly selected TaAP2 genes were validated through qPCR analyses. Overexpression of TaAP2-10-5D, the most likely homolog of Arabidopsis ANT gene, increased organ sizes in Arabidopsis. Our results extend our knowledge of the AP2 gene family in wheat, and contribute to further functional characterization of AP2s during wheat development with the ultimate goal of improving crop production.

Identification of TaPPH-7A haplotypes and development of a molecular marker associated with important agronomic traits in common wheat

BACKGROUND

Premature senescence of flag leaf severely affects wheat yield and quality. Chlorophyll (Chl) degradation is the most obvious symptom during leaf senescence and catalyzed by a series of enzymes. Pheophytin pheophorbide hydrolase (Pheophytinase, PPH) gene encodes a Chl degradation hydrolase.

RESULTS

In this study, the coding, genomic and promoter sequences of wheat TaPPH-A gene were cloned. The corresponding lengths were 1467 bp, 4479 bp and 3666 bp, respectively. Sequence structure analysis showed that TaPPH-A contained five exons and four introns. After the multiple sequences alignment of TaPPH-A genome from 36 accessions in a wheat diversity panel, four SNPs and one 2-bp InDel were observed, which formed two haplotypes, TaPPH-7A-1 and TaPPH-7A-2. Based on the SNP at 1299 bp (A/G), a molecular marker TaPPH-7A-dCAPS was developed to distinguish allelic variation (A/G). Using the molecular markers, 13 SSR, and 116 SNP markers, a linkage map of chromosome 7A were integrated. TaPPH-A was mapped on the chromosome region flanked by Xwmc9 (0.94 cM) and AX-95634545 (1.04 cM) on 7A in a DH population. Association analysis between TaPPH-7A allelic variation and agronomic traits found that TaPPH-7A was associated with TGW in 11 of 12 environments and Chl content at grain-filling stage under drought stress using Population 1 consisted of 323 accessions. The accessions possessed TaPPH-7A-1 (A) had higher TGW and Chl content than those possessed TaPPH-7A-2 (G), thus TaPPH-7A-1 (A) was a favorable allelic variation. By analyzing the frequency of favorable allelic variation TaPPH-7A-1 (A) in Population 2 with 157 landraces and Population 3 with 348 modern cultivars, we found it increased from pre-1950 (0) to 1960s (54.5%), then maintained a relatively stable level about 56% from 1960s to 1990s.

CONCLUSION

These results suggested the favorable allelic variation TaPPH-7A-1 (A) should be valuable in enhancing grain yield by improving the source (chlorophyll content) and sink (the developing grain) simultaneously. Furthermore, the newly developed molecular marker TaPPH-7A-dCAPS could be integrated into a breeding kit of screening high TGW wheat for marker-assisted selection.

EcoTILLING Reveals Natural Allelic Variations in Starch Synthesis Key Gene TaSSIV and Its Haplotypes Associated with Higher Thousand Grain Weight

Wheat is a staple food commodity grown worldwide, and wheat starch is a valuable source of energy and carbon that constitutes 80% of the grain weight. Manipulation of genes involved in starch synthesis significantly affects wheat grain weight and yield. TaSSIV plays an important role in starch synthesis and its main function is granule formation. To mine and stack more favorable alleles, single nucleotide polymorphisms (SNPs) of TaSSIV-A, B, and D were investigated across 362 wheat accessions by Ecotype-Targeting Induced Local Lesions IN Genome (EcoTILLING). As a result, a total of 38 SNPs in the amplified regions of three TaSSIV genes were identified, of which 10, 15, and 13 were in TaSSIV-A, B, and D, respectively. These 38 SNPs were evaluated by using KASP and six SNPs showed an allele frequency >5% whereas the rest were <5%, i.e., considered to be minor alleles. In the Chinese mini core collection, three haplotypes were detected for TaSSIV–A and three for TaSSIV–B. The results of an association study in the Chinese mini core collection with thousand grain weight (TGW) and spike length (SPL) showed that Hap-2-1A was significantly associated with TGW and Hap-3-1B with SPL. Allelic frequency and geographic distribution indicated that the favored haplotype (Hap-2-1A) has been positively selected in Chinese wheat breeding. These results suggested that the Kompetitive Allele Specific PCR (KASP) markers can be applied in starch improvement to ultimately improve wheat yield by marker assisted selection in wheat breeding.

A wheat protein kinase gene TaSnRK2.9-5A associated with yield contributing traits

We developed breeder-friendly high-throughput and cost-effective KASP marker for marker-assisted selection for grain yield related traits in wheat. Plant-specific protein kinase, SnRK2s, is a major family of signaling genes associated with metabolic regulations, nutrient utilization and response to external stimuli. In the present study, three copies of TaSnRK2.9 were isolated from chromosomes 5A, 5B and 5D of wheat (Triticum aestivum L.). The coding regions of TaSnRK2.9-5A, TaSnRK2.9-5B and promoter region of TaSnRK2.9-5D were investigated for sequence polymorphism. Single nucleotide polymorphisms (SNPs) were identified for TaSnRK2.9-5A, while no polymorphism was identified in TaSnRK2.9-5B and TaSnRK2.9-5D. The nucleotide sequence of TaSnRK2.9-5A consisted of 2180 bp having eight introns and nine exons. Three SNPs were identified at 308 nt, 698 nt and 1700 nt. For high-throughput genotyping, two kompetitive allele-specific PCR (KASP) markers were developed. Four haplotypes Hap-5A-1, Hap-5A-2, Hap-5A-3 and Hap-5A-4 were detected in wheat populations collected from China, Europe and Pakistan. Association analysis was performed with mixed linear model in TASSEL (v 5.0). The results indicated that Hap-5A-1/2 of TaSnRK2.9-5A were significantly associated with high thousand kernel weight, while Hap-5A-4 with high grains per spike. Overexpressing transgenic rice also showed higher grains per spike which is in accordance with association analysis results. Geographic distribution and allelic frequency indicted that the favored haplotypes were positively selected in Chinese (Hap-5A-1/2), Pakistani (Hap-5A-1), east European (Hap-5A-1) and west European (Hap-5A-4) wheat breeding. The results suggest that the developed KASP markers can be utilized in yield improvement by marker-assisted selection in wheat breeding.

Transcriptomic de novo analysis of pitaya (Hylocereus polyrhizus) canker disease caused by Neoscytalidium dimidiatum

BACKGROUND

Canker disease caused by Neoscytalidium dimidiatum is the most serious disease that attacks the pitaya industry. One pathogenic fungus, referred to as ND8, was isolated from the wild-type red-fleshed pitaya (Hylocereus polyrhizus) of Hainan Province. In the early stages of this disease, stems show little spots and a loss of green color. These spots then gradually spread until the stems became rotten due to infection by various strains. Canker disease caused by Neoscytalidium dimidiatum poses a significant threat to pitaya commercial plantations with the growth of stems and the yields, quality of pitaya fruits. However, a lack of transcriptomic and genomic information hinders our understanding of the molecular mechanisms underlying the pitaya defense response.

RESULTS

We investigated the host responses of red-fleshed pitaya (H. polyrhizus) cultivars against N. dimidiatum using Illumina RNA-Seq technology. Significant expression profiles of 23 defense-related genes were further analyzed by qRT-PCR. The total read length based on RNA-Seq was 25,010,007; mean length was 744, the N50 was 1206, and the guanine-cytosine content was 44.48%. Our investigation evaluated 33,584 unigenes, of which 6209 (18.49%) and 27,375 (81.51%) were contigs and singlets, respectively. These unigenes shared a similarity of 16.62% with Vitis vinifera, 7.48% with Theobroma cacao, 6.6% with Nelumbo nucifera and 5.35% with Jatropha curcas. The assembled unigenes were annotated into non-redundant (NR, 25161 unigenes), Kyoto Encyclopedia of Genes and Genomes (KEGG, 17895 unigenes), Clusters of Orthologous Groups (COG, 10475 unigenes), InterPro (19,045 unigenes), and Swiss-Prot public protein databases (16,458 unigenes). In addition, 24 differentially expressed genes, which were mainly associated with plant pathology pathways, were analyzed in-depth.

CONCLUSIONS

This study provides a basis for further in-depth research on the protein function of the annotated unigene assembly with cDNA sequences.

Pleiotropic effects of the wheat domestication gene Q on yield and grain morphology

Transformation from q to Q during wheat domestication functioned outside the boundary of threshability to increase yield, grains m^-2, grain weight and roundness, but to reduce grains per spike/spikelet. Mutation of the Q gene, well-known affecting wheat spike structure, represents a key domestication step in the formation of today's free-threshing, economically important wheats. In a previous study, multiple yield components and spike characteristics were associated with the Q gene interval in the bread wheat 'Forno' × European spelt 'Oberkulmer' recombinant inbred line population. Here, we reported that this interval was also associated with grain yield, grains m^-2, grain morphology, and spike dry weight at anthesis. To clarify the roles of Q in agronomic trait performance, a functional marker for the Q gene was developed. Analysis of allelic effects showed that the bread wheat Q allele conferred free-threshing habit, soft glumes, and short and compact spikes compared with q. In addition, the Q allele contributed to higher grain yield, more grains m^-2, and higher thousand grain weight, whereas q contributed to more grains per spike/spikelet likely resulting from increased preanthesis spike growth. For grain morphology, the Q allele was associated with reduced ratio of grain length to height, indicating a rounder grain. These results are supported by analysis of four Q mutant lines in the Chinese Spring background. Therefore, the transition from q to Q during wheat domestication had profound effects on grain yield and grain shape evolution as well, being a consequence of pleiotropy.

The genome sequence of 'Kurodagosun', a major carrot variety in Japan and China, reveals insights into biological research and carrot breeding

Carrot (Daucus carota L.) is one of the most economically important root vegetables in the world, providing numerous nutrients for human health. China is the largest country of carrot production in the world, and 'Kurodagosun' has been a major carrot variety in China. Carrot material used in this study was the inbred line 'DC-27', which was derived by forced selfing from 'Kurodagosun'. To understand the genetic system and plant-specific genes of 'Kurodagosun', we report the draft genome sequence of carrot 'DC-27' assembled using a combination of Roche454 and HiSeq 2000 sequencing technologies to achieve 32-fold genome coverage. A total of 31,891 predicted genes were identified. These assembled sequences provide candidate genes involved in biological processes including stress response and carotenoid biosynthesis. Genomic sequences corresponding to 371.6 Mb was less than 473 Mb, which is the estimated genome size. The availability of a draft sequence of the 'DC-27' genome advances knowledge on the biological research and breeding of carrot, as well as other Apiaceae plants. The 'DC-27' genome sequence data also provide a new resource to explore the evolution of other higher plants.

Functional Analysis and Marker Development of TaCRT-D Gene in Common Wheat (Triticum aestivum L.)

Calreticulin (CRT), an endoplasmic reticulum (ER)-localized Ca2+-binding/buffering protein, is highly conserved and extensively expressed in animal and plant cells. To understand the function of CRTs in wheat (Triticum aestivum L.), particularly their roles in stress tolerance, we cloned the full-length genomic sequence of the TaCRT-D isoform from D genome of common hexaploid wheat, and characterized its function by transgenic Arabidopsis system. TaCRT-D exhibited different expression patterns in wheat seedling under different abiotic stresses. Transgenic Arabidopsis plants overexpressing ORF of TaCRT-D displayed more tolerance to drought, cold, salt, mannitol, and other abiotic stresses at both seed germination and seedling stages, compared with the wild-type controls. Furthermore, DNA polymorphism analysis and gene mapping were employed to develop the functional markers of this gene for marker-assistant selection in wheat breeding program. One SNP, S440 (T→C) was detected at the TaCRT-D locus by genotyping a wheat recombinant inbred line (RIL) population (114 lines) developed from Opata 85 × W7984. The TaCRT-D was then fine mapped between markers Xgwm645 and Xgwm664 on chromosome 3DL, corresponding to genetic distances of 3.5 and 4.4 cM, respectively, using the RIL population and Chinese Spring nulli-tetrasomic lines. Finally, the genome-specific and allele-specific markers were developed for the TaCRT-D gene. These findings indicate that TaCRT-D function importantly in plant stress responses, providing a gene target for genetic engineering to increase plant stress tolerance and the functional markers of TaCRT-D for marker-assistant selection in wheat breeding.

Genome-wide association study for 13 agronomic traits reveals distribution of superior alleles in bread wheat from the Yellow and Huai Valley of China

Bread wheat is a leading cereal crop worldwide. Limited amount of superior allele loci restricted the progress of molecular improvement in wheat breeding. Here, we revealed new allelic variation distribution for 13 yield-related traits in series of genome-wide association studies (GWAS) using the wheat 90K genotyping assay, characterized in 163 bread wheat cultivars. Agronomic traits were investigated in 14 environments at three locations over 3 years. After filtering SNP data sets, GWAS using 20 689 high-quality SNPs associated 1769 significant loci that explained, on average, ~20% of the phenotypic variation, both detected already reported loci and new promising genomic regions. Of these, repetitive and pleiotropic SNPs on chromosomes 6AS, 6AL, 6BS, 5BL and 7AS were significantly linked to thousand kernel weight, for example BS00021705_51 on 6BS and wsnp_Ex_c32624_41252144 on 6AS, with phenotypic variation explained (PVE) of ~24%, consistently identified in 12 and 13 of the 14 environments, respectively. Kernel length-related SNPs were mainly identified on chromosomes 7BS, 6AS, 5AL and 5BL. Plant height-related SNPs on chromosomes 4DS, 6DL, 2DS and 1BL were, respectively, identified in more than 11 environments, with averaged PVE of ~55%. Four SNPs were confirmed to be important genetic loci in two RIL populations. Based on repetivity and PVE, a total of 41 SNP loci possibly played the key role in modulating yield-related traits of the cultivars surveyed. Distribution of superior alleles at the 41 SNP loci indicated that superior alleles were getting popular with time and modern cultivars had integrated many superior alleles, especially for peduncle length- and plant height-related superior alleles. However, there were still 19 SNP loci showing less than percentages of 50% in modern cultivars, suggesting they should be paid more attention to improve yield-related traits of cultivars in the Yellow and Huai wheat region. This study could provide useful information for dissection of yield-related traits and valuable genetic loci for marker-assisted selection in Chinese wheat breeding programme.

Elite Haplotypes of a Protein Kinase Gene TaSnRK2.3 Associated with Important Agronomic Traits in Common Wheat

Plant-specific protein kinase SnRK2s play crucial roles in response to various environmental stimuli. TaSnRK2.3, a SnRK2 member, was involved in the response to multiple abiotic stresses in wheat. To facilitate the use of TaSnRK2.3 in wheat breeding, the three genomic sequences of TaSnRK2.3, originating from the A, B, and D genomes of hexaploid wheat, were obtained. Sequence polymorphism assays showing 4 and 10 variations were detected at TaSnRK2.3-1A and at TaSnRK2.3-1B, respectively, yet no variation was identified at TaSnRK2.3-1D. Three haplotypes for A genome, and two main haplotypes for B genome of TaSnRK2.3 were identified in 32 genotypes. Functional markers (2.3AM1, 2.3AM2, 2.3BM1, 2.3BM2) were successfully developed to distinguish different haplotypes. Association analysis was performed with the general linear model in TASSEL 2.1. The results showed that both TaSnRK2.3-1A and TaSnRK2.3-1B were significantly associated with plant height (PH), length of peduncle and penultimate node, as well as 1,000-grain weight (TGW) under different environments. Additionally, TaSnRK2.3-1B was significantly associated with stem water-soluble carbohydrates at flowering and mid-grain filling stages. Hap-1A-1 had higher TGW and lower PH; Hap-1B-1 had higher TGW and stem water-soluble carbohydrates, as well as lower PH, thus the two haplotypes were considered as elite haplotypes. Geographic distribution and allelic frequencies indicated that the two preferred haplotypes Hap-1A-1 and Hap-1B-1 were positively selected in the process of Chinese wheat breeding. These results could be valuable for genetic improvement and germplasm enhancement using molecular marker assisted selection in wheat breeding.