In silico polymorphism analysis for the development of simple sequence repeat and transposon markers and construction of linkage map in cultivated peanut

Identification of candidate genes associated with peanut pod length by combined analysis of QTL-seq and RNA-seq

Pod length (PL) is one of the major traits determining pod size and yield of peanut. Discovering the quantitative trait loci (QTL) and identifying candidate genes associated with PL are essential for breeding high-yield peanut. In this study, quantitative trait loci sequencing (QTL-seq) was performed using the F2 population constructed by a short-pod variety Tifrunner (Tif) and a long-pod line Lps, and a 0.77 Mb genomic region on chromosome 07 was identified as the candidate region for PL. Then, the candidate region was narrowed to a 265.93 kb region by traditional QTL approach. RNA-seq analysis showed that there were four differentially expressed genes (DEGs) in the candidate region, among which Arahy.PF2L6F (AhCDC48) and Arahy.P4LK2T (AhTAA1) were speculated to be PL-related candidate genes. These results were informative for the elucidation of the underlying regulatory mechanism in peanut pod length and would facilitate further identification of valuable target genes.

Designing future peanut: the power of genomics-assisted breeding

KEY MESSAGE

Integrating GAB methods with high-throughput phenotyping, genome editing, and speed breeding hold great potential in designing future smart peanut cultivars to meet market and food supply demands. Cultivated peanut (Arachis hypogaea L.), a legume crop greatly valued for its nourishing food, cooking oil, and fodder, is extensively grown worldwide. Despite decades of classical breeding efforts, the actual on-farm yield of peanut remains below its potential productivity due to the complicated interplay of genotype, environment, and management factors, as well as their intricate interactions. Integrating modern genomics tools into crop breeding is necessary to fast-track breeding efficiency and rapid progress. When combined with speed breeding methods, this integration can substantially accelerate the breeding process, leading to faster access of improved varieties to farmers. Availability of high-quality reference genomes for wild diploid progenitors and cultivated peanuts has accelerated the process of gene/quantitative locus discovery, developing markers and genotyping assays as well as a few molecular breeding products with improved resistance and oil quality. The use of new breeding tools, e.g., genomic selection, haplotype-based breeding, speed breeding, high-throughput phenotyping, and genome editing, is probable to boost genetic gains in peanut. Moreover, renewed attention to efficient selection and exploitation of targeted genetic resources is also needed to design high-quality and high-yielding peanut cultivars with main adaptation attributes. In this context, the combination of genomics-assisted breeding (GAB), genome editing, and speed breeding hold great potential in designing future improved peanut cultivars to meet market and food supply demands.

Molecular tagging of seed size using MITE markers in an induced large seed mutant with higher cotyledon cell size in groundnut

A large seed mutant, TG 89 having a 76.7% increment in hundred kernel weight in comparison to its parent TG 26, was isolated from an electron beam-induced mutagenized population. Studies based on environmental scanning electron microscopy of both parent and mutant revealed that the mutant seed cotyledon had significantly bigger cell size than parent. A mapping population with 122 F2 plants derived from the mutant and a distant normal seed genotype (ICGV 15007) was utilized to map the QTL associated with higher HKW. Bulk segregant analysis revealed putative association of three markers with this mutant large seed trait. Further, genotyping of F2 individuals with polymorphic markers detected 14 linkage groups with a map distance of 1053 cM. QTL analysis revealed a significant additive major QTL for the mutant large seed trait on linkage group A05 explaining 12.7% phenotypic variation for the seed size. This QTL was located between flanking markers AhTE333 and AhTE810 having a map interval of 4.7 cM which corresponds to 90.65 to 107.24 Mbp in A05 chromosome, respectively. Within this genomic fragment, an ortholog of the BIG SEEDS 1 gene was found at 102,476,137 bp. Real-time PCR revealed down-regulation of this BIG SEEDS 1 gene in the mutant indicating a loss of function mutation giving rise to a large seed phenotype. This QTL was validated in 11 advanced breeding lines having large seed size from this mutant but with varied genetic backgrounds. This validation showcased a highly promising selection accuracy of 90.9% for the marker-assisted selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03909-0.

Identification of QTLs associated with very-long chain fatty acid (VLCFA) content via linkage mapping and BSA-seq in peanut

KEY MESSAGE

Three major QTLs qA01, qB04.1 and qB05 for VLCFA content and their corresponding allele-specific markers will benefit peanut low VLCFA breeding, and a candidate gene Arahy.IF1JV3 was predicted. Peanut is a globally significant oilseed crop worldwide, and contains a high content (20%) of saturated fatty acid (SFA) in its seeds. As high level SFA intake in human dietary may increase the cardiovascular disease risk, reducing the SFA content in peanut is crucial for improving its nutritional quality. Half of the SFAs in peanut are very long-chain fatty acids (VLCFA), so reducing the VLCFA content is a feasible strategy to decrease the total SFA content. Luoaowan with extremely low VLCFA (4.80%) was crossed with Jihua16 (8.00%) to construct an F2:4 population. Three major QTLs including qA01, qB04.1 and qB05 for VLCFA content were detected with 4.43 ~ 14.32% phenotypic variation explained through linkage mapping. Meanwhile, three genomic regions on chromosomes B03, B04 and B05 were identified via BSA-seq approach. Two co-localized intervals on chromosomes B04 (100.10 ~ 103.97 Mb) and B05 (6.39 ~ 10.90 Mb) were identified. With markers developed based on SNP/InDel variations in qA01 between the two parents, the remaining interval was refined to 103.58 ~ 111.14 Mb. A candidate gene Arahy.IF1JV3 encoding a β-ketoacyl-CoA synthase was found in qA01, and its expression level in Luoaowan was significantly lower than that in Jihua16. Allele-specific markers targeting qA01, qB04.1 and qB05 were developed and validated in F4 population, and an elite line with high oleic, low VLCFA (5.05%) and low SFA (11.48%) contents was selected. This study initially revealed the genetic mechanism of VLCFA content, built a marker-assisted selection system for low VLCFA breeding, and provided an effective method to decrease the SFA content in peanut.

Mapping Quantitative Trait Loci (QTLs) for Hundred-Pod and Hundred-Seed Weight under Seven Environments in a Recombinant Inbred Line Population of Cultivated Peanut (Arachis hypogaea L.)

The cultivated peanut (Arachis hypogaea L.) is a significant oil and cash crop globally. Hundred-pod and -seed weight are important components for peanut yield. To unravel the genetic basis of hundred-pod weight (HPW) and hundred-seed weight (HSW), in the current study, a recombinant inbred line (RIL) population with 188 individuals was developed from a cross between JH5 (JH5, large pod and seed weight) and M130 (small pod and seed weight), and was utilized to identify QTLs for HPW and HSW. An integrated genetic linkage map was constructed by using SSR, AhTE, SRAP, TRAP and SNP markers. This map consisted of 3130 genetic markers, which were assigned to 20 chromosomes, and covered 1998.95 cM with an average distance 0.64 cM. On this basis, 31 QTLs for HPW and HSW were located on seven chromosomes, with each QTL accounting for 3.7-10.8% of phenotypic variance explained (PVE). Among these, seven QTLs were detected under multiple environments, and two major QTLs were found on B04 and B08. Notably, a QTL hotspot on chromosome A08 contained seven QTLs over a 2.74 cM genetic interval with an 0.36 Mb physical map, including 18 candidate genes. Of these, Arahy.D52S1Z, Arahy.IBM9RL, Arahy.W18Y25, Arahy.CPLC2W and Arahy.14EF4H might play a role in modulating peanut pod and seed weight. These findings could facilitate further research into the genetic mechanisms influencing pod and seed weight in cultivated peanut.

BSA‑seq and genetic mapping identified candidate genes for branching habit in peanut

The candidate gene AhLBA1 controlling lateral branch angel of peanut was fine-mapped to a 136.65-kb physical region on chromosome 15 using the BSA-seq and QTL mapping. Lateral branch angel (LBA) is an important plant architecture trait of peanut, which plays key role in lodging, peg soil penetration and pod yield. However, there are few reports of fine mapping and quantitative trait loci (QTLs)/cloned genes for LBA in peanut. In this project, a mapping population was constructed using a spreading variety Tifrunner and the erect variety Fuhuasheng. Through bulked segregant analysis sequencing (BSA-seq), a major gene related to LBA, named as AhLBA1, was preliminarily mapped at the region of Chr.15: 150-160 Mb. Then, using traditional QTL approach, AhLBA1 was narrowed to a 1.12 cM region, corresponding to a 136.65-kb physical interval of the reference genome. Of the nine genes housed in this region, three of them were involved in hormone metabolism and regulation, including one "F-box protein" and two "2-oxoglutarate (2OG) and Fe(II)-dependent oxygenase (2OG oxygenase)" encoding genes. In addition, we found that the level of some classes of cytokinin (CK), auxin and ethylene showed significant differences between spreading and erect peanuts at the junction of main stem and lateral branch. These findings will aid further elucidation of the genetic mechanism of LBA in peanut and facilitating marker-assisted selection (MAS) in the future breeding program.

Multi-Omics Profiling Identifies Candidate Genes Controlling Seed Size in Peanut

Seed size is the major yield component and a key target trait that is selected during peanut breeding. However, the mechanisms that regulate peanut seed size are unknown. Two peanut mutants with bigger seed size were isolated in this study by ⁶⁰Co treatment of a common peanut landrace, Huayu 22, and were designated as the "big seed" mutant lines (hybs). The length and weight of the seed in hybs were about 118% and 170% of those in wild-type (WT), respectively. We adopted a multi-omics approach to identify the genomic locus underlying the hybs mutants. We performed whole genome sequencing (WGS) of WT and hybs mutants and identified thousands of large-effect variants (SNPs and indels) that occurred in about four hundred genes in hybs mutants. Seeds from both WT and hybs lines were sampled 20 days after flowering (DAF) and were used for RNA-Seq analysis; the results revealed about one thousand highly differentially expressed genes (DEGs) in hybs compared to WT. Using a method that combined large-effect variants with DEGs, we identified 45 potential candidate genes that shared gene product mutations and expression level changes in hybs compared to WT. Among the genes, two candidate genes encoding cytochrome P450 superfamily protein and NAC transcription factors may be associated with the increased seed size in hybs. The present findings provide new information on the identification and functional research into candidate genes responsible for the seed size phenotype in peanut.

Comparative repeatome analysis reveals new evidence on genome evolution in wild diploid Arachis (Fabaceae) species

Opposing changes in the abundance of satellite DNA and long terminal repeat (LTR) retroelements are the main contributors to the variation in genome size and heterochromatin amount in Arachis diploids. The South American genus Arachis (Fabaceae) comprises 83 species organized in nine taxonomic sections. Among them, section Arachis is characterized by species with a wide genome and karyotype diversity. Such diversity is determined mainly by the amount and composition of repetitive DNA. Here we performed computational analysis on low coverage genome sequencing to infer the dynamics of changes in major repeat families that led to the differentiation of genomes in diploid species (x = 10) of genus Arachis, focusing on section Arachis. Estimated repeat content ranged from 62.50 to 71.68% of the genomes. Species with different genome composition tended to have different landscapes of repeated sequences. Athila family retrotransposons were the most abundant and variable lineage among Arachis repeatomes, with peaks of transpositional activity inferred at different times in the evolution of the species. Satellite DNAs (satDNAs) were less abundant, but differentially represented among species. High rates of evolution of an AT-rich superfamily of satDNAs led to the differential accumulation of heterochromatin in Arachis genomes. The relationship between genome size variation and the repetitive content is complex. However, largest genomes presented a higher accumulation of LTR elements and lower contents of satDNAs. In contrast, species with lowest genome sizes tended to accumulate satDNAs in detriment of LTR elements. Phylogenetic analysis based on repetitive DNA supported the genome arrangement of section Arachis. Altogether, our results provide the most comprehensive picture on the repeatome dynamics that led to the genome differentiation of Arachis species.

Identification and characterization of transposable element AhMITE1 in the genomes of cultivated and two wild peanuts

Background

The cultivated peanut (Arachis hypogaea L., AABB) is an allotetraploid hybrid between two diploid peanuts, A. duranensis (AA genome) and A. ipaensis (BB genome). Miniature inverted-repeat transposable elements (MITEs), some of which are known as active nonautonomous DNA transposons with high copy numbers, play important roles in genome evolution and diversification. AhMITE1, a member of the MITE family of transposons, but information on the peanut genomes is still limited. Here, we analyzed AhMITE1, AuMITE1 and ApMITE1 in the cultivated (A. hypogaea) and two wild peanut (A. duranensis and A. ipaensis) genomes.

Results

The cultivated and the two wild peanut genomes harbored 142, 14 and 21 AhMITE1, AuMITE1 and ApMITE1 family members, respectively. These three family members exhibited highly conserved TIR sequences, and insertions preferentially occurred within 2 kb upstream and downstream of gene-coding and AT-rich regions. Phylogenetic and pairwise nucleotide diversity analysis showed that AhMITE1 and ApMITE1 family members have undergone one round of amplification bursts during the evolution of the peanut genome. PCR analyses were performed in 23 peanut varieties and demonstrated that AhMITE1 is an active transposon and that hybridization or chemical mutagenesis can promote the mobilization of AhMITE1.

Conclusions

AhMITE1, AuMITE1 and ApMITE1 family members were identified based on local BLAST search with MAK between the cultivated and the two wild peanut genomes. The phylogenetic, nucleotide diversity and variation copy numbers of AhMITE1, AuMITE1 and ApMITE1 members provides opportunities for investigating their roles during peanut evolution. These findings will contribute to knowledge on diversity of AhMITE1, provide information about the potential impact on the gene expression and promote the development of DNA markers in peanut.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08732-0.

Detection of a major QTL and development of KASP markers for seed weight by combining QTL-seq, QTL-mapping and RNA-seq in peanut

Combining QTL-seq, QTL-mapping and RNA-seq identified a major QTL and candidate genes, which contributed to the development of KASP markers and understanding of molecular mechanisms associated with seed weight in peanut. Seed weight, as an important component of seed yield, is a significant target of peanut breeding. However, relatively little is known about the quantitative trait loci (QTLs) and candidate genes associated with seed weight in peanut. In this study, three major QTLs on chromosomes A05, B02, and B06 were determined by applying the QTL-seq approach in a recombinant inbred line (RIL) population. Based on conventional QTL-mapping, these three QTL regions were successfully narrowed down through newly developed single nucleotide polymorphism (SNP) and simple sequence repeat markers. Among these three QTL regions, qSWB06.3 exhibited stable expression, contributing mainly to phenotypic variance across environments. Furthermore, differentially expressed genes (DEGs) were identified at the three seed developmental stages between the two parents of the RIL population. It was found that the DEGs were widely distributed in the ubiquitin-proteasome pathway, the serine/threonine-protein pathway, signal transduction of hormones and transcription factors. Notably, DEGs at the early stage were mostly involved in regulating cell division, whereas DEGs at the middle and late stages were primarily involved in cell expansion during seed development. The expression patterns of candidate genes related to seed weight in qSWB06.3 were investigated using quantitative real-time PCR. In addition, the allelic diversity of qSWB06.3 was investigated in peanut germplasm accessions. The marker Ah011475 has higher efficiency for discriminating accessions with different seed weights, and it would be useful as a diagnostic marker in marker-assisted breeding. This study provided insights into the genetic and molecular mechanisms of seed weight in peanut.

Insights into the Genomic Architecture of Seed and Pod Quality Traits in the U.S. Peanut Mini-Core Diversity Panel

Traits such as seed weight, shelling percent, percent sound mature kernels, and seed dormancy determines the quality of peanut seed. Few QTL (quantitative trait loci) studies using biparental mapping populations have identified QTL for seed dormancy and seed grade traits. Here, we report a genome-wide association study (GWAS) to detect marker–trait associations for seed germination, dormancy, and seed grading traits in peanut. A total of 120 accessions from the U.S. peanut mini-core collection were evaluated for seed quality traits and genotyped using Axiom SNP (single nucleotide polymorphism) array for peanut. We observed significant variation in seed quality traits in different accessions and different botanical varieties. Through GWAS, we were able to identify multiple regions associated with sound mature kernels, seed weight, shelling percent, seed germination, and dormancy. Some of the genomic regions that were SNP associated with these traits aligned with previously known QTLs. For instance, QTL for seed dormancy has been reported on chromosome A05, and we also found SNP on the same chromosome associated with seed dormancy, explaining around 20% of phenotypic variation. In addition, we found novel genomic regions associated with seed grading, seed germination, and dormancy traits. SNP markers associated with seed quality and dormancy identified here can accelerate the selection process. Further, exploring the function of candidate genes identified in the vicinity of the associated marker will assist in understanding the complex genetic network that governs seed quality.

Fine-mapping and gene candidate analysis for AhRt1, a major dominant locus responsible for testa color in cultivated peanut

AhRt1 controlling red testa color in peanut was fine-mapped to an interval of 580 kb on chromosome A03, and one gene encoding bHLH transcriptional factor was identified as the putative candidate gene. Peanut with red testa has higher nutritional and economic value than the traditional pink testa varieties. Identification of genes controlling red testa color will accelerate the breeding program and facilitate uncovering the genetic mechanism. In this study, in order to identify gene underlying the red testa color in peanut, a F₂ population derived from a cross between a pink testa peanut variety "Fuhua 8" and a red testa variety "Quanhonghua 1" was constructed. The genetic analysis for the F₂ population revealed that the red testa color was controlled by one single dominant locus. This locus, named as AhRt1 (Arachis hypogaea Red Testa 1), was preliminary identified in chromosome A03 by BSA-sequencing analysis. Using a segregation mapping population, AhRt1 was fine-mapped to a 580-kb genomic region by substitution mapping strategy. Gene candidate analysis suggested that one predicted gene encoding bHLH transcriptional factor may be the possible candidate gene for AhRt1. A diagnostic marker closely linked to candidate gene has been developed for validating the fine-mapping result in different populations and peanut germplasm. Our findings will benefit the breeding program for developing new varieties with red testa color and laid foundation for map-based cloning gene responsible for red testa in peanut.

Dissection of the Genetic Basis of Yield-Related Traits in the Chinese Peanut Mini-Core Collection Through Genome-Wide Association Studies

Peanut is an important legume crop worldwide. To uncover the genetic basis of yield features and assist breeding in the future, we conducted genome-wide association studies (GWAS) for six yield-related traits of the Chinese peanut mini-core collection. The seed (pod) size and weight of the population were investigated under four different environments, and these traits showed highly positive correlations in pairwise combinations. We sequenced the Chinese peanut mini-core collection using genotyping-by-sequencing approach and identified 105,814 high-quality single-nucleotide polymorphisms (SNPs). The population structure analysis showed essentially subspecies patterns in groups and obvious geographical distribution patterns in subgroups. A total of 79 significantly associated loci (P < 4.73 × 10^-7) were detected for the six yield-related traits through GWAS. Of these, 31 associations were consistently detected in multiple environments, and 15 loci were commonly detected to be associated with multiple traits. Two major loci located on chromosomal pseudomolecules A06 and A02 showed pleiotropic effects on yield-related traits, explaining ∼20% phenotypic variations across environments. The two genomic regions were found 46 putative candidate genes based on gene annotation and expression profile. The diagnostic marker for the yield-related traits from non-synonymous SNP (Aradu-A06-107901527) was successfully validated, achieving a high correlation between nucleotide polymorphism and phenotypic variation. This study provided insights into the genetic basis of yield-related traits in peanut and verified one diagnostic marker to facilitate marker-assisted selection for developing high-yield peanut varieties.

Comprehensive Transcriptome Analyses Reveal Candidate Genes for Variation in Seed Size/Weight During Peanut (Arachis hypogaea L.) Domestication

Seed size/weight, a key domestication trait, is also an important selection target during peanut breeding. However, the mechanisms that regulate peanut seed development are unknown. We re-sequenced 12 RNA samples from developing seeds of two cultivated peanut accessions (Lines 8106 and 8107) and wild Arachis monticola at 15, 30, 45, and 60 days past flowering (DPF). Transcriptome analyses showed that ∼36,000 gene loci were expressed in each of the 12 RNA samples, with nearly half exhibiting moderate (2 ≤ FPKM < 10) expression levels. Of these genes, 12.2% (4,523) were specifically expressed during seed development, mainly at 15 DPF. Also, ∼12,000 genes showed significant differential expression at 30, 45, and/or 60 DPF within each of the three peanut accessions, accounting for 31.8-34.1% of the total expressed genes. Using a method that combined comprehensive transcriptome analysis and previously mapped QTLs, we identified several candidate genes that encode transcription factor TGA7, topless-related protein 2, IAA-amino acid hydrolase ILR1-like 5, and putative pentatricopeptide repeat-containing (PPR) protein. Based on sequence variations identified in these genes, SNP markers were developed and used to genotype both 30 peanut landraces and a genetic segregated population, implying that EVM0025654 encoding a PPR protein may be associated with the increased seed size/weight of the cultivated accessions in comparison with the allotetraploid wild peanut. Our results provide additional knowledge for the identification and functional research into candidate genes responsible for the seed size/weight phenotype in peanut.

Fine-Mapping of a Wild Genomic Region Involved in Pod and Seed Size Reduction on Chromosome A07 in Peanut (Arachis hypogaea L.)

Fruit and seed size are important yield component traits that have been selected during crop domestication. In previous studies, Advanced Backcross Quantitative Trait Loci (AB-QTL) and Chromosome Segment Substitution Line (CSSL) populations were developed in peanut by crossing the cultivated variety Fleur11 and a synthetic wild allotetraploid (Arachis ipaensis × Arachis duranensis)^4x. In the AB-QTL population, a major QTL for pod and seed size was detected in a ~5 Mb interval in the proximal region of chromosome A07. In the CSSL population, the line 12CS_091, which carries the QTL region and that produces smaller pods and seeds than Fleur11, was identified. In this study, we used a two-step strategy to fine-map the seed size QTL region on chromosome A07. We developed new SSR and SNP markers, as well as near-isogenic lines (NILs) in the target QTL region. We first located the QTL in ~1 Mb region between two SSR markers, thanks to the genotyping of a large F₂ population of 2172 individuals and a single marker analysis approach. We then used nine new SNP markers evenly distributed in the refined QTL region to genotype 490 F₃ plants derived from 88 F₂, and we selected 10 NILs. The phenotyping of the NILs and marker/trait association allowed us to narrowing down the QTL region to a 168.37 kb chromosome segment, between the SNPs Aradu_A07_1148327 and Aradu_A07_1316694. This region contains 22 predicted genes. Among these genes, Aradu.DN3DB and Aradu.RLZ61, which encode a transcriptional regulator STERILE APETALA-like (SAP) and an F-box SNEEZY (SNE), respectively, were of particular interest. The function of these genes in regulating the variation of fruit and seed size is discussed. This study will contribute to a better knowledge of genes that have been targeted during peanut domestication.

Pod and Seed Trait QTL Identification To Assist Breeding for Peanut Market Preferences

Although seed and pod traits are important for peanut breeding, little is known about the inheritance of these traits. A recombinant inbred line (RIL) population of 156 lines from a cross of Tifrunner x NC 3033 was genotyped with the Axiom_Arachis1 SNP array and SSRs to generate a genetic map composed of 1524 markers in 29 linkage groups (LG). The genetic positions of markers were compared with their physical positions on the peanut genome to confirm the validity of the linkage map and explore the distribution of recombination and potential chromosomal rearrangements. This linkage map was then used to identify Quantitative Trait Loci (QTL) for seed and pod traits that were phenotyped over three consecutive years for the purpose of developing trait-associated markers for breeding. Forty-nine QTL were identified in 14 LG for seed size index, kernel percentage, seed weight, pod weight, single-kernel, double-kernel, pod area and pod density. Twenty QTL demonstrated phenotypic variance explained (PVE) greater than 10% and eight more than 20%. Of note, seven of the eight major QTL for pod area, pod weight and seed weight (PVE >20% variance) were attributed to NC 3033 and located in a single linkage group, LG B06_1. In contrast, the most consistent QTL for kernel percentage were located on A07/B07 and derived from Tifrunner.

High-density SNP map facilitates fine mapping of QTLs and candidate genes discovery for Aspergillus flavus resistance in peanut (Arachis hypogaea)

Two novel resistant QTLs mapped and candidate genes identified for Aspergillus flavus resistance in cultivated peanut using SLAF-seq. Aflatoxin contamination in peanuts caused by Aspergillus flavus is a serious food safety issue for human health around the world. Host plant resistance to fungal infection and reduction in aflatoxin are crucial for mitigating this problem. Identification of the resistance-linked markers can be used in marker-assisted breeding for varietal development. Here we report construction of two high-density genetic linkage maps with 1975 SNP loci and 5022 SNP loci, respectively. Two consistent quantitative trait loci (QTL) were identified as qRAF-3-1 and qRAF-14-1, which located on chromosomes A03 and B04, respectively. QTL qRAF-3-1 was mapped within 1.67 cM and had more than 19% phenotypic variance explained (PVE), while qRAF-14-1 was located within 1.34 cM with 5.15% PVE. While comparing with the reference genome, the mapped QTLs, qRAF-3-1 and qRAF-14-1, were located within a physical distance of 1.44 Megabase pair (Mbp) and 2.22 Mbp, harboring 67 and 137 genes, respectively. Among the identified candidate genes, six genes with the same function were found within both QTLs regions. In addition, putative disease resistance RPP13-like protein 1 (RPP13), lipoxygenase (Lox), WRKY transcription factor (WRKY) and cytochrome P450 71B34 genes were also identified. Using microarray analysis, genes responded to A. flavus infection included coding for RPP13, pentatricopeptide repeat-containing-like protein, and Lox which may be possible candidate genes for resistance to A. flavus. The QTLs and candidate genes will further facilitate marker development and validation of genes for deployment in the molecular breeding programs against A. flavus in peanuts.

Dissection of the genetic basis of oil content in Chinese peanut cultivars through association mapping

Background

Peanut is one of the primary sources for vegetable oil worldwide, and enhancing oil content is the main objective in several peanut breeding programs of the world. Tightly linked markers are required for faster development of high oil content peanut varieties through genomics-assisted breeding (GAB), and association mapping is one of the promising approaches for discovery of such associated markers.

Results

An association mapping panel consisting of 292 peanut varieties extensively distributed in China was phenotyped for oil content and genotyped with 583 polymorphic SSR markers. These markers amplified 3663 alleles with an average of 6.28 alleles per locus. The structure, phylogenetic relationship, and principal component analysis (PCA) indicated two subgroups majorly differentiating based on geographic regions. Genome-wide association analysis identified 12 associated markers including one (AGGS1014_2) highly stable association controlling up to 9.94% phenotypic variance explained (PVE) across multiple environments. Interestingly, the frequency of the favorable alleles for 12 associated markers showed a geographic difference. Two associated markers (AGGS1014_2 and AHGS0798) with 6.90–9.94% PVE were verified to enhance oil content in an independent RIL population and also indicated selection during the breeding program.

Conclusion

This study provided insights into the genetic basis of oil content in peanut and verified highly associated two SSR markers to facilitate marker-assisted selection for developing high-oil content breeding peanut varieties.

Miniature inverted-repeat transposable elements (MITEs), derived insertional polymorphism as a tool of marker systems for molecular plant breeding

Plant molecular breeding is expected to give significant gains in cultivar development through development and utilization of suitable molecular marker systems for genetic diversity analysis, rapid DNA fingerprinting, identification of true hybrids, trait mapping and marker-assisted selection. Transposable elements (TEs) are the most abundant component in a genome and being used as genetic markers in the plant molecular breeding. Here, we review on the high copious transposable element belonging to class-II DNA TEs called "miniature inverted-repeat transposable elements" (MITEs). MITEs are ubiquitous, short and non-autonomous DNA transposable elements which have a tendency to insert into genes and genic regions have paved a way for the development of functional DNA marker systems in plant genomes. This review summarises the characteristics of MITEs, principles and methodologies for development of MITEs based DNA markers, bioinformatics tools and resources for plant MITE discovery and their utilization in crop improvement.

Developing Transposable Element Marker System for Molecular Breeding

Transposable element (TE) marker system was developed considering the useful properties of the transposable elements such as their large number in the animal and plant genomes, high rate of insertion polymorphism, and ease of detection. Various methods have been employed for developing a large number of TE markers in several crop plants for genomics studies. Here we describe some of these methods including the recent whole genome search. We also review the application of TE markers in molecular breeding.

QTL identification for seed weight and size based on a high-density SLAF-seq genetic map in peanut (Arachis hypogaea L.)

BACKGROUND

The cultivated peanut is an important oil and cash crop grown worldwide. To meet the growing demand for peanut production each year, genetic studies and enhanced selection efficiency are essential, including linkage mapping, genome-wide association study, bulked-segregant analysis and marker-assisted selection. Specific locus amplified fragment sequencing (SLAF-seq) is a powerful tool for high density genetic map (HDGM) construction and quantitative trait loci (QTLs) mapping. In this study, a HDGM was constructed using SLAF-seq leading to identification of QTL for seed weight and size in peanut.

RESULTS

A recombinant inbred line (RIL) population was advanced from a cross between a cultivar 'Huayu36' and a germplasm line '6-13' with contrasting seed weight, size and shape. Based on the cultivated peanut genome, a HDGM was constructed with 3866 loci consisting of SLAF-seq and simple sequence repeat (SSR) markers distributed on 20 linkage groups (LGs) covering a total map distance of 1266.87 cM. Phenotypic data of four seed related traits were obtained in four environments, which mostly displayed normal distribution with varied levels of correlation. A total of 27 QTLs for 100 seed weight (100SW), seed length (SL), seed width (SW) and length to width ratio (L/W) were identified on 8 chromosomes, with LOD values of 3.16-31.55 and explaining phenotypic variance (PVE) from 0.74 to 83.23%. Two stable QTL regions were identified on chromosomes 2 and 16, and gene content within these regions provided valuable information for further functional analysis of yield component traits.

CONCLUSIONS

This study represents a new HDGM based on the cultivated peanut genome using SLAF-seq and SSRs. QTL mapping of four seed related traits revealed two stable QTL regions on chromosomes 2 and 16, which not only facilitate fine mapping and cloning these genes, but also provide opportunity for molecular breeding of new peanut cultivars with improved seed weight and size.

Construction of a genetic linkage map of Lentinula edodes based on SSR, SRAP and TRAP markers

Genetic mapping is a basic tool for eukaryotic genomic research. It allows the localization of genes or quantitative trait loci (QTLs) and map-based cloning. In this study, we constructed a linkage map based on DNA samples from a commercial strain L808, including two parental monokaryons and 93 single spore isolates considered with segregating to 1:1:1:1 at four mating types (A₁B₁, A₁B₂, A₂B₁ and A₂B₂). Using Simple Sequence Repeats (SSR), Sequence Related Amplified Polymorphism (SRAP), Target Region Amplified Polymorphism (TRAP) molecular markers, 182 molecular markers and two mating factors were located on 11 linkage groups (LGs). The total length of the map was 948.083 centimorgan (cM), with an average marker interval distance of 4.817 cM. Only two gaps spanning more than 20 cM was observed. The probability of 20 cM, 10 cM, 5 cM genetic distance cover one marker was 99.68%, 94.36%, 76.43% in our genetic linkage map, respectively. This is the first linkage map of Lentinula edodes using SSR markers, which provides essential information for quantitative trait analyses and improvement of genome assembly.

Identification of major consensus QTLs for seed size and minor QTLs for pod traits in cultivated groundnut (Arachis hypogaea L.)

Hundred kernel weight is an important indicator for large-seeded genotype selection. A recombinant inbred line population was used to decipher the genetic architecture of seed size and three pod traits in cultivated groundnut based on the phenotypic data from six and three environments, respectively. The study revealed a consensus major QTL for HKW in B07 group that explained 10.5-23.9% phenotypic variation due to seed size. Further, two other minor QTLs were identified in B03 and B08 group for the seed size. Two minor QTLs for pod beak were positioned in B03 and A08. A minor QTL for pod reticulation was also mapped in the same map interval with the pod beak QTL in A08. Another minor QTL for pod constriction was co-mapped with the minor QTL for HKW in B08. The other minor QTL for pod constriction was placed in the neighboring map interval with the consensus QTL for seed size in B07 that suggests linkage of pod constriction with large seed trait. Analysis of the flanking markers profile in 71 cultivated groundnut genotypes revealed a strong association of pPGPseq_2E06 marker with large seed trait.

Genetics, genomics and breeding of groundnut (Arachis hypogaea L.)

Groundnut is an important food and oil crop in the semiarid tropics, contributing to household food consumption and cash income. In Asia and Africa, yields are low attributed to various production constraints. This review paper highlights advances in genetics, genomics and breeding to improve the productivity of groundnut. Genetic studies concerning inheritance, genetic variability and heritability, combining ability and trait correlations have provided a better understanding of the crop's genetics to develop appropriate breeding strategies for target traits. Several improved lines and sources of variability have been identified or developed for various economically important traits through conventional breeding. Significant advances have also been made in groundnut genomics including genome sequencing, marker development and genetic and trait mapping. These advances have led to a better understanding of the groundnut genome, discovery of genes/variants for traits of interest and integration of marker-assisted breeding for selected traits. The integration of genomic tools into the breeding process accompanied with increased precision of yield trialing and phenotyping will increase the efficiency and enhance the genetic gain for release of improved groundnut varieties.

Discovery of genomic regions and candidate genes controlling shelling percentage using QTL-seq approach in cultivated peanut (Arachis hypogaea L.)

Cultivated peanut (Arachis hypogaea L.) is an important grain legume providing high-quality cooking oil, rich proteins and other nutrients. Shelling percentage (SP) is the 2nd most important agronomic trait after pod yield and this trait significantly affects the economic value of peanut in the market. Deployment of diagnostic markers through genomics-assisted breeding (GAB) can accelerate the process of developing improved varieties with enhanced SP. In this context, we deployed the QTL-seq approach to identify genomic regions and candidate genes controlling SP in a recombinant inbred line population (Yuanza 9102 × Xuzhou 68-4). Four libraries (two parents and two extreme bulks) were constructed and sequenced, generating 456.89-790.32 million reads and achieving 91.85%-93.18% genome coverage and 14.04-21.37 mean read depth. Comprehensive analysis of two sets of data (Yuanza 9102/two bulks and Xuzhou 68-4/two bulks) using the QTL-seq pipeline resulted in discovery of two overlapped genomic regions (2.75 Mb on A09 and 1.1 Mb on B02). Nine candidate genes affected by 10 SNPs with non-synonymous effects or in UTRs were identified in these regions for SP. Cost-effective KASP (Kompetitive Allele-Specific PCR) markers were developed for one SNP from A09 and three SNPs from B02 chromosome. Genotyping of the mapping population with these newly developed KASP markers confirmed the major control and stable expressions of these genomic regions across five environments. The identified candidate genomic regions and genes for SP further provide opportunity for gene cloning and deployment of diagnostic markers in molecular breeding for achieving high SP in improved varieties.

Heterochromatin evolution in Arachis investigated through genome-wide analysis of repetitive DNA

The most conspicuous difference among chromosomes and genomes in Arachis species, the patterns of heterochromatin, was mainly modeled by differential amplification of different members of one superfamily of satellite DNAs. Divergence in repetitive DNA is a primary driving force for genome and chromosome evolution. Section Arachis is karyotypically diverse and has six different genomes. Arachis glandulifera (D genome) has the most asymmetric karyotype and the highest reproductive isolation compared to the well-known A and B genome species. These features make A. glandulifera an interesting model species for studying the main repetitive components that accompanied the genome and chromosome diversification in the section. Here, we performed a genome-wide analysis of repetitive sequences in A. glandulifera and investigated the chromosome distribution of the identified satellite DNA sequences (satDNAs). LTR retroelements, mainly the Ty3-gypsy families "Fidel/Feral" and "Pipoka/Pipa", were the most represented. Comparative analyses with the A and B genomes showed that many of the previously described transposable elements (TEs) were differently represented in the D genome, and that this variation accompanied changes in DNA content. In addition, four major satDNAs were characterized. Agla_CL8sat was the major component of pericentromeric heterochromatin, while Agla_CL39sat, Agla_CL69sat, and Agla_CL122sat were found in heterochromatic and/or euchromatic regions. Even though Agla_CL8sat belong to a different family than that of the major satDNA (ATR-2) found in the heterochromatin of the A, K, and F genomes, both satDNAs are members of the same superfamily. This finding suggests that closely related satDNAs of an ancestral library were differentially amplified leading to the major changes in the heterochromatin patterns that accompanied the karyotype and genome differentiation in Arachis.

Construction of High-Density Genetic Map and Mapping Quantitative Trait Loci for Growth Habit-Related Traits of Peanut (Arachis hypogaea L.)

Plant growth habit is an important and complex agronomic trait and is associated with yield, disease resistance, and mechanized harvesting in peanuts. There are at least two distinct growth habits (erect and prostrate) and several intermediate forms existing in the peanut germplasm. A recombinant inbred line population containing 188 individuals was developed from a cross of "Jihua 5" and "M130" for genetically dissecting the architecture of the growth habit. A new high-density genetic linkage map was constructed by using specific locus amplified fragment sequencing technology. The map contains 2,808 single-nucleotide polymorphism markers distributed on 20 linkage groups with a total length of 1,308.20 cM and an average inter-marker distance of 0.47 cM. The quantitative trait locus (QTL) analysis of the growth habit-related traits was conducted based on phenotyping data from seven environments. A total of 39 QTLs for growth habit-related traits was detected on 10 chromosomes explaining 4.55-27.74% of the phenotypic variance, in which 6 QTLs were for lateral branch angle, 8 QTLs were for extent radius, 7 QTLs were for the index of plant type, 11 QTLs were for main stem height, and 7 QTLs were for lateral branch length. Among these QTLs, 12 were co-localized on chromosome B05 spanning an approximately 0.17 Mb physical interval in comparison with the allotetraploid reference genome of "Tifrunner." Analysis of the co-localized genome region has shown that the putative genes are involved in light and hormones and will facilitate peanut growth habit molecular breeding and study of peanut domestication.

A Major and Stable QTL for Bacterial Wilt Resistance on Chromosome B02 Identified Using a High-Density SNP-Based Genetic Linkage Map in Cultivated Peanut Yuanza 9102 Derived Population

Bacterial wilt (BW) is one of the important diseases limiting the production of peanut (Arachis hypogaea L.) worldwide. The sufficient precise information on the quantitative trait loci (QTL) for BW resistance is essential for facilitating gene mining and applying in molecular breeding. Cultivar Yuanza 9102 is BW resistant, bred from wide cross between cultivated peanut Baisha 1016 and a wild diploid peanut species A. chacoense with BW resistance. In this study, we aim to map the major QTLs related to BW-resistance in Yuanza 9102. A high density SNP-based genetic linkage map was constructed through double-digest restriction-site-associated DNA sequencing (ddRADseq) technique based on Yuanza 9102 derived recombinant inbred lines (RILs) population. The map contained 2,187 SNP markers distributed on 20 linkage groups (LGs) spanning 1566.10 cM, and showed good synteny with AA genome from A. duranensis and BB genome from A. ipaensis. Phenotypic frequencies of BW resistance among RIL population showed two-peak distribution in four environments. Four QTLs explaining 5.49 to 23.22% phenotypic variance were identified to be all located on chromosome B02. The major QTL, qBWB02.1 (12.17–23.33% phenotypic variation explained), was detected in three environments showing consistent and stable expression. Furthermore, there was positive additive effect among these major and minor QTLs. The major QTL region was mapped to a region covering 2.3 Mb of the pseudomolecule B02 of A. ipaensis which resides in 21 nucleotide-binding site -leucine-rich repeat (NBS-LRR) encoding genes. The result of the major stable QTL (qBWB02.1) not only offers good foundation for discovery of BW resistant gene but also provide opportunity for deployment of the QTL in marker-assisted breeding in peanut.

High-Density Genetic Map Construction and Identification of QTLs Controlling Oleic and Linoleic Acid in Peanut using SLAF-seq and SSRs

The cultivated peanut, A. hypogaea L., is an important oil and food crop globally.High-density genetic linkage mapping is a valuable and effective method for exploring complex quantitative traits. In this context, a recombinant inbred line (RIL) of 146 lines was developed by crossing Huayu28 and P76. We developed 433,679 high-quality SLAFs, of which 29,075 were polymorphic. 4,817 SLAFs were encoded and grouped into different segregation patterns. A high-resolution genetic map containing 2,334 markers (68 SSRs and 2,266 SNPs) on 20 linkage groups (LGs) spanning 2586.37 cM was constructed for peanut. The average distance between adjacent markers was 2.25 cM. Based on phenotyping in seven environments, QTLs for oleic acid (C18:1), linoleic acid (C18:2) and the ratio of oleic acid to linoleic acid (O/L) were identified and positioned on linkage groups A03, A04, A09, B09 and B10. Marker2575339 and Marker2379598 in B09 were associated with C18:1, C18:2 and O/L in seven environments, Marker4391589 and Marker4463600 in A09 were associated with C18:1, C18:2 and O/L in six environments. This map exhibits high resolution and accuracy, which will facilitate QTL discovery for essential agronomic traits in peanut.

Chromosomes A07 and A05 associated with stable and major QTLs for pod weight and size in cultivated peanut (Arachis hypogaea L.)

Co-localized intervals and candidate genes were identified for major and stable QTLs controlling pod weight and size on chromosomes A07 and A05 in an RIL population across four environments. Cultivated peanut (Arachis hypogaea L.) is an important legume crops grown in > 100 countries. Hundred-pod weight (HPW) is an important yield trait in peanut, but its underlying genetic mechanism was not well studied. In this study, a mapping population (Xuhua 13 × Zhonghua 6) with 187 recombinant inbred lines (RILs) was developed to map quantitative trait loci (QTLs) for HPW together with pod length (PL) and pod width (PW) by both unconditional and conditional QTL analyses. A genetic map covering 1756.48 cM was constructed with 817 markers. Additive effects, epistatic interactions, and genotype-by-environment interactions were analyzed using the phenotyping data generated across four environments. Twelve additive QTLs were identified on chromosomes A05, A07, and A08 by unconditional analysis, and five of them (qPLA07, qPLA05.1, qPWA07, qHPWA07.1, and qHPWA05.2) showed major and stable expressions in all environments. Conditional QTL mapping found that PL had stronger influences on HPW than PW. Notably, qHPWA07.1, qPLA07, and qPWA07 that explained 17.93-43.63% of the phenotypic variations of the three traits were co-localized in a 5 cM interval (1.48 Mb in physical map) on chromosome A07 with 147 candidate genes related to catalytic activity and metabolic process. In addition, qHPWA05.2 and qPLA05.1 were co-localized with minor QTL qPWA05.2 to a 1.3 cM genetic interval (280 kb in physical map) on chromosome A05 with 12 candidate genes. This study provides a comprehensive characterization of the genetic components controlling pod weight and size as well as candidate QTLs and genes for improving pod yield in future peanut breeding.

Development of AhMITE1 markers through genome-wide analysis in peanut (Arachis hypogaea L.)

OBJECTIVE

In peanut, the DNA polymorphism is very low despite enormous phenotypic variations. This limits the use of genomics-assisted breeding to enhance peanut productivity. This study aimed to develop and validate new AhMITE1 and cleaved amplified polymorphic sequences (CAPS) markers.

RESULTS

In total, 2957 new AhMITE1 markers were developed in addition to identifying 465 already reported markers from the whole genome re-sequencing data (WGRS) of 33 diverse genotypes of peanut. The B sub-genome (1620) showed more number of markers than the A sub-genome (1337). Distribution also varied among the chromosomes of both the sub-genomes. Further, 52.6% of the markers were from genic regions; where 31.0% were from intronic regions and 5.2% were from exonic regions. Of the 343 randomly selected markers, 82.2% showed amplification validation, with up to 35.5% polymorphism. From the SNPs on the A03, B01, B02 and B03 chromosomes, 11,730 snip-SNPs (potential CAPS sites) were identified, and 500 CAPS markers were developed from chromosome A03. Of these markers, 30.0% showed validation and high polymorphism. This study demonstrated the potential of the WGRS data to develop AhMITE1 and CAPS markers, which showed high level of validation and polymorphism. These marker resources will be useful for various genetic studies and mapping in peanut.

Stable QTLs for Plant Height on Chromosome A09 Identified From Two Mapping Populations in Peanut (Arachis hypogaea L.)

The peanut (Arachis hypogaea L.) is an important grain legume extensively cultivated worldwide, supplying edible oil and protein for human consumption. As in many other crops, plant height is a crucial factor in determining peanut architecture traits and has a unique effect on resistance to lodging and efficiency of mechanized harvesting as well as yield. Currently, the genetic basis underlying plant height remains unclear in peanut, which have hampered marker-assisted selection in breeding. In this study, we conducted a quantitative trait locus (QTL) analysis for peanut plant height by using two recombinant inbred line (RIL) populations including "Yuanza 9102 × Xuzhou 68-4 (YX)" and "Xuhua 13 × Zhonghua 6 (XZ)". In the YX population, 38 QTLs including 10 major QTLs from 9 chromosomes were detected in 4 environments, and 8 consensus QTLs integrated by meta-analysis expressed stably across multiple environments. In the XZ population, 3 major QTLs and seven minor QTLs from 6 chromosomes were detected across 3 environments. Generally, most major QTLs from the two populations were located on pseudomolecule chromosome 9 of Arachis duranesis (A09), indicating there would be key genes on A09 controlling plant height. Further analysis revealed that qPHA09.1a from the XZ population and one consensus QTL, cqPHA09.d from the YX population were co-localized in a reliable 3.4 Mb physical interval on A09, which harbored 161 genes including transcription factors and enzymes related to signaling transduction and cell wall formation. The major and stable QTLs identified in this study may be useful for further gene cloning and identification of molecular markers applicable for breeding.

Development of a High-Density Genetic Map Based on Specific Length Amplified Fragment Sequencing and Its Application in Quantitative Trait Loci Analysis for Yield-Related Traits in Cultivated Peanut

High-density genetic maps (HDGMs) are very useful for genomic studies and quantitative trait loci (QTL) mapping. However, the low frequency of DNA polymorphisms in peanut has limited the quantity of available markers and hindered the construction of a HDGM. This study generated a peanut genetic map with the highest number of high-quality SNPs based on specific locus amplified fragment sequencing (SLAF-seq) technology and a newly constructed RIL population ("ZH16" × "sd-H1"). The constructed HDGM included 3,630 SNP markers belonging to 2,636 bins on 20 linkage groups (LGs), and it covers 2,098.14 cM in length, with an average marker distance of 0.58 cM. This HDGM was applied for the following collinear comparison, scaffold anchoring and analysis of genomic characterization including recombination rates and segregation distortion in peanut. For QTL mapping of investigated 14 yield-related traits, a total of 62 QTLs were detected on 12 chromosomes across 3 environments, and the co-localization of QTLs was observed for these traits which were significantly correlated on phenotype. Two stable co-located QTLs for seed- and pod-related traits were significantly identified in the chromosomal end of B06 and B07, respectively. The construction of HDGM and QTL analysis for yield-related traits in this study provide useful information for fine mapping and functional analysis of genes as well as molecular marker-assisted breeding.

Mapping of important taxonomic and productivity traits using genic and non-genic transposable element markers in peanut (Arachis hypogaea L.)

A mapping population of recombinant inbred lines (RILs) derived from TMV 2 and its mutant, TMV 2-NLM was employed for mapping important taxonomic and productivity traits using genic and non-genic transposable element markers in peanut. Single nucleotide polymorphism and copy number variation using RAD-Sequencing data indicated very limited polymorphism between TMV 2 and TMV 2-NLM. But phenotypically they differed significantly for many taxonomic and productivity traits. Also, the RIL population showed significant variation for a few additional agronomic traits. A genetic linkage map of 1,205.66 cM was constructed using 91 genic and non-genic Arachis hypogaea transposable element (AhTE) markers. Using single marker analysis and QTL analysis, the markers with high phenotypic variance explained (PVE) were identified for branching pattern (32.3%), number of primary and secondary branches (19.9% and 28.4%, respectively), protein content (26.4%), days to 50% flowering (22.0%), content of oleic acid (15.1%), test weight (13.6%) and pod width (12.0%). Three genic markers (AhTE0357, AhTE0391, AhTE0025) with Arachis hypogaea miniature inverted-repeat transposable element (AhMITE1) activity in the genes Araip.TG1BL (B02 chromosome), Aradu.7N61X (A09 chromosome) and Aradu.7065G (A07 chromosome), respectively showed strong linkage with these taxonomic, productivity and quality traits. Since TMV 2 and TMV 2-NLM differed subtly at DNA level, the background noise in detecting the marker-trait associations was minimum; therefore, the markers identified in this study for the taxonomic and productivity traits may be significant and useful in peanut molecular breeding.

Development of SSR markers and identification of major quantitative trait loci controlling shelling percentage in cultivated peanut (Arachis hypogaea L.)

A total of 204,439 SSR markers were developed in diploid genomes, and 25 QTLs for shelling percentage were identified in a RIL population across 4 years including five consistent QTLs. Cultivated peanut (Arachis hypogaea L.) is an important grain legume providing edible oil and protein for human nutrition. Genome sequences of its diploid ancestors, Arachis duranensis and A. ipaensis, were reported, but their SSRs have not been well exploited and utilized hitherto. Shelling percentage is an important economic trait and its improvement has been one of the major objectives in peanut breeding programs. In this study, the genome sequences of A. duranensis and A. ipaensis were used to develop SSR markers, and a mapping population (Yuanza 9102 × Xuzhou 68-4) with 195 recombinant inbred lines was used to map QTLs controlling shelling percentage. The numbers of newly developed SSR markers were 84,383 and 120,056 in the A. duranensis and A. ipaensis genomes, respectively. Genotyping of the mapping population was conducted with both newly developed and previously reported markers. QTL analysis using the phenotyping data generated in Wuhan across four consecutive years and genotyping data of 830 mapped loci identified 25 QTLs with 4.46-17.01% of phenotypic variance explained in the four environments. Meta-analysis revealed five consistent QTLs that could be detected in at least two environments. Notably, the consistent QTL cqSPA09 was detected in all four environments and explained 10.47-17.01% of the phenotypic variance. The segregation in the progeny of a residual heterozygous line confirmed that the cpSPA09 locus had additive effect in increasing shelling percentage. These consistent and major QTL regions provide opportunity not only for further gene discovery, but also for the development of functional markers for breeding.

Co-localization of major quantitative trait loci for pod size and weight to a 3.7 cM interval on chromosome A05 in cultivated peanut (Arachis hypogaea L.)

Background

Cultivated peanut (Arachis hypogaea L.), an important source of edible oil and protein, is widely grown in tropical and subtropical areas of the world. Genetic improvement of yield-related traits is essential for improving yield potential of new peanut varieties. Genomics-assisted breeding (GAB) can accelerate the process of genetic improvement but requires linked markers for the traits of interest. In this context, we developed a recombinant inbred line (RIL) mapping population (Yuanza 9102 × Xuzhou 68-4) with 195 individuals and used to map quantitative trait loci (QTLs) associated with three important pod features, namely pod length, pod width and hundred-pod weight.

Results

QTL analysis using the phenotyping data generated across four environments in two locations and genotyping data on 743 mapped loci identified 15 QTLs for pod length, 11 QTLs for pod width and 16 QTLs for hundred-pod weight. The phenotypic variation explained (PVE) ranged from 3.68 to 27.84%. Thirteen QTLs were consistently detected in at least two environments and three QTLs (qPLA05.7, qPLA09.3 and qHPWA05.6) were detected in all four environments indicating their consistent and stable expression. Three major QTLs, detected in at least three environments, were found to be co-localized to a 3.7 cM interval on chromosome A05, and they were qPLA05.7 for pod length (16.89–27.84% PVE), qPWA05.5 for pod width (13.73–14.12% PVE), and qHPWA05.6 for hundred-pod weight (13.75–26.82% PVE). This 3.7 cM linkage interval corresponds to ~2.47 Mb genomic region of the pseudomolecule A05 of A. duranensis, including 114 annotated genes related to catalytic activity and metabolic process.

Conclusions

This study identified three major consistent and stable QTLs for pod size and weight which were co-localized in a 3.7 cM interval on chromosome A05. These QTL regions not only offer further investigation for gene discovery and development of functional markers but also provide opportunity for deployment of these QTLs in GAB for improving yield in peanut.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3456-x) contains supplementary material, which is available to authorized users.

Genetic mapping of yield traits using RIL population derived from Fuchuan Dahuasheng and ICG6375 of peanut (Arachis hypogaea L.)

The genetic architecture determinants of yield traits in peanut (Arachis hypogaea L.) are poorly understood. In the present study, an effort was made to map quantitative trait loci (QTLs) for yield traits using recombinant inbred lines (RIL). A genetic linkage map was constructed containing 609 loci, covering a total of 1557.48 cM with an average distance of 2.56 cM between adjacent markers. The present map exhibited good collinearity with the physical map of diploid species of Arachis. Ninety-two repeatable QTLs were identified for 11 traits including height of main stem, total branching number, and nine pod- and seed-related traits. Of the 92 QTLs, 15 QTLs were expressed across three environments and 65 QTLs were newly identified. Twelve QTLs for the height of main stem and the pod- and seed-related traits explaining more than 10 % of phenotypic variation showed a great potential for marker-assisted selection in improving these traits. The trait-by-trait meta-analysis revealed 33 consensus QTLs. The consensus QTLs and other QTLs were further integrated into 29 pleiotropic unique QTLs with the confidence interval of 1.86 cM on average. The significant co-localization of QTLs was consistent with the significant phenotypic correlations among these traits. The complexity of the genetic architecture of yield traits was demonstrated. The present QTLs for pod- and seed-related traits could be the most fundamental genetic factors contributing to the yield traits in peanut. The results provide a good foundation for fine mapping, cloning and designing molecular breeding of favorable genes in peanut.

Genetic Variation and Association Mapping of Seed-Related Traits in Cultivated Peanut (Arachis hypogaea L.) Using Single-Locus Simple Sequence Repeat Markers

Cultivated peanut (Arachis hypogaea L.) is an allotetraploid (AABB, 2n = 4x = 40), valued for its edible oil and digestible protein. Seed size and weight are important agronomical traits significantly influence the yield and nutritional composition of peanut. However, the genetic basis of seed-related traits remains ambiguous. Association mapping is a powerful approach for quickly and efficiently exploring the genetic basis of important traits in plants. In this study, a total of 104 peanut accessions were used to identify molecular markers associated with seed-related traits using 554 single-locus simple sequence repeat (SSR) markers. Most of the accessions had no or weak relationship in the peanut panel. The linkage disequilibrium (LD) decayed with the genetic distance of 1cM at the genome level and the LD of B subgenome decayed faster than that of the A subgenome. Large phenotypic variation was observed for four seed-related traits in the association panel. Using mixed linear model with population structure and kinship, a total of 30 significant SSR markers were detected to be associated with four seed-related traits (P < 1.81 × 10^-3) in different environments, which explained 11.22-32.30% of the phenotypic variation for each trait. The marker AHGA44686 was simultaneously and repeatedly associated with seed length and hundred-seed weight in multiple environments with large phenotypic variance (26.23 ∼ 32.30%). The favorable alleles of associated markers for each seed-related trait and the optimal combination of favorable alleles of associated markers were identified to significantly enhance trait performance, revealing a potential of utilization of these associated markers in peanut breeding program.

Development and deployment of a high-density linkage map identified quantitative trait loci for plant height in peanut (Arachis hypogaea L.)

Plant height is one of the most important architecture traits in crop plants. In peanut, the genetic basis of plant height remains ambiguous. In this context, we genotyped a recombinant inbred line (RIL) population with 140 individuals developed from a cross between two peanut varieties varying in plant height, Zhonghua 10 and ICG 12625. Genotyping data was generated for 1,175 SSR and 42 transposon polymorphic markers and a high-density genetic linkage map was constructed with 1,219 mapped loci covering total map length of 2,038.75 cM i.e., accounted for nearly 80% of the peanut genome. Quantitative trait locus (QTL) analysis using genotyping and phenotyping data for three environments identified 8 negative-effect QTLs and 10 positive-effect QTLs for plant height. Among these QTLs, 8 QTLs had a large contribution to plant height that explained ≥10% phenotypic variation. Two major-effect consensus QTLs namely cqPHA4a and cqPHA4b were identified with stable performance across three environments. Further, the allelic recombination of detected QTLs proved the existence of the phenomenon of transgressive segregation for plant height in the RIL population. Therefore, this study not only successfully reported a high-density genetic linkage map of peanut and identified genomic region controlling plant height but also opens opportunities for further gene discovery and molecular breeding for plant height in peanut.

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.

Genomic survey sequencing for development and validation of single-locus SSR markers in peanut (Arachis hypogaea L.)

BACKGROUND

Single-locus markers have many advantages compared with multi-locus markers in genetic and breeding studies because their alleles can be assigned to particular genomic loci in diversity analyses. However, there is little research on single-locus SSR markers in peanut. Through the de novo assembly of DNA sequencing reads of A. hypogaea, we developed single-locus SSR markers in a genomic survey for better application in genetic and breeding studies of peanut.

RESULTS

In this study, DNA libraries with four different insert sizes were used for sequencing with 150 bp paired-end reads. Approximately 237 gigabases of clean data containing 1,675,631,984 reads were obtained after filtering. These reads were assembled into 2,102,446 contigs with an N50 length of 1,782 bp, and the contigs were further assembled into 1,176,527 scaffolds with an N50 of 3,920 bp. The total length of the assembled scaffold sequences was 2.0 Gbp, and 134,652 single-locus SSRs were identified from 375,180 SSRs. Among these developed single-locus SSRs, trinucleotide motifs were the most abundant, followed by tetra-, di-, mono-, penta- and hexanucleotide motifs. The most common motif repeats for the various types of single-locus SSRs have a tendency to be A/T rich. A total of 1,790 developed in silico single-locus SSR markers were chosen and used in PCR experiments to confirm amplification patterns. Of them, 1,637 markers that produced single amplicons in twelve inbred lines were considered putative single-locus markers, and 290 (17.7 %) showed polymorphisms. A further F2 population study showed that the segregation ratios of the 97 developed SSR markers, which showed polymorphisms between the parents, were consistent with the Mendelian inheritance law for single loci (1:2:1). Finally, 89 markers were assigned to an A. hypogaea linkage map. A subset of 100 single-locus SSR markers was shown to be highly stable and universal in a collection of 96 peanut accessions. A neighbor-joining tree of this natural population showed that genotypes have obviously correlation with botanical varieties.

CONCLUSIONS

We have shown that the detection of single-locus SSR markers from a de novo genomic assembly of a combination of different-insert-size libraries is highly efficient. This is the first report of the development of genome-wide single-locus markers for A. hypogaea, and the markers developed in this study will be useful for gene tagging, sequence scaffold assignment, linkage map construction, diversity analysis, variety identification and association mapping in peanut.

Characterization and Transferable Utility of Microsatellite Markers in the Wild and Cultivated Arachis Species

Microsatellite or simple sequence repeat (SSR) is one of the most widely distributed molecular markers that have been widely utilized to assess genetic diversity and genetic mapping for important traits in plants. However, the understanding of microsatellite characteristics in Arachis species and the currently available amount of high-quality SSR markers remain limited. In this study, we identified 16,435 genome survey sequences SSRs (GSS-SSRs) and 40,199 expressed sequence tag SSRs (EST-SSRs) in Arachis hypogaea and its wild relative species using the publicly available sequence data. The GSS-SSRs had a density of 159.9–239.8 SSRs/Mb for wild Arachis and 1,015.8 SSR/Mb for cultivated Arachis, whereas the EST-SSRs had the density of 173.5–384.4 SSR/Mb and 250.9 SSRs/Mb for wild and cultivated Arachis, respectively. The trinucleotide SSRs were predominant across Arachis species, except that the dinucleotide accounted for most in A. hypogaea GSSs. From Arachis GSS-SSR and EST-SSR sequences, we developed 2,589 novel SSR markers that showed a high polymorphism in six diverse A. hypogaea accessions. A genetic linkage map that contained 540 novel SSR loci and 105 anchor SSR loci was constructed by case of a recombinant inbred lines F₆ population. A subset of 82 randomly selected SSR markers were used to screen 39 wild and 22 cultivated Arachis accessions, which revealed a high transferability of the novel SSRs across Arachis species. Our results provided informative clues to investigate microsatellite patterns across A. hypogaea and its wild relative species and potentially facilitate the germplasm evaluation and gene mapping in Arachis species.

The genome sequences of Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanut

Cultivated peanut (Arachis hypogaea) is an allotetraploid with closely related subgenomes of a total size of ∼2.7 Gb. This makes the assembly of chromosomal pseudomolecules very challenging. As a foundation to understanding the genome of cultivated peanut, we report the genome sequences of its diploid ancestors (Arachis duranensis and Arachis ipaensis). We show that these genomes are similar to cultivated peanut's A and B subgenomes and use them to identify candidate disease resistance genes, to guide tetraploid transcript assemblies and to detect genetic exchange between cultivated peanut's subgenomes. On the basis of remarkably high DNA identity of the A. ipaensis genome and the B subgenome of cultivated peanut and biogeographic evidence, we conclude that A. ipaensis may be a direct descendant of the same population that contributed the B subgenome to cultivated peanut.

Quantitative trait locus analysis for pod- and kernel-related traits in the cultivated peanut (Arachis hypogaea L.)

BACKGROUND

The cultivated peanut (Arachis hypogaea L.) is an important oil and food crop in the world. Pod- and kernel-related traits are direct factors involved in determining the yield of the peanut. However, the genetic basis underlying pod- and kernel-related traits in the peanut remained largely unknown, which hampered the improvement of peanut through marker-assisted selection. To understand the genetic basis underlying pod- and kernel-related traits in the peanut and provide more useful information for marker-assisted breeding, we conducted quantitative trait locus (QTL) analysis for pod length and width and seed length and width by use of two F2:3 populations derived from cultivar Fuchuan Dahuasheng × ICG 6375 (FI population) and cultivar Xuhua 13 × cultivar Zhonghua 6 (XZ population) in this study.

RESULTS

Two genetic maps containing 347 and 228 polymorphic markers were constructed for FI and XZ populations respectively. In total, 39 QTLs explaining 1.25-26.11% of the phenotypic variations were detected in two populations. For the FI population, 26 QTLs were detected between the two environments, among which twelve were not mapped before. For the XZ population, thirteen QTLs were detected, among which eight were not reported before. One QTL for pod width was repeatedly mapped between the two populations.

CONCLUSION

The QTL analyses for pod length and width and seed length and width were conducted in this study using two mapping populations. Novel QTLs were identified, which included two for pod length, four for pod width, five for seed length and one for seed width in the FI population, and three for pod length, three for pod width and two for seed width in the XZ population. Our results will be helpful for improving pod- and seed-related traits in peanuts through marker-assisted selection.

Exciting journey of 10 years from genomes to fields and markets: Some success stories of genomics-assisted breeding in chickpea, pigeonpea and groundnut

Legume crops such as chickpea, pigeonpea and groundnut, mostly grown in marginal environments, are the major source of nutrition and protein to the human population in Asia and Sub-Saharan Africa. These crops, however, have a low productivity, mainly due to their exposure to several biotic and abiotic stresses in the marginal environments. Until 2005, these crops had limited genomics resources and molecular breeding was very challenging. During the last decade (2005-2015), ICRISAT led demand-driven innovations in genome science and translated the massive genome information in breeding. For instance, large-scale genomic resources including draft genome assemblies, comprehensive genetic and physical maps, thousands of SSR markers, millions of SNPs, several high-throughput as well as low cost marker genotyping platforms have been developed in these crops. After mapping several breeding related traits, several success stories of translational genomics have become available in these legumes. These include development of superior lines with enhanced drought tolerance in chickpea, enhanced and pyramided resistance to Fusarium wilt and Ascochyta blight in chickpea, enhanced resistance to leaf rust in groundnut, improved oil quality in groundnut and utilization of markers for assessing purity of hybrids/parental lines in pigeonpea. Some of these stories together with future prospects have been discussed.

PSR: polymorphic SSR retrieval

BACKGROUND

With the advent of high-throughput sequencing technologies large-scale identification of microsatellites became affordable and was especially directed to non-model species. By contrast, few efforts have been published toward the automatic identification of polymorphic microsatellites by exploiting sequence redundancy. Few tools for genotyping microsatellite repeats have been implemented so far that are able to manage huge amount of sequence data and handle the SAM/BAM file format. Most of them have been developed for and tested on human or model organisms with high quality reference genomes.

RESULTS

In this note we describe polymorphic SSR retrieval (PSR), a read counter and simple sequence repeat (SSR) length polymorphism detection tool. It is written in Perl and was developed to identify length polymorphisms in perfect microsatellites exploiting next generation sequencing (NGS) data. PSR has been developed bearing in mind plant non-model species for which de novo transcriptome assembly is generally the first sequence resource available to be used for SSR-mining. PSR is divided into two modules: the read-counting module (PSR_read_retrieval) identifies all the reads that cover the full-length of perfect microsatellites; the comparative module (PSR_poly_finder) detects both heterozygous and homozygous alleles at each microsatellite locus across all genotypes under investigation. Two threshold values to call a length polymorphism and reduce the number of false positives can be defined by the user: the minimum number of reads overlapping the repetitive stretch and the minimum read depth. The first parameter determines if the microsatellite-containing sequence must be processed or not, while the second one is decisive for the identification of minor alleles. PSR was tested on two different case studies. The first study aims at the identification of polymorphic SSRs in a set of de novo assembled transcripts defined by RNA-sequencing of two different plant genotypes. The second research activity aims to investigate sequence variations within a collection of newly sequenced chloroplast genomes. In both the cases PSR results are in agreement with those obtained by capillary gel separation.

CONCLUSION

PSR has been specifically developed from the need to automate the gene-based and genome-wide identification of polymorphic microsatellites from NGS data. It overcomes the limits related to the existing and time-consuming efforts based on tools developed in the pre-NGS era.

Quantitative trait locus analysis of agronomic and quality-related traits in cultivated peanut (Arachis hypogaea L.)

SSR-based QTL mapping provides useful information for map-based cloning of major QTLs and can be used to improve the agronomic and quality traits in cultivated peanut by marker-assisted selection. Cultivated peanut (Arachis hypogaea L.) is an allotetraploid species (AABB, 2n = 4× = 40), valued for its edible oil and digestible protein. Linkage mapping has been successfully conducted for most crops, and it has been applied to detect the quantitative trait loci (QTLs) of biotic and abiotic traits in peanut. However, the genetic basis of agronomic and quality-related traits remains unclear. In this study, high levels of phenotypic variation, broad-sense heritability and significant correlations were observed for agronomic and quality-related traits in an F 2:3 population. A genetic linkage map was constructed for cultivated peanut containing 470 simple sequence repeat (SSR) loci, with a total length of 1877.3 cM and average distance of 4.0 cM between flanking markers. For 10 agronomic traits, 24 QTLs were identified and each QTL explained 1.69-18.70 % of the phenotypic variance. For 8 quality-related traits, 12 QTLs were identified that explained 1.72-20.20 % of the phenotypic variance. Several QTLs for multiple traits were overlapped, reflecting the phenotypic correlation between these traits. The majority of QTLs exhibited obvious dominance or over-dominance effects on agronomic and quality traits, highlighting the importance of heterosis for breeding. A comparative analysis revealed genomic duplication and arrangement of peanut genome, which aids the assembly of scaffolds in genomic sequencing of Arachis hypogaea. Our QTL analysis results enabled us to clearly understand the genetic base of agronomic and quality traits in cultivated peanut, further accelerating the progress of map-based cloning of major QTLs and marker-assisted selection in future breeding.

Generation of peanut drought tolerant plants by pingyangmycin-mediated in vitro mutagenesis and hydroxyproline-resistance screening

In order to enlarge the potential resources of drought-tolerant peanuts, we conducted in vitro mutagenesis with Pingyangmycin (PYM) as the mutagen as well as directed screening on a medium supplemented with Hydroxyproline (HYP). After being extracted from mature seeds (cv. Huayu 20), the embryonic leaflets were cultured on somatic embryogenesis-induction medium with 4 mg/L PYM and the generated embryos were successively transferred to a germination medium with 4 and then 8 mmol/L HYP to screen HYP-tolerant plantlets. After that, these plantlets were grafted and transplanted to the experimental field. In the next generation, all seeds were sown in the field, and phenotype variation and trait segregation can be observed in most of the offspring (M₂ generation). The M₃ generation individuals were subjected to drought stress at the seedling stages. The activities of SOD and POD were substantially increased in eight offspring of 11 HYP-tolerant, regenerated plants than in their mutagenic parents. To determine the correlation between mutant phenotypes and genomic modification, we carried out a comparison of the DNA polymorphisms between the mutagenic parents and 13 M₃ generation individuals from different HYP-tolerant, regenerated plants with SSR primers. Results showed that most mutants and parent plants had signs of polymorphisms. Under drought stress, some M₃ generation individuals of 10 original HYP-tolerant, regenerated plants produced more pods than the mutagenic parent; twenty individuals among them produced >60 g pods/plant. M₄-generation seeds were tested for quality characteristics by Near Infrared Spectroscopy (NIS) and nine individuals with higher protein content (>30%) and 21 individuals with higher oil content (>58%) were screened. We concluded that the use of PYM-based in vitro mutagenesis in combination with directed screening with HYP is effective for the creation of potential drought-tolerant mutants of peanut.