Identification of QTNs, QTN-by-environment interactions for plant height and ear height in maize multi-environment GWAS.

Plant height (PH) and ear height (EH) are important traits associated with biomass, lodging resistance, and grain yield in maize. There were strong effects of genotype x environment interaction (GEI) on plant height and ear height of maize. In this study, 203 maize inbred lines were grown at five locations across China's Spring and Summer corn belts, and plant height (PH) and ear height (EH) phenotype data were collected and grouped using GGE biplot. Five locations fell into two distinct groups (or mega environments) that coincide with two corn ecological zones called Summer Corn Belt and Spring Corn Belt. In total, 73,174 SNPs collected using GBS sequencing platform were used as genotype data and a recently released multi-environment GWAS software package IIIVmrMLM was employed to identify QTNs and QTN x environment (corn belt) interaction (QEIs); 12 and 11 statistically significant QEIs for PH and EH were detected respectively and their phenotypic effects were further partitioned into Add*E and Dom*E components. There were 28 and 25 corn-belt-specific QTNs for PH and EH identified, respectively. The result shows that there are a large number of genetic loci underlying the PH and EH GEIs and IIIVmrMLM is a powerful tool in discovering QTNs that have significant QTN-by-Environment interaction. PH and EH candidate genes were annotated based on transcriptomic analysis and haplotype analysis. EH related-QEI S10_135 (Zm00001d025947, saur76, small auxin up RNA76) and PH related-QEI S4_4 (Zm00001d049692, mads32, encoding MADS-transcription factor 32), and corn-belt specific QTNs including S10_4 (Zm00001d023333, sdg127, set domain gene127) and S7_1 (Zm00001d018614, GLR3.4, and glutamate receptor 3.4 or Zm00001d018616, DDRGK domain-containing protein) were reported, and the relationship among GEIs, QEIs and phenotypic plasticity and their biological and breeding implications were discussed.

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.

[Cloning and functional verification of PhAEP gene, a key enzyme for biosynthesis of heterophyllin A in Pseudostellaria heterophylla].

This paper aimed to study the role of asparagine endopeptidase(AEP) gene in the biosynthesis mechanism of cyclic peptide compounds in Pseudostellaria heterophylla. The transcriptome database of P. heterophylla was systematically mined and screened, and an AEP gene, tentatively named PhAEP, was successfully cloned. The heterologous function verification by Nicotiana benthamiana showed that the expression of the gene played a role in the biosynthesis of heterophyllin A in P. heterophylla. Bioinformatics analysis showed that the cDNA of PhAEP was 1 488 bp in length, encoding 495 amino acids with a molecular weight of 54.72 kDa. The phylogenetic tree showed that the amino acid sequence encoded by PhAEP was highly similar to that of Butelase-1 in Clitoria ternatea, reaching 80%. The sequence homology and cyclase active site analysis revealed that the PhAEP enzyme may specifically hydrolyse the C-terminal Asn/Asp(Asx) site of the core peptide in the HA linear precursor peptide of P. heterophylla, thereby participating in the ring formation of the linear precursor peptide. The results of real-time quantitative polymerase chain reaction(RT-qPCR) showed that the expression level of PhAEP was the highest in fruits, followed by in roots, and the lowest in leaves. The heterophyllin A of P. heterophylla was detected in N. benthamiana that co-expressed PrePhHA and PhAEP genes instantaneously. In this study, the PhAEP gene, a key enzyme in the biosynthesis of heterophyllin A in P. heterophylla, has been successfully cloned, which lays a foundation for further analysis of the molecular mechanism of PhAEP enzyme in the biosynthesis of heterophyllin A in P. heterophylla and has important significance for the study of synthetic biology of cyclic peptide compounds in P. heterophylla.

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.

AabHLH113 integrates jasmonic acid and abscisic acid signaling to positively regulate artemisinin biosynthesis in Artemisia annua.

Artemisinin, a sesquiterpene lactone isolated from Artemisia annua, is in huge market demand due to its efficient antimalarial action, especially after the COVID-19 pandemic. Many researchers have elucidated that phytohormones jasmonic acid (JA) and abscisic acid (ABA) positively regulate artemisinin biosynthesis via types of transcription factors (TFs). However, the crosstalk between JA and ABA in regulating artemisinin biosynthesis remains unclear. Here, we identified a novel ABA- and JA-induced bHLH TF, AabHLH113, which positively regulated artemisinin biosynthesis by directly binding to the promoters of artemisinin biosynthetic genes, DBR2 and ALDH1. The contents of artemisinin and dihydroartemisinic acid increased by 1.71- to 2.06-fold and 1.47- to 2.23-fold, respectively, in AabHLH1113 overexpressed A. annua, whereas they decreased by 14-36% and 26-53%, respectively, in RNAi-AabHLH113 plants. Furthermore, we demonstrated that AabZIP1 and AabHLH112, which, respectively, participate in ABA and JA signaling pathway to regulate artemisinin biosynthesis, directly bind to and activate the promoter of AabHLH113. Collectively, we revealed a complex network in which AabHLH113 plays a key interrelational role to integrate ABA- and JA-mediated regulation of artemisinin biosynthesis.

AabHLH112, a bHLH transcription factor, positively regulates sesquiterpenes biosynthesis in Artemisia annua.

The bHLH transcription factors play important roles in the regulation of plant growth, development, and secondary metabolism. ß-Caryophyllene, epi-cedrol, and ß-farnesene, three kinds of sesquiterpenes mainly found in plants, are widely used as spice in the food industry and biological pesticides in agricultural production. Furthermore, they also have a significant value in the pharmaceutical industry. However, there is currently a lack of knowledge on the function of bHLH family TFs in ß-caryophyllene, epi-cedrol, and ß-farnesene biosynthesis. Here, we found that AabHLH112 transcription factor had a novel function to positively regulate ß-carophyllene, epi-cedrol, and ß-farnesene biosynthesis in Artemisia annua. Exogenous MeJA enhanced the expression of AabHLH112 and genes of ß-caryophyllene synthase (CPS), epi-cedrol synthase (ECS), and ß-farnesene synthase (BFS), as well as sesquiterpenes content. Dual-LUC assay showed the activation of AaCPS, AaECS, and AaBFS promoters were enhanced by AabHLH112. Yeast one-hybrid assay showed AabHLH112 could bind to the G-box (CANNTG) cis-element in promoters of both AaCPS and AaECS. In addition, overexpression of AabHLH112 in A. annua significantly elevated the expression levels of AaCPS, AaECS, and AaBFS as well as the contents of ß-caryophyllene, epi-cedrol, and ß-farnesene, while suppressing AabHLH112 expression by RNAi reduced the expression of the three genes and the contents of the three sesquiterpenes. These results suggested that AabHLH112 is a positive regulator of ß-caryophyllene, epi-cedrol, and ß-farnesene biosynthesis in A. annua.

Molecular insights into AabZIP1-mediated regulation on artemisinin biosynthesis and drought tolerance in Artemisia annua.

Artemisia annua is the main natural source of artemisinin production. In A. annua, extended drought stress severely reduces its biomass and artemisinin production while short-term water-withholding or abscisic acid (ABA) treatment can increase artemisinin biosynthesis. ABA-responsive transcription factor AabZIP1 and JA signaling AaMYC2 have been shown in separate studies to promote artemisinin production by targeting several artemisinin biosynthesis genes. Here, we found AabZIP1 promote the expression of multiple artemisinin biosynthesis genes including AaDBR2 and AaALDH1, which AabZIP1 does not directly activate. Subsequently, it was found that AabZIP1 up-regulates AaMYC2 expression through direct binding to its promoter, and that AaMYC2 binds to the promoter of AaALDH1 to activate its transcription. In addition, AabZIP1 directly transactivates wax biosynthesis genes AaCER1 and AaCYP86A1. The biosynthesis of artemisinin and cuticular wax and the tolerance of drought stress were significantly increased by AabZIP1 overexpression, whereas they were significantly decreased in RNAi-AabZIP1 plants. Collectively, we have uncovered the AabZIP1-AaMYC2 transcriptional module as a point of cross-talk between ABA and JA signaling in artemisinin biosynthesis, which may have general implications. We have also identified AabZIP1 as a promising candidate gene for the development of A. annua plants with high artemisinin content and drought tolerance in metabolic engineering breeding.

The RppC-AvrRppC NLR-effector interaction mediates the resistance to southern corn rust in maize.

Southern corn rust (SCR), caused by the fungal pathogen Puccinia polysora, is a major threat to maize production worldwide. Efficient breeding and deployment of resistant hybrids are key to achieving durable control of SCR. Here, we report the molecular cloning and characterization of RppC, which encodes an NLR-type immune receptor and is responsible for a major SCR resistance quantitative trait locus. Furthermore, we identified the corresponding avirulence effector, AvrRppC, which is secreted by P. polysora and triggers RppC-mediated resistance. Allelic variation of AvrRppC directly determines the effectiveness of RppC-mediated resistance, indicating that monitoring of AvrRppC variants in the field can guide the rational deployment of RppC-containing hybrids in maize production. Currently, RppC is the most frequently deployed SCR resistance gene in China, and a better understanding of its mode of action is critical for extending its durability.

Genetic variation and population structure in China summer maize germplasm.

Maize (Zea mays L.) germplasm in China Summer maize ecological region (CSM) or central corn-belt of China is diverse but has not been systematically characterized at molecular level. In this study, genetic variation, genome diversity, linkage disequilibrium patterns, population structure, and characteristics of different heterotic groups were studied using 525,141 SNPs obtained by Genotyping-By-Sequencing (GBS) for 490 inbred lines collected from researchers at CSM region. The SNP density is lower near centromere, but higher near telomere region of maize chromosome, the degree of linkage disequilibrium (r2) vary at different chromosome regions. Majority of the inbred lines (66.05%) show pairwise relative kinship near zero, indicating a large genetic diversity in the CSM breeding germplasm. Using 4849 tagSNPs derived from 3618 haplotype blocks, the 490 inbred lines were delineated into 3 supergroups, 6 groups, and 10 subgroups using ADMIXTURE software. A procedure of assigning inbred lines into heterotic groups using genomic data and tag-SNPs was developed and validated. Genome differentiation among different subgroups measured by Fst, and the genetic diversity within each subgroup measured by GD are both large. The share of heterotic groups that have significant North American germplasm contribution: P, SS, IDT, and X, accounts about 54% of the CSM breeding germplasm collection and has increased significantly in the last two decades. Two predominant types of heterotic pattern in CSM region are: M-Reid group × TSPT group, and X subgroup × Local subgroups.

Effect of marker segregation distortion on high density linkage map construction and QTL mapping in Soybean (Glycine max L.).

Marker segregation distortion is a natural phenomenon. Severely distorted markers are usually excluded in the construction of linkage maps. We investigated the effect of marker segregation distortion on linkage map construction and quantitative trait locus (QTL) mapping. A total of 519 recombinant inbred lines of soybean from orthogonal and reciprocal crosses between LSZZH and NN493-1 were genotyped by specific length amplified fragment markers and seed linoleic acid content was measured in three environments. As a result, twenty linkage groups were constructed with 11,846 markers, including 1513 (12.77%) significantly distorted markers, on 20 chromosomes, and the map length was 2475.86 cM with an average marker-interval of 0.21 cM. The inclusion of distorted markers in the analysis was shown to not only improve the grouping of the markers from the same chromosomes, and the consistency of linkage maps with genome, but also increase genome coverage by markers. Combining genotypic data from both orthogonal and reciprocal crosses decreased the proportion of distorted markers and then improved the quality of linkage maps. Validation of the linkage maps was confirmed by the high collinearity between positions of markers in the soybean reference genome and in linkage maps and by the high consistency of 24 QTL regions in this study compared with the previously reported QTLs and lipid metabolism related genes. Additionally, linkage maps that include distorted markers could add more information to the outputs from QTL mapping. These results provide important information for linkage mapping, gene cloning and marker-assisted selection in soybean.

An assessment of transgenomics as a tool for identifying genes involved in the evolutionary differentiation of closely related plant species.

• Transgenomics is the process of introducing genomic clones from a donor species into a recipient species and then screening the resultant transgenic lines for phenotypes of interest. This method might allow us to find genes involved in the evolution of phenotypic differences between species as well as genes that have the potential to contribute to reproductive isolation: potential speciation genes. • More than 1100 20-kbp genomic clones from Leavenworthia alabamica were moved into Arabidopsis thaliana by transformation. After screening a single primary transformant for each line, clones associated with mutant phenotypes were tested for repeatability and co-segregation. • We found 84 clones with possible phenotypic effects, of which eight were repeatedly associated with the same phenotype. One clone, 11_11B, co-segregated with a short fruit phenotype. Further study showed that 11_11B affects seed development, with as much as one-third of the seeds aborted in some fruit. • Transgenomics is a viable strategy for discovering genes of evolutionary interest. We identify methods to reduce false positives and false negatives in the future. 11_11B can be viewed as a potential speciation gene, illustrating the value of transgenomics for studying the molecular basis of reproductive isolation.

Performance assessment of kernel density clustering for gene expression profile data.

Kernel density smoothing techniques have been used in classification or supervised learning of gene expression profile (GEP) data, but their applications to clustering or unsupervised learning of those data have not been explored and assessed. Here we report a kernel density clustering method for analysing GEP data and compare its performance with the three most widely-used clustering methods: hierarchical clustering, K-means clustering, and multivariate mixture model-based clustering. Using several methods to measure agreement, between-cluster isolation, and withincluster coherence, such as the Adjusted Rand Index, the Pseudo F test, the r(2) test, and the profile plot, we have assessed the effectiveness of kernel density clustering for recovering clusters, and its robustness against noise on clustering both simulated and real GEP data. Our results show that the kernel density clustering method has excellent performance in recovering clusters from simulated data and in grouping large real expression profile data sets into compact and well-isolated clusters, and that it is the most robust clustering method for analysing noisy expression profile data compared to the other three methods assessed.