POWR1 is a domestication gene pleiotropically regulating seed quality and yield in soybean.

Seed protein, oil content and yield are highly correlated agronomically important traits that essentially account for the economic value of soybean. The underlying molecular mechanisms and selection of these correlated seed traits during soybean domestication are, however, less known. Here, we demonstrate that a CCT gene, POWR1, underlies a large-effect protein/oil QTL. A causative TE insertion truncates its CCT domain and substantially increases seed oil content, weight, and yield while decreasing protein content. POWR1 pleiotropically controls these traits likely through regulating seed nutrient transport and lipid metabolism genes. POWR1 is also a domestication gene. We hypothesize that the TE insertion allele is exclusively fixed in cultivated soybean due to selection for larger seeds during domestication, which significantly contributes to shaping soybean with increased yield/seed weight/oil but reduced protein content. This study provides insights into soybean domestication and is significant in improving seed quality and yield in soybean and other crop species.

Selection of GmSWEET39 for oil and protein improvement in soybean.

Soybean [Glycine max (L.) Merr.] was domesticated from wild soybean (G. soja Sieb. and Zucc.) and has been further improved as a dual-use seed crop to provide highly valuable oil and protein for food, feed, and industrial applications. However, the underlying genetic and molecular basis remains less understood. Having combined high-confidence bi-parental linkage mapping with high-resolution association analysis based on 631 whole sequenced genomes, we mapped major soybean protein and oil QTLs on chromosome15 to a sugar transporter gene (GmSWEET39). A two-nucleotide CC deletion truncating C-terminus of GmSWEET39 was strongly associated with high seed oil and low seed protein, suggesting its pleiotropic effect on protein and oil content. GmSWEET39 was predominantly expressed in parenchyma and integument of the seed coat, and likely regulates oil and protein accumulation by affecting sugar delivery from maternal seed coat to the filial embryo. We demonstrated that GmSWEET39 has a dual function for both oil and protein improvement and undergoes two different paths of artificial selection. A CC deletion (CC-) haplotype H1 has been intensively selected during domestication and extensively used in soybean improvement worldwide. H1 is fixed in North American soybean cultivars. The protein-favored (CC+) haplotype H3 still undergoes ongoing selection, reflecting its sustainable role for soybean protein improvement. The comprehensive knowledge on the molecular basis underlying the major QTL and GmSWEET39 haplotypes associated with soybean improvement would be valuable to design new strategies for soybean seed quality improvement using molecular breeding and biotechnological approaches.

Genetic Architecture of Soybean Yield and Agronomic Traits.

Soybean is the world's leading source of vegetable protein and demand for its seed continues to grow. Breeders have successfully increased soybean yield, but the genetic architecture of yield and key agronomic traits is poorly understood. We developed a 40-mating soybean nested association mapping (NAM) population of 5,600 inbred lines that were characterized by single nucleotide polymorphism (SNP) markers and six agronomic traits in field trials in 22 environments. Analysis of the yield, agronomic, and SNP data revealed 23 significant marker-trait associations for yield, 19 for maturity, 15 for plant height, 17 for plant lodging, and 29 for seed mass. A higher frequency of estimated positive yield alleles was evident from elite founder parents than from exotic founders, although unique desirable alleles from the exotic group were identified, demonstrating the value of expanding the genetic base of US soybean breeding.

Author Correction: Dynamic Changes of Genome-Wide DNA Methylation during Soybean Seed Development.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Dynamic Changes of Genome-Wide DNA Methylation during Soybean Seed Development.

Seed development is programmed by expression of many genes in plants. Seed maturation is an important developmental process to soybean seed quality and yield. DNA methylation is a major epigenetic modification regulating gene expression. However, little is known about the dynamic nature of DNA methylation and its effects on gene expression during plant development. Through whole-genome bisulfite sequencing, we showed that DNA methylation went through dynamic changes during seed maturation. An average of 66% CG, 45% CHG and 9% CHH contexts was methylated in cotyledons. CHH methylation levels in cotyledons changed greatly from 6% at the early stage to 11% at the late stage. Transcribed genes were approximately two-fold more likely to be differentially methylated than non-transcribed genes. We identified 40, 66 and 2136 genes containing differentially methylated regions (DMRs) with negative correlation between their expression and methylation in the CG, CHG and CHH contexts, respectively. The majority of the DMR genes in the CHH context were transcriptionally down-regulated as seeds mature: 99% of them during early maturation were down-regulated, and preferentially associated with DNA replication and cell division. The results provide novel insights into the dynamic nature of DNA methylation and its relationship with gene regulation in seed development.

Identification and characterization of large DNA deletions affecting oil quality traits in soybean seeds through transcriptome sequencing analysis.

KEY MESSAGE: Identification and characterization of a 254-kb genomic deletion on a duplicated chromosome segment that resulted in a low level of palmitic acid in soybean seeds using transcriptome sequencing. A large number of soybean genotypes varying in seed oil composition and content have been identified. Understanding the molecular mechanisms underlying these variations is important for breeders to effectively utilize them as a genetic resource. Through design and application of a bioinformatics approach, we identified nine co-regulated gene clusters by comparing seed transcriptomes of nine soybean genotypes varying in oil composition and content. We demonstrated that four gene clusters in the genotypes M23, Jack and N0304-303-3 coincided with large-scale genome rearrangements. The co-regulated gene clusters in M23 and Jack mapped to a previously described 164-kb deletion and a copy number amplification of the Rhg1 locus, respectively. The coordinately down-regulated gene clusters in N0304-303-3 were caused by a 254-kb deletion containing 19 genes including a fatty acyl-ACP thioesterase B gene (FATB1a). This deletion was associated with reduced palmitic acid content in seeds and was the molecular cause of a previously reported nonfunctional FATB1a allele, fap nc . The M23 and N0304-304-3 deletions were located in duplicated genome segments retained from the Glycine-specific whole genome duplication that occurred 13 million years ago. The homoeologous genes in these duplicated regions shared a strong similarity in both their encoded protein sequences and transcript accumulation levels, suggesting that they may have conserved and important functions in seeds. The functional conservation of homoeologous genes may result in genetic redundancy and gene dosage effects for their associated seed traits, explaining why the large deletion did not cause lethal effects or completely eliminate palmitic acid in N0304-303-3.

Identification and characterization of transcript polymorphisms in soybean lines varying in oil composition and content.

BACKGROUND: Variation in seed oil composition and content among soybean varieties is largely attributed to differences in transcript sequences and/or transcript accumulation of oil production related genes in seeds. Discovery and analysis of sequence and expression variations in these genes will accelerate soybean oil quality improvement. RESULTS: In an effort to identify these variations, we sequenced the transcriptomes of soybean seeds from nine lines varying in oil composition and/or total oil content. Our results showed that 69,338 distinct transcripts from 32,885 annotated genes were expressed in seeds. A total of 8,037 transcript expression polymorphisms and 50,485 transcript sequence polymorphisms (48,792 SNPs and 1,693 small Indels) were identified among the lines. Effects of the transcript polymorphisms on their encoded protein sequences and functions were predicted. The studies also provided independent evidence that the lack of FAD2-1A gene activity and a non-synonymous SNP in the coding sequence of FAB2C caused elevated oleic acid and stearic acid levels in soybean lines M23 and FAM94-41, respectively. CONCLUSIONS: As a proof-of-concept, we developed an integrated RNA-seq and bioinformatics approach to identify and functionally annotate transcript polymorphisms, and demonstrated its high effectiveness for discovery of genetic and transcript variations that result in altered oil quality traits. The collection of transcript polymorphisms coupled with their predicted functional effects will be a valuable asset for further discovery of genes, gene variants, and functional markers to improve soybean oil quality.

Systems and evolutionary characterization of microRNAs and their underlying regulatory networks in soybean cotyledons.

MicroRNAs (miRNAs) are an emerging class of small RNAs regulating a wide range of biological processes. Soybean cotyledons evolved as sink tissues to synthesize and store seed reserves which directly affect soybean seed yield and quality. However, little is known about miRNAs and their regulatory networks in soybean cotyledons. We sequenced 292 million small RNA reads expressed in soybean cotyledons, and discovered 130 novel miRNA genes and 72 novel miRNA families. The cotyledon miRNAs arose at various stages of land plant evolution. Evolutionary analysis of the miRNA genes in duplicated genome segments from the recent Glycine whole genome duplication revealed that the majority of novel soybean cotyledon miRNAs were young, and likely arose after the duplication event 13 million years ago. We revealed the evolutionary pathway of a soybean cotyledon miRNA family (soy-miR15/49) that evolved from a neutral invertase gene through an inverted duplication and a series of DNA amplification and deletion events. A total of 304 miRNA genes were expressed in soybean cotyledons. The miRNAs were predicted to target 1910 genes, and form complex miRNA networks regulating a wide range of biological pathways in cotyledons. The comprehensive characterization of the miRNAs and their underlying regulatory networks at gene, pathway and system levels provides a foundation for further studies of miRNAs in cotyledons.

Paramutagenicity of a p1 epiallele in maize.

Complex silencing mechanisms in plants and other kingdoms target transposons, repeat sequences, invasive viral nucleic acids and transgenes, but also endogenous genes and genes involved in paramutation. Paramutation occurs in a heterozygote when a transcriptionally active allele heritably adopts the epigenetic state of a transcriptionally and/or post-transcriptionally repressed allele. P1-rr and its silenced epiallele P1-pr, which encode a Myb-like transcription factor mediating pigmentation in floral organs of Zea mays, differ in their cytosine methylation pattern and chromatin structure at a complex enhancer site. Here, we tested whether P1-pr is able to heritably silence its transcriptionally active P1-rr allele in a heterozygote and whether DNA methylation is associated with the establishment and maintenance of P1-rr silencing. We found that P1-pr participates in paramutation as the repressing allele and P1-rr as the sensitive allele. Silencing of P1-rr is highly variable compared to the inducing P1-pr resulting in a wide range of gene expression. Whereas cytosine methylation at P1-rr is negatively correlated with transcription and pigment levels after segregation of P1-pr, methylation lags behind the establishment of the repressed p1 gene expression. We propose a model in which P1-pr paramutation is triggered by changing epigenetic states of transposons immediately adjacent to a P1-rr enhancer sequence. Considering the vast amount of transposable elements in the maize genome close to regulatory elements of genes, numerous loci could undergo paramutation-induced allele silencing, which could also have a significant impact on breeding agronomically important traits.

Epiallele biogenesis in maize.

We have correlated cytosine methylation of two epialleles, P1-rr and P1-pr, with variation in gene expression and therefore phenotype. The p1 gene in maize encodes a transcription factor that controls phlobaphene pigment accumulation in floral tissues. While cytosine methylation was assayed in various regions spanning 17 kb, the only difference in DNA methylation pattern between the expressed P1-rr allele and the silenced P1-pr allele was detected in a region that consists of a complex arrangement of transposons and adjacent repeats. This region, which comprises the distal enhancer element of P1-rr, is hypermethylated in P1-pr compared to P1-rr. Based on other precedents, we hypothesize that DNA methylation spreads from the transposable elements into the flanking P1-rr enhancer, thereby transcriptionally silencing the gene. Interestingly, P1-pr is reactivated in mutants of the dominant epigenetic modifier Ufo1. DNA methylation in the distal enhancer sequence is significantly reduced, which inversely correlates with increased transcript levels and pigmentation in P1-pr Ufo1 plants. If in general DNA methylation spreads from transposons into adjacent sequences containing regulatory elements for neighboring genes, the corresponding genes could be silenced by chance. Given the large amount of transposable elements in the maize genome, epialleles may be far more frequent than previously estimated.

Divergence of gene regulation through chromosomal rearrangements.

BACKGROUND: The molecular mechanisms that modify genome structures to give birth and death to alleles are still not well understood. To investigate the causative chromosomal rearrangements, we took advantage of the allelic diversity of the duplicated p1 and p2 genes in maize. Both genes encode a transcription factor involved in maysin synthesis, which confers resistance to corn earworm. However, p1 also controls accumulation of reddish pigments in floral tissues and has therefore acquired a new function after gene duplication. p1 alleles vary in their tissue-specific expression, which is indicated in their allele designation: the first suffix refers to red or white pericarp pigmentation and the second to red or white glume pigmentation. RESULTS: Comparing chromosomal regions comprising p1-ww[4Co63], P1-rw1077 and P1-rr4B2 alleles with that of the reference genome, P1-wr[B73], enabled us to reconstruct additive events of transposition, chromosome breaks and repairs, and recombination that resulted in phenotypic variation and chimeric regulatory signals. The p1-ww[4Co63] null allele is probably derived from P1-wr[B73] by unequal crossover between large flanking sequences. A transposon insertion in a P1-wr-like allele and NHEJ (non-homologous end-joining) could have resulted in the formation of the P1-rw1077 allele. A second NHEJ event, followed by unequal crossover, probably led to the duplication of an enhancer region, creating the P1-rr4B2 allele. Moreover, a rather dynamic picture emerged in the use of polyadenylation signals by different p1 alleles. Interestingly, p1 alleles can be placed on both sides of a large retrotransposon cluster through recombination, while functional p2 alleles have only been found proximal to the cluster. CONCLUSIONS: Allelic diversity of the p locus exemplifies how gene duplications promote phenotypic variability through composite regulatory signals. Transposition events increase the level of genomic complexity based not only on insertions but also on excisions that cause DNA double-strand breaks and trigger illegitimate recombination.

Non-Mendelian regulation and allelic variation of methionine-rich delta-zein genes in maize.

Sufficient methionine levels in the seed are critical for the supply of a balanced diet for feed and food. Currently, animal feed is supplemented with chemically synthesized methionine, which could be completely replaced with naturally synthesized methionine. However, insufficient levels of methionine are due to alleles of two genes in the maize genome that are expressed during seed development, which have a high percentage of methionine codons, ranging from 23 to 28%, while free methionine is very low. The two genes, dzs10 and dzs18, belong to the prolamin gene family that arose during the evolution of the grasses and were duplicated during a whole genome duplication event. We have found several dzs10 and dzs18 null alleles caused either by transposon insertion or frame shift mutations. Maize seeds with null mutations of both genes have a normal phenotype in contrast to other prolamin genes, explaining the accumulation of methionine deficiency in normal breeding efforts. Moreover, the trans-regulation of these genes deviates from Mendelian inheritance. One allele of the regulatory locus dzr1 is inherited in a parent-of-origin fashion, while another allele appears to prevent Mendelian segregation of the high-methionine phenotype in backcrosses.

Change of gene structure and function by non-homologous end-joining, homologous recombination, and transposition of DNA.

An important objective in genome research is to relate genome structure to gene function. Sequence comparisons among orthologous and paralogous genes and their allelic variants can reveal sequences of functional significance. Here, we describe a 379-kb region on chromosome 1 of maize that enables us to reconstruct chromosome breakage, transposition, non-homologous end-joining, and homologous recombination events. Such a high-density composition of various mechanisms in a small chromosomal interval exemplifies the evolution of gene regulation and allelic diversity in general. It also illustrates the evolutionary pace of changes in plants, where many of the above mechanisms are of somatic origin. In contrast to animals, somatic alterations can easily be transmitted through meiosis because the germline in plants is contiguous to somatic tissue, permitting the recovery of such chromosomal rearrangements. The analyzed region contains the P1-wr allele, a variant of the genetically well-defined p1 gene, which encodes a Myb-like transcriptional activator in maize. The P1-wr allele consists of eleven nearly perfect P1-wr 12-kb repeats that are arranged in a tandem head-to-tail array. Although a technical challenge to sequence such a structure by shotgun sequencing, we overcame this problem by subcloning each repeat and ordering them based on nucleotide variations. These polymorphisms were also critical for recombination and expression analysis in presence and absence of the trans-acting epigenetic factor Ufo1. Interestingly, chimeras of the p1 and p2 genes, p2/p1 and p1/p2, are framing the P1-wr cluster. Reconstruction of sequence amplification steps at the p locus showed the evolution from a single Myb-homolog to the multi-gene P1-wr cluster. It also demonstrates how non-homologous end-joining can create novel gene fusions. Comparisons to orthologous regions in sorghum and rice also indicate a greater instability of the maize genome, probably due to diploidization following allotetraploidization.