An InDel-based linkage map of hot pepper (Capsicum annuum)

Development of SNP and InDel markers by genome resequencing and transcriptome sequencing in radish (Raphanus sativus L.)

BACKGROUND

Single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) are the most abundant genetic variations and widely distribute across the genomes in plant. Development of SNP and InDel markers is a valuable tool for genetics and genomic research in radish (Raphanus sativus L.).

RESULTS

In this study, a total of 366,679 single nucleotide polymorphisms (SNPs) and 97,973 insertion-deletion (InDel) markers were identified based on genome resequencing between 'YZH' and 'XHT'. In all, 53,343 SNPs and 4,257 InDels were detected in two cultivars by transcriptome sequencing. Among the InDel variations, 85 genomic and 15 transcriptomic InDels were newly developed and validated PCR. The 100 polymorphic InDels markers generated 207 alleles among 200 Chinese radish germplasm, with an average 2.07 of the number of alleles (Na) and with an average 0.33 of the polymorphism information content (PIC). Population structure and phylogenetic relationship revealed that the radish cultivars from northern China were clustered together and the southwest China cultivars were clustered together. RNA-Seq analysis revealed that 11,003 differentially expressed genes (DEGs) were identified between the two cultivars, of which 5,020 were upregulated and 5,983 were downregulated. In total, 145 flowering time-related DGEs were detected, most of which were involved in flowering time integrator, circadian clock/photoperiod autonomous, and vernalization pathways. In flowering time-related DGEs region, 150 transcriptomic SNPs and 9 InDels were obtained.

CONCLUSIONS

The large amount of SNPs and InDels identified in this study will provide a valuable marker resource for radish genetic and genomic studies. The SNPs and InDels within flowering time-related DGEs provide fundamental insight into for dissecting molecular mechanism of bolting and flowering in radish.

A Genomic Evaluation of Six Selected Inbred Lines of the Naturalized Plants of Milk Thistle (Silybum marianum L. Gaertn.) in Korea

Milk thistle (Silybum marianum) belongs to the Asteraceae family and is a medicinal plant native to the Mediterranean Basin. Silymarin in achene is a widely used herbal product for chronic liver disease. There is growing interest in natural medicine using milk thistle in Korea, but the raw material completely relies on imports. Despite its economic importance, phenotypic evaluations of native resources of milk thistle in Korea have not been carried out. In addition, genomic research and molecular marker development are very limited in milk thistle. In this study, we evaluated 220 milk thistle resources consisting of 172 accessions collected from the domestic market, and 48 accessions isolated from 6 accessions distributed by the National Agrobiodiversity Center in Korea. Six plant characteristics (height, seed weight, number of flowers, seed weight per flower, spine length, and color at harvest) were measured, and six samples (M01-M06) were selected to represent the genetic diversity of the population for genomic research. To develop PCR-based and co-dominant insertion/deletion (InDel) markers, we performed genome-wide InDel detection by comparing the whole-genome resequencing data of the six selected accessions with the reference genome sequence (GCA_001541825). As a result, 177 InDel markers with high distinguishability and reproducibility were selected from the 30,845 InDel variants. Unknowingly imported alien plant resources could easily be genetically mixed, and jeopardized seed purity can cause continuous difficulties in the development of high value-added agricultural platforms utilizing natural products. The selected plant materials and 177 validated InDel markers developed via whole-genome resequencing analysis could be valuable resources for breeding, conservation, and ecological studies of natives to Korea, along with acceleration of Silybum marianum industrialization.

Development of Omni InDel and supporting database for maize

Insertions-deletions (InDels) are the second most abundant molecular marker in the genome and have been widely used in molecular biology research along with simple sequence repeats (SSR) and single-nucleotide polymorphisms (SNP). However, InDel variant mining and marker development usually focuses on a single type of dimorphic InDel, which does not reflect the overall InDel diversity across the genome. Here, we developed Omni InDels for maize, soybean, and rice based on sequencing data and genome assembly that included InDel variants with base lengths from 1 bp to several Mb, and we conducted a detailed classification of Omni InDels. Moreover, we screened a set of InDels that are easily detected and typed (Perfect InDels) from the Omni InDels, verified the site authenticity using 3,587 germplasm resources from 11 groups, and analyzed the germplasm resources. Furthermore, we developed a Multi-InDel set based on the Omni InDels; each Multi-InDel contains multiple InDels, which greatly increases site polymorphism, they can be detected in multiple platforms such as fluorescent capillary electrophoresis and sequencing. Finally, we developed an online database website to make Omni InDels easy to use and share and developed a visual browsing function called "Variant viewer" for all Omni InDel sites to better display the variant distribution.

InDelGT: An integrated pipeline for extracting indel genotypes for genetic mapping in a hybrid population using next-generation sequencing data

Premise

Although several software packages are available for genotyping insertion/deletion (indel) polymorphisms in genomes using next-generation sequencing data, simultaneously calling indel genotypes across many individuals for use in genetic mapping remains challenging.

Methods and Results

We present an integrated pipeline, InDelGT, for the extraction of indel genotypes from a segregating population such as backcross or F2 lines, or from an F1 cross between outbred species. The InDelGT algorithm is implemented in three steps: generating an indel catalog, calling indel genotypes, and analyzing indel segregation. We demonstrated the use of the pipeline with an example data set from an F1 hybrid population of Populus and successfully constructed the two parental genetic linkage maps.

Conclusions

InDelGT is a practical tool that can quickly genotype a large number of indel markers within a population following Mendelian segregation. The InDelGT pipeline is freely available on GitHub (https://github.com/tongchf/InDelGT).

Genome-wide identification and development of InDel markers in tobacco (Nicotiana tabacum L.) using RAD-seq

Insertions and deletions (InDels) can be used as molecular markers in genetic studies and marker-assisted selection breeding. However, genetic improvement in tobacco has been hindered by limited genetic diversity information and relatedness within available germplasm. A Chinese tobacco variety, Yueyan-98, was resequenced using restriction-site associated DNA (RAD-seq) approach to develop InDel markers. In total, 32,884 InDel loci were detected between Yueyan-98 and the K326 reference sequence [18,598 (56.55%) deletions and 14,288 (43.45%) insertions], ranging from 1 to 62 bp in length. Of the 6,733 InDels (> 4 bp) that were suitable for polyacrylamide gel electrophoresis, 150 were randomly selected. These 150 InDels were unevenly distributed on 23 chromosomes, and the highest numbers of InDels were observed on chromosomes Nt05, Nt13, and Nt23. The average density of adjacent InDels was 19.36 Mb. Thirty-seven InDels were located in genic regions. Polymerase chain reaction (PCR)-based markers were developed to validate polymorphism; 113 (79.80%) of the 150 InDel markers showed polymorphism and were further used for genetic diversity analysis of 50 tobacco accessions (13 from China, 1 from Mexico, and 36 from the USA). The average expected heterozygosity (He) and polymorphism information content (PIC) values were 0.28 ± 0.16 and 0.38 ± 0.10, respectively. The average Shannon diversity index (I) was 0.34 ± 0.18, with genetic diversity ranging from 0.13-0.57. The 50 accessions were classified into two groups with a genetic similarity coefficient of 0.68. Principal coordinate analysis (PCoA) and population structure analysis showed similar results and divided the population into two groups unrelated to their geographical origins. AMOVA showed 4% variance among the population and the remaining 96% within the population, suggesting low genetic differentiation between two subpopulations. Furthermore, 10 InDels (19 alleles) were significantly identified for tobacco plant height using GLM+Q model at P < 0.005. Among these, three markers (Nt-I-26, Nt-I-41, and Nt-I-44) were detected in at least two environments, with phenotypic variance explained (PVE) ranging from 14.03 to 32.68%. The polymorphic InDel markers developed can be used for hybrid identification, genetic diversity, genetic linkage map construction, gene mapping, and MAS breeding programs of tobacco.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01187-3.

Genome-wide development of insertion-deletion (InDel) markers for Cannabis and its uses in genetic structure analysis of Chinese germplasm and sex-linked marker identification

Background

Cannabis sativa L., a dioecious plant derived from China, demonstrates important medicinal properties and economic value worldwide. Cannabis properties have been usually harnessed depending on the sex of the plant. To analyse the genetic structure of Chinese Cannabis and identify sex-linked makers, genome-wide insertion-deletion (InDel) markers were designed and used.

Results

In this study, a genome-wide analysis of insertion-deletion (InDel) polymorphisms was performed based on the recent genome sequences. In total, 47,558 InDels were detected between the two varieties, and the length of InDels ranged from 4 bp to 87 bp. The most common InDels were tetranucleotides, followed by pentanucleotides. Chromosome 5 exhibited the highest number of InDels among the Cannabis chromosomes, while chromosome 10 exhibited the lowest number. Additionally, 31,802 non-redundant InDel markers were designed, and 84 primers evenly distributed in the Cannabis genome were chosen for polymorphism analysis. A total of 38 primers exhibited polymorphisms among three accessions, and of the polymorphism primers, 14 biallelic primers were further used to analyse the genetic structure. A total of 39 fragments were detected, and the PIC value ranged from 0.1209 to 0.6351. According to the InDel markers and the flowering time, the 115 Chinese germplasms were divided into two subgroups, mainly composed of cultivars obtained from the northernmost and southernmost regions, respectively. Additional two markers, “Cs-I1–10” and “Cs-I1–15”, were found to amplify two bands (398 bp and 251 bp; 293 bp and 141 bp) in the male plants, while 389-bp or 293-bp bands were amplified in female plants. Using the two markers, the feminized and dioecious varieties could also be distinguished.

Conclusion

Based on the findings obtained herein, we believe that this study will facilitate the genetic improvement and germplasm conservation of Cannabis in China, and the sex-linked InDel markers will provide accurate sex identification strategies for Cannabis breeding and production.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07883-w.

Development and validation of genome-wide InDel markers with high levels of polymorphism in bitter gourd (Momordica charantia)

Background

The preferred choice for molecular marker development is identifying existing variation in populations through DNA sequencing. With the genome resources currently available for bitter gourd (Momordica charantia), it is now possible to detect genome-wide insertion-deletion (InDel) polymorphisms among bitter gourd populations, which guides the efficient development of InDel markers.

Results

Here, using bioinformatics technology, we detected 389,487 InDels from 61 Chinese bitter gourd accessions with an average density of approximately 1298 InDels/Mb. Then we developed a total of 2502 unique InDel primer pairs with a polymorphism information content (PIC) ≥0.6 distributed across the whole genome. Amplification of InDels in two bitter gourd lines ‘47–2–1-1-3’ and ‘04–17,’ indicated that the InDel markers were reliable and accurate. To highlight their utilization, the InDel markers were employed to construct a genetic map using 113 ‘47–2–1-1-3’ × ‘04–17’ F₂ individuals. This InDel genetic map of bitter gourd consisted of 164 new InDel markers distributed on 15 linkage groups with a coverage of approximately half of the genome.

Conclusions

This is the first report on the development of genome-wide InDel markers for bitter gourd. The validation of the amplification and genetic map construction suggests that these unique InDel markers may enhance the efficiency of genetic studies and marker-assisted selection for bitter gourd.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07499-0.

Fine mapping of the major anthracnose resistance QTL AnRGO5 in Capsicum chinense 'PBC932'

BACKGROUND

Colletotrichum species are the causal agents of anthracnose, a major disease affecting the yield and quality of pepper (Capsicum spp.). Colletotrichum scovillei is widespread in China, has strong pathogenicity and drug resistance, and causes anthracnose disease in pepper fruits that severely reduces production. Previously, an anti-anthracnose locus AnR_GO5 was mapped to the P5 chromosome on the basis of analyses of fruit at the green mature stage. The aim of this study was to narrow down the interval of this locus and identify the gene responsible for conferring resistance.

RESULTS

On the basis of results of re-sequencing of Capsicum chinense 'PBC932' and C. annuum '77013', we developed Kompetitive allele-specific PCR (KASPar) markers and insertion-deletion (InDel) markers linked to AnR_GO5 at the green mature fruit stage and used them to construct a genetic linkage map (42 markers, 24.4 cM in length). Using data obtained in phenotypic and genotypic analyses of BC₄S₁, BC₄S₂, and BC₄S₃ populations, AnR_GO5 was located between the markers P5in-2266-404 and P5in-2268-978 within a physical distance of 164 kb. This region contained five genes, including CA05g17730. CA05g17730 encodes 'R1C-3-like' putative late blight resistance protein homologs. The transcript level of CA05g17730 differed between 'PBC932' and '77013'. The structure of the CA05g17730 gene also differed between 'PBC932' and '77013'.

CONCLUSIONS

We narrowed down the QTL interval to a region containing five genes. These results will be useful for further research on the mechanisms of resistance to anthracnose, and for marker assisted selection for anthracnose-resistant capsicum lines.

Molecular mapping of dominant gene responsible for leaf curl virus resistance in chilli pepper (Capsicum annuum L.)

A resistant source (S-343) having monogenic dominant resistance to chilli leaf curl virus disease (ChiLCVD) has been identified at Punjab Agricultural University (PAU), Ludhiana. The F₂ mapping population of 204 plants was derived from the cross MS-341 (susceptible) × S-343 (resistant) to identify the linked marker with the disease-resistant gene. Out of the 685 single-sequence repeats (SSRs) used, only 160 primers showed parental polymorphism. These 160 polymorphic primers were used for bulk segregant analysis and only eight SSR primers were able to differentiate the resistant and susceptible bulks. The linkage analysis revealed that the two markers CA 516044 and PAU-LC-343-1 were found linked with the disease-resistant gene covering a total distance of 15.7 centimorgan (cM). The two primers CA 516044 and PAU-LC-343-1 were found located on chromosome 6 of the pepper genome at a genetic distance of 6.8 cM and 8.9 cM, respectively, from the resistant gene. The validation of linked markers was performed using 26 resistant and susceptible genotypes developed at PAU, Ludhiana by former researchers. The validation of the primers revealed that there was a correlation between phenotypic and genotypic data of the used genotypes, and these markers can be used for the marker-assisted breeding procedures for transferring ChiLCVD resistance until the gene-based markers will be developed. The markers described in this study are the first-ever molecular markers identified as linked to the ChiLCVD-resistant gene.

Fine mapping of restorer-of-fertility gene based on high-density genetic mapping and collinearity analysis in pepper (Capsicum annuum L.)

The pepper restorer-of-fertility (CaRf) gene was fine mapped based on conjoint analysis of recombinants and collinearity between the two pepper reference genomes. Capana06g003028, encoding an Rf-like PPR protein, was proposed as the strongest candidate for pepper CaRf based on sequence comparison and expression analysis. The cytoplasmic male sterility (CMS)/restorer-of-fertility (Rf) system not only provides an excellent model to dissect genetic interactions between the mitochondria and nucleus but also plays a vital role in high-efficiency hybrid seed production in crops including pepper (Capsicum spp.). Although two important CMS candidate genes, orf507 and Ψatp6-2, have previously been suggested, the pepper Rf gene (CaRf) has not yet been isolated. In this study, a high-density genetic map comprising 7566 SNP markers in 1944 bins was first constructed with the array genotyping results from 317 F₂ individuals. Based on this map, the CaRf gene was preliminarily mapped to a region of 1.15 Mb in length at the end of chromosome P6. Then, by means of a conjoint analysis of recombinants and collinearity between the two pepper reference genomes, an important candidate interval with 270.10 kb in length was delimited for CaRf. Finally, Capana06g003028, which encodes an Rf-like PPR protein, was proposed as the strongest candidate for CaRf based on sequence analysis and expression characteristics in sterile and fertile plants. The high-density genetic map and fine mapping results provided here lay a foundation for the application of molecular breeding, as well as cloning and functional analysis of CaRf, in pepper.

Target sequencing reveals genetic diversity, population structure, core-SNP markers, and fruit shape-associated loci in pepper varieties

BACKGROUND

The widely cultivated pepper (Capsicum spp.) is one of the most diverse vegetables; however, little research has focused on characterizing the genetic diversity and relatedness of commercial varieties grown in China. In this study, a panel of 92 perfect single-nucleotide polymorphisms (SNPs) was identified using re-sequencing data from 35 different C. annuum lines. Based on this panel, a Target SNP-seq genotyping method was designed, which combined multiplex amplification of perfect SNPs with Illumina sequencing, to detect polymorphisms across 271 commercial pepper varieties.

RESULTS

The perfect SNPs panel had a high discriminating capacity due to the average value of polymorphism information content, observed heterozygosity, expected heterozygosity, and minor allele frequency, which were 0.31, 0.28, 0.4, and 0.31, respectively. Notably, the studied pepper varieties were morphologically categorized based on fruit shape as blocky-, long horn-, short horn-, and linear-fruited. The long horn-fruited population exhibited the most genetic diversity followed by the short horn-, linear-, and blocky-fruited populations. A set of 35 core SNPs were then used as kompetitive allele-specific PCR (KASPar) markers, another robust genotyping technique for variety identification. Analysis of genetic relatedness using principal component analysis and phylogenetic tree construction indicated that the four fruit shape populations clustered separately with limited overlaps. Based on STRUCTURE clustering, it was possible to divide the varieties into five subpopulations, which correlated with fruit shape. Further, the subpopulations were statistically different according to a randomization test and F_st statistics. Nine loci, located on chromosomes 1, 2, 3, 4, 6, and 12, were identified to be significantly associated with the fruit shape index (p < 0.0001).

CONCLUSIONS

Target SNP-seq developed in this study appears as an efficient power tool to detect the genetic diversity, population relatedness and molecular breeding in pepper. Moreover, this study demonstrates that the genetic structure of Chinese pepper varieties is significantly influenced by breeding programs focused on fruit shape.

InDel Marker Based Estimation of Multi-Gene Allele Contribution and Genetic Variations for Grain Size and Weight in Rice (Oryza sativa L.)

The market success of any rice cultivar is exceedingly dependent on its grain appearance, as well as its grain yield, which define its demand by consumers as well as growers. The present study was undertaken to explore the contribution of nine major genes, qPE9~1, GW2, SLG7, GW5, GS3, GS7, GW8, GS5, and GS2, in regulating four size and weight related traits, i.e., grain length (GL), grain width (GW), grain thickness (GT), and thousand grain weight (TGW) in 204 diverse rice germplasms using Insertion/Deletion (InDel) markers. The studied germplasm displayed wide-ranging variability in the four studied traits. Except for three genes, all six genes showed considerable association with these traits with varying strengths. Whole germplasm of 204 genotypes could be categorized into three major clusters with different grain sizes and weights that could be utilized in rice breeding programs where grain appearance and weight are under consideration. The study revealed that TGW was 24.9% influenced by GL, 37.4% influenced by GW, and 49.1% influenced by GT. Hence, assuming the trend of trait selection, i.e., GT > GW > GL, for improving TGW in the rice yield enhancement programs. The InDel markers successfully identified a total of 38 alleles, out of which 27 alleles were major and were found in more than 20 genotypes. GL was associated with four genes (GS3, GS7, GW8, and GS2). GT was also found to be regulated by four different genes (GS3, GS7, GW8, and GS2) out of the nine studied genes. GW was found to be under the control of three studied genes (GW5, GW8, and GS2), whereas TGW was found to be under the influence of four genes (SLG7, GW5, GW8, and GS5) in the germplasm under study. The Unweighted Pair Group Method with Arithmetic means (UPGMA) tree based on the studied InDel marker loci segregated the whole germplasm into three distinct clusters with dissimilar grain sizes and weights. A two-dimensional scatter plot constructed using Principal Coordinate Analysis (PCoA) based on InDel markers further separated the 204 rice germplasms into four sub-populations with prominent demarcations of extra-long, long, medium, and short grain type germplasms that can be utilized in breeding programs accordingly. The present study could help rice breeders to select a suitable InDel marker and in formulation of breeding strategies for improving grain appearance, as well as weight, to develop rice varieties to compete international market demands with higher yield returns. This study also confirms the efficient application of InDel markers in studying diverse types of rice germplasm, allelic frequencies, multiple-gene allele contributions, marker-trait associations, and genetic variations that can be explored further.

InDel markers: An extended marker resource for molecular breeding in chickpea

Chickpea is one of the most important food legumes that holds the key to meet rising global food and nutritional demand. In order to deploy molecular breeding approaches in crop improvement programs, user friendly and cost effective marker resources remain prerequisite. The advent of next generation sequencing (NGS) technology has resulted in the generation of several thousands of markers as part of several large scale genome sequencing and re-sequencing initiatives. Very recently, PCR based Insertion-deletions (InDels) are becoming a popular gel based genotyping solution because of their co-dominant, inexpensive, and highly polymorphic nature. With an objective to expand marker resources for genomics assisted breeding (GAB) in chickpea, whole genome re-sequencing data generated on five parental lines of one interspecific (ICC 4958 × PI 489777) and two intra-specific (ICC 283 × ICC 8261 and ICC 4958 × ICC 1882) mapping populations, were used for identification of InDels. A total of 231,658 InDels were identified using Dindel software with default parameters. Further, a total of 8,307 InDels with ≥20 bp size were selected for development of gel based markers, of which primers could be designed for 7,523 (90.56%) markers. On average, markers appeared at a frequency of 1,038 InDels/LG with a maximum number of markers on CaLG04 (1,952 InDels) and minimum on CaLG08 (360 InDels). In order to validate these InDels, a total of 423 primer pairs were randomly selected and tested on the selected parental lines. A high amplification rate of 80% was observed ranging from 46.06 to 58.01% polymorphism rate across parents on 3% agarose gel. This study clearly reflects the usefulness of available sequence data for the development of genome-wide InDels in chickpea that can further contribute and accelerate a wide range of genetic and molecular breeding activities in chickpea.

Development and Application of InDel Markers for Capsicum spp. Based on Whole-Genome Re-Sequencing

Genome-wide identification of Insertion/Deletion polymorphisms (InDels) in Capsicum spp. was performed through comparing whole-genome re-sequencing data from two Capsicum accessions, C. annuum cv. G29 and C. frutescens cv. PBC688, with the reference genome sequence of C. annuum cv. CM334. In total, we identified 1,664,770 InDels between CM334 and PBC688, 533,523 between CM334 and G29, and 1,651,856 between PBC688 and G29. From these InDels, 1605 markers of 3-49 bp in length difference between PBC688 and G29 were selected for experimental validation: 1262 (78.6%) showed polymorphisms, 90 (5.6%) failed to amplify, and 298 (18.6%) were monomorphic. For further validation of these InDels, 288 markers were screened across five accessions representing five domesticated species. Of these assayed markers, 194 (67.4%) were polymorphic, 87 (30.2%) monomorphic and 7 (2.4%) failed. We developed three interspecific InDels, which associated with three genes and showed specific amplification in five domesticated species and clearly differentiated the interspecific hybrids. Thus, our novel PCR-based InDel markers provide high application value in germplasm classification, genetic research and marker-assisted breeding in Capsicum species.

Genome-wide identification and expression analyses of WRKY transcription factor family members from chickpea (Cicer arietinum L.) reveal their role in abiotic stress-responses

BACKGROUND

WRKY proteins play a vital role in the regulation of several imperative plant metabolic processes and pathways, especially under biotic and abiotic stresses. Although WRKY genes have been characterized in various major crop plants, their identification and characterization in pulse legumes is still in its infancy. Chickpea (Cicer arietinum L.) is the most important pulse legume grown in arid and semi-arid tropics.

OBJECTIVE

In silico identification and characterization of WRKY transcription factor-encoding genes in chickpea genome.

METHODS

For this purpose, a systematic genome-wide analysis was carried out to identify the non-redundant WRKY transcription factors in the chickpea genome.

RESULTS

We have computationally identified 70 WRKY-encoding non-redundant genes which were randomly distributed on all the chickpea chromosomes except chromosome 8. The evolutionary phylogenetic analysis classified the WRKY proteins into three major groups (I, II and III) and seven sub-groups (IN, IC, IIa, IIb, IIc, IId and IIe). The gene structure analysis revealed the presence of 2-7 introns among the family members. Along with the presence of absolutely conserved signatory WRKY domain, 19 different domains were also found to be conserved in a group-specific manner. Insights of gene duplication analysis revealed the predominant role of segmental duplications for the expansion of WRKY genes in chickpea. Purifying selection seems to be operated during the evolution and expansion of paralogous WRKY genes. The transcriptome data-based in silico expression analysis revealed the differential expression of CarWRKY genes in root and shoot tissues under salt, drought, and cold stress conditions. Moreover, some of these genes showed identical expression pattern under these stresses, revealing the possibility of involvement of these genes in conserved abiotic stress-response pathways.

CONCLUSION

This genome-wide computational analysis will serve as a base to accelerate the functional characterization of WRKY TFs especially under biotic and abiotic stresses.

Construction of a high-density genetic map and identification of QTLs for cucumber mosaic virus resistance in pepper (Capsicum annuum L.) using specific length amplified fragment sequencing (SLAF-seq)

Pepper (Capsicum) is one of the most important vegetable and spice crops. Aphid-transmitted cucumber mosaic virus (CMV) causes significant damage to pepper crops across the world. The genetic basis of CMV resistance in pepper is complex, and the mechanisms underlying resistance remain largely unknown. Here, we employed a SLAF-seq approach to generate a high-density genetic map of pepper. The map spanned 1,785.46 cM, containing 12,727 markers on 12 chromosomes, with a mean marker distance of 0.16 cM between adjacent markers. We used this map and the interval mapping (IM) and multiple QTL mapping (MQM) procedures to detect genetic regions associated with quantitative trait for CMV resistance. Three QTLs, qcmv11.1, qcmv11.2 and qcmv12.1, conferred resistance to CMV and showed trait variation of 10.2%, 19.2% and 7.3% respectively. Our results will help to develop markers linked to CMV-resistant QTLs to improve pepper resistance to CMV.

Sex-specific effects of cytotoxic chemotherapy agents cyclophosphamide and mitomycin C on gene expression, oxidative DNA damage, and epigenetic alterations in the prefrontal cortex and hippocampus - an aging connection

Recent research shows that chemotherapy agents can be more toxic to healthy brain cells than to the target cancer cells. They cause a range of side effects, including memory loss and cognitive dysfunction that can persist long after the completion of treatment. This condition is known as chemo brain. The molecular and cellular mechanisms of chemo brain remain obscure. Here, we analyzed the effects of two cytotoxic chemotherapy drugs—cyclophosphamide (CPP) and mitomycin C (MMC) - on transcriptomic and epigenetic changes in the murine prefrontal cortex (PFC) and hippocampal regions. We for the first time showed that CPP and MMC treatments led to profound sex- and brain region-specific alterations in gene expression profiles. Gene expression changes were most prominent in the PFC tissues of female mice 3 weeks after MMC treatment, and the gene expression response was much greater for MCC than CPP exposure. MMC exposure resulted in oxidative DNA damage, evidenced by accumulation of 8-oxo-2′-deoxyguanosine (8-oxodG) and a decrease in the level of 8-oxodG repair protein OGG1 in the PFC of female animals 3 weeks after treatment. MMC treatment decreased global DNA methylation and increased DNA hydroxymethylation in the PFC tissues of female mice. The majority of the changes induced by chemotherapy in the PFC tissues of female mice resembled those that occur during the brain's aging processes. Therefore, our study suggests a link between chemotherapy-induced chemo brain and brain aging, and provides an important roadmap for future analysis.

Abnormal circadian oscillation of hippocampal MAPK activity and power spectrums in NF1 mutant mice

Studies have implied that the circadian oscillation of mitogen-activated protein kinase (MAPK) signal pathways is crucial for hippocampus-dependent memory. NF1 mouse models (Nf1 heterozygous null mutants; Nf1 +/-) displayed enhanced MAPK activity in the hippocampus and resulted in memory deficits. We assumed a link between MAPK pathways and hippocampal rhythmic oscillations, which have never been explored in Nf1 +/- mice. We demonstrated that the level of extracellular signal-regulated kinases 1 and 2 (ERK1/2) phosphorylation in Nf1 +/- mice were significantly higher at nighttime than at daytime. Moreover, the in vivo recording revealed that for the Nf1 +/- group, the power spectral density of theta rhythm significantly decreased and the firing rates of pyramidal neurons increased. Our results indicated that the hippocampal MAPK oscillation and theta rhythmic oscillations in Nf1 +/- mice were disturbed and hinted about a possible mechanism for the brain dysfunction in Nf1 +/- mice.

Development of a SNP array and its application to genetic mapping and diversity assessment in pepper (Capsicum spp.)

The development and application of single nucleotide polymorphisms (SNPs) is in its infancy for pepper. Here, a set of 15,000 SNPs were chosen from the resequencing data to develop an array for pepper with 12,720 loci being ultimately synthesized. Of these, 8,199 (~64.46%) SNPs were found to be scorable and covered ~81.18% of the whole genome. With this array, a high-density interspecific genetic map with 5,569 SNPs was constructed using 297 F2 individuals, and genetic diversity of a panel of 399 pepper elite/landrace lines was successfully characterized. Based on the genetic map, one major QTL, named Up12.1, was detected for the fruit orientation trait. A total of 65 protein-coding genes were predicted within this QTL region based on the current annotation of the Zunla-1 genome. In summary, the thousands of well-validated SNP markers, high-density genetic map and genetic diversity information will be useful for molecular genetics and innovative breeding in pepper. Furthermore, the mapping results lay foundation for isolating the genes underlying variation in fruit orientation of Capsicum.

A comprehensive characterization of simple sequence repeats in pepper genomes provides valuable resources for marker development in Capsicum

The sequences of the full set of pepper genomes including nuclear, mitochondrial and chloroplast are now available for use. However, the overall of simple sequence repeats (SSR) distribution in these genomes and their practical implications for molecular marker development in Capsicum have not yet been described. Here, an average of 868,047.50, 45.50 and 30.00 SSR loci were identified in the nuclear, mitochondrial and chloroplast genomes of pepper, respectively. Subsequently, systematic comparisons of various species, genome types, motif lengths, repeat numbers and classified types were executed and discussed. In addition, a local database composed of 113,500 in silico unique SSR primer pairs was built using a homemade bioinformatics workflow. As a pilot study, 65 polymorphic markers were validated among a wide collection of 21 Capsicum genotypes with allele number and polymorphic information content value per marker raging from 2 to 6 and 0.05 to 0.64, respectively. Finally, a comparison of the clustering results with those of a previous study indicated the usability of the newly developed SSR markers. In summary, this first report on the comprehensive characterization of SSR motifs in pepper genomes and the very large set of SSR primer pairs will benefit various genetic studies in Capsicum.

New Insights on Eggplant/Tomato/Pepper Synteny and Identification of Eggplant and Pepper Orthologous QTL

Eggplant, pepper, and tomato are the most exploited berry-producing vegetables within the Solanaceae family. Their genomes differ in size, but each has 12 chromosomes which have undergone rearrangements causing a redistribution of loci. The genome sequences of all three species are available but differ in coverage, assembly quality and percentage of anchorage. Determining their syntenic relationship and QTL orthology will contribute to exploit genomic resources and genetic data for key agronomic traits. The syntenic analysis between tomato and pepper based on the alignment of 34,727 tomato CDS to the pepper genome sequence, identified 19,734 unique hits. The resulting synteny map confirmed the 14 inversions and 10 translocations previously documented, but also highlighted 3 new translocations and 4 major new inversions. Furthermore, each of the 12 chromosomes exhibited a number of rearrangements involving small regions of 0.5-0.7 Mbp. Due to high fragmentation of the publicly available eggplant genome sequence, physical localization of most eggplant QTL was not possible, thus, we compared the organization of the eggplant genetic map with the genome sequence of both tomato and pepper. The eggplant/tomato syntenic map confirmed all the 10 translocations but only 9 of the 14 known inversions; on the other hand, a newly detected inversion was recognized while another one was not confirmed. The eggplant/pepper syntenic map confirmed 10 translocations and 8 inversions already detected and suggested a putative new translocation. In order to perform the assessment of eggplant and pepper QTL orthology, the eggplant and pepper sequence-based markers located in their respective genetic map were aligned onto the pepper genome. GBrowse in pepper was used as reference platform for QTL positioning. A set of 151 pepper QTL were located as well as 212 eggplant QTL, including 76 major QTL (PVE ≥ 10%) affecting key agronomic traits. Most were confirmed to cluster in orthologous chromosomal regions. Our results highlight that the availability of genome sequences for an increasing number of crop species and the development of "ultra-dense" physical maps provide new and key tools for detailed syntenic and orthology studies between related plant species.

Development of INDEL Markers for Genetic Mapping Based on Whole Genome Resequencing in Soybean

Soybean [Glycine max (L.) Merrill] is an important crop worldwide. In this study, a Chinese local soybean cultivar, Hedou 12, was resequenced by next generation sequencing technology to develop INsertion/DELetion (INDEL) markers for genetic mapping. 49,276 INDEL polymorphisms and 242,059 single nucleotide polymorphisms were detected between Hedou 12 and the Williams 82 reference sequence. Of these, 243 candidate INDEL markers ranging from 5-50 bp in length were chosen for validation, and 165 (68%) of them revealed polymorphisms between Hedou 12 and Williams 82. The validated INDEL markers were also tested in 12 other soybean cultivars. The number of polymorphisms in the pairwise comparisons of 14 soybean cultivars varied from 27 to 165. To test the utility of these INDEL markers, they were used to perform genetic mapping of a crinkly leaf mutant, and the CRINKLY LEAF locus was successfully mapped to a 360 kb region on chromosome 7. This research shows that high-throughput sequencing technologies can facilitate the development of genome-wide molecular markers for genetic mapping in soybean.

Construction of an interspecific genetic map based on InDel and SSR for mapping the QTLs affecting the initiation of flower primordia in pepper (Capsicum spp.)

Re-sequencing permits the mining of genome-wide variations on a large scale and provides excellent resources for the research community. To accelerate the development and application of molecular markers and identify the QTLs affecting the flowering time-related trait in pepper, a total of 1,038 pairs of InDel and 674 SSR primers from different sources were used for genetic mapping using the F2 population (n = 154) derived from a cross between BA3 (C. annuum) and YNXML (C. frutescens). Of these, a total of 224 simple PCR-based markers, including 129 InDels and 95 SSRs, were validated and integrated into a map, which was designated as the BY map. The BY map consisted of 13 linkage groups (LGs) and spanned a total genetic distance of 1,249.77 cM with an average marker distance of 5.60 cM. Comparative analysis of the genetic and physical map based on the anchored markers showed that the BY map covered nearly the whole pepper genome. Based on the BY map, one major and five minor QTLs affecting the number of leaves on the primary axis (Nle) were detected on chromosomes P2, P7, P10 and P11 in 2012. The major QTL on P2 was confirmed based on another subset of the same F2 population (n = 147) in 2014 with selective genotyping of markers from the BY map. With the accomplishment of pepper whole genome sequencing and annotations (release 2.0), 153 candidate genes were predicted to embed in the Nle2.2 region, of which 12 important flowering related genes were obtained. The InDel/SSR-based interspecific genetic map, QTLs and candidate genes obtained by the present study will be useful for the downstream isolation of flowering time-related gene and other genetic applications for pepper.