Pericentromeric effects shape the patterns of divergence, retention, and expression of duplicated genes in the paleopolyploid soybean

The evolutionary forces that govern the divergence and retention of duplicated genes in polyploids are poorly understood. In this study, we first investigated the rates of nonsynonymous substitution (Ka) and the rates of synonymous substitution (Ks) for a nearly complete set of genes in the paleopolyploid soybean (Glycine max) by comparing the orthologs between soybean and its progenitor species Glycine soja and then compared the patterns of gene divergence and expression between pericentromeric regions and chromosomal arms in different gene categories. Our results reveal strong associations between duplication status and Ka and gene expression levels and overall low Ks and low levels of gene expression in pericentromeric regions. It is theorized that deleterious mutations can easily accumulate in recombination-suppressed regions, because of Hill-Robertson effects. Intriguingly, the genes in pericentromeric regions-the cold spots for meiotic recombination in soybean-showed significantly lower Ka and higher levels of expression than their homoeologs in chromosomal arms. This asymmetric evolution of two members of individual whole genome duplication (WGD)-derived gene pairs, echoing the biased accumulation of singletons in pericentromeric regions, suggests that distinct genomic features between the two distinct chromatin types are important determinants shaping the patterns of divergence and retention of WGD-derived genes.     

A Comprehensive Analysis of the Cupin Gene Family in Soybean (Glycine max)

Cupin superfamily of proteins, including germin and germin-like proteins (GLPs) from higher plants, is known to play crucial roles in plant development and defense. To date, no systematic analysis has been conducted in soybean (Glycine max) incorporating genome organization, gene structure, expression compendium. In this study, 69 putative Cupin genes were identified from the whole-genome of soybean, which were non-randomly distributed on 17 of the 20 chromosomes. These Gmcupin proteins were phylogenetically clustered into ten distinct subgroups among which the gene structures were highly conserved. Eighteen pairs (52.2%) of duplicate paralogous genes were preferentially retained in duplicated regions of the soybean genome. The distributions of GmCupin genes implied that long segmental duplications contributed significantly to the expansion of the GmCupin gene family. According to the RNA-seq data analysis, most of the Gmcupins were differentially expressed in tissue-specific expression pattern and the expression of some duplicate genes were partially redundant while others showed functional diversity, suggesting the Gmcupins have been retained by substantial subfunctionalization during soybean evolutionary processes. Selective analysis based on single nucleotide polymorphisms (SNPs) in cultivated and wild soybeans revealed sixteen Gmcupins had selected site(s), with all SNPs in Gmcupin10.3 and Gmcupin07.2 genes were selected sites, which implied these genes may have undergone strong selection effects during soybean domestication. Taken together, our results contribute to the functional characterization of Gmcupin genes in soybean.     

Genomic,molecular evolution,and expression analysis of NOX genes in soybean (Glycine max)

Reactive oxygen species (ROS) are versatile signaling molecules in sensing stresses and play critical roles in signaling and development. Plasma membrane NADPH oxidases (NOXs) are key producers of ROS, and play important roles in the regulation of plant-pathogen interactions. Here, we performed a comprehensive analysis of the NOX gene family in the soybean genome (Glycine max) and 17 NOX (GmNOX) genes were identified. Structural analysis revealed that the GmNOX proteins in soybean were as conserved as those in other plants. 8 duplicated gene pairs were formed by a Glycine-specific whole-genome duplication (WGD) event approximately 13 million years ago (Mya). The Ka/Ks ratios of GmNOX genes ranged from 0.04 to 0.28, suggesting that the GmNOX family had undergone purifying selection in soybean. Gene expression patterns showed different expression of these duplicate genes, suggesting that the GmNOXs were retained by substantial subfunctionalization during the soybean evolutionary processes. Subsequently, the expression of GmNOXs in response to drought and phytohormones were characterized via qPCR. Importantly, four GmNOXs showed strong expression in nodules, pointing to their probable involvement in nodulation. Thus, our results shed light on the evolutionary history of this family in soybean and contribute to the functional characterization of GmNOX genes in soybean.     

Nucleotide diversity of the upstream region of the putative MADS-box gene controlling soybean maturity

MADS-box genes are involved in plant reproductive development. However, the role of gene nucleotide diversity in soybean flowering and maturity remains unknown. Therefore, in this study, the distribution of DNA polymorphisms in the putative MADS-box gene located near the quantitative trait loci (QTL) for flowering time and maturity was targeted for association analysis using Glycine max (cultivated soybean) and Glycine soja (wild soybean). Sixteen single nucleotide polymorphisms identified in the upstream region of the putative MADS-box gene around QTL Pod mat 13-7 and Fflr 4-2 on chromosome 7 were found to be highly associated with maturity in soybean. The genetic diversity between cultivated soybeans and the wild relative was comparable, although the early maturity group (EMG) was less diverse than the late maturity group (LMG) of the cultivated soybean. Population size changes of the MADS-box gene in this soybean germplasm appeared to result from non-random selection. A selective pressure seemed to act on this gene in the EMG, while the LMG and G. soja were in genetic equilibrium. Neutrality tests and the constructed neighbor-joining tree indicate that the EMG of G. max has experienced strong artificial selection for its domestication and genetic improvement.     

Evolutionary history of mitogen-activated protein kinase (MAPK) genes in Lotus,Medicago, and Phaseolus

Mitogen-Activated Protein Kinase (MAPK) genes encode proteins that mediate various signaling pathways associated with biotic and abiotic stress responses in eukaryotes. The MAPK genes form a 3-tier signal transduction cascade between cellular stimuli and physiological responses. Recent identification of soybean MAPKs and availability of genome sequences from other legume species allowed us to identify their MAPK genes. The main objectives of this study were to identify MAPKs in 3 legume species, Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, and to assess their phylogenetic relationships. We used approaches in comparative genomics for MAPK gene identification and named the newly identified genes following Arabidopsis MAPK nomenclature model. We identified 19, 18, and 15 MAPKs and 7, 4, and 9 MAPKKs in the genome of Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, respectively. Within clade placement of MAPKs and MAPKKs in the 3 legume species were consistent with those in soybean and Arabidopsis. Among 5 clades of MAPKs, 4 founder clades were consistent to MAPKs of other plant species and orthologs of MAPK genes in the fifth clade-"Clade E" were consistent with those in soybean. Our results also indicated that some gene duplication events might have occurred prior to eudicot-monocot divergence. Highly diversified MAPKs in soybean relative to those in 3 other legume species are attributable to the polyploidization events in soybean. The identification of the MAPK genes in the legume species is important for the legume crop improvement; and evolutionary relationships and functional divergence of these gene members provide insights into plant genome evolution.     

Comparison of the starch synthesis genes between maize and rice: copies, chromosome location and expression divergence

Gene duplication and divergence are important evolutionary processes. It has been suggested that a whole genome duplication (WGD) event occurred in the Gramineae, predating its divergence, and a second WGD occurred in maize during its evolution. In this study we compared the fate of the genes involved in the core pathway of starch biosynthesis following the ancient and second WGDs in maize and rice. In total, thirty starch synthesis genes were detected in the maize genome, which covered all the starch synthesis gene families encoded by 27 genes in rice. All of these genes, except ZmGBSSIIb and ZmBEIII, are anchored within large-scale synteny blocks of rice and maize chromosomes. Previous findings and our results indicate that two of the current copies of many starch synthesis genes (including AGPL, AGPS, GBSS, SSII, SSIII, and BEII) probably arose from the ancient WGD in the Gramineae and are still present in the maize and rice genome. Furthermore, two copies of at least six genes (AGPS1, SSIIb, SSIIIb, GBSSII, BEI, and ISA3) appear to have been retained in the maize genome after its second WGD, although complete coding regions were only detected among the duplicate sets of AGPS1, SSIIb, and SSIIIb. The expression patterns of the remaining duplicate sets of starch synthesis genes (AGPL1/2, AGPS1/2, SSIIa/b, SSIIIa/b, GBSSI/II, and BEIIa/b) differ in their expression and could be classified into two groups in maize. The first group is mainly expressed in the endosperm, whereas the second is expressed in other organs and the early endosperm development. The four duplicate sets of ZmGBSSII, ZmSSIIb, ZmSSIIIb and AGPS1, which arose from the second WGD diverged in gene structure and/or expression patterns in maize. These results indicated that some duplicated starch synthesis genes were remained, whereas others diverged in gene structure and/or expression pattern in maize. For most of the duplicated genes, one of the copies has disappeared in the maize genome after the WGD and the subsequent “diploidization”. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The Role of cis Regulatory Evolution in Maize Domestication

Gene expression differences between divergent lineages caused by modification of cis regulatory elements are thought to be important in evolution. We assayed genome-wide cis and trans regulatory differences between maize and its wild progenitor, teosinte, using deep RNA sequencing in F₁ hybrid and parent inbred lines for three tissue types (ear, leaf and stem). Pervasive regulatory variation was observed with approximately 70% of ∼17,000 genes showing evidence of regulatory divergence between maize and teosinte. However, many fewer genes (1,079 genes) show consistent cis differences with all sampled maize and teosinte lines. For ∼70% of these 1,079 genes, the cis differences are specific to a single tissue. The number of genes with cis regulatory differences is greatest for ear tissue, which underwent a drastic transformation in form during domestication. As expected from the domestication bottleneck, maize possesses less cis regulatory variation than teosinte with this deficit greatest for genes showing maize-teosinte cis regulatory divergence, suggesting selection on cis regulatory differences during domestication. Consistent with selection on cis regulatory elements, genes with cis effects correlated strongly with genes under positive selection during maize domestication and improvement, while genes with trans regulatory effects did not. We observed a directional bias such that genes with cis differences showed higher expression of the maize allele more often than the teosinte allele, suggesting domestication favored up-regulation of gene expression. Finally, this work documents the cis and trans regulatory changes between maize and teosinte in over 17,000 genes for three tissues.     

Global investigation of the co‐evolution of MIRNA genes and microRNA targets during soybean domestication

Although the selection of coding genes during plant domestication has been well studied, the evolution of MIRNA genes (MIRs) and the interaction between microRNAs (miRNAs) and their targets in this process are poorly understood. Here, we present a genome‐wide survey of the selection of MIRs and miRNA targets during soybean domestication and improvement. Our results suggest that, overall, MIRs have higher evolutionary rates than miRNA targets. Nonetheless, they do demonstrate certain similar evolutionary patterns during soybean domestication: MIRs and miRNA targets with high expression and duplication status, and with greater numbers of partners, exhibit lower nucleotide divergence than their counterparts without these characteristics, suggesting that expression level, duplication status, and miRNA–target interaction are essential for evolution of MIRs and miRNA targets. Further investigation revealed that miRNA–target pairs that are subjected to strong purifying selection have greater similarities than those that exhibited genetic diversity. Moreover, mediated by domestication and improvement, the similarities of a large number of miRNA–target pairs in cultivated soybean populations were increased compared to those in wild soybeans, whereas a small number of miRNA–target pairs exhibited decreased similarity, which may be associated with the adoption of particular domestication traits. Taken together, our results shed light on the co‐evolution of MIRs and miRNA targets during soybean domestication.     

Ribosomal RNA genes in soybean and common bean: chromosomal organization, expression, and evolution

Ribosomal RNA (5S and 45S) genes were investigated by FISH in two related legumes: soybean [Glycine max (L.) Merr.] and common bean (Phaseolis vulgaris L.). These species are both members of the same tribe (Phaseoleae), but common bean is diploid while soybean is a tetraploid which has undergone diploidization. In contrast to ploidy expectations, soybean had only one 5S and one 45S rDNA locus whereas common bean had more than two 5S rDNA loci and two 45S rDNA loci. Double hybridization experiments with differentially labelled probes indicated that the soybean 45S and 5S rDNA loci are located on different chromosomes and in their distal regions. Likewise, the common bean 45S and 5S rDNA loci were on unique chromosomes, though two of the 5S rDNA loci were on the same chromosome. FISH analysis of interphase nuclei revealed the spatial arrangement of rDNA loci and suggested expression patterns. In both species, we observed one or more 5S rDNA hybridization sites and two 45S rDNA hybridization sites associated with the nucleolar periphery. The 45S rDNA hybridization patterns frequently exhibited gene puffs as de-condensed chromatin strings within the nucleoli. The other condensed rDNA sites (both 5S and 45S) were spatially distant from the nucleolus in nucleoplasmic regions containing heterochromatin. The distribution of rDNA between the nucleoplasm and the nucleoli is consistent with differential gene expression between homologous alleles and among homoeologous loci.     

The IQD Gene Family in Soybean: Structure,Phylogeny, Evolution and Expression

Members of the plant-specific IQ67-domain (IQD) protein family are involved in plant development and the basal defense response. Although systematic characterization of this family has been carried out in Arabidopsis, tomato (Solanum lycopersicum), Brachypodium distachyon and rice (Oryza sativa), systematic analysis and expression profiling of this gene family in soybean (Glycine max) have not previously been reported. In this study, we identified and structurally characterized IQD genes in the soybean genome. A complete set of 67 soybean IQD genes (GmIQD1–67) was identified using Blast search tools, and the genes were clustered into four subfamilies (IQD I–IV) based on phylogeny. These soybean IQD genes are distributed unevenly across all 20 chromosomes, with 30 segmental duplication events, suggesting that segmental duplication has played a major role in the expansion of the soybean IQD gene family. Analysis of the Ka/Ks ratios showed that the duplicated genes of the GmIQD family primarily underwent purifying selection. Microsynteny was detected in most pairs: genes in clade 1–3 might be present in genome regions that were inverted, expanded or contracted after the divergence; most gene pairs in clade 4 showed high conservation with little rearrangement among these gene-residing regions. Of the soybean IQD genes examined, six were most highly expressed in young leaves, six in flowers, one in roots and two in nodules. Our qRT-PCR analysis of 24 soybean IQD III genes confirmed that these genes are regulated by MeJA stress. Our findings present a comprehensive overview of the soybean IQD gene family and provide insights into the evolution of this family. In addition, this work lays a solid foundation for further experiments aimed at determining the biological functions of soybean IQD genes in growth and development.     

SATB1 Mediates Long-Range Chromatin Interactions: A Dual Regulator of Anti-Apoptotic BCL2 and Pro-Apoptotic NOXA Genes

Aberrant expression of special AT-rich binding protein 1 (SATB1), a global genomic organizer, has been associated with various cancers, which raises the question of how higher-order chromatin structure contributes to carcinogenesis. Disruption of apoptosis is one of the hallmarks of cancer. We previously demonstrated that SATB1 mediated specific long-range chromosomal interactions between the mbr enhancer located within 3’-UTR of the BCL2 gene and the promoter to regulate BCL2 expression during early apoptosis. In the present study, we used chromosome conformation capture (3C) assays and molecular analyses to further investigate the function of the SATB1-mediated higher-order chromatin structure in co-regulation of the anti-apoptotic BCL2 gene and the pro-apoptotic NOXA gene located 3.4Mb downstream on Chromosome 18. We demonstrated that the mbr enhancer spatially juxtaposed the promoters of BCL2 and NOXA genes through SATB1-mediated chromatin-loop in Jurkat cells. Decreased SATB1 levels switched the mbr-BCL2 loop to mbr-NOXA loop, and thus changed expression of these two genes. The SATB1-mediated dynamic switch of the chromatin loop structures was essential for the cooperative expression of the BCL2 and NOXA genes in apoptosis. Notably, the role of SATB1 was specific, since inhibition of SATB1 degradation by caspase-6 inhibitor or caspase-6-resistant SATB1 mutant reversed expression of BCL-2 and NOXA in response to apoptotic stimulation. This study reveals the critical role of SATB1-organized higher-order chromatin structure in regulating the dynamic equilibrium of apoptosis-controlling genes with antagonistic functions and suggests that aberrant SATB1 expression might contribute to cancer development by disrupting the co-regulated genes in apoptosis pathways.     

Structure and chromosomal arrangement of leghemoglobin genes in kidney bean suggest divergence in soybean leghemoglobin gene loci following tetraploidization

We have determined the structure of one of the leghemoglobin (Lb) genes of Phaseolus vulgaris (kidney bean) and deduced the chromosomal arrangement of leghemoglobin genes by genomic hybridizations with Lb and two other sequences, each specific to the 5' or 3' region of the soybean leghemoglobin loci. By comparing this organization with two other species of legumes, Glycine max (soybean) and G. soja (wild soybean), a phylogeny of leghemoglobin gene loci was traced. The intragenic structure of the kidney bean leghemoglobin gene shows the same intron/exon arrangement as that of soybean leghemoglobin genes and extensive sequence homologies in both coding as well as 5' and 3' non-coding regions. The presence in the kidney bean genome of four leghemoglobin genes suggests that tandem duplications of a single primordial plant globin gene had occurred to generate four leghemoglobin genes in an `Lb-locus' before Glycine and Phaseolus species diverged. Chromosome duplication by tetraploidization in Glycine generated two loci containing four genes each. A large deletion in one of the two four-gene loci in soybean resulted in the generation of the Lbc₂ locus containing two leghemoglobin genes, one functional and another pseudo (LbΨ₂). This pseudogene, unlike that present on the main locus, is represented by only two and a half exons and appears to be truncated. The two other truncated genes (LbT₁ and LbT₂) were probably generated similarly in the genome of Glycine spp. following tetraploidization before the divergence of G. max and G. soja.     

Characterization of soybean β-expansin genes and their expression responses to symbiosis,nutrient deficiency,and hormone treatment

Expansins are plant cell wall-loosening proteins encoded by a superfamily of genes including α-expansin, β-expansin, expansin-like A, and expansin-like B proteins. They play a variety of biological roles during plant growth and development. Expansin genes have been reported in many plant species, and results primarily from graminaceous members indicate that β-expansins are more abundant in monocots than in dicots. Soybean [Glycine max (L.) Merr] is an important legume crop. This work identified nine β-expansin gene family members in soybean (GmEXPBs) that were divided into two distinct classes based on phylogeny and gene structure, with divergence between the two groups occurring more in introns than in exons. A total of 887 hormone-responsive and environmental stress-related putative cis-elements from 188 families were found in the 2-kb upstream region of GmEXPBs. Variations in number and type of cis-elements associated with each gene indicate that the function of these genes is differentially regulated by these signals. Expression analysis confirmed that the family members were ubiquitously, yet differentially expressed in soybean. Responsiveness to nutrient deficiency stresses and regulation by auxin (indole-3-acetic acid) and cytokinin (6-benzylaminopurine) varied among GmEXPBs. In addition, most β-expansin genes were associated with symbiosis of soybean inoculated with Rhizobium or abuscular mycorrhizal fungi (AMF). Taken together, these results systematically investigate the characteristics of the entire GmEXPB family in soybean and comprise the first report analyzing the relationship of GmEXPBs with rhizobial or AMF symbiosis. This information is a valuable step in the process of understanding the expansin protein functions in soybean and opens avenues for continued researches.     

A single origin and moderate bottleneck during domestication of soybean (Glycine max): implications from microsatellites and nucleotide sequences

Background and Aims

It is essential to illuminate the evolutionary history of crop domestication in order to understand further the origin and development of modern cultivation and agronomy; however, despite being one of the most important crops, the domestication origin and bottleneck of soybean (Glycine max) are poorly understood. In the present study, microsatellites and nucleotide sequences were employed to elucidate the domestication genetics of soybean.

Methods

The genomes of 79 landrace soybeans (endemic cultivated soybeans) and 231 wild soybeans (G. soja) that represented the species-wide distribution of wild soybean in East Asia were scanned with 56 microsatellites to identify the genetic structure and domestication origin of soybean. To understand better the domestication bottleneck, four nucleotide sequences were selected to simulate the domestication bottleneck.

Key Results

Model-based analysis revealed that most of the landrace genotypes were assigned to the inferred wild soybean cluster of south China, South Korea and Japan. Phylogeny for wild and landrace soybeans showed that all landrace soybeans formed a single cluster supporting a monophyletic origin of all the cultivars. The populations of the nearest branches which were basal to the cultivar lineage were wild soybeans from south China. The coalescent simulation detected a bottleneck severity of K′ = 2 during soybean domestication, which could be explained by a foundation population of 6000 individuals if domestication duration lasted 3000 years.

Conclusions

As a result of integrating geographic distribution with microsatellite genotype assignment and phylogeny between landrace and wild soybeans, a single origin of soybean in south China is proposed. The coalescent simulation revealed a moderate genetic bottleneck with an effective wild soybean population used for domestication estimated to be ≈2 % of the total number of ancestral wild soybeans. Wild soybeans in Asia, especially in south China contain tremendous genetic resources for cultivar improvement.