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.

BackgroundR genes are a key component of genetic interactions between plants and biotrophic bacteria and are known to regulate resistance against bacterial invasion. The most common R proteins contain a nucleotide-binding site and a leucine-rich repeat (NBS-LRR) domain. Some NBS-LRR genes in the soybean genome have also been reported to function in disease resistance. In this study, the number of NBS-LRR genes was found to correlate with the number of disease resistance quantitative trait loci (QTL) that flank these genes in each chromosome. NBS-LRR genes co-localized with disease resistance QTL. The study also addressed the functional redundancy of disease resistance on recently duplicated regions that harbor NBS-LRR genes and NBS-LRR gene expression in the bacterial leaf pustule (BLP)-induced soybean transcriptome.ResultsA total of 319 genes were determined to be putative NBS-LRR genes in the soybean genome. The number of NBS-LRR genes on each chromosome was highly correlated with the number of disease resistance QTL in the 2-Mb flanking regions of NBS-LRR genes. In addition, the recently duplicated regions contained duplicated NBS-LRR genes and duplicated disease resistance QTL, and possessed either an uneven or even number of NBS-LRR genes on each side. The significant difference in NBS-LRR gene expression between a resistant near-isogenic line (NIL) and a susceptible NIL after inoculation of Xanthomonas axonopodis pv. glycines supports the conjecture that NBS-LRR genes have disease resistance functions in the soybean genome.ConclusionsThe number of NBS-LRR genes and disease resistance QTL in the 2-Mb flanking regions of each chromosome was significantly correlated, and several recently duplicated regions that contain NBS-LRR genes harbored disease resistance QTL for both sides. In addition, NBS-LRR gene expression was significantly different between the BLP-resistant NIL and the BLP-susceptible NIL in response to bacterial infection. From these observations, NBS-LRR genes are suggested to contribute to disease resistance in soybean. Moreover, we propose models for how NBS-LRR genes were duplicated, and apply Ks values for each NBS-LRR gene cluster.

Association analysis studies can be used to test for associations between molecular markers and quantitative trait loci (QTL). In this study, a genome-wide scan was performed using 150 simple sequence repeat (SSR) markers to identify QTL associated with seed protein content in soybean. The initial mapping population consisted of two subpopulations of 48 germplasm accessions each, with high or low protein levels based on data from the USDA's Germplasm Resources Information Network website. Intrachromosomal LD extended up to 50 cM with r 2 > 0.1 and 10 cM with r 2 > 0.2 across the accessions. An association map consisting of 150 markers was constructed on the basis of differences in allele frequency distributions between the two subpopulations. Eleven putative QTL were identified on the basis of highly significant markers. Nine of these are in regions where protein QTL have been mapped, but the genomic regions containing Satt431 on LG J and Satt551 on LG M have not been reported in previous linkage mapping studies. Furthermore, these new putative protein QTL do not map near any QTL known to affect maturity. Since biased population structure was known to exist in the original association analysis population, association analyses were also conducted on two similar but independent confirmation populations. Satt431 and Satt551 were also significant in those analyses. These results suggest that our association analysis approach could be a useful alternative to linkage mapping for the identification of unreported regions of the soybean genome containing putative QTL.

To survive in various conditions of CO 2 availability, Chlamydomonas reinhardtii shows adaptive changes, such as induction of a CO 2 -concentrating mechanism, changes in cell organization, and induction of several genes, including a periplasmic carbonic anhydrase (pCA1) encoded by Cah1. Among a collection of insertionally generated mutants, a mutant has been isolated that showed no pCA1 protein and no Cah1 mRNA. This mutant strain, designated cah1-1, has been confirmed to have a disruption in the Cah1 gene caused by a single Arg7 insert. The most interesting feature of cah1-1 is its lack of any significant growth phenotype. There is no major difference in growth or photosynthesis between the wild type and cah1-1 over a pH range from 5.0 to 9.0 even though this mutant apparently lacks Cah1 expression in air. Although the presence of pCA1 apparently gives some minor benefit at very low CO 2 concentrations, the characteristics of this Cah1 null mutant demonstrate that pCA1 is not essential for function of the CO 2 -concentrating mechanism or for growth of C. reinhardtii at limiting CO 2 concentrations.Aquatic and soil-borne photosynthetic organisms, including Chlamydomonas reinhardtii, live in quite variable conditions of CO 2 availability. To survive in limiting CO 2 conditions, C. reinhardtii and other microalgae show adaptive changes, such as induction of a CCM (Badger et al

Bacterial leaf pustule (BLP) disease is caused by Xanthomonas axonopodis pv. glycines (Xag). To investigate the plant basal defence mechanisms induced in response to Xag, differential gene expression in near-isogenic lines (NILs) of BLP-susceptible and BLP-resistant soybean was analysed by RNA-Seq. Of a total of 46 367 genes that were mapped to soybean genome reference sequences, 1978 and 783 genes were found to be up- and down-regulated, respectively, in the BLP-resistant NIL relative to the BLP-susceptible NIL at 0, 6, and 12h after inoculation (hai). Clustering analysis revealed that these genes could be grouped into 10 clusters with different expression patterns. Functional annotation based on gene ontology (GO) categories was carried out. Among the putative soybean defence response genes identified (GO:0006952), 134 exhibited significant differences in expression between the BLP-resistant and -susceptible NILs. In particular, pathogen-associated molecular pattern (PAMP) and damage-associated molecular pattern (DAMP) receptors and the genes induced by these receptors were highly expressed at 0 hai in the BLP-resistant NIL. Additionally, pathogenesis-related (PR)-1 and -14 were highly expressed at 0 hai, and PR-3, -6, and -12 were highly expressed at 12 hai. There were also significant differences in the expression of the core JA-signalling components MYC2 and JASMONATE ZIM-motif. These results indicate that powerful basal defence mechanisms involved in the recognition of PAMPs or DAMPs and a high level of accumulation of defence-related gene products may contribute to BLP resistance in soybean.

Key message Genomic regions associated with seed protein, oil and amino acid contents were identified by genome-wide association analyses. Geographic distributions of haplotypes indicate scope of improvement of these traits. Abstract Soybean [ Glycine max (L.) Merr.] protein and oil are used worldwide in feed, food and industrial materials. Increasing seed protein and oil contents is important; however, protein content is generally negatively correlated with oil content. We conducted a genome-wide association study using phenotypic data collected from five environments for 621 accessions in maturity groups I–IV and 34,014 markers to identify quantitative trait loci (QTL) for seed content of protein, oil and several essential amino acids. Three and five genomic regions were associated with seed protein and oil contents, respectively. One, three, one and four genomic regions were associated with cysteine, methionine, lysine and threonine content (g kg −1 crude protein), respectively. As previously shown, QTL on chromosomes 15 and 20 were associated with seed protein and oil contents, with both exhibiting opposite effects on the two traits, and the chromosome 20 QTL having the most significant effect. A multi-trait mixed model identified trait-specific QTL. A QTL on chromosome 5 increased oil with no effect on protein content, and a QTL on chromosome 10 increased protein content with little effect on oil content. The chromosome 10 QTL co-localized with maturity gene E2 / GmGIa . Identification of trait-specific QTL indicates feasibility to reduce the negative correlation between protein and oil contents. Haplotype blocks were defined at the QTL identified on chromosomes 5, 10, 15 and 20. Frequencies of positive effect haplotypes varied across maturity groups and geographic regions, providing guidance on which alleles have potential to contribute to soybean improvement for specific regions. Electronic supplementary material The online version of this article (10.1007/s00122-019-03304-5) contains supplementary material, which is available to authorized users.

One of the most notable contrasts between the photorespiratory pathway of higher plants and that of many of the green algae including Chlamydomonas reinhardtii lies in the enzymes that serve for oxidation of glycolate to glyoxylate. The gene disrupted by insertional mutagenesis in a high-CO2-requiring mutant, HCR89, of C. reinhardtii was determined to encode glycolate dehydrogenase (EC 1.1.99.14), which serves as the counterpart of glycolate oxidase (EC 1.1.3.15) in classical higher plant photorespiration. Neither glycolate nor D-lactate oxidation from the membrane fraction of HCR89 was detected. Excretion of over-accumulated glycolate into media due to the absence of glycolate dehydrogenase activity was observed for HCR89 under both high- and low-CO2 conditions. Chlamydomonas glycolate dehydrogenase, CrGDH, with a molecular mass of 118 851 Da, comprises a relatively hydrophobic N-terminal region, a FAD-containing domain homologous to the D subunit of the glycolate oxidase complex from Escherischia coli, and an iron–sulfur cluster containing domain homologous to the C subunit of anaerobic glycerol-3-phosphate dehydrogenase complex from Escherichia coli. The second Cys residue in the second iron–sulfur cluster motif of CrGDH is replaced by Asp, as CxxDxxCxxxCP, indicating the second iron–sulfur cluster coordinates most likely 3Fe–4S instead of 4Fe–4S. The membrane association of the glycolate dehydrogenase activity agrees with three predicted transmembrane regions on the iron–sulfur domain.Key words: algae, Chlamydomonas, CO2, glycolate, lactate, mitochondria, photorespiration, photosynthesis.

Two soybean recombinant inbred line populations, Jinpumkong 2 x SS2-2 (J x S) and Iksannamulkong x SS2-2 (I x S) showed population-specific quantitative trait loci (QTLs) for days to flowering (DF) and days to maturity (DM) and these were closely correlated within population. In the present study, we identified QTLs for six yield-related traits with simple sequence repeat markers, and biological correlations between flowering traits and yield-related traits. The yield-related traits included plant height (PH), node numbers of main stem (NNMS), pod numbers per plant (PNPP), seed numbers per pod (SNPP), 100-seed weight (SW), and seed yield per plant (SYPP). Eighteen QTLs for six yield-related traits were detected on nine chromosomes (Chrs), containing four QTLs for PH, two for NNMS, two for PNPP, three for SNPP, five for SW, and two for SYPP. Two highly significant QTLs for PH and NNMS were identified on Chr 6 (LG C2) in both populations where the major flowering gene, E1, and two DF and DM QTLs were located. One other PNPP QTL was also located on this region, explaining 12.9% of phenotypic variation. Other QTLs for yield-related traits showed population-specificity. Two significant SYPP QTLs potentially related with QTLs for SNPP and PNPP were found on the same loci of Chrs 8 (Satt390) and 10 (Sat_108). Also, highly significant positive phenotypic correlations (P < 0.01) were found between DF with PH, NNMS, PNPP, and SYPP in both populations, while flowering was negatively correlated with SNPP and SW in the J x S (P < 0.05) and I x S (P < 0.01) populations. Similar results were also shown between DM and yield-related traits, except for one SW. These QTLs identified may be useful for marker-assisted selection by soybean breeders.