Transcriptome Analysis of Differentially Expressed Genes Associated with Salt Stress in Cowpea (Vigna unguiculata L.) during the Early Vegetative Stage.

Cowpea (Vigna unguiculata (L.), 2n = 22) is a tropical crop grown in arid and semiarid regions that is tolerant to abiotic stresses such as heat and drought. However, in these regions, salt in the soil is generally not eluted by rainwater, leading to salt stress for a variety of plant species. This study was conducted to identify genes related to salt stress using the comparative transcriptome analysis of cowpea germplasms with contrasting salt tolerance. Using the Illumina Novaseq 6000 platform, 1.1 billion high-quality short reads, with a total length of over 98.6 billion bp, were obtained from four cowpea germplasms. Of the differentially expressed genes identified for each salt tolerance type following RNA sequencing, 27 were shown to exhibit significant expression levels. These candidate genes were subsequently narrowed down using reference-sequencing analysis, and two salt stress-related genes (Vigun_02G076100 and Vigun_08G125100) with single-nucleotide polymorphism (SNP) variation were selected. Of the five SNPs identified in Vigun_02G076100, one that caused significant amino acid variation was identified, while all nucleotide variations in Vigun_08G125100 was classified as missing in the salt-resistant germplasms. The candidate genes and their variation, identified in this study provide, useful information for the development of molecular markers for cowpea breeding programs.

Rice protein-binding microarrays: a tool to detect cis-acting elements near promoter regions in rice.

MAIN CONCLUSION: The present study showed that a rice (Oryza sativa)-specific protein-binding microarray (RPBM) can be applied to analyze DNA-binding motifs with a TF where binding is evaluated in extended natural promoter regions. The analysis may facilitate identifying TFs and their downstream genes and constructing gene networks through cis-elements. Transcription factors (TFs) regulate gene expression at the transcriptional level by binding a specific DNA sequence. Thus, predicting the DNA-binding motifs of TFs is one of the most important areas in the functional analysis of TFs in the postgenomic era. Although many methods have been developed to address this challenge, many TFs still have unknown DNA-binding motifs. In this study, we designed RPBM with 40-bp probes and 20-bp of overlap, yielding 49 probes spanning the 1-kb upstream region before the translation start site of each gene in the entire genome. To confirm the efficiency of RPBM technology, we selected two previously studied TFs, OsWOX13 and OsSMF1, and an uncharacterized TF, OsWRKY34. We identified the ATTGATTG and CCACGTCA DNA-binding sequences of OsWOX13 and OsSMF1, respectively. In total, 635 and 932 putative feature genes were identified for OsWOX13 and OsSMF1, respectively. We discovered the CGTTGACTTT DNA-binding sequence and 195 putative feature genes of OsWRKY34. RPBM could be applicable in the analysis of DNA-binding motifs for TFs where binding is evaluated in the promoter and 5' upstream CDS regions. The analysis may facilitate identifying TFs and their downstream genes and constructing gene networks through cis-elements.

Genome-wide analysis of epigenetic and transcriptional changes associated with heterosis in pigeonpea.

Hybrids are extensively used in agriculture to deliver an increase in yield, yet the molecular basis of heterosis is not well understood. Global DNA methylation analysis, transcriptome analysis and small RNA profiling were aimed to understand the epigenetic effect of the changes in gene expression level in the two hybrids and their parental lines. Increased DNA methylation was observed in both the hybrids as compared to their parents. This increased DNA methylation in hybrids showed that majority of the 24-nt siRNA clusters had higher expression in hybrids than the parents. Transcriptome analysis revealed that various phytohormones (auxin and salicylic acid) responsive hybrid-MPV DEGs were significantly altered in both the hybrids in comparison to MPV. DEGs associated with plant immunity and growth were overexpressed whereas DEGs associated with basal defence level were repressed. This antagonistic patterns of gene expression might contribute to the greater growth of the hybrids. It was also noticed that some common as well as unique changes in the regulatory pathways were associated with heterotic growth in both the hybrids. Approximately 70% and 67% of down-regulated hybrid-MPV DEGs were found to be differentially methylated in ICPH 2671 and ICPH 2740 hybrid, respectively. This reflected the association of epigenetic regulation in altered gene expressions. Our findings also revealed that miRNAs might play important roles in hybrid vigour in both the hybrids by regulating their target genes, especially in controlling plant growth and development, defence and stress response pathways. The above finding provides an insight into the molecular mechanism of pigeonpea heterosis.

QTL Analysis of Resistance to High-Intensity UV-B Irradiation in Soybean (Glycine max [L.] Merr.).

High-intensity ultraviolet-B (UV-B) irradiation is a complex abiotic stressor resulting in excessive light exposure, heat, and dehydration, thereby affecting crop yields. In the present study, we identified quantitative trait loci (QTLs) for resistance to high-intensity UV-B irradiation in soybean (Glycine max [L.]). We used a genotyping-by-sequencing approach using an F6 recombinant inbred line (RIL) population derived from a cross between Cheongja 3 (UV-B sensitive) and Buseok (UV-B resistant). We evaluated the degree of leaf damage by high-intensity UV-B radiation in the RIL population and identified four QTLs, UVBR12-1, 6-1, 10-1, and 14-1, for UV-B stress resistance, together explaining 20% of the observed phenotypic variation. The genomic regions containing UVBR12-1 and UVBR6-1 and their syntenic blocks included other known biotic and abiotic stress-related QTLs. The QTL with the highest logarithm of odds (LOD) score of 3.76 was UVBR12-1 on Chromosome 12, containing two genes encoding spectrin beta chain, brain (SPTBN, Glyma.12g088600) and bZIP transcription factor21/TGACG motif-binding 9 (bZIP TF21/TGA9, Glyma.12g088700). Their amino acid sequences did not differ between the mapping parents, but both genes were significantly upregulated by UV-B stress in Buseok but not in Cheongja 3. Among five genes in UVBR6-1 on Chromosome 6, Glyma.06g319700 (encoding a leucine-rich repeat family protein) had two nonsynonymous single nucleotide polymorphisms differentiating the parental lines. Our findings offer powerful genetic resources for efficient and precise breeding programs aimed at developing resistant soybean cultivars to multiple stresses. Furthermore, functional validation of the candidate genes will improve our understanding of UV-B stress defense mechanisms.

Genetic and epigenetic divergence of duplicate genes in two legume species.

Soybean (Glycine max) and common bean (Phaseolus vulgaris) share a polyploidy event ~59 MYA, followed by a Glycine-specific whole genome duplication (WGD) ~8-13 MYA. Duplicated genes were classified into five categories: singletons, dispersed, proximal, tandem, or WGD/segmental and found strong correlations between gene category and functional annotation. Photosynthesis and transcriptional regulation-related Gene Ontology terms were significantly over-represented in singletons and WGD genes, respectively, aligning with the gene balance hypothesis. We found that the divergence of gene expression and DNA methylation between WGD-derived paralogs increased with age and that WGD genes, initially retained via dosage constraints, subsequently underwent expression divergence, associated with other factors such as DNA methylation. Genes derived from different modes of duplication differed in breadth, level, and specificity of expression in both species. Orthologous genes and ungrouped genes (genes not in an ortholog group) differed in expression patterns. The protein divergence rates of WGD paralog pairs containing an ungrouped gene were higher than those for which both copies had orthologs. We propose that many ungrouped genes are derived from divergent and redundant gene copies, concordant with the neofunctionalization hypothesis. Tandemly duplicated genes were distinct from WGD-derived genes, indicating that mode of duplication contributes to the evolutionary fate of duplicated genes.

Divergent cytosine DNA methylation patterns in single-cell, soybean root hairs.

Chromatin modifications, such as cytosine methylation of DNA, play a significant role in mediating gene expression in plants, which affects growth, development, and cell differentiation. As root hairs are single-cell extensions of the root epidermis and the primary organs for water uptake and nutrients, we sought to use root hairs as a single-cell model system to measure the impact of environmental stress. We measured changes in cytosine DNA methylation in single-cell root hairs as compared with multicellular stripped roots, as well as in response to heat stress. Differentially methylated regions (DMRs) in each methylation context showed very distinct methylation patterns between cell types and in response to heat stress. Intriguingly, at normal temperature, root hairs were more hypermethylated than were stripped roots. However, in response to heat stress, both root hairs and stripped roots showed hypomethylation in each context, especially in the CHH context. Moreover, expression analysis of mRNA from similar tissues and treatments identified some associations between DMRs, genes and transposons. Taken together, the data indicate that changes in DNA methylation are directly or indirectly associated with expression of genes and transposons within the context of either specific tissues/cells or stress (heat).

A new approach for annotation of transposable elements using small RNA mapping.

Transposable elements (TEs) are mobile genomic DNA sequences found in most organisms. They so densely populate the genomes of many eukaryotic species that they are often the major constituents. With the rapid generation of many plant genome sequencing projects over the past few decades, there is an urgent need for improved TE annotation as a prerequisite for genome-wide studies. Analogous to the use of RNA-seq for gene annotation, we propose a new method for de novo TE annotation that uses as a guide 24 nt-siRNAs that are a part of TE silencing pathways. We use this new approach, called TASR (for Transposon Annotation using Small RNAs), for de novo annotation of TEs in Arabidopsis, rice and soybean and demonstrate that this strategy can be successfully applied for de novo TE annotation in plants.Executable PERL is available for download from: http://tasr-pipeline.sourceforge.net/.

A Comparative Epigenomic Analysis of Polyploidy-Derived Genes in Soybean and Common Bean.

Soybean (Glycine max) and common bean (Phaseolus vulgaris) share a paleopolyploidy (whole-genome duplication [WGD]) event, approximately 56.5 million years ago, followed by a genus Glycine-specific polyploidy, approximately 10 million years ago. Cytosine methylation is an epigenetic mark that plays an important role in the regulation of genes and transposable elements (TEs); however, the role of DNA methylation in the fate/evolution of genes following polyploidy and speciation has not been fully explored. Whole-genome bisulfite sequencing was used to produce nucleotide resolution methylomes for soybean and common bean. We found that, in soybean, CG body-methylated genes were abundant in WGD genes, which were, on average, more highly expressed than single-copy genes and had slower evolutionary rates than unmethylated genes, suggesting that WGD genes evolve more slowly than single-copy genes. CG body-methylated genes were also enriched in shared single-copy genes (single copy in both species) that may be responsible for the broad and high expression patterns of this class of genes. In addition, diverged methylation patterns in non-CG contexts between paralogs were due mostly to TEs in or near genes, suggesting a role for TEs and non-CG methylation in regulating gene expression post polyploidy. Reference methylomes for both soybean and common bean were constructed, providing resources for investigating epigenetic variation in legume crops. Also, the analysis of methylation patterns of duplicated and single-copy genes has provided insights into the functional consequences of polyploidy and epigenetic regulation in plant genomes.

Transcriptomic changes due to water deficit define a general soybean response and accession-specific pathways for drought avoidance.

BACKGROUND: Among abiotic stresses, drought is the most common reducer of crop yields. The slow-wilting soybean genotype PI 416937 is somewhat robust to water deficit and has been used previously to map the trait in a bi-parental population. Since drought stress response is a complex biological process, whole genome transcriptome analysis was performed to obtain a deeper understanding of the drought response in soybean. RESULTS: Contrasting data from PI 416937 and the cultivar 'Benning', we developed a classification system to identify genes that were either responding to water-deficit in both genotypes or that had a genotype x environment (GxE) response. In spite of very different wilting phenotypes, 90% of classifiable genes had either constant expression in both genotypes (33%) or very similar response profiles (E genes, 57%). By further classifying E genes based on expression profiles, we were able to discern the functional specificity of transcriptional responses at particular stages of water-deficit, noting both the well-known reduction in photosynthesis genes as well as the less understood up-regulation of the protein transport pathway. Two percent of classifiable genes had a well-defined GxE response, many of which are located within slow-wilting QTLs. We consider these strong candidates for possible causal genes underlying PI 416937's unique drought avoidance strategy. CONCLUSIONS: There is a general and functionally significant transcriptional response to water deficit that involves not only known pathways, such as down-regulation of photosynthesis, but also up-regulation of protein transport and chromatin remodeling. Genes that show a genotypic difference are more likely to show an environmental response than genes that are constant between genotypes. In this study, at least five genes that clearly exhibited a genotype x environment response fell within known QTL and are very good candidates for further research into slow-wilting.

Chromosome-level genome assembly of milk thistle (Silybum marianum (L.) Gaertn.).

Silybum marianum (L.) Gaertn., commonly known as milk thistle, is a medicinal plant belonging to the Asteraceae family. This plant has been recognized for its medicinal properties for over 2,000 years. However, the genome of this plant remains largely undiscovered, having no reference genome at a chromosomal level. Here, we assembled the chromosome-level genome of S. marianum, allowing for the annotation of 53,552 genes and the identification of transposable elements comprising 58% of the genome. The genome assembly from this study showed 99.1% completeness as determined by BUSCO assessment, while the previous assembly (ASM154182v1) showed 36.7%. Functional annotation of the predicted genes showed 50,329 genes (94% of total genes) with known protein functions in public databases. Comparative genome analysis among Asteraceae plants revealed a striking conservation of collinearity between S. marianum and C. cardunculus. The genomic information generated from this study will be a valuable resource for milk thistle breeding and for use by the larger research community.

Vertically Stackable Ovonic Threshold Switch Oscillator Using Atomic Layer Deposited Ge0.6Se0.4 Film for High-Density Artificial Neural Networks.

Nanodevice oscillators (nano-oscillators) have received considerable attention to implement in neuromorphic computing as hardware because they can significantly improve the device integration density and energy efficiency compared to complementary metal oxide semiconductor circuit-based oscillators. This work demonstrates vertically stackable nano-oscillators using an ovonic threshold switch (OTS) for high-density neuromorphic hardware. A vertically stackable Ge0.6Se0.4 OTS-oscillator (VOTS-OSC) is fabricated with a vertical crossbar array structure by growing Ge0.6Se0.4 film conformally on a contact hole structure using atomic layer deposition. The VOTS-OSC can be vertically integrated onto peripheral circuits without causing thermal damage because the fabrication temperature is <400 °C. The fabricated device exhibits oscillation characteristics, which can serve as leaky integrate-and-fire neurons in spiking neural networks (SNNs) and coupled oscillators in oscillatory neural networks (ONNs). For practical applications, pattern recognition and vertex coloring are demonstrated with SNNs and ONNs, respectively, using semiempirical simulations. This structure increases the oscillator integration density significantly, enabling complex tasks with a large number of oscillators. Moreover, it can enhance the computational speed of neural networks due to its rapid switching speed.

Whole-genome sequencing and intensive analysis of the undomesticated soybean (Glycine soja Sieb. and Zucc.) genome.

The genome of soybean (Glycine max), a commercially important crop, has recently been sequenced and is one of six crop species to have been sequenced. Here we report the genome sequence of G. soja, the undomesticated ancestor of G. max (in particular, G. soja var. IT182932). The 48.8-Gb Illumina Genome Analyzer (Illumina-GA) short DNA reads were aligned to the G. max reference genome and a consensus was determined for G. soja. This consensus sequence spanned 915.4 Mb, representing a coverage of 97.65% of the G. max published genome sequence and an average mapping depth of 43-fold. The nucleotide sequence of the G. soja genome, which contains 2.5 Mb of substituted bases and 406 kb of small insertions/deletions relative to G. max, is ∼0.31% different from that of G. max. In addition to the mapped 915.4-Mb consensus sequence, 32.4 Mb of large deletions and 8.3 Mb of novel sequence contigs in the G. soja genome were also detected. Nucleotide variants of G. soja versus G. max confirmed by Roche Genome Sequencer FLX sequencing showed a 99.99% concordance in single-nucleotide polymorphism and a 98.82% agreement in insertion/deletion calls on Illumina-GA reads. Data presented in this study suggest that the G. soja/G. max complex may be at least 0.27 million y old, appearing before the relatively recent event of domestication (6,000∼9,000 y ago). This suggests that soybean domestication is complicated and that more in-depth study of population genetics is needed. In any case, genome comparison of domesticated and undomesticated forms of soybean can facilitate its improvement.

A nationwide epidemiological study of nocturnal enuresis in Korean adolescents and adults: population based cross sectional study.

We performed a nationwide epidemiological study to evaluate the prevalence and characteristics of nocturnal enuresis (NE) in Korean adolescents and adults. A questionnaire was sent via e-mail to 51,073 people aged 16-40 yr by stratified sampling according to age, sex, and region among a 200,000 internet survey panel pool. The questionnaire included following information; presence or absence of NE, frequency of NE, possible risk factors for NE, self-esteem scale score and depression score results, and measures for the treatment of NE. Among the 2,117 responders, 54 (2.6%) had NE (≥1 enuretic episode within 6 months). Of 54 bedwetters, 9.3% wet ≥1 night per week and 20.5% wet ≥1 per month. The prevalence rates remained relatively stable with no apparent trend of reduction with age. The presence of sleep disturbance, family history, urgency, or urge incontinence increased the probability of NE episode significantly. The self-esteem score was lower (P=0.053) and the depression scale score was higher (P=0.003) in bedwetters compared with non-bedwetters. Overall 2.6% of Korean aged 16-40 yr have NE. The higher rate of urgency and urge incontinence in adolescent and adult enuretics suggests that bladder function has an important role in adolescent and adult NE.

Impact of operation voltage and NH3 annealing on the fatigue characteristics of ferroelectric AlScN thin films grown by sputtering.

This work investigates the impact of the magnitude of cycling voltage on the fatigue characteristics of 40 nm-thick AlScN ferroelectric thin film. The fatigue rate and the rejuvenation of remanent polarization vary with the cycling voltage. The primary fatigue mechanism is identified to be the interfacial layer formation and domain wall pinning at high and low cycling voltages, respectively. Additionally, annealing the film under the NH3 atmosphere decreases the fatigue rate and improves endurance by eliminating impurities in the film. The amount of trapped charges at the interface also decreases after NH3 annealing, leading to a reduction in leakage current. Furthermore, the ferroelectric performance of the AlScN film is not degraded after the thermal annealing at 900 °C under the NH3 environment, suggesting its robustness against the severe thermal budget. It is concluded that NH3 annealing is a promising method to address the reliability issue of the AlScN film.

Molecular sequence variations of the lipoxygenase-2 gene in soybean.

Soybean lipoxygenase genes comprise a multi-gene family, with the seed lipoxygenase isozymes LOX1, LOX2, and LOX3 present in soybean seeds. Among these, the LOX2 isozyme is primarily responsible for the "beany" flavor of most soybean seeds. The variety, Jinpumkong 2, having null alleles (lx1, lx2, and lx3) lacks the three seed lipoxygenases; so, sequence variations between the lipoxygenase-2 genes of Pureunkong (Lx2) and Jinpumkong 2 (lx2) cultivars were examined. One indel, four single nucleotide polymorphisms (SNPs), a 175-bp fragment in the 5'-flanking sequence, and a missense mutation within the coding region were found in Jinpumkong 2. The distribution of the sequence variations was investigated among 90 recombinant inbred lines (RILs) derived from a cross of Pureunkong × Jinpumkong 2 and in 480 germplasm accessions with various origins and maturity groups. Evidence for a genetic bottleneck was observed: the 175-bp fragment was rare in Glycine max, but present in the majority of the G. soja accessions. Furthermore, the 175-bp fragment was not detected in the 5' upstream region of the Lx2 gene on chromosome (Chr) 13 in Williams 82; instead, a similar 175-bp fragment was positioned in the homeologous region on Chr 15. The findings indicated that the novel fragment identified was originally present in the Lx2 region prior to the recent genome duplication in soybean, but became rare in the G. max gene pool. The missense mutation of the conserved histidine residue of the lx2 allele was developed into a single nucleotide-amplified polymorphism (SNAP) marker. The missense mutation showed a perfect correlation with the LOX2-lacking phenotype, so the SNAP marker is expected to facilitate breeding of soybean cultivars which lack the LOX2 isozyme.

Evaluation of maximum voided volume in Korean children by use of a 48-h frequency volume chart.

UNLABELLED: Study Type - Diagnostic (validating cohort). Level of Evidence 2a. What's known on the subject? and What does the study add? The relationship between the maximum voided volume followed a linear curve. The formula presented, bladder capacity (mL) = 12 ×[age (years) + 11], is thought to be a reasonable one for Korean children. Korean children have a smaller bladder capacity than that reported in previous Western studies. OBJECTIVE: • To develop practical guidelines for the prediction of normal bladder capacity in Korean children measured by a frequency volume chart (FVC), maximum voided volume (MVV) is an important factor in the diagnosis of children with abnormal voiding function. SUBJECTS AND METHODS: • In all, 298 children, aged 3-13 years, with no history of voiding disorders volunteered for the study. The MVV was determined in 219 subjects by use of a completely recorded FVC. RESULTS: • Linear regression analysis was used to define the exact relationship between age and bladder capacity. An approximate formula related age to bladder capacity as follows: bladder capacity (mL) = 12 ×[age (years) + 11]. CONCLUSIONS: • The relationship between the MVV measured by a FVC by age (3-13 years) of Korean children followed a linear curve. • When applied to normal voiding patterns, the formula presented appears to be a reasonable one for Korean children.

Sequence variants of Toll-like receptor 4 (TLR4) and the risk of prostate cancer in Korean men.

PURPOSE: Chronic inflammation has been considered a potential risk factor for prostate cancer. Toll-like receptors (TLRs) are important in the innate immune response to pathogens and in cross talk between innate immunity and adaptive immunity. In this study, sequence variants in the TLR4 gene were investigated to determine whether they were associated with prostate cancer risk in a Korean cohort. METHODS: An association study between 11 single-nucleotide polymorphisms (SNPs) of the TLR4 gene and prostate cancer was performed in 463 Korean male subjects including 240 prostate cancer patients and 223 healthy controls. SNPs were genotyped using the TaqMan assay, and their association with the risk of prostate cancer was evaluated using logistic regression analysis. RESULTS: The statistical analysis revealed that one SNP at the 3'UTR (rs11536889) showed significant association with the risk of prostate cancer (P (corr) = 0.005, OR = 1.81). One common haplotype (ht2) was also significantly associated with the risk of prostate cancer (P (corr) = 0.009, OR = 1.77). However, further analysis showed no association between any of the SNPs and prostate cancer prognostic factors such as the Gleason score or tumor stage. CONCLUSIONS: The findings of this study suggest that polymorphisms of the TLR4 gene might be associated with the risk of prostate cancer in Korean men.

Harnessing epigenetic variability for crop improvement: current status and future prospects.

BACKGROUND: The epigenetic mechanisms play critical roles in a vast diversity of biological processes of plants, including development and response to environmental challenges. Particularly, DNA methylation is a stable epigenetic signature that supplements the genetics-based view of complex life phenomena. In crop breeding, the decrease in genetic diversity due to artificial selection of conventional breeding methods has been a long-standing concern. Therefore, the epigenetic diversity has been proposed as a new resource for future crop breeding, which will be hereinafter referred to as epibreeding. DISCUSSION: The induction of methylome changes has been performed in plants by several methods including chemical drugs treatment and tissue culture. Target-specific epigenetic engineering has been also attempted by exogenous RNAi mediated by virus-induced gene silencing and grafting. Importantly, the new and innovative techniques including the CRISPR-Cas9 system have recently been adopted in epigenetic engineering of plant genomes, facilitating the efforts for epibreeding. CONCLUSION: In this review, we introduce several examples of natural and induced epigenetic changes impacting on agronomic traits and discuss the methods for generating epigenomic diversity and site-specific epigenetic engineering.

A Cold-Shock Protein from the South Pole-Dwelling Soil Bacterium Arthrobacter sp. Confers Cold Tolerance to Rice.

Low temperature is a critical environmental factor restricting the physiology of organisms across kingdoms. In prokaryotes, cold shock induces the expression of various genes and proteins involved in cellular processes. Here, a cold-shock protein (ArCspA) from the South Pole-dwelling soil bacterium Arthrobacter sp. A2-5 was introduced into rice, a monocot model plant species. Four-week-old 35S:ArCspA transgenic rice plants grown in a cold chamber at 4 °C survived for 6 days. Cold stress significantly decreased the chlorophyll content in WT plants after 4 days compared with that in 35S:ArCspA transgenic plants. RNA-seq analysis was performed on WT and 35S:ArCspA transgenic rice with/without cold stress. GO terms such as "response to stress (GO:0006950)", "response to cold (GO:0009409)", and "response to heat (GO:0009408)" were significantly enriched among the upregulated genes in the 35S:ArCspA transgenic rice under normal conditions, even without cold-stress treatment. The expression of five cold stress-related genes, Rab16B (Os11g0454200), Rab21 (Os11g0454300), LEA22 (Os01g0702500), ABI5 (Os01 g0859300), and MAPK5 (Os03g0285800), was significantly upregulated in the transgenic rice compared with the WT rice. These results indicate that the ArCspA gene might be involved in the induction of cold-responsive genes and provide cold tolerance.

Genome-Wide Association Study for Ultraviolet-B Resistance in Soybean (Glycine max L.).

The depletion of the stratospheric ozone layer is a major environmental issue and has increased the dosage of ultraviolet-B (UV-B) radiation reaching the Earth's surface. Organisms are negatively affected by enhanced UV-B radiation, and especially in crop plants this may lead to severe yield losses. Soybean (Glycine max L.), a major legume crop, is sensitive to UV-B radiation, and therefore, it is required to breed the UV-B-resistant soybean cultivar. In this study, 688 soybean germplasms were phenotyped for two categories, Damage of Leaf Chlorosis (DLC) and Damage of Leaf Shape (DLS), after supplementary UV-B irradiation for 14 days. About 5% of the germplasms showed strong UV-B resistance, and GCS731 was the most resistant genotype. Their phenotypic distributions showed similar patterns to the normal, suggesting UV-B resistance as a quantitative trait governed by polygenes. A total of 688 soybean germplasms were genotyped using the Axiom® Soya 180K SNP array, and a genome-wide association study (GWAS) was conducted to identify SNPs significantly associated with the two traits, DLC and DLS. Five peaks on chromosomes 2, 6, 10, and 11 were significantly associated with either DLC or DLS, and the five adjacent genes were selected as candidate genes responsible for UV-B resistance. Among those candidate genes, Glyma.02g017500 and Glyma.06g103200 encode cryptochrome (CRY) and cryptochrome 1 (CRY1), respectively, and are known to play a role in DNA repair during photoreactivation. Real-time quantitative RT-PCR (qRT-PCR) results revealed that CRY1 was expressed significantly higher in the UV-B-resistant soybean compared to the susceptible soybean after 6 h of UV-B irradiation. This study is the first GWAS report on UV-B resistance in soybean, and the results will provide valuable information for breeding UV-B-resistant soybeans in preparation for climate change.