Calcium-permeable AMPA receptors govern PV neuron feature selectivity.

The brain helps us survive by forming internal representations of the external world1,2. Excitatory cortical neurons are often precisely tuned to specific external stimuli3,4. However, inhibitory neurons, such as parvalbumin-positive (PV) interneurons, are generally less selective5. PV interneurons differ from excitatory neurons in their neurotransmitter receptor subtypes, including AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptors (AMPARs)6,7. Excitatory neurons express calcium-impermeable AMPARs that contain the GluA2 subunit (encoded by GRIA2), whereas PV interneurons express receptors that lack the GluA2 subunit and are calcium-permeable (CP-AMPARs). Here we demonstrate a causal relationship between CP-AMPAR expression and the low feature selectivity of PV interneurons. We find low expression stoichiometry of GRIA2 mRNA relative to other subunits in PV interneurons that is conserved across ferrets, rodents, marmosets and humans, and causes abundant CP-AMPAR expression. Replacing CP-AMPARs in PV interneurons with calcium-impermeable AMPARs increased their orientation selectivity in the visual cortex. Manipulations to induce sparse CP-AMPAR expression demonstrated that this increase was cell-autonomous and could occur with changes beyond development. Notably, excitatory-PV interneuron connectivity rates and unitary synaptic strength were unaltered by CP-AMPAR removal, which suggested that the selectivity of PV interneurons can be altered without markedly changing connectivity. In Gria2-knockout mice, in which all AMPARs are calcium-permeable, excitatory neurons showed significantly degraded orientation selectivity, which suggested that CP-AMPARs are sufficient to drive lower selectivity regardless of cell type. Moreover, hippocampal PV interneurons, which usually exhibit low spatial tuning, became more spatially selective after removing CP-AMPARs, which indicated that CP-AMPARs suppress the feature selectivity of PV interneurons independent of modality. These results reveal a new role of CP-AMPARs in maintaining low-selectivity sensory representation in PV interneurons and implicate a conserved molecular mechanism that distinguishes this cell type in the neocortex.

Massively parallel dissection of RNA in RNA-protein interactions in vivo.

Many of the biological functions performed by RNA are mediated by RNA-binding proteins (RBPs), and understanding the molecular basis of these interactions is fundamental to biology. Here, we present massively parallel RNA assay combined with immunoprecipitation (MPRNA-IP) for in vivo high-throughput dissection of RNA-protein interactions and describe statistical models for identifying RNA domains and parsing the structural contributions of RNA. By using custom pools of tens of thousands of RNA sequences containing systematically designed truncations and mutations, MPRNA-IP is able to identify RNA domains, sequences, and secondary structures necessary and sufficient for protein binding in a single experiment. We show that this approach is successful for multiple RNAs of interest, including the long noncoding RNA NORAD, bacteriophage MS2 RNA, and human telomerase RNA, and we use it to interrogate the hitherto unknown sequence or structural RNA-binding preferences of the DNA-looping factor CTCF. By integrating systematic mutation analysis with crosslinking immunoprecipitation, MPRNA-IP provides a novel high-throughput way to elucidate RNA-based mechanisms behind RNA-protein interactions in vivo.

Evolutionary divergence of Firre localization and expression.

Long noncoding RNAs (lncRNAs) are rapidly evolving and thus typically poorly conserved in their sequences. How these sequence differences affect the characteristics and potential functions of lncRNAs with shared synteny remains unclear. Here we show that the syntenically conserved lncRNA Firre displays distinct expression and localization patterns in human and mouse. Single molecule RNA FISH reveals that in a range of cell lines, mouse Firre (mFirre) is predominantly nuclear, while human FIRRE (hFIRRE) is distributed between the cytoplasm and nucleus. This localization pattern is maintained in human/mouse hybrid cells expressing both human and mouse Firre, implying that the localization of the lncRNA is species autonomous. We find that the majority of hFIRRE transcripts in the cytoplasm are comprised of isoforms that are enriched in RRD repeats. We furthermore determine that in various tissues, mFirre is more highly expressed than its human counterpart. Our data illustrate that the rapid evolution of syntenic lncRNAs can lead to variations in lncRNA localization and abundance, which in turn may result in disparate lncRNA functions even in closely related species.

Characterizing the nuclear and cytoplasmic transcriptomes in developing and mature human cortex uncovers new insight into psychiatric disease gene regulation.

Transcriptome compartmentalization by the nuclear membrane provides both stochastic and functional buffering of transcript activity in the cytoplasm, and has recently been implicated in neurodegenerative disease processes. Although many mechanisms regulating transcript compartmentalization are also prevalent in brain development, the extent to which subcellular localization differs as the brain matures has yet to be addressed. To characterize the nuclear and cytoplasmic transcriptomes during brain development, we sequenced both RNA fractions from homogenate prenatal and adult human postmortem cortex using poly(A)+ and Ribo-Zero library preparation methods. We find that while many genes are differentially expressed by fraction and developmental expression changes are similarly detectable in nuclear and cytoplasmic RNA, the compartmented transcriptomes become more distinct as the brain matures, perhaps reflecting increased utilization of nuclear retention as a regulatory strategy in adult brain. We examined potential mechanisms of this developmental divergence including alternative splicing, RNA editing, nuclear pore composition, RNA-binding protein motif enrichment, and RNA secondary structure. Intron retention is associated with greater nuclear abundance in a subset of transcripts, as is enrichment for several splicing factor binding motifs. Finally, we examined disease association with fraction-regulated gene sets and found nuclear-enriched genes were also preferentially enriched in gene sets associated with neurodevelopmental psychiatric disorders. These results suggest that although gene-level expression is globally comparable between fractions, nuclear retention of transcripts may play an underappreciated role in developmental regulation of gene expression in brain, particularly in genes whose dysregulation is related to neuropsychiatric disorders.

Genome-wide investigation of intragenic DNA methylation identifies ZMIZ1 gene as a prognostic marker in glioblastoma and multiple cancer types.

DNA methylation has long been recognized as a tumor-promoting factor when aberrantly regulated in the promoter region of genes. However, the effect of intragenic DNA methylation remains poorly understood on the clinical aspects of cancer. Here, we first evaluated the significance of intragenic DNA methylation for survival outcomes of cancer patients in a genome-wide manner. Glioblastoma patients with hypermethylated intragenic regions exhibited better survival than hypomethylated patients. Enrichment analyses of intragenic DNA methylation profiles with epigenetic signatures prioritized the intragenic DNA methylation of ZMIZ1 as a possible glioblastoma prognostic marker that is independent of MGMT methylation in IDH1 wild-type patients. This intragenic region harbored molecular signatures of alternative transcription across many cell types. Furthermore, we found that the intragenic region of ZMIZ1 can serve as a molecular marker in multiple cancers including astrocytomas, bladder cancer and renal cell carcinoma according to DNA methylation status. Finally, in vitro and in vivo experiments uncovered the role of ZMIZ1 as a driver of tumor cell migration. Altogether, our results identify ZMIZ1 as a prognostic marker in cancer and highlight the clinical significance of intragenic methylation in cancer.

Restoration of miR-29b exerts anti-cancer effects on glioblastoma.

BACKGROUND: Glioblastoma multiforme (GBM) is known as one of the most fatal forms of cancer. MicroRNAs have been widely implicated in the regulation of mammalian development and pathogenesis. The brain-enriched miR-29 subfamilies are known to be exclusively expressed in the developing brain, and they are aberrantly down-regulated in GBM. This study aims to elucidate the role of miR-29b in GBM development and the feasibility of therapeutic targeting using conjugated nanoparticles. METHODS: After confirmation of miR-29b expression levels in GBM tissues by analysis of open source data, the anticancer effect of miR-29b was tested by the introduction of syn-hsa-miR-29b-3p in the A172 GBM cell line. In vitro studies of cell viability and apoptosis and ex vivo study using GBM tissue slice cultures from 3 patients and nanoparticle delivery of miR-29b were performed. RESULTS: We discovered an increase in apoptotic cell populations with the introduction of miR-29b in the GBM cell line. An established human-derived GBM tissue slice culture system confirmed the anticancer effect of miR-29b-conjugated nanoparticles. Using PCR array, we found that exogenous miR-29b inhibits the expression of COL1A2, COL3A1, COL4A1, ELN, ITGA11, MMP24, and SPARC, which mediates an anticancer effect. CONCLUSIONS: miR-29b may serve as a putative therapeutic molecule when its expression is restored using a nanoparticle delivery system in GBM.

Variation block-based genomics method for crop plants.

BACKGROUND: In contrast with wild species, cultivated crop genomes consist of reshuffled recombination blocks, which occurred by crossing and selection processes. Accordingly, recombination block-based genomics analysis can be an effective approach for the screening of target loci for agricultural traits. RESULTS: We propose the variation block method, which is a three-step process for recombination block detection and comparison. The first step is to detect variations by comparing the short-read DNA sequences of the cultivar to the reference genome of the target crop. Next, sequence blocks with variation patterns are examined and defined. The boundaries between the variation-containing sequence blocks are regarded as recombination sites. All the assumed recombination sites in the cultivar set are used to split the genomes, and the resulting sequence regions are termed variation blocks. Finally, the genomes are compared using the variation blocks. The variation block method identified recurring recombination blocks accurately and successfully represented block-level diversities in the publicly available genomes of 31 soybean and 23 rice accessions. The practicality of this approach was demonstrated by the identification of a putative locus determining soybean hilum color. CONCLUSIONS: We suggest that the variation block method is an efficient genomics method for the recombination block-level comparison of crop genomes. We expect that this method will facilitate the development of crop genomics by bringing genomics technologies to the field of crop breeding.

A major and stable QTL associated with seed weight in soybean across multiple environments and genetic backgrounds.

KEY MESSAGE: We detected a QTL for single seed weight in soybean that was stable across multiple environments and genetic backgrounds with the use of two recombinant inbred line populations. Single seed weight (SSW) in soybean is a key determinant of both seed yield and the quality of soy food products, and it exhibits wide variation. SSW is under genetic control, but the molecular mechanisms of such control remain unclear. We have now investigated quantitative trait loci (QTLs) for SSW in soybean and have identified such a QTL that is stable across multiple environments and genetic backgrounds. Two populations of 225 and 250 recombinant inbred lines were developed from crosses between Japanese and US cultivars of soybean that differ in SSW by a factor of ~2, and these populations were grown in at least three different environments. A whole-genome panel comprising 304 simple sequence repeat (SSR) loci was applied to mapping in each population. We identified 15 significant QTLs for SSW dispersed among 11 chromosomes in the two populations. One QTL located between Sat_284 and Sat_292 on chromosome 17 was detected (3.6 < LOD < 14.1) in both populations grown in all environments. This QTL, tentatively designated qSw17-1, accounted for 9.4-20.9 % of phenotypic variation in SSW, with a dominant allele being associated with increased SSW. Given its substantial effect on SSW, qSw17-1 is an attractive target for positional cloning, and SSR markers closely associated with this locus may prove useful for marker-assisted selection for SSW control in soybean.

Cognitive subdomain responses to galantamine in Alzheimer's disease.

We investigated the effects of galantamine on cognitive subdomains in Alzheimer's disease (AD). Sixty-six patients with mild-to-moderate AD received open-label galantamine for 52 weeks. Cognitive function was measured using the Korean version of the Alzheimer's Disease Assessment Scale-cognitive subscale (ADAS-cog-K). Response to galantamine was defined as "improvement or no deterioration" on the total scores of the ADAS-cog-K at 26 weeks. In the overall intent-to-treat sample, we found less cognitive decline during 26 and 52 weeks than the expected untreated course as predicted by Stern's equation. The operationally defined response rate at 26 weeks was 66.7%. The responders differed significantly from the nonresponders only in the memory and language domains but not in the domains of praxis or frontal/executive function or in secondary outcome measures of neuropsychiatric symptoms and activities of daily living. The subdomain analysis revealed an effect of galantamine on preservation of memory that was not apparent in the overall analysis. Failure to achieve responder status by 26 weeks was associated with no further possibility of response.

Evaluation of the RNA-dependence of PRC2 binding to chromatin in human pluripotent stem cells.

Polycomb Repressive Complex 2 (PRC2), an important histone modifier and epigenetic repressor, has been known to interact with RNA for almost two decades. In our previous publication (Long, Hwang et al. 2020), we presented data supporting the functional importance of RNA interaction in maintaining PRC2 occupancy on chromatin, using comprehensive approaches including an RNA-binding mutant of PRC2 and an rChIP-seq assay. Recently, concerns have been expressed regarding whether the RNA-binding mutant has impaired histone methyltransferase activity and whether the rChIP-seq assay can potentially generate artifacts. Here we provide new data that support a number of our original findings. First, we found the RNA-binding mutant to be fully capable of maintaining H3K27me3 levels in human induced pluripotent stem cells. The mutant had reduced methyltransferase activity in vitro, but only on some substrates at early time points. Second, we found that our rChIP-seq method gave consistent data across antibodies and cell lines. Third, we further optimized rChIP-seq by using lower concentrations of RNase A and incorporating a catalytically inactive mutant RNase A as a control, as well as using an alternative RNase (RNase T1). The EZH2 rChIP-seq results using the optimized protocols supported our original finding that RNA interaction contributes to the chromatin occupancy of PRC2.

Identification of QTL controlling post-flowering period in soybean.

The length of the reproductive period affects the grain yield of soybean (Glycine max [L.] Merr), and genetic control of the period might contribute to yield improvement. To detect genetic factor(s) controlling the reproductive period, a population of recombinant inbred lines (RILs) was developed from a cross between Japanese landrace 'Ippon-Sangoh' and, Japanese cultivar 'Fukuyutaka' which differ in their duration from flowering to maturation (DFM) relative to the difference in the duration from sowing to flowering (DSF). In the RIL population, the DFM correlated poorly (r = -0.16 to 0.34) with the DSF in all field trials over 3 years. Two stable QTLs for the DFM on chromosomes (Chr-) 10 and 11 as well as two stable QTLs for the DSF on Chr-10 and -16 were identified. The QTL on Chr-11 for the reproductive period (designated as qDfm1; quantitative trait locus for duration from flowering to maturation 1) affected all three trials, and the difference in the DFM between the Fukuyutaka and Ippon-Sangoh was mainly accounted for qDfm1, in which the Fukuyutaka allele promoted a longer period. qDfm1 affected predominantly the reproductive period, and thus it might be possible to alter the period with little influence on the vegetative period.

Genome-wide perturbations of Alu expression and Alu-associated post-transcriptional regulations distinguish oligodendroglioma from other gliomas.

Alu is a primate-specific repeat element in the human genome and has been increasingly appreciated as a regulatory element in many biological processes. But the appreciation of Alu has been limited in tumorigenesis, especially for brain tumor. To investigate the relevance of Alu to the gliomagenesis, we studied Alu element-associated post-transcriptional processes and the RNA expression of the element by performing RNA-seq for a total of 41 pairs of neurotypical and diverse glioma brain tissues. We find that A-to-I editing and circular RNA levels, as well as Alu RNA expression, are decreased overall in gliomas, compared to normal tissue. Interestingly, grade 2 oligodendrogliomas are least affected in A-to-I editing and circular RNA levels among gliomas, whereas they have a higher proportion of down-regulated Alu subfamilies, compared to the other gliomas. These findings collectively imply a unique pattern of Alu-associated transcriptomes in grade 2 oligodendroglioma, providing an insight to gliomagenesis from the perspective of an evolutionary genetic element.

RNA is essential for PRC2 chromatin occupancy and function in human pluripotent stem cells.

Many chromatin-binding proteins and protein complexes that regulate transcription also bind RNA. One of these, Polycomb repressive complex 2 (PRC2), deposits the H3K27me3 mark of facultative heterochromatin and is required for stem cell differentiation. PRC2 binds RNAs broadly in vivo and in vitro. Yet, the biological importance of this RNA binding remains unsettled. Here, we tackle this question in human induced pluripotent stem cells by using multiple complementary approaches. Perturbation of RNA-PRC2 interaction by RNase A, by a chemical inhibitor of transcription or by an RNA-binding-defective mutant all disrupted PRC2 chromatin occupancy and localization genome wide. The physiological relevance of PRC2-RNA interactions is further underscored by a cardiomyocyte differentiation defect upon genetic disruption. We conclude that PRC2 requires RNA binding for chromatin localization in human pluripotent stem cells and in turn for defining cellular state.

The Sox2 transcription factor binds RNA.

Certain transcription factors are proposed to form functional interactions with RNA to facilitate proper regulation of gene expression. Sox2, a transcription factor critical for maintenance of pluripotency and neurogenesis, has been found associated with several lncRNAs, although it is unknown whether these interactions are direct or via other proteins. Here we demonstrate that human Sox2 interacts directly with one of these lncRNAs with high affinity through its HMG DNA-binding domain in vitro. These interactions are primarily with double-stranded RNA in a non-sequence specific fashion, mediated by a similar but not identical interaction surface. We further determined that Sox2 directly binds RNA in mouse embryonic stem cells by UV-cross-linked immunoprecipitation of Sox2 and more than a thousand Sox2-RNA interactions in vivo were identified using fRIP-seq. Together, these data reveal that Sox2 employs a high-affinity/low-specificity paradigm for RNA binding in vitro and in vivo.

The Tug1 lncRNA locus is essential for male fertility.

BACKGROUND: Several long noncoding RNAs (lncRNAs) have been shown to function as components of molecular machines that play fundamental roles in biology. While the number of annotated lncRNAs in mammalian genomes has greatly expanded, studying lncRNA function has been a challenge due to their diverse biological roles and because lncRNA loci can contain multiple molecular modes that may exert function. RESULTS: We previously generated and characterized a cohort of 20 lncRNA loci knockout mice. Here, we extend this initial study and provide a more detailed analysis of the highly conserved lncRNA locus, taurine-upregulated gene 1 (Tug1). We report that Tug1-knockout male mice are sterile with underlying defects including a low number of sperm and abnormal sperm morphology. Because lncRNA loci can contain multiple modes of action, we wanted to determine which, if any, potential elements contained in the Tug1 genomic region have any activity. Using engineered mouse models and cell-based assays, we provide evidence that the Tug1 locus harbors two distinct noncoding regulatory activities, as a cis-DNA repressor that regulates neighboring genes and as a lncRNA that can regulate genes by a trans-based function. We also show that Tug1 contains an evolutionary conserved open reading frame that when overexpressed produces a stable protein which impacts mitochondrial membrane potential, suggesting a potential third coding function. CONCLUSIONS: Our results reveal an essential role for the Tug1 locus in male fertility and uncover evidence for distinct molecular modes in the Tug1 locus, thus highlighting the complexity present at lncRNA loci.

Identification of differentially expressed subnetworks based on multivariate ANOVA.

BACKGROUND: Since high-throughput protein-protein interaction (PPI) data has recently become available for humans, there has been a growing interest in combining PPI data with other genome-wide data. In particular, the identification of phenotype-related PPI subnetworks using gene expression data has been of great concern. Successful integration for the identification of significant subnetworks requires the use of a search algorithm with a proper scoring method. Here we propose a multivariate analysis of variance (MANOVA)-based scoring method with a greedy search for identifying differentially expressed PPI subnetworks. RESULTS: Given the MANOVA-based scoring method, we performed a greedy search to identify the subnetworks with the maximum scores in the PPI network. Our approach was successfully applied to human microarray datasets. Each identified subnetwork was annotated with the Gene Ontology (GO) term, resulting in the phenotype-related functional pathway or complex. We also compared these results with those of other scoring methods such as t statistic- and mutual information-based scoring methods. The MANOVA-based method produced subnetworks with a larger number of proteins than the other methods. Furthermore, the subnetworks identified by the MANOVA-based method tended to consist of highly correlated proteins. CONCLUSION: This article proposes a MANOVA-based scoring method to combine PPI data with expression data using a greedy search. This method is recommended for the highly sensitive detection of large subnetworks.

Integrative analysis of DNA methylation suggests down-regulation of oncogenic pathways and reduced somatic mutation rates in survival outliers of glioblastoma.

The study of survival outliers of glioblastoma can provide important clues on gliomagenesis as well as on the ways to alter clinical course of this almost uniformly lethal cancer type. However, there has been little consensus on genetic and epigenetic signatures of the long-term survival outliers of glioblastoma. Although the two classical molecular markers of glioblastoma including isocitrate dehydrogenase 1 or 2 (IDH1/2) mutation and O6-methylguanine DNA methyltransferase (MGMT) promoter methylation are associated with overall survival rate of glioblastoma patients, they are not specific to the survival outliers. In this study, we compared the two groups of survival outliers of glioblastoma with IDH wild-type, consisting of the glioblastoma patients who lived longer than 3 years (n = 17) and the patients who lived less than 1 year (n = 12) in terms of genome-wide DNA methylation profile. Statistical analyses were performed to identify differentially methylated sites between the two groups. Functional implication of DNA methylation patterns specific to long-term survivors of glioblastoma were investigated by comprehensive enrichment analyses with genomic and epigenomic features. We found that the genome of long-term survivors of glioblastoma is differentially methylated relative to short-term survivor patients depending on CpG density: hypermethylation near CpG islands (CGIs) and hypomethylation far from CGIs. Interestingly, these two patterns are associated with distinct oncogenic aspects in gliomagenesis. In the long-term survival glioblastoma-specific sites distant from CGI, somatic mutations of glioblastoma are enriched with higher DNA methylation, suggesting that the hypomethylation in long-term survival glioblastoma can contribute to reduce the rate of somatic mutation. On the other hand, the hypermethylation near CGIs associates with transcriptional downregulation of genes involved in cancer progression pathways. Using independent cohorts of IDH1/2- wild type glioblastoma, we also showed that these two patterns of DNA methylation can be used as molecular markers of long-term survival glioblastoma. Our results provide extended understanding of DNA methylation, especially of DNA hypomethylation, in cancer genome and reveal clinical importance of DNA methylation pattern as prognostic markers of glioblastoma.

The Firre locus produces a trans-acting RNA molecule that functions in hematopoiesis.

RNA has been classically known to play central roles in biology, including maintaining telomeres, protein synthesis, and in sex chromosome compensation. While thousands of long noncoding RNAs (lncRNAs) have been identified, attributing RNA-based roles to lncRNA loci requires assessing whether phenotype(s) could be due to DNA regulatory elements, transcription, or the lncRNA. Here, we use the conserved X chromosome lncRNA locus Firre, as a model to discriminate between DNA- and RNA-mediated effects in vivo. We demonstrate that (i) Firre mutant mice have cell-specific hematopoietic phenotypes, and (ii) upon exposure to lipopolysaccharide, mice overexpressing Firre exhibit increased levels of pro-inflammatory cytokines and impaired survival. (iii) Deletion of Firre does not result in changes in local gene expression, but rather in changes on autosomes that can be rescued by expression of transgenic Firre RNA. Together, our results provide genetic evidence that the Firre locus produces a trans-acting lncRNA that has physiological roles in hematopoiesis.

Influence of heat and moisture treatment on carotenoids, phenolic content, and antioxidant capacity of orange maize flour.

The aim of this work was to study the effect of heat and moisture treatment (HMT) on carotenoids, phenolic content and antioxidant capacity of ground, orange maize. Total carotenoid content (TCC) of untreated sample (53.39 mg/kg) was 2.2 times higher than measured in treated orange maize f (24.61 mg/kg). The rates of degradation with HMT were in the following order: ß-carotene > ß-cryptoxanthin > zeaxanthin > lutein. There was a significant interaction between longer heating time and higher moisture content on carotenoid degradation (p < .05). Total phenolic content (TPC) in raw sample (1664.74 mg/kg) was two-fold higher than in treated orange maize (827.89 mg/kg). Ferulic acid was the most abundant and stable phenolic acid in raw and treated orange maize. The antioxidant capacity of orange maize was higher in methanol than in butanol extracts. The highest correlation (0.924) was observed between TPC and ABTS+ scavenging capacity of methanol extracts.

Barcode System for Genetic Identification of Soybean [Glycine max (L.) Merrill] Cultivars Using InDel Markers Specific to Dense Variation Blocks.

For genetic identification of soybean [Glycine max (L.) Merrill] cultivars, insertions/deletions (InDel) markers have been preferred currently because they are easy to use, co-dominant and relatively abundant. Despite their biological importance, the investigation of InDels with proven quality and reproducibility has been limited. In this study, we described soybean barcode system approach based on InDel makers, each of which is specific to a dense variation block (dVB) with non-random recombination due to many variations. Firstly, 2,274 VBs were mined by analyzing whole genome data in six soybean cultivars (Backun, Sinpaldal 2, Shingi, Daepoong, Hwangkeum, and Williams 82) for transferability to dVB-specific InDel markers. Secondly, 73,327 putative InDels in the dVB regions were identified for the development of soybean barcode system. Among them, 202 dVB-specific InDels from all soybean cultivars were selected by gel electrophoresis, which were converted as 2D barcode types according to comparing amplicon polymorphisms in the five cultivars to the reference cultivar. Finally, the polymorphism of the markers were assessed in 147 soybean cultivars, and the soybean barcode system that allows a clear distinction among soybean cultivars is also detailed. In addition, the changing of the dVBs in a chromosomal level can be quickly identified due to investigation of the reshuffling pattern of the soybean cultivars with 27 maker sets. Especially, a backcross-inbred offspring, "Singang" and a recurrent parent, "Sowon" were identified by using the 27 InDel markers. These results indicate that the soybean barcode system enables not only the minimal use of molecular markers but also comparing the data from different sources due to no need of exploiting allele binning in new varieties.