Dynamic rearrangements determine genome organization and useful traits in soybean

Evolutionary fates of gene-body methylation and its divergent association with gene expression in pigeonpea

Pigeonpea (Cajanus cajan L. Huth) is an agronomically important legume cultivated worldwide. In this study, we extensively analyzed gene-body methylation (GbM) patterns in pigeonpea. We found a bimodal distribution of CG and CHG methylation patterns. GbM features- slow evolution rate and increased length remained conserved. Genes with moderate CG body methylation showed highest expression where as highly-methylated genes showed lowest expression. Transposable element (TE)-related genes were methylated in multiple contexts and hence classified as C-methylated genes. A low expression among C-methylated genes was associated with transposons insertion in gene-body and upstream regulatory regions. The CG methylation patterns were found to be conserved in orthologs compared with non-CG methylation. By comparing methylation patterns between differentially methylated regions (DMRs) of the three genotypes, we found that variably methylated marks are less likely to target evolutionary conserved sequences. Finally, our analysis showed enrichment of nitrogen-related genes in GbM orthologs of legumes, which could be promising candidates for generating epialleles for crop improvement.

Altered chromatin architecture and gene expression during polyploidization and domestication of soybean

Polyploidy or whole-genome duplication (WGD) is widespread in plants and is a key driver of evolution and speciation, accompanied by rapid and dynamic changes in genomic structure and gene expression. The 3D structure of the genome is intricately linked to gene expression, but its role in transcription regulation following polyploidy and domestication remains unclear. Here, we generated high-resolution (∼2 kb) Hi-C maps for cultivated soybean (Glycine max), wild soybean (Glycine soja), and common bean (Phaseolus vulgaris). We found polyploidization in soybean may induce architecture changes of topologically associating domains and subsequent diploidization led to chromatin topology alteration around chromosome-rearrangement sites. Compared with single-copy and small-scale duplicated genes, WGD genes displayed more long-range chromosomal interactions and were coupled with higher levels of gene expression and chromatin accessibilities but void of DNA methylation. Interestingly, chromatin loop reorganization was involved in expression divergence of the genes during soybean domestication. Genes with chromatin loops were under stronger artificial selection than genes without loops. These findings provide insights into the roles of dynamic chromatin structures on gene expression during polyploidization, diploidization, and domestication of soybean.

Identification of Genomic Insertion and Flanking Sequence of G2-EPSPS and GAT Transgenes in Soybean Using Whole Genome Sequencing Method

Molecular characterization of sequence flanking exogenous fragment insertion is essential for safety assessment and labeling of genetically modified organism (GMO). In this study, the T-DNA insertion sites and flanking sequences were identified in two newly developed transgenic glyphosate-tolerant soybeans GE-J16 and ZH10-6 based on whole genome sequencing (WGS) method. More than 22.4 Gb sequence data (∼21 × coverage) for each line was generated on Illumina HiSeq 2500 platform. The junction reads mapped to boundaries of T-DNA and flanking sequences in these two events were identified by comparing all sequencing reads with soybean reference genome and sequence of transgenic vector. The putative insertion loci and flanking sequences were further confirmed by PCR amplification, Sanger sequencing, and co-segregation analysis. All these analyses supported that exogenous T-DNA fragments were integrated in positions of Chr19: 50543767-50543792 and Chr17: 7980527-7980541 in these two transgenic lines. Identification of genomic insertion sites of G2-EPSPS and GAT transgenes will facilitate the utilization of their glyphosate-tolerant traits in soybean breeding program. These results also demonstrated that WGS was a cost-effective and rapid method for identifying sites of T-DNA insertions and flanking sequences in soybean.

Of the Nine Cytidine Deaminase-Like Genes in Arabidopsis, Eight Are Pseudogenes and Only One Is Required to Maintain Pyrimidine Homeostasis in Vivo

CYTIDINE DEAMINASE (CDA) catalyzes the deamination of cytidine to uridine and ammonia in the catabolic route of C nucleotides. The Arabidopsis (Arabidopsis thaliana) CDA gene family comprises nine members, one of which (AtCDA) was shown previously in vitro to encode an active CDA. A possible role in C-to-U RNA editing or in antiviral defense has been discussed for other members. A comprehensive bioinformatic analysis of plant CDA sequences, combined with biochemical functionality tests, strongly suggests that all Arabidopsis CDA family members except AtCDA are pseudogenes and that most plants only require a single CDA gene. Soybean (Glycine max) possesses three CDA genes, but only two encode functional enzymes and just one has very high catalytic efficiency. AtCDA and soybean CDAs are located in the cytosol. The functionality of AtCDA in vivo was demonstrated with loss-of-function mutants accumulating high amounts of cytidine but also CMP, cytosine, and some uridine in seeds. Cytidine hydrolysis in cda mutants is likely caused by NUCLEOSIDE HYDROLASE1 (NSH1) because cytosine accumulation is strongly reduced in a cda nsh1 double mutant. Altered responses of the cda mutants to fluorocytidine and fluorouridine indicate that a dual specific nucleoside kinase is involved in cytidine as well as uridine salvage. CDA mutants display a reduction in rosette size and have fewer leaves compared with the wild type, which is probably not caused by defective pyrimidine catabolism but by the accumulation of pyrimidine catabolism intermediates reaching toxic concentrations.

Genome-wide identification and evolutionary analysis of leucine-rich repeat receptor-like protein kinase genes in soybean

BACKGROUND

Leucine-rich repeat receptor-like kinases (LRR-RLKs) constitute the largest subfamily of receptor-like kinases in plant. A number of reports have demonstrated that plant LRR-RLKs play important roles in growth, development, differentiation, and stress responses. However, no comprehensive analysis of this gene family has been carried out in legume species.

RESULTS

Based on the principles of sequence similarity and domain conservation, a total of 467 LRR-RLK genes were identified in soybean genome. The GmLRR-RLKs are non-randomly distributed across all 20 chromosomes of soybean and about 73.3 % of them are located in segmental duplicated regions. The analysis of synonymous substitutions for putative paralogous gene pairs indicated that most of these gene pairs resulted from segmental duplications in soybean genome. Furthermore, the exon/intron organization, motif composition and arrangements were considerably conserved among members of the same groups or subgroups in the constructed phylogenetic tree. The close phylogenetic relationship between soybean LRR-RLK genes with identified Arabidopsis genes in the same group also provided insight into their putative functions. Expression profiling analysis of GmLRR-RLKs suggested that they appeared to be differentially expressed among different tissues and some of duplicated genes exhibited divergent expression patterns. In addition, artificial selected GmLRR-RLKs were also identified by comparing the SNPs between wild and cultivated soybeans and 17 genes were detected in regions previously reported to contain domestication-related QTLs.

CONCLUSIONS

Comprehensive and evolutionary analysis of soybean LRR-RLK gene family was performed at whole genome level. The data provides valuable tools in future efforts to identify functional divergence of this gene family and gene diversity among different genotypes in legume species.

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.

Genome duplication and gene loss affect the evolution of heat shock transcription factor genes in legumes

Whole-genome duplication events (polyploidy events) and gene loss events have played important roles in the evolution of legumes. Here we show that the vast majority of Hsf gene duplications resulted from whole genome duplication events rather than tandem duplication, and significant differences in gene retention exist between species. By searching for intraspecies gene colinearity (microsynteny) and dating the age distributions of duplicated genes, we found that genome duplications accounted for 42 of 46 Hsf-containing segments in Glycine max, while paired segments were rarely identified in Lotus japonicas, Medicago truncatula and Cajanus cajan. However, by comparing interspecies microsynteny, we determined that the great majority of Hsf-containing segments in Lotus japonicas, Medicago truncatula and Cajanus cajan show extensive conservation with the duplicated regions of Glycine max. These segments formed 17 groups of orthologous segments. These results suggest that these regions shared ancient genome duplication with Hsf genes in Glycine max, but more than half of the copies of these genes were lost. On the other hand, the Glycine max Hsf gene family retained approximately 75% and 84% of duplicated genes produced from the ancient genome duplication and recent Glycine-specific genome duplication, respectively. Continuous purifying selection has played a key role in the maintenance of Hsf genes in Glycine max. Expression analysis of the Hsf genes in Lotus japonicus revealed their putative involvement in multiple tissue-/developmental stages and responses to various abiotic stimuli. This study traces the evolution of Hsf genes in legume species and demonstrates that the rates of gene gain and loss are far from equilibrium in different species.

The ureide-degrading reactions of purine ring catabolism employ three amidohydrolases and one aminohydrolase in Arabidopsis, soybean, and rice

Several ureides are intermediates of purine base catabolism, releasing nitrogen from the purine nucleotides for reassimilation into amino acids. In some legumes like soybean (Glycine max), ureides are used for nodule-to-shoot translocation of fixed nitrogen. Four enzymes of Arabidopsis (Arabidopsis thaliana), (1) allantoinase, (2) allantoate amidohydrolase (AAH), (3) ureidoglycine aminohydrolase, and (4) ureidoglycolate amidohydrolase (UAH), catalyze the complete hydrolysis of the ureide allantoin in vitro. However, the metabolic route in vivo remains controversial. Here, in growth and metabolite analyses of Arabidopsis mutants, we demonstrate that these enzymes are required for allantoin degradation in vivo. Orthologous enzymes are present in soybean, encoded by one to four gene copies. All isoenzymes are active in vitro, while some may be inefficiently translated in vivo. Surprisingly, transcript and protein amounts are not significantly regulated by nitrogen fixation or leaf ureide content. A requirement for soybean AAH and UAH for ureide catabolism in leaves has been demonstrated by the use of virus-induced gene silencing. Functional AAH, ureidoglycine aminohydrolase, and UAH are also present in rice (Oryza sativa), and orthologous genes occur in all other plant genomes sequenced to date, indicating that the amidohydrolase route of ureide degradation is universal in plants, including mosses (e.g. Physcomitrella patens) and algae (e.g. Chlamydomomas reinhardtii).

A comparison of the molecular organization of genomic regions associated with resistance to common bacterial blight in two Phaseolus vulgaris genotypes

Resistance to common bacterial blight, caused by Xanthomonas axonopodis pv. phaseoli, in Phaseolus vulgaris is conditioned by several loci on different chromosomes. Previous studies with OAC-Rex, a CBB-resistant, white bean variety of Mesoamerican origin, identified two resistance loci associated with the molecular markers Pv-CTT001 and SU91, on chromosome 4 and 8, respectively. Resistance to CBB is assumed to be derived from an interspecific cross with Phaseolus acutifolius in the pedigree of OAC-Rex. Our current whole genome sequencing effort with OAC-Rex provided the opportunity to compare its genome in the regions associated with CBB resistance with the v1.0 release of the P. vulgaris line G19833, which is a large seeded bean of Andean origin, and (assumed to be) CBB susceptible. In addition, the genomic regions containing SAP6, a marker associated with P. vulgaris-derived CBB-resistance on chromosome 10, were compared. These analyses indicated that gene content was highly conserved between G19833 and OAC-Rex across the regions examined (>80%). However, fifty-nine genes unique to OAC Rex were identified, with resistance gene homologues making up the largest category (10 genes identified). Two unique genes in OAC-Rex located within the SU91 resistance QTL have homology to P. acutifolius ESTs and may be potential sources of CBB resistance. As the genomic sequence assembly of OAC-Rex is completed, we expect that further comparisons between it and the G19833 genome will lead to a greater understanding of CBB resistance in bean.

Gene divergence of homeologous regions associated with a major seed protein content QTL in soybean

Understanding several modes of duplication contributing on the present genome structure is getting an attention because it could be related to numerous agronomically important traits. Since soybean serves as a rich protein source for animal feeds and human consumption, breeding efforts in soybean have been directed toward enhancing seed protein content. The publicly available soybean sequences and its genomically featured elements facilitate comprehending of quantitative trait loci (QTL) for seed protein content in concordance with homeologous regions in soybean genome. Although parts of chromosome (Chr) 20 and Chr 10 showed synteny, QTLs for seed protein content present only on Chr 20. Using comparative analysis of gene contents in recently duplicated genomic regions harboring QTL for protein/oil content on Chrs 20 and 10, a total of 27 genes are present in duplicated regions of both Chrs. Notably, 4 tandem duplicates of the putative homeobox protein 22 (HB22) are present only on Chr 20 and this Medicago truncatula homolog expressed in endosperm at seed filling stage. These tandem duplicates could contribute on the protein/oil QTL of Chr 20. Our study suggests that non-shared gene contents within the duplicated genomic regions might lead to absence/presence of QTL related to protein/oil content.

Open or close the gate - stomata action under the control of phytohormones in drought stress conditions

Two highly specialized cells, the guard cells that surround the stomatal pore, are able to integrate environmental and endogenous signals in order to control the stomatal aperture and thereby the gas exchange. The uptake of CO2 is associated with a loss of water by leaves. Control of the size of the stomatal aperture optimizes the efficiency of water use through dynamic changes in the turgor of the guard cells. The opening and closing of stomata is regulated by the integration of environmental signals and endogenous hormonal stimuli. The various different factors to which the guard cells respond translates into the complexity of the network of signaling pathways that control stomatal movements. The perception of an abiotic stress triggers the activation of signal transduction cascades that interact with or are activated by phytohormones. Among these, abscisic acid (ABA), is the best-known stress hormone that closes the stomata, although other phytohormones, such as jasmonic acid, brassinosteroids, cytokinins, or ethylene are also involved in the stomatal response to stresses. As a part of the drought response, ABA may interact with jasmonic acid and nitric oxide in order to stimulate stomatal closure. In addition, the regulation of gene expression in response to ABA involves genes that are related to ethylene, cytokinins, and auxin signaling. In this paper, recent findings on phytohormone crosstalk, changes in signaling pathways including the expression of specific genes and their impact on modulating stress response through the closing or opening of stomata, together with the highlights of gaps that need to be elucidated in the signaling network of stomatal regulation, are reviewed.

Comparative sequence analysis of nitrogen fixation-related genes in six legumes

Legumes play an important role as food and forage crops in international agriculture especially in developing countries. Legumes have a unique biological process called nitrogen fixation (NF) by which they convert atmospheric nitrogen to ammonia. Although legume genomes have undergone polyploidization, duplication and divergence, NF-related genes, because of their essential functional role for legumes, might have remained conserved. To understand the relationship of divergence and evolutionary processes in legumes, this study analyzes orthologs and paralogs for selected 20 NF-related genes by using comparative genomic approaches in six legumes i.e., Medicago truncatula (Mt), Cicer arietinum, Lotus japonicus, Cajanus cajan (Cc), Phaseolus vulgaris (Pv), and Glycine max (Gm). Subsequently, sequence distances, numbers of synonymous substitutions per synonymous site (Ks) and non-synonymous substitutions per non-synonymous site (Ka) between orthologs and paralogs were calculated and compared across legumes. These analyses suggest the closest relationship between Gm and Cc and the highest distance between Mt and Pv in six legumes. Ks proportional plots clearly showed ancient genome duplication in all legumes, whole genome duplication event in Gm and also speciation pattern in different legumes. This study also reports some interesting observations e.g., no peak at Ks 0.4 in Gm-Gm, location of two independent genes next to each other in Mt and low Ks values for outparalogs for three genes as compared to other 12 genes. In summary, this study underlines the importance of NF-related genes and provides important insights in genome organization and evolutionary aspects of six legume species analyzed.

Divergence of flowering genes in soybean

Soybean genome sequences were blasted with Arabidopsis thaliana regulatory genes involved in photoperioddependent flowering. This approach enabled the identification of 118 genes involved in the flowering pathway. Two genome sequences of cultivated (Williams 82) and wild (IT182932) soybeans were employed to survey functional DNA variations in the flowering-related homologs. Forty genes exhibiting nonsynonymous substitutions between G. max and G. soja were catalogued. In addition, 22 genes were found to co-localize with QTLs for six traits including flowering time, first flower, pod maturity, beginning of pod, reproductive period, and seed filling period. Among the genes overlapping the QTL regions, two LHY/CCA1 genes, GI and SFR6 contained amino acid changes. The recently duplicated sequence regions of the soybean genome were used as additional criteria for the speculation of the putative function of the homologs. Two duplicated regions showed redundancy of both flowering-related genes and QTLs. ID 12398025, which contains the homeologous regions between chr 7 and chr 16, was redundant for the LHY/CCA1 and SPA1 homologs and the QTLs. Retaining of the CRY1 gene and the pod maturity QTLs were observed in the duplicated region of ID 23546507 on chr 4 and chr 6. Functional DNA variation of the LHY/CCA1 gene (Glyma07g05410) was present in a counterpart of the duplicated region on chr 7, while the gene (Glyma16g01980) present in the other portion of the duplicated region on chr 16 did not show a functional sequence change. The gene list catalogued in this study provides primary insight for understanding the regulation of flowering time and maturity in soybean.

Genome-wide mapping of NBS-LRR genes and their association with disease resistance in soybean

BACKGROUND

R 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.

RESULTS

A 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.

CONCLUSIONS

The 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.

Four allantoinase genes are expressed in nitrogen-fixing soybean

Soybean (Glycine max L. [Merr]) plants export nitrogen from the nodules as ureides during symbiotic biological nitrogen fixation. Ureides also play a role as nitrogen storage compounds in the seeds and are broken down in germinating seedlings. In this work we identified four soybean genes encoding allantoinase (E.C. 3.5.2.5), an enzyme involved in both ureide production in nodules and ureide catabolism in leaves and other sink tissues. We examined ureide content, allantoinase enzyme activity and expression of these genes, which we term GmALN1 through GmALN4, in germinating seedlings and in vegetative tissues from 45 day old soybean plants. GmALN1 and GmALN2 transcripts were measured in all tissues, but similar levels of expression of GmALN3 and GmALN4 was only observed in nodules. The soybean allantoinase genes seem to have arisen through tandem gene duplication followed by a whole genome duplication. We looked for evidence of the tandem duplication in common bean (Phaseolus vulgaris L.) and present evidence that it occured sometime in the bean lineage before these two species diverged, but before soybean became a tetraploid.

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.

Genome-enabled insights into legume biology

Legumes are the third-largest family of angiosperms, the second-most-important crop family, and a key source of biological nitrogen in agriculture. Recently, the genome sequences of Glycine max (soybean), Medicago truncatula, and Lotus japonicus were substantially completed. Comparisons among legume genomes reveal a key role for duplication, especially a whole-genome duplication event approximately 58 Mya that is shared by most agriculturally important legumes. A second and more recent genome duplication occurred only in the lineage leading to soybean. Outcomes of genome duplication, including gene fractionation and sub- and neofunctionalization, have played key roles in shaping legume genomes and in the evolution of legume-specific traits. Analysis of legume genome sequences also enables the discovery of legume-specific gene families and provides a framework for genome-wide association mapping that will target phenotypes of special importance in legumes. Translating genomic resources from sequenced species to less studied but still important "orphan" legumes will enhance prospects for world food production.

RNA-Seq analysis of a soybean near-isogenic line carrying bacterial leaf pustule-resistant and -susceptible alleles

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.

Differential expression of cysteine desulfurases in soybean

BACKGROUND

Iron-sulfur [Fe-S] clusters are prosthetic groups required to sustain fundamental life processes including electron transfer, metabolic reactions, sensing, signaling, gene regulation and stabilization of protein structures. In plants, the biogenesis of Fe-S protein is compartmentalized and adapted to specific needs of the cell. Many environmental factors affect plant development and limit productivity and geographical distribution. The impact of these limiting factors is particularly relevant for major crops, such as soybean, which has worldwide economic importance.

RESULTS

Here we analyze the transcriptional profile of the soybean cysteine desulfurases NFS1, NFS2 and ISD11 genes, involved in the biogenesis of [Fe-S] clusters, by quantitative RT-PCR. NFS1, ISD11 and NFS2 encoding two mitochondrial and one plastid located proteins, respectively, are duplicated and showed distinct transcript levels considering tissue and stress response. NFS1 and ISD11 are highly expressed in roots, whereas NFS2 showed no differential expression in tissues. Cold-treated plants showed a decrease in NFS2 and ISD11 transcript levels in roots, and an increased expression of NFS1 and ISD11 genes in leaves. Plants treated with salicylic acid exhibited increased NFS1 transcript levels in roots but lower levels in leaves. In silico analysis of promoter regions indicated the presence of different cis-elements in cysteine desulfurase genes, in good agreement with differential expression of each locus. Our data also showed that increasing of transcript levels of mitochondrial genes, NFS1/ISD11, are associated with higher activities of aldehyde oxidase and xanthine dehydrogenase, two cytosolic Fe-S proteins.

CONCLUSIONS

Our results suggest a relationship between gene expression pattern, biochemical effects, and transcription factor binding sites in promoter regions of cysteine desulfurase genes. Moreover, data show proportionality between NFS1 and ISD11 genes expression.

Changes in twelve homoeologous genomic regions in soybean following three rounds of polyploidy

With the advent of high-throughput sequencing, the availability of genomic sequence for comparative genomics is increasing exponentially. Numerous completed plant genome sequences enable characterization of patterns of the retention and evolution of genes within gene families due to multiple polyploidy events, gene loss and fractionation, and differential evolutionary pressures over time and across different gene families. In this report, we trace the changes that have occurred in 12 surviving homoeologous genomic regions from three rounds of polyploidy that contributed to the current Glycine max genome: a genome triplication before the origin of the rosids (~130 to 240 million years ago), a genome duplication early in the legumes (~58 million years ago), and a duplication in the Glycine lineage (~13 million years ago). Patterns of gene retention following the genome triplication event generally support predictions of the Gene Balance Hypothesis. Finally, we find that genes in networks with a high level of connectivity are more strongly conserved than those with low connectivity and that the enrichment of these highly connected genes in the 12 highly conserved homoeologous segments may in part explain their retention over more than 100 million years and repeated polyploidy events.

QTL identification of flowering time at three different latitudes reveals homeologous genomic regions that control flowering in soybean

Since the genetic control of flowering time is very important in photoperiod-sensitive soybean (Glycine max (L.) Merr.), genes affecting flowering under different environment conditions have been identified and described. The objectives were to identify quantitative trait loci (QTLs) for flowering time in different latitudinal and climatic regions, and to understand how chromosomal rearrangement and genome organization contribute to flowering time in soybean. Recombinant inbred lines from a cross between late-flowering 'Jinpumkong 2' and early-flowering 'SS2-2' were used to evaluate the phenotypic data for days to flowering (DF) collected from Kamphaeng Saen, Thailand (14°01'N), Suwon, Korea (37°15'N), and Longjing, China (42°46'N). A weakly positive phenotypic correlation (r = 0.36) was found between DF in Korea and Thailand; however, a strong correlation (r = 0.74) was shown between Korea and China. After 178 simple sequence repeat (SSR) markers were placed on a genetic map spanning 2,551.7 cM, four independent DF QTLs were identified on different chromosomes (Chrs). Among them, three QTLs on Chrs 9, 13 and 16 were either Thailand- or Korea-specific. The DF QTL on Chr 6 was identified in both Korea and China, suggesting it is less environment-sensitive. Comparative analysis of four DF QTL regions revealed a syntenic relationship between two QTLs on Chrs 6 and 13. All five duplicated gene pairs clustered in the homeologous genomic regions were found to be involved in the flowering. Identification and comparative analysis of multiple DF QTLs from different environments will facilitate the significant improvement in soybean breeding programs with respect to control of flowering time.

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.

The Protein Disulfide Isomerase gene family in bread wheat (T. aestivum L.)

BACKGROUND

The Protein Disulfide Isomerase (PDI) gene family encodes several PDI and PDI-like proteins containing thioredoxin domains and controlling diversified metabolic functions, including disulfide bond formation and isomerisation during protein folding. Genomic, cDNA and promoter sequences of the three homologous wheat genes encoding the "typical" PDI had been cloned and characterized in a previous work. The purpose of present research was the cloning and characterization of the complete set of genes encoding PDI and PDI like proteins in bread wheat (Triticum aestivum cv Chinese Spring) and the comparison of their sequence, structure and expression with homologous genes from other plant species.

RESULTS

Eight new non-homologous wheat genes were cloned and characterized. The nine PDI and PDI-like sequences of wheat were located in chromosome regions syntenic to those in rice and assigned to eight plant phylogenetic groups. The nine wheat genes differed in their sequences, genomic organization as well as in the domain composition and architecture of their deduced proteins; conversely each of them showed high structural conservation with genes from other plant species in the same phylogenetic group. The extensive quantitative RT-PCR analysis of the nine genes in a set of 23 wheat samples, including tissues and developmental stages, showed their constitutive, even though highly variable expression.

CONCLUSIONS

The nine wheat genes showed high diversity, while the members of each phylogenetic group were highly conserved even between taxonomically distant plant species like the moss Physcomitrella patens. Although constitutively expressed the nine wheat genes were characterized by different expression profiles reflecting their different genomic organization, protein domain architecture and probably promoter sequences; the high conservation among species indicated the ancient origin and diversification of the still evolving gene family. The comprehensive structural and expression characterization of the complete set of PDI and PDI-like wheat genes represents a basis for the functional characterization of this gene family in the hexaploid context of bread wheat.

Fine mapping of a resistance gene to bacterial leaf pustule in soybean

Soybean bacterial leaf pustule (BLP) is a prevalent disease caused by Xanthomonas axonopodis pv. glycines. Fine mapping of the BLP resistant gene, rxp, is needed to select BLP resistant soybean cultivars by marker-assisted selection (MAS). We used a total of 227 recombinant inbred lines (RILs) derived from a cross between 'Taekwangkong' (BLP susceptible) and 'Danbaekkong' (BLP resistant) for rxp fine mapping and two different sets of near isogenic lines (NILs) from Hwangkeumkong x SS2-2 and Taekwangkong x SS2-2 were used for confirmation. Using sequences between Satt372 and Satt486 flanking rxp from soybean genome sequences, eight simple sequence repeats (SSR) and two single nucleotide polymorphism (SNP) markers were newly developed in a 6.2-cM interval. Linkage mapping with the RILs and NILs allowed us to map the rxp region with high resolution. The genetic order of all markers was completely consistent with their physical order. QTL analysis by comparison of the BLP phenotyping data with all markers showed rxp was located between SNUSSR17_9 and SNUSNP17_12. Gene annotation analysis of the 33 kb region between SNUSSR17_9 and SNUSNP17_12 suggested three predicted genes, two of which could be candidate genes of BLP resistance: membrane protein and zinc finger protein. Candidate genes showed high similarity with their paralogous genes, which were located on the duplicated regions obtaining BLP resistance QTLs. High-resolution map in rxp region with eight SSR and two SNP markers will be useful for not only MAS of BLP resistance but also characterization of rxp.