Liying Yan

Diversity characterization and association analysis of agronomic traits in a Chinese peanut (Arachis hypogaea L.) mini-core collection

Association mapping is a powerful approach for exploring the molecular basis of phenotypic variations in plants. A peanut (Arachis hypogaea L.) mini-core collection in China comprising 298 accessions was genotyped using 109 simple sequence repeat (SSR) markers, which identified 554 SSR alleles and phenotyped for 15 agronomic traits in three different environments, exhibiting abundant genetic and phenotypic diversity within the panel. A model-based structure analysis assigned all accessions to three groups. Most of the accessions had the relative kinship of less than 0.05, indicating that there were no or weak relationships between accessions of the mini-core collection. For 15 agronomic traits in the peanut panel, generally the Q + K model exhibited the best performance to eliminate the false associated positives compared to the Q model and the general linear model-simple model. In total, 89 SSR alleles were identified to be associated with 15 agronomic traits of three environments by the Q + K model-based association analysis. Of these, eight alleles were repeatedly detected in two or three environments, and 15 alleles were commonly detected to be associated with multiple agronomic traits. Simple sequence repeat allelic effects confirmed significant differences between different genotypes of these repeatedly detected markers. Our results demonstrate the great potential of integrating the association analysis and marker-assisted breeding by utilizing the peanut mini-core collection.

Genetic Diversity and Population Structure of the Major Peanut (Arachis hypogaea L.) Cultivars Grown in China by SSR Markers

One hundred and forty-six highly polymorphic simple sequence repeat (SSR) markers were used to assess the genetic diversity and population structure of 196 peanut (Arachis Hypogaea L.) cultivars which had been extensively planted in different regions in China. These SSR markers amplified 440 polymorphic bands with an average of 2.99, and the average gene diversity index was 0.11. Eighty-six rare alleles with a frequency of less than 1% were identified in these cultivars. The largest Fst or genetic distance was found between the cultivars that adapted to the south regions and those to the north regions in China. A neighbor-joining tree of cultivars adapted to different ecological regions was constructed based on pairwise Nei’s genetic distances, which showed a significant difference between cultivars from the south and the north regions. A model-based population structure analysis divided these peanut cultivars into five subpopulations (P1a, P1b, P2, P3a and P3b). P1a and P1b included most the cultivars from the southern provinces including Guangdong, Guangxi and Fujian. P2 population consisted of the cultivars from Hubei province and parts from Shandong and Henan. P3a and P3b had cultivars from the northern provinces including Shandong, Anhui, Henan, Hebei, Jiangsu and the Yangtze River region including Sichuan province. The cluster analysis, PCoA and PCA based on the marker genotypes, revealed five distinct clusters for the entire population that were related to their germplasm regions. The results indicated that there were obvious genetic variations between cultivars from the south and the north, and there were distinct genetic differentiation among individual cultivars from the south and the north. Taken together, these results provided a molecular basis for understanding genetic diversity of Chinese peanut cultivars.

Deep sequencing analysis of transcriptomes in Aspergillus flavus in response to resveratrol

BackgroundResveratrol has been reported as a natural phytoalexin that inhibits infection or the growth of certain fungi including Aspergillus flavus. Our previous research revealed that aflatoxin production in A. flavus was reduced in medium with resveratrol. To understand the molecular mechanism of the A. flavus response to resveratrol treatment, the high-throughput paired-end RNA-Seq was applied to analyze the transcriptomic profiles of A. flavus.ResultsIn total, 366 and 87 genes of A. flavus were significantly up- and down- regulated, respectively, when the fungus was treated with resveratrol. Gene Ontology (GO) functional enrichment analysis revealed that 48 significantly differentially expressed genes were involved in 6 different terms. Most genes in the aflatoxin biosynthetic pathway genes cluster (#54) did not show a significant change when A. flavus was treated with resveratrol, but 23 of the 30 genes in the #54 cluster were down-regulated. The transcription of aflA and aflB was significantly suppressed under resveratrol treatment, resulting in an insufficient amount of the starter unit hexanoate for aflatoxin biosynthesis. In addition, resveratrol significantly increased the activity of antioxidative enzymes that destroy radicals, leading to decreased aflatoxin production. Moreover, stuA, fluG, flbC, and others genes involved in mycelial and conidial development were down-regulated, which disrupted the cell’s orderly differentiation and blocked conidia formation and mycelia development. The transcripts of laeA and veA were slightly inhibited by resveratrol, which may partly decrease aflatoxin production and depress conidia formation.ConclusionsResveratrol can affect the expression of A. flavus genes that are related to developmental and secondary metabolic processes, resulting in decreased aflatoxin production and conidia formation and could also cause abnormal mycelia development. These results provide insight into the transcriptome of A. flavus in response to resveratrol and a new clew for further study in regulation of aflatoxin biosynthesis in A. flavus.

Genome-wide analysis of the basic leucine zipper (bZIP) transcription factor gene family in six legume genomes.

BackgroundPlant bZIP proteins characteristically harbor a highly conserved bZIP domain with two structural features: a DNA-binding basic region and a leucine (Leu) zipper dimerization region. They have been shown to be diverse transcriptional regulators, playing crucial roles in plant development, physiological processes, and biotic/abiotic stress responses. Despite the availability of six completely sequenced legume genomes, a comprehensive investigation of bZIP family members in legumes has yet to be presented.ResultsIn this study, we identified 428 bZIP genes encoding 585 distinct proteins in six legumes, Glycine max, Medicago truncatula, Phaseolus vulgaris, Cicer arietinum, Cajanus cajan, and Lotus japonicus. The legume bZIP genes were categorized into 11 groups according to their phylogenetic relationships with genes from Arabidopsis. Four kinds of intron patterns (a–d) within the basic and hinge regions were defined and additional conserved motifs were identified, both presenting high group specificity and supporting the group classification. We predicted the DNA-binding patterns and the dimerization properties, based on the characteristic features in the basic and hinge regions and the Leu zipper, respectively, which indicated that some highly conserved amino acid residues existed across each major group. The chromosome distribution and analysis for WGD-derived duplicated blocks revealed that the legume bZIP genes have expanded mainly by segmental duplication rather than tandem duplication. Expression data further revealed that the legume bZIP genes were expressed constitutively or in an organ-specific, development-dependent manner playing roles in multiple seed developmental stages and tissues. We also detected several key legume bZIP genes involved in drought- and salt-responses by comparing fold changes of expression values in drought-stressed or salt-stressed roots and leaves.ConclusionsIn summary, this genome-wide identification, characterization and expression analysis of legume bZIP genes provides valuable information for understanding the molecular functions and evolution of the legume bZIP transcription factor family, and highlights potential legume bZIP genes involved in regulating tissue development and abiotic stress responses.

Comparison of Genetic Diversity Based on SSR Markers Between Peanut Mini Core Collections from China and ICRISAT

Mini core collection plays an important role in evaluating genetic resources of peanut (Arachis hypogaea L.). This study aimed at comparing the genetic diversities of a domestic and exotic peanut mini core collections and providing basic data for germplasm introduction and peanut breeding. The exotic mini core collection, which was introduced from the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT, Patancheru, India), consisted of 168 accessions, whereas the Chinese mini core collection consisted of 298 accessions. A total of 26 polymorphic simple sequence repeat (SSR) markers were used to evaluate the genetic diversity of these genotypes and their similarities. The similarity coefficient ranged from 0.49 to 0.99 among the total 466 accessions and the largest genetic distance was between “L2 Gangguo” from the Chinese core collection and “ICG12625” from the ICRISAT core collection. Among the 6 botanical types of peanut, the fastigiata type had the largest genetic diversity indices in both Chinese and ICRISAT mini core collections, which were 1.11 and 0.97, respectively. The hypogaea type had the second largest diversity with genetic diversity indices of 0.89 and 0.88 for Chinese and ICRISAT mini core collections, respectively. The genetic distance between the 2 mini core collections was larger than that within a mini core collection. Particularly, the aequatoriana genotype ICG12625 from ICRISAT had the largest differences to the Chinese genotypes. According to similarity coefficient and genetic diversity index, the genetic diversity was greater among genotypes from the Chinese mini core collection than that from the ICRISAT mini core collection

Identification of ERF genes in peanuts and functional analysis of AhERF008 and AhERF019 in abiotic stress response

Ethylene-responsive factor (ERF) play an important role in regulating gene expression in plant development and response to stresses. In peanuts (Arachis hypogaea L.), which produce flowers aerially and pods underground, only a few ERF genes have been identified so far. This study identifies 63 ERF unigenes from 247,313 peanut EST sequences available in the NCBI database. The phylogeny, gene structures, and putative conserved motifs in the peanut ERF proteins were analysed. Comparative analysis revealed the absence of two subgroups (A1 and A3) of the ERF family in peanuts; only 10 subgroups were identified in peanuts compared to 12 subgroups in Arabidopsis and soybeans. AP2/ERF domains were found to be conserved among peanuts, Arabidopsis, and soybeans. Outside the AP2/ERF domain, many soybean-specific conserved motifs were also detected in peanuts. The expression analysis of ERF family genes representing each clade revealed differential expression patterns in response to biotic and abiotic stresses. Overexpression of AhERF008 influenced the root gravity of Arabidopsis, whereas overexpression of AhERF019 enhanced tolerance to drought, heat, and salt stresses in Arabidopsis. The information generated in this study will be helpful to further investigate the function of ERFs in plant development and stress response.

Quantitative Trait Locus Analysis of Late Leaf Spot Resistance and Plant-Type-Related Traits in Cultivated Peanut (Arachis hypogaea L.) under Multi-Environments.

Late leaf spot (LLS) is one of the most serious foliar diseases affecting peanut worldwide leading to huge yield loss. To understand the genetic basis of LLS and assist breeding in the future, we conducted quantitative trait locus (QTL) analysis for LLS and three plant-type-related traits including height of main stem (HMS), length of the longest branch (LLB) and total number of branches (TNB). Significant negative correlations were observed between LLS and the plant-type-related traits in multi-environments of a RIL population from the cross Zhonghua 5 and ICGV 86699. A total of 20 QTLs were identified for LLS, of which two QTLs were identified in multi-environments and six QTLs with phenotypic variation explained (PVE) more than 10%. Ten, seven, fifteen QTLs were identified for HMS, LLB and TNB, respectively. Of these, one, one, two consensus QTLs and three, two, three major QTLs were detected for HMS, LLB and TNB, respectively. Of all 52 unconditional QTLs for LLS and plant-type-related traits, 10 QTLs were clustered in five genetic regions, of which three clusters including five robust major QTLs overlapped between LLS and one of the plant-type-related traits, providing evidence that the correlation could be genetically constrained. On the other hand, conditional mapping revealed different numbers and different extent of additive effects of QTLs for LLS conditioned on three plant-type-related traits (HMS, LLB and TNB), which improved our understanding of interrelationship between LLS and plant-type-related traits at the QTL level. Furthermore, two QTLs, qLLSB6-7 and qLLSB1 for LLS resistance, were identified residing in two clusters of NB-LRR—encoding genes. This study not only provided new favorable QTLs for fine-mapping, but also suggested that the relationship between LLS and plant-type-related traits of HMS, LLB and TNB should be considered while breeding for improved LLS resistance in peanut.

Dynamics in the resistant and susceptible peanut (Arachis hypogaea L.) root transcriptome on infection with the Ralstonia solanacearum.

BackgroundBacterial wilt caused by Ralstonia solanacearum is a serious soil-borne disease of peanut (Arachis hypogaea L). The molecular basis of peanut response to R. solanacearum remains unknown. To understand the resistance mechanism behind peanut resistance to R. solanacearum, we used RNA-Seq to perform global transcriptome profiling on the roots of peanut resistant (R) and susceptible (S) genotypes under R. solanacearum infection.ResultsA total of 4.95 x 108 raw sequence reads were generated and subsequently assembled into 271, 790 unigenes with an average length of 890 bp and a N50 of 1, 665 bp. 179, 641 unigenes could be annotated by public protein databases. The pairwise transcriptome comparsions of time course (6, 12, 24, 48 and 72 h post inoculation) were conducted 1) between inoculated and control samples of each genotype, 2) between inoculated samples of R and S genotypes. The linear dynamics of transcriptome profile was observed between adjacent samples for each genotype, two genotypes shared similar transcriptome pattern at early time points with most significant up regulation at 12 hour, and samples from R genotype at 24 h and S genotype at 48 h showed similar transcriptome pattern, significant differences of transcriptional profile were observed in pairwise comparisons between R and S genotypes. KEGG analysis showed that the primary metabolisms were inhibited in both genotypes and stronger inhibition in R genotype post inoculation. The defense related genes (R gene, LRR-RLK, cell wall genes, etc.) generally showed a genotype-specific down regulation and different expression between both genotypes.ConclusionThis transcriptome profiling provided the largest data set that explores the dynamic in crosstalk between peanut and R. solanacearum. The results suggested that the down-regulation of primary metabolism is contributed to the resistance difference between R and S genotypes. The genotype-specific expression pattern of defense related DEGs also contributed to the resistance difference between R and S genotype. This study will strongly contribute to better understand the molecular interaction between plant and R. solanacearum.

Functional Genomic Analysis of Aspergillus flavus Interacting with Resistant and Susceptible Peanut

In the Aspergillus flavus (A. flavus)–peanut pathosystem, development and metabolism of the fungus directly influence aflatoxin contamination. To comprehensively understand the molecular mechanism of A. flavus interaction with peanut, RNA-seq was used for global transcriptome profiling of A. flavus during interaction with resistant and susceptible peanut genotypes. In total, 67.46 Gb of high-quality bases were generated for A. flavus-resistant (af_R) and -susceptible peanut (af_S) at one (T1), three (T2) and seven (T3) days post-inoculation. The uniquely mapped reads to A. flavus reference genome in the libraries of af_R and af_S at T2 and T3 were subjected to further analysis, with more than 72% of all obtained genes expressed in the eight libraries. Comparison of expression levels both af_R vs. af_S and T2 vs. T3 uncovered 1926 differentially expressed genes (DEGs). DEGs associated with mycelial growth, conidial development and aflatoxin biosynthesis were up-regulated in af_S compared with af_R, implying that A. flavus mycelia more easily penetrate and produce much more aflatoxin in susceptible than in resistant peanut. Our results serve as a foundation for understanding the molecular mechanisms of aflatoxin production differences between A. flavus-R and -S peanut, and offer new clues to manage aflatoxin contamination in crops.

Transcriptome Analysis of a New Peanut Seed Coat Mutant for the Physiological Regulatory Mechanism Involved in Seed Coat Cracking and Pigmentation.

Seed-coat cracking and undesirable color of seed coat highly affects external appearance and commercial value of peanuts (Arachis hypogaea L.). With an objective to find genetic solution to the above problems, a peanut mutant with cracking and brown colored seed coat (testa) was identified from an EMS treated mutant population and designated as “peanut seed coat crack and brown color mutant line (pscb).” The seed coat weight of the mutant was almost twice of the wild type, and the germination time was significantly shorter than wild type. Further, the mutant had lower level of lignin, anthocyanin, proanthocyanidin content, and highly increased level of melanin content as compared to wild type. Using RNA-Seq, we examined the seed coat transcriptome in three stages of seed development in the wild type and the pscb mutant. The RNA-Seq analysis revealed presence of highly differentially expressed phenylpropanoid and flavonoid pathway genes in all the three seed development stages, especially at 40 days after flowering (DAF40). Also, the expression of polyphenol oxidases and peroxidase were found to be activated significantly especially in the late seed developmental stage. The genome-wide comparative study of the expression profiles revealed 62 differentially expressed genes common across all the three stages. By analyzing the expression patterns and the sequences of the common differentially expressed genes of the three stages, three candidate genes namely c36498_g1 (CCoAOMT1), c40902_g2 (kinesin), and c33560_g1 (MYB3) were identified responsible for seed-coat cracking and brown color phenotype. Therefore, this study not only provided candidate genes but also provided greater insights and molecular genetic control of peanut seed-coat cracking and color variation. The information generated in this study will facilitate further identification of causal gene and diagnostic markers for breeding improved peanut varieties with smooth and desirable seed coat color.