Genetic variation and population structure in China summer maize germplasm

Using ddRADseq to assess the genetic diversity of in-farm and gene bank cacao resources in the Baracoa region, eastern Cuba, for use and conservation purposes

The Baracoa region, eastern Cuba, hosts around 80 % of the country cacao (Theobroma cacao L.) plantations. Cacao plants in farms are diverse in origin and propagation, with grafted and hybrid plants being the more common ones. Less frequent are plants from cuttings, TSH progeny, and traditional Cuban cacao. A national cacao gene bank is also present in Baracoa, with 282 accessions either prospected in Cuba or introduced from other countries. A breeding program associated with the gene bank started in the 1990s based on agro-morphological descriptors. The genetic diversity of cacao resources in Baracoa has been poorly described, except for traditional Cuban cacao, affecting the proper development of the breeding program and the cacao planting policies in the region. To assess the population structure and genetic diversity of cacao resources in Baracoa region, we genotyped plants from both cacao gene bank (CG) and cacao farms (CF) applying a new ddRADseq protocol for cacao. After data processing, two SNPs datasets containing 11,425 and 6,481 high-quality SNPs were generated with 238 CG and 135 CF plants, respectively. SNPs were unevenly distributed along the 10 cacao chromosomes and laid mainly in noncoding regions of the genome. Population structure analysis with these SNP datasets identified seven and four genetic groups in CG and CF samples, respectively. Clustering using UPGMA and principal component analysis mostly agree with population structure results. Amelonado was the predominant cacao ancestry, accounting for 49.22 % (CG) and 57.73 % (CF) of the total. Criollo, Contamana, Iquitos, and Nanay ancestries were detected in both CG and CF samples, while Nacional and Marañon backgrounds were only identified in CG. Genetic differentiation among CG (FST ranging from 0.071 to 0.407) was higher than among CF genetic groups (FST: 0.093-0.282). Genetic diversity parameters showed similar values for CG and CF samples. The CG and CF genetic groups with the lowest genetic diversity parameters had the highest proportion of Amelonado ancestry. These results should contribute to reinforcing the ongoing breeding program and updating the planting policies on cacao farms, with an impact on the social and economic life of the region.

Genome-wide association study, population structure, and genetic diversity of the tea plant in Guizhou Plateau

BACKGROUND

Guizhou Plateau, as one of the original centers of tea plant, has a profound multi-ethnic cultural heritage and abundant tea germplasm resources. However, the impact of indigenous community factors on the genetic diversity, population structure and geographical distribution of tea plant is still unclear.

RESULTS

Using the genotyping-by-sequencing (GBS) approach, we collected 415 tea plant accessions from the study sites, estimated genetic diversity, developed a core collection, and conducted a genome-wide association study (GWAS) based on 99,363 high-quality single-nucleotide polymorphisms (SNPs). A total of 415 tea accessions were clustered into six populations (GP01, GP02, GP03, GP04, GP05 and GP06), and the results showed that GP04 and GP05 had the highest and lowest genetic diversity (Pi = 0.214 and Pi = 0.145, respectively). Moreover, 136 tea accessions (33%) were selected to construct the core set that can represent the genetic diversity of the whole collection. By analyzing seven significant SNP markers associated with the traits such as the germination period of one bud and two leaves (OTL) and the germination period of one bud and three leaves (OtL), four candidate genes possibly related to OTL and OtL were identified.

CONCLUSIONS

This study revealed the impact of indigenous communities on the population structure of 415 tea accessions, indicating the importance of cultural practices for protection and utilization of tea plant genetic resources. Four potential candidate genes associated with the OTL and OtL of tea plant were also identified, which will facilitate genetic research, germplasm conservation, and breeding.

Identification of southern corn rust resistance QTNs in Chinese summer maize germplasm via multi-locus GWAS and post-GWAS analysis

Southern corn rust (SCR) caused by Puccinia polysora Underw is a major disease leading to severe yield losses in China Summer Corn Belt. Using six multi-locus GWAS methods, we identified a set of SCR resistance QTNs from a diversity panel of 140 inbred lines collected from China Summer Corn Belt. Thirteen QTNs on chromosomes 1, 2, 4, 5, 6, and 8 were grouped into three types of allele effects and their associations with SCR phenotypes were verified by post-GWAS case-control sampling, allele/haplotype effect analysis. Relative resistance (RRR) and relative susceptibility (RRs) catering to its inbred carrier were estimated from single QTN and QTN-QTN combos and epistatitic effects were estimated for QTN-QTN combos. By transcriptomic annotation, a set of candidate genes were predicted to be involved in transcriptional regulation (S5_145, Zm00001d01613, transcription factor GTE4), phosphorylation (S8_123, Zm00001d010672, Pgk2- phosphoglycerate kinase 2), and temperature stress response (S6_164a/S6_164b, Zm00001d038806, hsp101, and S5_211, Zm00001d017978, cellulase25). The breeding implications of the above findings were discussed.

Nitrogen reduction combined with ETc irrigation maintained summer maize yield and increased water and nitrogen use efficiency

Introduction

High rainfall and excessive urea application are counterproductive to summer maize growth requirements and lower grain yield and water/nitrogen (N) use efficiency. The objective of this study was to determine whether ETc irrigation based on summer maize demand and reduced nitrogen rate in the Huang Huai Hai Plain increased water and nitrogen use efficiency without sacrificing yield.

Methods

To achieve this, we conducted an experiment with four irrigation levels [ambient rainfall (I0) and 50% (I1), 75% (I2), and 100% (I3) of actual crop evapotranspiration (ET_c)] and four nitrogen rates [no nitrogen fertilizer (N0), recommended nitrogen rate of urea (NU), recommended nitrogen rate of blending controlled-release urea with conventional urea fertilizer (BCRF) (NC), and reduced nitrogen rate of BCRF (NR)] in 2016-2018.

Results

The results show that reduced irrigation and nitrogen rate reduced Fv/Fm, ¹³C-photosynthate, and nitrogen accumulation both in the kernel and plant. I3NC and I3NU accumulated higher ¹³C-photosynthate, nitrogen, and dry matter. However, ¹³C-photosynthate and nitrogen distribution to the kernel was decreased from I2 to I3 and was higher in BCRF than in urea. I2NC and I2NR promoted their distribution to the kernel, resulting in a higher harvest index. Compared with I3NU, I2NR increased root length density by 32.8% on average, maintaining considerable leaf Fv/Fm and obtaining similar kernel number and kernel weight. The higher root length density of I2NR of 40-60 cm promoted ¹³C-photosynthate and nitrogen distribution to the kernel and increased the harvest index. As a result, the water use efficiency (WUE) and nitrogen agronomic use efficiency (NAUE) in I2NR increased by 20.5%-31.9% and 11.0%-38.0% than that in I3NU, respectively.

Discussion

Therefore, 75%ET_c deficit irrigation and BCRF fertilizer with 80% nitrogen rate improved root length density, maintained leaf Fv/Fm in the milking stage, promoted 13C-photosynthate, and distributed nitrogen to the kernel, ultimately providing a higher WUE and NAUE without significantly reducing grain yield.

Population structure analysis to explore genetic diversity and geographical distribution characteristics of wild tea plant in Guizhou Plateau

BACKGROUND

Tea, the second largest consumer beverage in the world after water, is widely cultivated in tropical and subtropical areas. However, the effect of environmental factors on the distribution of wild tea plants is unclear.

RESULTS

A total of 159 wild tea plants were collected from different altitudes and geological types of the Guizhou Plateau. Using the genotyping-by-sequencing method, 98,241 high-quality single nucleotide polymorphisms were identified. Genetic diversity, population structure analysis, principal component analysis, phylogenetic analysis, and linkage disequilibrium were performed. The genetic diversity of the wild tea plant population from the Silicate Rock Classes of Camellia gymnogyna was higher than that from the Carbonate Rock Classes of Camellia tachangensis. In addition, the genetic diversity of wild tea plants from the second altitude gradient was significantly higher than that of wild tea plants from the third and first altitude gradients. Two inferred pure groups (GP01 and GP02) and one inferred admixture group (GP03) were identified by population structure analysis and were verified by principal component and phylogenetic analyses. The highest differentiation coefficients were determined for GP01 vs. GP02, while the lowest differentiation coefficients were determined for GP01 vs. GP03.

CONCLUSIONS

This study revealed the genetic diversity and geographical distribution characteristics of wild tea plants in the Guizhou Plateau. There are significant differences in genetic diversity and evolutionary direction between Camellia tachangensis with Carbonate Rock Classes at the first altitude gradient and Camellia gymnogyna with Silicate Rock Classes at the third altitude gradient. Geological environment, soil mineral element content, soil pH, and altitude markedly contributed to the genetic differentiation between Camellia tachangensis and Camellia gymnogyna.

Modeling the influence of phenotypic plasticity on maize hybrid performance

Phenotypic plasticity, the ability of an individual to alter its phenotype in response to changes in the environment, has been proposed as a target for breeding crop varieties with high environmental fitness. Here, we used phenotypic and genotypic data from multiple maize (Zea mays L.) populations to mathematically model phenotypic plasticity in response to the environment (PPRE) in inbred and hybrid lines. PPRE can be simply described by a linear model in which the two main parameters, intercept a and slope b, reflect two classes of genes responsive to endogenous (class A) and exogenous (class B) signals that coordinate plant development. Together, class A and class B genes contribute to the phenotypic plasticity of an individual in response to the environment. We also made connections between phenotypic plasticity and hybrid performance or general combining ability (GCA) of yield using 30 F₁ hybrid populations generated by crossing the same maternal line with 30 paternal lines from different maize heterotic groups. We show that the parameters a and b from two given parental lines must be concordant to reach an ideal GCA of F₁ yield. We hypothesize that coordinated regulation of the two classes of genes in the F₁ hybrid genome is the basis for high GCA. Based on this theory, we built a series of predictive models to evaluate GCA in silico between parental lines of different heterotic groups.

Genetic analysis of global faba bean diversity, agronomic traits and selection signatures

KEY MESSAGE

We identified marker-trait associations for key faba bean agronomic traits and genomic signatures of selection within a global germplasm collection. Faba bean (Vicia faba L.) is a high-protein grain legume crop with great potential for sustainable protein production. However, little is known about the genetics underlying trait diversity. In this study, we used 21,345 high-quality SNP markers to genetically characterize 2678 faba bean genotypes. We performed genome-wide association studies of key agronomic traits using a seven-parent-MAGIC population and detected 238 significant marker-trait associations linked to 12 traits of agronomic importance. Sixty-five of these were stable across multiple environments. Using a non-redundant diversity panel of 685 accessions from 52 countries, we identified three subpopulations differentiated by geographical origin and 33 genomic regions subjected to strong diversifying selection between subpopulations. We found that SNP markers associated with the differentiation of northern and southern accessions explained a significant proportion of agronomic trait variance in the seven-parent-MAGIC population, suggesting that some of these traits were targets of selection during breeding. Our findings point to genomic regions associated with important agronomic traits and selection, facilitating faba bean genomics-based breeding.

Pattern and variation in simple sequence repeat (SSR) at different genomic regions and its implications to maize evolution and breeding

BACKGROUND

Repetitive DNA sequences accounts for over 80% of maize genome. Although simple sequence repeats (SSRs) account for only 0.03% of the genome, they have been widely used in maize genetic research and breeding as highly informative codominant DNA markers. The genome-wide distribution and polymorphism of SSRs are not well studied due to the lack of high-quality genome DNA sequence data.

RESULTS

In this study, using data from high-quality de novo-sequenced maize genomes of five representative maize inbred lines, we revealed that SSRs were more densely present in telomeric region than centromeric region, and were more abundant in genic sequences than intergenic sequences. On genic sequences, tri- and hexanucleotide motifs were more abundant in CDS sequence and some mono- and dinucleotide motifs were more abundant in UTR sequences. Median length and chromosomal density of SSRs were both narrowly range-bound, with median length of 14-18 bp and genome-wide average density of 3355.77 bp/Mbp. LTR-RTs of < 0.4 Mya had higher SSR density (4498-4992 bp/Mbp). The genome-specific and motif-specific SSR polymorphism were studied. Their potential breeding applications were discussed.

CONCLUSIONS

We found that the median length of SSR sequences of different SSR motifs was nearly constant. SSR density in genic regions was much higher than intergenic regions. In addition, SSR density at LTR-RTs of different evolutionary ages varied in a narrow range. The SSRs and their LTR-RT carriers evolved at an equal rate. All these observations indicated that SSR length and density were under control of yet unknown evolutionary forces. The chromosome region-specific and motif-specific SSR polymorphisms we observed supported the notion that SSR polymorphism was invaluable genome resource for developing highly informative genome and gene markers in maize genetic research and molecular breeding.

Effects of post-silking low temperature on the starch and protein metabolism, endogenous hormone contents, and quality of grains in waxy maize

Waxy maize has many excellent characteristics in food and nonfood industries. However, post-silking low temperature (LT) has severe limitations on its grain yield and quality. In this study, field and pot trials were conducted to investigate the effects of post-silking LT on the physiological, biochemical, and functional characteristics of two waxy maize grains. The field and pot trials were performed with sowing date and artificial climate chamber, respectively, for LT treatment from silking stage to maturity. Results in pot trial were used to explain and validate the findings in field trial. Compared with the ambient treatment, the LT treatment significantly reduced kernel weight during the grain filling stage (P < 0.05). LT treatment in both environments resulted in an average decrease in dry weight of SYN5 and YN7 at maturity by 36.6% and 42.8%, respectively. Enzymatic activities related to starch and protein biosynthesis decreased under the LT treatment during the filling stage, accompanied by a decrease in the accumulation amounts and contents of soluble sugar and starch, and a decrease in protein accumulation amount. Meanwhile, the contents of abscisic acid, indole-3-acetic acid, and gibberellin 3 in grains decreased under the LT treatment during the filling stage. Peak, trough, breakdown, final, and setback viscosities of grains decreased by LT. LT treatment decreased the gelatinization enthalpy of grains and increased the retrogradation percentage. In conclusion, post-silking LT stress altered the content of grain components by inhibiting the production of phytohormones and down-regulating the enzymatic activities involved in starch and protein metabolism, which resulted in the deterioration of grain pasting and thermal properties.

Transcriptomic and alternative splicing analyses provide insights into the roles of exogenous salicylic acid ameliorating waxy maize seedling growth under heat stress

BACKGROUND

Salicylic acid (SA) is a phytohormone which works to regulate the abiotic stress response of plants. However, the molecular mechanism by which SA mediates heat tolerance in waxy maize (Zea mays L. sinsensis Kulesh) remains unknown.

RESULTS

Two varieties of waxy maize seedlings, heat-tolerant 'Yunuo7' (Y7) and heat-sensitive 'Suyunuo5' (S5), were pretreated with SA prior to heat stress (HTS). After treatment, physiological and transcriptomic changes were analyzed. Compared with HTS, the exogenous application of SA enhanced the shoot dry weight, the activities of antioxidant enzymes (e.g., SOD, POD, CAT and APX), and the concentration of endogenous phytohormones (e.g., SA, ABA, IAA, GA3), while decreased the MDA content. Transcriptome analysis showed that the number of differentially expressed genes (DEGs) identified in the control (CK) vs HTS and HTS vs HTS + SA comparisons were more in S5 than in Y7. HTS induced the downregulation of genes involved in photosynthesis and the upregulation of genes encoding heat shock transcription factors (HSFs) and heat shock proteins (HSPs). Compared with HTS, SA pretreatment reversed the expression of 5 photosynthesis-related genes, 26 phytohormone-related genes, and all genes encoding HSFs and HSPs in S5. Furthermore, the number of alternative splicing (AS) events increased under HTS treatment for both varieties, while decreased under SA pretreatment of S5. Differentially spliced genes (DSGs) showed little overlap with DEGs, and DEGs and DSGs differed significantly in functional enrichment.

CONCLUSIONS

Physiological and transcriptional together indicated that HTS and SA pretreatment had a greater effect on S5 than Y7. Additionally, it appears that transcriptional regulation and AS work synergistically to enhance thermotolerance in heat-sensitive waxy maize. Our study revealed the regulatory effects and underlying molecular mechanisms of SA on waxy maize seedling under HTS.

Population Structure Analysis and Association Mapping for Turcicum Leaf Blight Resistance in Tropical Maize Using SSR Markers

Maize is an important cereal crop in the world for feed, food, fodder, and raw materials of industries. Turcicum leaf blight (TLB) is a major foliar disease that can cause more than 50% yield losses in maize. Considering this, the molecular diversity, population structure, and genome-wide association study (GWAS) for TLB resistance were studied in 288 diverse inbred lines genotyped using 89 polymorphic simple sequence repeats (SSR) markers. These lines werescreened for TLB disease at two hot-spot locations under artificially inoculated conditions. The average percent disease incidence (PDI) calculated for each genotype ranged from 17 (UMI 1201) to 78% (IML 12-22) with an overall mean of 40%. The numbers of alleles detected at a locus ranged from twoto nine, with a total of 388 alleles. The polymorphic information content (PIC) of each marker ranged between 0.04 and 0.86. Out of 89 markers, 47 markers were highly polymorphic (PIC ≥ 0.60). This indicated that the SSR markers used were very informative and suitable for genetic diversity, population structure, and marker-trait association studies.The overall observed homozygosity for highly polymorphic markers was 0.98, which indicated that lines used were genetically pure. Neighbor-joining clustering, factorial analysis, and population structure studies clustered the 288 lines into 3–5 groups. The patterns of grouping were in agreement with the origin and pedigree records of the genotypesto a greater extent.A total of 94.10% lines were successfully assigned to one or another group at a membership probability of ≥0.60. An analysis of molecular variance (AMOVA) revealed highly significant differences among populations and within individuals. Linkage disequilibrium for r² and D′ between loci ranged from 0 to 0.77 and 0 to 1, respectively. A marker trait association analysis carried out using a general linear model (GLM) and mixed linear model (MLM), identified 15 SSRs markers significantly associated with TLB resistance.These 15 markers were located on almost all chromosomes (Chr) except 7, 8, and 9. The phenotypic variation explained by these loci ranged from 6% (umc1367) to 26% (nc130, phi085). Maximum 7 associated markers were located together on Chr 2 and 5. The selected regions identified on Chr 2 and 5 corroborated the previous studies carried out in the Indian maize germplasm. Further, 11 candidate genes were identified to be associated with significant markers. The identified sources for TLB resistance and associated markers may be utilized in molecular breeding for the development of suitable genotypes.

Genetic Diversity, Population Structure and Linkage Disequilibrium Analyses in Tropical Maize Using Genotyping by Sequencing

Several maize breeding programs in India have developed numerous inbred lines but the lines have not been characterized using high-density molecular markers. Here, we studied the molecular diversity, population structure, and linkage disequilibrium (LD) patterns in a panel of 314 tropical normal corn, two sweet corn, and six popcorn inbred lines developed by 17 research centers in India, and 62 normal corn from the International Maize and Wheat Improvement Center (CIMMYT). The 384 inbred lines were genotyped with 60,227 polymorphic single nucleotide polymorphisms (SNPs). Most of the pair-wise relative kinship coefficients (58.5%) were equal or close to 0, which suggests the lack of redundancy in the genomic composition in the majority of inbred lines. Genetic distance among most pairs of lines (98.3%) varied from 0.20 to 0.34 as compared with just 1.7% of the pairs of lines that differed by <0.20, which suggests greater genetic variation even among sister lines. The overall average of 17% heterogeneity was observed in the panel indicated the need for further inbreeding in the high heterogeneous genotypes. The mean nucleotide diversity and frequency of polymorphic sites observed in the panel were 0.28 and 0.02, respectively. The model-based population structure, principal component analysis, and phylogenetic analysis revealed three to six groups with no clear patterns of clustering by centers-wise breeding lines, types of corn, kernel characteristics, maturity, plant height, and ear placement. However, genotypes were grouped partially based on their source germplasm from where they derived.

Population structure analysis to explore genetic diversity and geographical distribution characteristics of cultivated-type tea plant in Guizhou Plateau

BACKGROUND

Tea plants originated in southwestern China. Guizhou Plateau is an original center of tea plants, and is rich in germplasm resources. However, the genetic diversity, population structure and distribution characteristics of cultivated-type tea plants in the region are unknown. In this study, we explored the genetic diversity and geographical distribution of cultivated-type tea accessions in Guizhou Plateau.

RESULTS

We used 112,072 high-quality genotyping-by-sequencing to analyze the genetic diversity, principal components, phylogeny, population structure, and linkage disequilibrium, and develop a core collection of 253 cultivated-type tea plant accessions from Guizhou Plateau. The results showed Genetic diversity of the cultivated-type tea accessions of the Pearl River Basin was significantly higher than that of the cultivated-type tea accessions of the Yangtze River Basin. Three inferred pure groups (CG-1, CG-2 and CG-3) and one inferred admixture group (CG-4), were identified by a population structure analysis, and verified by principal component and phylogenetic analyses. The highest genetic distance and differentiation coefficients were determined for CG-2 vs CG-3. The lower genetic distance and differentiation coefficients were determined for CG-4 vs CG-2 and CG-4 vs CG-3, respectively. We developed a core set and a primary set. The primary and core sets contained 77.0 and 33.6% of all individuals in the initial set, respectively. The primary set may serve as the primary population in genome-wide association studies, while the core collection may serve as the core population in multiple treatment setting studies.

CONCLUSIONS

The present study demonstrated the genetic diversity and geographical distribution characteristics of cultivated-type tea plants in Guizhou Plateau. Significant differences in genetic diversity and evolutionary direction were detected between the ancient landraces of the Pearl River Basin and the those of the Yangtze River Basin. Major rivers and ancient hubs were largely responsible for the genetic exchange between the Pearl River Basin and the Yangtze River Basin ancient landraces as well as the formation of the ancient hubs evolutionary group. Genetic diversity, population structure and core collection elucidated by this study will facilitate further genetic studies, germplasm protection, and breeding of tea plants.