Arabidopsis Transcription Regulatory Factor Domain/Domain Interaction Analysis Tool-Liquid/Liquid Phase Separation, Oligomerization, GO Analysis: A Toolkit for Interaction Data-Based Domain Analysis.

Although a large number of databases are available for regulatory elements, a bottleneck has been created by the lack of bioinformatics tools to predict the interaction modes of regulatory elements. To reduce this gap, we developed the Arabidopsis Transcription Regulatory Factor Domain/Domain Interaction Analysis Tool-liquid/liquid phase separation (LLPS), oligomerization, GO analysis (ART FOUNDATION-LOG), a useful toolkit for protein-nucleic acid interaction (PNI) and protein-protein interaction (PPI) analysis based on domain-domain interactions (DDIs). LLPS, protein oligomerization, the structural properties of protein domains, and protein modifications are major components in the orchestration of the spatiotemporal dynamics of PPIs and PNIs. Our goal is to integrate PPI/PNI information into the development of a prediction model for identifying important genetic variants in peaches. Our program unified interdatabase relational keys based on protein domains to facilitate inference from the model species. A key advantage of this program lies in the integrated information of related features, such as protein oligomerization, LOG analysis, structural characterizations of domains (e.g., domain linkers, intrinsically disordered regions, DDIs, domain-motif (peptide) interactions, beta sheets, and transmembrane helices), and post-translational modification. We provided simple tests to demonstrate how to use this program, which can be applied to other eukaryotic organisms.

Chromosomal features revealed by comparison of genetic maps of Glycine max and Glycine soja.

Recombination is a crucial component of evolution and breeding. New combinations of variation on chromosomes are shaped by recombination. Recombination is also involved in chromosomal rearrangements. However, recombination rates vary tremendously among chromosome segments. Genome-wide genetic maps are one of the best tools to study variation of recombination. Here, we describe high density genetic maps of Glycine max and Glycine soja constructed from four segregating populations. The maps were used to identify chromosomal rearrangements and find the highly predictable pattern of cross-overs on the broad scale in soybean. Markers on these genetic maps were used to evaluate assembly quality of the current soybean reference genome sequence. We find a strong inversion candidate larger than 3 Mb based on patterns of cross-overs. We also identify quantitative trait loci (QTL) that control number of cross-overs. This study provides fundamental insights relevant to practical strategy for breeding programs and for pan-genome researches.

Dissection of soybean populations according to selection signatures based on whole-genome sequences.

BACKGROUND: Domestication and improvement processes, accompanied by selections and adaptations, have generated genome-wide divergence and stratification in soybean populations. Simultaneously, soybean populations, which comprise diverse subpopulations, have developed their own adaptive characteristics enhancing fitness, resistance, agronomic traits, and morphological features. The genetic traits underlying these characteristics play a fundamental role in improving other soybean populations. RESULTS: This study focused on identifying the selection signatures and adaptive characteristics in soybean populations. A core set of 245 accessions (112 wild-type, 79 landrace, and 54 improvement soybeans) selected from 4,234 soybean accessions was re-sequenced. Their genomic architectures were examined according to the domestication and improvement, and accessions were then classified into 3 wild-type, 2 landrace, and 2 improvement subgroups based on various population analyses. Selection and gene set enrichment analyses revealed that the landrace subgroups have selection signals for soybean-cyst nematode HG type 0 and seed development with germination, and that the improvement subgroups have selection signals for plant development with viability and seed development with embryo development, respectively. The adaptive characteristic for soybean-cyst nematode was partially underpinned by multiple resistance accessions, and the characteristics related to seed development were supported by our phenotypic findings for seed weights. Furthermore, their adaptive characteristics were also confirmed as genome-based evidence, and unique genomic regions that exhibit distinct selection and selective sweep patterns were revealed for 13 candidate genes. CONCLUSIONS: Although our findings require further biological validation, they provide valuable information about soybean breeding strategies and present new options for breeders seeking donor lines to improve soybean populations.

Korean soybean core collection: Genotypic and phenotypic diversity population structure and genome-wide association study.

A core collection is a subset that represents genetic diversity of the total collection. Soybean (Glycine max (L.) Merr.) is one of major food and feed crops. It is the world's most cultivated annual herbaceous legume. Constructing a core collection for soybean could play a pivotal role in conserving and utilizing its genetic variability for research and breeding programs. To construct and evaluate a Korean soybean core collection, genotypic and phenotypic data as well as population structure, were analyzed. The Korean soybean core collection consisted of 430 accessions selected from 2,872 collections based on Affymetrix Axiom® 180k SoyaSNP array data. The core collection represented 99% of genotypic diversity of the total collection. Analysis of population structure clustered the core collection into five subpopulations. Accessions from South Korea and North Korea were distributed across five subpopulations. Analysis of molecular variance indicated that only 2.01% of genetic variation could be explained by geographic origins while 16.18% of genetic variation was accounted for by subpopulations. Genome-wide association study (GWAS) for days to flowering, flower color, pubescent color, and growth habit confirmed that the core collection had the same genetic diversity for tested traits as the total collection. The Korean soybean core collection was constructed based on genotypic information of the 180k SNP data. Size and phenotypic diversity of the core collection accounted for approximately 14.9% and 18.1% of the total collection, respectively. GWAS of core and total collections successfully confirmed loci associated with tested traits. Consequently, the present study showed that the Korean soybean core collection could provide fundamental and practical material and information for both soybean genetic research and breeding programs.

Genetic diversity patterns and domestication origin of soybean.

KEY MESSAGE: Genotyping data of a comprehensive Korean soybean collection obtained using a large SNP array were used to clarify global distribution patterns of soybean and address the evolutionary history of soybean. Understanding diversity and evolution of a crop is an essential step to implement a strategy to expand its germplasm base for crop improvement research. Accessions intensively collected from Korea, which is a small but central region in the distribution geography of soybean, were genotyped to provide sufficient data to underpin population genetic questions. After removing natural hybrids and duplicated or redundant accessions, we obtained a non-redundant set comprising 1957 domesticated and 1079 wild accessions to perform population structure analyses. Our analysis demonstrates that while wild soybean germplasm will require additional sampling from diverse indigenous areas to expand the germplasm base, the current domesticated soybean germplasm is saturated in terms of genetic diversity. We then showed that our genome-wide polymorphism map enabled us to detect genetic loci underlying flower color, seed-coat color, and domestication syndrome. A representative soybean set consisting of 194 accessions was divided into one domesticated subpopulation and four wild subpopulations that could be traced back to their geographic collection areas. Population genomics analyses suggested that the monophyletic group of domesticated soybeans was likely originated at a Japanese region. The results were further substantiated by a phylogenetic tree constructed from domestication-associated single nucleotide polymorphisms identified in this study.

Molecular breeding of a novel orange-brown tomato fruit with enhanced beta-carotene and chlorophyll accumulation.

BACKGROUND: Tomatoes provide a significant dietary source of the carotenoids, lycopene and ß-carotene. During ripening, carotenoid accumulation determines the fruit colors while chlorophyll degradation. These traits have been, and continue to be, a significant focus for plant breeding efforts. Previous work has found strong evidence for a relationship between CYC-B gene expression and the orange color of fleshy fruit. Other work has identified a point mutation in SGR that impedes chlorophyll degradation and causes brown flesh color to be retained in some tomato varieties. METHODS: We crossed two inbred lines, KNY2 (orange) and KNB1 (brown) and evaluated the relationship between these genes for their effect on fruit color. Phenotypes of F2 generation plants were analyzed and a novel 'orange-brown' fruit color was identified. RESULTS: We confirm two SNPs, one in CYC-B and another in SGR gene sequence, associated with segregation of 'orange-brown' fruit color in F2 generation. The carotenoid and chlorophyll content of a fleshy fruit was assessed across the different phenotypes and showed a strong correlation with expression pattern of carotenoid biosynthesis genes and SGR function. The orange-brown fruit has high ß-carotene and chlorophyll. Our results provide valuable information for breeders to develop tomato fruit of a novel color using molecular markers.

Identification of haplotypes at the Rsv4 genomic region in soybean associated with durable resistance to soybean mosaic virus.

KEY MESSAGE: Discovery of new germplasm sources and identification of haplotypes for the durable Soybean mosaic virus resistance gene, Rsv 4, provide novel resources for map-based cloning and genetic improvement efforts in soybean. The Soybean mosaic virus (SMV) resistance locus Rsv4 is of interest because it provides a durable type of resistance in soybean [Glycine max (L.) Merr.]. To better understand its molecular basis, we used a population of 309 BC3F2 individuals to fine-map Rsv4 to a ~120 kb interval and leveraged this genetic information in a second study to identify accessions 'Haman' and 'Ilpumgeomjeong' as new sources of Rsv4. These two accessions along with three other Rsv4 and 14 rsv4 accessions were used to examine the patterns of nucleotide diversity at the Rsv4 region based on high-depth resequencing data. Through a targeted association analysis of these 19 accessions within the ~120 kb interval, a cluster of four intergenic single-nucleotide polymorphisms (SNPs) was found to perfectly associate with SMV resistance. Interestingly, this ~120 kb interval did not contain any genes similar to previously characterized dominant disease resistance genes. Therefore, a haplotype analysis was used to further resolve the association signal to a ~94 kb region, which also resulted in the identification of at least two Rsv4 haplotypes. A haplotype phylogenetic analysis of this region suggests that the Rsv4 locus in G. max is recently introgressed from G. soja. This integrated study provides a strong foundation for efforts focused on the cloning of this durable virus resistance gene and marker-assisted selection of Rsv4-mediated SMV resistance in soybean breeding programs.

Development, validation and genetic analysis of a large soybean SNP genotyping array.

Cultivated soybean (Glycine max) suffers from a narrow germplasm relative to other crop species, probably because of under-use of wild soybean (Glycine soja) as a breeding resource. Use of a single nucleotide polymorphism (SNP) genotyping array is a promising method for dissecting cultivated and wild germplasms to identify important adaptive genes through high-density genetic mapping and genome-wide association studies. Here we describe a large soybean SNP array for use in diversity analyses, linkage mapping and genome-wide association analyses. More than four million high-quality SNPs identified from high-depth genome re-sequencing of 16 soybean accessions and low-depth genome re-sequencing of 31 soybean accessions were used to select 180,961 SNPs for creation of the Axiom(®) SoyaSNP array. Validation analysis for a set of 222 diverse soybean lines showed that 170,223 markers were of good quality for genotyping. Phylogenetic and allele frequency analyses of the validation set data indicated that accessions showing an intermediate morphology between cultivated and wild soybeans collected in Korea were natural hybrids. More than 90 unanchored scaffolds in the current soybean reference sequence were assigned to chromosomes using this array. Finally, dense average spacing and preferential distribution of the SNPs in gene-rich chromosomal regions suggest that this array may be suitable for genome-wide association studies of soybean germplasm. Taken together, these results suggest that use of this array may be a powerful method for soybean genetic analyses relating to many aspects of soybean breeding.

Population structure and domestication revealed by high-depth resequencing of Korean cultivated and wild soybean genomes.

Despite the importance of soybean as a major crop, genome-wide variation and evolution of cultivated soybeans are largely unknown. Here, we catalogued genome variation in an annual soybean population by high-depth resequencing of 10 cultivated and 6 wild accessions and obtained 3.87 million high-quality single-nucleotide polymorphisms (SNPs) after excluding the sites with missing data in any accession. Nuclear genome phylogeny supported a single origin for the cultivated soybeans. We identified 10-fold longer linkage disequilibrium (LD) in the wild soybean relative to wild maize and rice. Despite the small population size, the long LD and large SNP data allowed us to identify 206 candidate domestication regions with significantly lower diversity in the cultivated, but not in the wild, soybeans. Some of the genes in these candidate regions were associated with soybean homologues of canonical domestication genes. However, several examples, which are likely specific to soybean or eudicot crop plants, were also observed. Consequently, the variation data identified in this study should be valuable for breeding and for identifying agronomically important genes in soybeans. However, the long LD of wild soybeans may hinder pinpointing causal gene(s) in the candidate regions.

Dynamic genetic features of chromosomes revealed by comparison of soybean genetic and sequence-based physical maps.

Despite the intensive soybean [Glycine max (L.) Merrill] genome studies, the high chromosome number (20) of the soybean plant relative to many other major crops has hindered the development of a high-resolution genomewide genetic map derived from a single population. Here, we report such a map, which was constructed in an F15 population derived from a cross between G. max and G. soja lines using indel polymorphisms detected via a G. soja genome resequencing. By targeting novel indel markers to marker-poor regions, all marker intervals were reduced to under 6 cM on a genome scale. Comparison of the Williams 82 soybean reference genome sequence and our genetic map indicated that marker orders of 26 regions were discrepant with each other. In addition, our comparison showed seven misplaced and two absent markers in the current Williams 82 assembly and six markers placed on the scaffolds that were not incorporated into the pseudomolecules. Then, we showed that, by determining the missing sequences located at the presumed beginning points of the five major discordant segments, these observed discordant regions are mostly errors in the Williams 82 assembly. Distributions of the recombination rates along the chromosomes were similar to those of other organisms. Genotyping of indel markers and genome resequencing of the two parental lines suggested that some marker-poor chromosomal regions may represent introgression regions, which appear to be prevalent in soybean. Given the even and dense distribution of markers, our genetic map can serve as a bridge between genomics research and breeding programs.

Ln is a key regulator of leaflet shape and number of seeds per pod in soybean.

Narrow leaflet soybean (Glycine max) varieties tend to have more seeds per pod than broad leaflet varieties. Narrow leaflet in soybean is conferred by a single recessive gene, ln. Here, we show that the transition from broad (Ln) to narrow leaflet (ln) is associated with an amino acid substitution in the EAR motif encoded by a gene (designated Gm-JAGGED1) homologous to Arabidopsis JAGGED (JAG) that regulates lateral organ development and the variant exerts a pleiotropic effect on fruit patterning. The genomic region that regulates both the traits was mapped to a 12.6-kb region containing only one gene, Gm-JAG1. Introducing the Gm-JAG1 allele into a loss-of-function Arabidopsis jagged mutant partially restored the wild-type JAG phenotypes, including leaf shape, flower opening, and fruit shape, but the Gm-jag1 (ln) and EAR-deleted Gm-JAG1 alleles in the jagged mutant did not result in an apparent phenotypic change. These observations indicate that despite some degree of functional change of Gm-JAG1 due to the divergence from Arabidopsis JAG, Gm-JAG1 complemented the functions of JAG in Arabidopsis thaliana. However, the Gm-JAG1 homoeolog, Gm-JAG2, appears to be sub- or neofunctionalized, as revealed by the differential expression of the two genes in multiple plant tissues, a complementation test, and an allelic analysis at both loci.

Glycine max non-nodulation locus rj1: a recombinogenic region encompassing a SNP in a lysine motif receptor-like kinase (GmNFR1α).

The rj1 mutation of soybean is a simple recessive allele in a single line that arose as a spontaneous mutation in a population; it exhibits non-nodulation with virtually all Bradyrhizobium and Sinorhizobium strains. Here, we described fine genetic and physical mapping of the rj1 locus on soybean chromosome 2. The initial mapping of the rj1 locus using public markers indicated that A343.p2, a sequence-based marker that contains sequence similar to a part of the LjNFR1 gene regulating nodule formation as a member of lysin motif-type receptor-like kinase (LYK) family, maps very close to or cosegregates with the rj1 locus. The sequence of A343.p2 is 100% identical to parts of two BAC clone sequences (GM_WBb0002O19 and GM_WBb098N11) that contain three members of the LYK family. We analyzed the sequence contig (262 kbp) of the two BAC clones by resequencing and subsequent fine genetic and physical mapping. The results indicated that rj1 is located in a gene-rich region with a recombination rate of 120 kbp/cM: several fold higher than the genome average. Among the LYK genes, NFR1α is most likely the gene encoded at the Rj1 locus. The non-nodulating rj1 allele was created by a single base-pair deletion that results in a premature stop codon. Taken together, the fine genetic and physical mapping of the Rj1-residing chromosomal region, combined with the unexpected observation of a putative recombination hotspot, allowed us to demonstrate that the Rj1 locus most likely encodes the NFR1α gene.

Fine genetic mapping of the genomic region controlling leaflet shape and number of seeds per pod in the soybean.

Narrow leaflet cultivars tend to have more seeds per pod than broad leaflet cultivars in soybean [Glycine max (L.) Merr.], which suggests that the leaflet-shape trait locus is tightly linked to or cosegregates with the trait locus controlling the number of seeds per pod (NSPP). Here, we attempted to further elucidate the relationship between leaflet shape and NSPP. A BC(3)F(2) population from a cross between the 'Sowon' (narrow leaflets and high NSPP) and 'V94-5152' (broad leaflets and low NSPP) variants was used. The results of the molecular genetic analyses indicated that, although the NSPP characteristic, in particular, the occurrence of 4-seeded pods, is governed by additional modifying genes that are likely present in Sowon, the two traits cosegregate in the BC(3)F(2) population. The mapping results generated using public markers demonstrated that the narrow leaflet-determining gene in Sowon is an allele of the previously highly studied ln gene on chromosome 20. A high-resolution map delimited the genomic region controlling both the leaflet shape and NSPP traits to a sequence length of 66 kb, corresponding to 0.7 cM. Among the three genes annotated in this 66 kb region, Glyma20g25000.1 appeared to be a good candidate for the Ln-encoding gene, owing to its 47.8% homology with the protein encoding for the JAGGED gene that regulates lateral organ development in Arabidopsis. Taken together, our results suggested that phenotypic variations for narrow leaflet and NSPP are predominantly from the pleiotropic effects of the ln gene. Thus, our results should provide a molecular framework for soybean breeding programs with the objective of improving soybean yield.

Genetic analysis of genes controlling natural variation of seed coat and flower colors in soybean.

Soybean exhibits natural variation in flower and seed coat colors via the deposition of various anthocyanin pigments in the respective tissues. Although pigmentation in seeds or flowers has been well dissected at molecular level in several plant species, the genes controlling natural variation in anthocyanin traits in the soybean are not completely understood. To evaluate the genetic correlation between genetic loci and genes, 8 enzyme-encoding gene families and a transcription factor were localized in a soybean genome-wide genetic map. Among the seed coat color-controlling loci, the genetic location of the gene encoding for W1 was substantiated in the context of the current soybean molecular genetic map and O was postulated to correspond to anthocyanidin reductase. Among the genetic loci that regulate flower pigmentation, the genetic locations of the genes encoding for W1, W4, and Wp were identified, W3 was mapped on soybean linkage group B2 (chromosome 14), and W2 was postulated to correspond to an MYB transcription factor. Correlation studies between the developed markers and 3 color-controlling loci provided important empirical data that should prove useful in the design of marker-assisted breeding schemes as well as future association studies involving soybean.

Genome structure in soybean revealed by a genomewide genetic map constructed from a single population.

A complete genetic linkage map of the soybean, in which sequence-based (SB) genetic markers are evenly distributed genomewide, was constructed from an F(12) population composed of 113 recombinant inbred lines derived from an interspecific cross involving Korean genotypes Hwangkeum and IT182932. Several approaches were employed for the development of 112 novel SB markers targeting both the gaps and the ends of the linkage groups (LGs). The resultant map harbored 20 well-resolved LGs presumed to correspond to the 20 pairs of soybean chromosomes. The map allowed us to identify the important chromosomal structures that were not observed in the integrated genetic maps, to identify the new potentially gene-rich regions, to detect segregation distortion regions within the whole genome, and to extend the ends of the LGs. The results will facilitate the further discovery of agronomically relevant genetic loci in the heretofore neglected chromosomal regions and should also provide some important links between the soybean genetic, physical, and genome sequence maps in the regions.