Methylation affects transposition and splicing of a large CACTA transposon from a MYB transcription factor regulating anthocyanin synthase genes in soybean seed coats

Identifying Candidate Genes Related to Soybean (Glycine max) Seed Coat Color via RNA-Seq and Coexpression Network Analysis

BACKGROUND

The quality of soybeans is reflected in the seed coat color, which indicates soybean quality and commercial value. Researchers have identified genes related to seed coat color in various plants. However, research on the regulation of genes related to seed coat color in soybeans is rare.

METHODS

In this study, four lines of seed coats with different colors (medium yellow 14, black, green, and brown) were selected from the F2:5 population, with Beinong 108 as the female parent and green bean as the male parent, and the dynamic changes in the anthocyanins in the seed coat were stained with 4-dimethylaminocinnamaldehyde (DMACA) during the grain maturation process (20 days from grain drum to seed harvest). Through RNA-seq of soybean lines with four different colored seed coats at 30 and 50 days after seeding, we can further understand the key pathways and gene regulation modules between soybean seed coats of different colors.

RESULTS

DMACA revealed that black seed coat soybeans produce anthocyanins first and have the deepest staining. Clustering and principal component analysis (PCA) of the RNA-seq data divided the eight samples into two groups, resulting in 16,456 DEGs, including 5359 TFs. GO and KEGG enrichment analyses revealed that the flavonoid biosynthesis, starch and sucrose metabolism, carotenoid biosynthesis, and circadian rhythm pathways were significantly enriched. We also conducted statistical and expression pattern analyses on the differentially expressed transcription factors. Based on weighted gene coexpression network analysis (WGCNA), we identified seven specific modules that were significantly related to the four soybean lines with different seed coat colors. The connectivity and functional annotation of genes within the modules were calculated, and 21 candidate genes related to soybean seed coat color were identified, including six transcription factor (TF) genes and three flavonoid pathway genes.

CONCLUSIONS

These findings provide a theoretical basis for an in-depth understanding of the molecular mechanisms underlying differences in soybean seed coat color and provide new genetic resources.

Tandemly duplicated MYB genes are functionally diverged in the regulation of anthocyanin biosynthesis in soybean

Gene duplications have long been recognized as a driving force in the evolution of genes, giving rise to novel functions. The soybean (Glycine max) genome is characterized by a large number of duplicated genes. However, the extent and mechanisms of functional divergence among these duplicated genes in soybean remain poorly understood. In this study, we revealed that 4 MYB genes (GmMYBA5, GmMYBA2, GmMYBA1, and Glyma.09g235000)-presumably generated by tandem duplication specifically in the Phaseoleae lineage-exhibited a stronger purifying selection in soybean compared to common bean (Phaseolus vulgaris). To gain insights into the diverse functions of these tandemly duplicated MYB genes in anthocyanin biosynthesis, we examined the expression, transcriptional activity, induced metabolites, and evolutionary history of these 4 MYB genes. Our data revealed that Glyma.09g235000 is a pseudogene, while the remaining 3 MYB genes exhibit strong transcriptional activation activity, promoting anthocyanin biosynthesis in different soybean tissues. GmMYBA5, GmMYBA2, and GmMYBA1 induced anthocyanin accumulation by upregulating the expression of anthocyanin pathway-related genes. Notably, GmMYBA5 showed a lower capacity for gene induction compared to GmMYBA2 and GmMYBA1. Metabolomics analysis further demonstrated that GmMYBA5 induced distinct anthocyanin accumulation in Nicotiana benthamiana leaves and soybean hairy roots compared to GmMYBA2 and GmMYBA1, suggesting their functional divergence leading to the accumulation of different metabolites accumulation following gene duplication. Together, our data provide evidence of functional divergence within the MYB gene cluster following tandem duplication, which sheds light on the potential evolutionary directions of gene duplications during legume evolution.

An omics strategy increasingly improves the discovery of genetic loci and genes for seed-coat color formation in soybean

The phenotypic color of seeds is a complex agronomic trait and has economic and biological significance. The genetic control and molecular regulation mechanisms have been extensively studied. Here, we used a multi-omics strategy to explore the color formation in soybean seeds at a big data scale. We identified 13 large quantitative trait loci (QTL) for color with bulk segregating analysis in recombinant inbreeding lines. GWAS analysis of colors and decomposed attributes in 763 germplasms revealed associated SNP sites perfectly falling in five major QTL, suggesting inherited regulation on color during natural selection. Further transcriptomics analysis before and after color accumulation revealed 182 differentially expression genes (DEGs) in the five QTL, including known genes CHS, MYB, and F3'H involved in pigment accumulation. More DEGs with consistently upregulation or downregulation were identified as shared regulatory genes for two or more color formations while some DEGs were only for a specific color formation. For example, five upregulated DEGs in QTL qSC-3 were in flavonoid biosynthesis responsible for black and brown seed. The DEG (Glyma.08G085400) was identified in the purple seed only, which encodes gibberellin 2-beta-dioxygenase in the metabolism of colorful terpenoids. The candidate genes are involved in flavonoid biosynthesis, transcription factor regulation, gibberellin and terpenoid metabolism, photosynthesis, ascorbate and aldarate metabolism, and lipid metabolism. Seven differentially expressed transcription factors were also speculated that may regulate color formation, including a known MYB. The finds expand QTL and gene candidates for color formation, which could guide to breed better cultivars with designed colors.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01414-z.

Long-read direct RNA sequencing reveals epigenetic regulation of chimeric gene-transposon transcripts in Arabidopsis thaliana

Transposable elements (TEs) are accumulated in both intergenic and intragenic regions in plant genomes. Intragenic TEs often act as regulatory elements of associated genes and are also co-transcribed with genes, generating chimeric TE-gene transcripts. Despite the potential impact on mRNA regulation and gene function, the prevalence and transcriptional regulation of TE-gene transcripts are poorly understood. By long-read direct RNA sequencing and a dedicated bioinformatics pipeline, ParasiTE, we investigated the transcription and RNA processing of TE-gene transcripts in Arabidopsis thaliana. We identified a global production of TE-gene transcripts in thousands of A. thaliana gene loci, with TE sequences often being associated with alternative transcription start sites or transcription termination sites. The epigenetic state of intragenic TEs affects RNAPII elongation and usage of alternative poly(A) signals within TE sequences, regulating alternative TE-gene isoform production. Co-transcription and inclusion of TE-derived sequences into gene transcripts impact regulation of RNA stability and environmental responses of some loci. Our study provides insights into TE-gene interactions that contributes to mRNA regulation, transcriptome diversity, and environmental responses in plants.

Morphological, Molecular, and Biochemical Characterization of a Unique Lentil (Lens culinaris Medik.) Genotype Showing Seed-Coat Color Anomalies Due to Altered Anthocyanin Pathway

This study reports the identification of a unique lentil (Lens culinaris Medik.) genotype L4717-NM, a natural mutant (NM) derived from a variety L4717, producing brown, black, and spotted seed-coat colored seeds in a single plant, generation after generation, in different frequencies. The genetic similarity of L4717 with that of L4717-NM expressing anomalous seed-coat color was established using 54 SSR markers. In addition, various biochemical parameters such as TPC (total phenolic content), TFC (total flavonoid content), DPPH (2,2-diphenyl-1-picrylhydrazyl), FRAP (ferric reducing antioxidant power), H2O2 (peroxide quantification), TCC (total carotenoids content), TAC (total anthocyanin content), and TAA (total ascorbic acid) were also studied in the seeds, sprouts, and seedlings of L4717, brown, black, and spotted seed-coat colored seeds. Stage-specific variations for the key biochemical parameters were recorded, and seedling stage was found the best for many parameters. Moreover, seeds with black seed coat showed better nutraceutical values for most of the studied traits. A highly significant (p ≤ 0.01) and positive correlation was observed between DPPH and TPC, TAA, TFC, etc., whereas, protein content showed a negative correlation with the other studied parameters. The seed coat is maternal tissue and we expect expression of seed-coat color as per the maternal genotype. However, such an anomalous seed-coat expression, which seems to probably be governed by some transposable element in the identified genotype, warrants more detailed studies involving exploitation of the anthocyanin pathway.

A pubescence color gene enhances tolerance to cold-induced seed cracking in yellow soybean

In yellow soybean, severe cold weather causes seed cracking on the dorsal side. Yellow soybeans carry the I or iⁱ allele of the I locus and have a yellow (I) or pigmented (iⁱ ) hilum. We previously isolated an additional allele, designated as Ic, of the I locus, and reported that yellow soybeans with the IcIc genotype may be tolerant to cold-induced seed cracking. The Ic allele by itself, however, does not confer high tolerance. The association of a pubescence color gene (T) with suppression of low-temperature-induced seed coat deterioration has been previously reported. In the present study, we tested whether T is effective for the suppression of cold-induced seed cracking using two pairs of near-isogenic lines for the T locus in the iⁱiⁱ or IcIc background. In both backgrounds, the cracked seed rate of the near-isogenic line with the TT genotype was significantly lower than that with the tt genotype, which indicates that T has an inhibitory effect on cold-induced seed cracking. Furthermore, we also showed that gene pyramiding of Ic and T can improve tolerance to cold-induced seed cracking. Our findings should aid the development of highly SC-tolerant cultivars in soybean breeding programs.

Dissecting the transcriptional regulation of proanthocyanidin and anthocyanin biosynthesis in soybean (Glycine max)

Proanthocyanidins (PAs), also known as condensed tannins, are plant natural products that are beneficial for human and livestock health. As one of the largest grown crops in the world, soybean (Glycine max) is widely used as human food and animal feed. Many cultivated soybeans with yellow seed coats lack PAs or anthocyanins, although some soybean cultivars have coloured seed coats that contain these compounds. Here, we analyse the transcriptional control of PA and anthocyanin biosynthesis in soybean. Ectopic expression of the transcription factors (TFs) GmTT2A, GmTT2B, GmMYB5A or R in soybean hairy roots induced the accumulation of PAs (primarily in phloem tissues) or anthocyanins and led to up-regulation of 1775, 856, 1411 and 1766 genes, respectively, several of which encode enzymes involved in PA biosynthesis. The genes regulated by GmTT2A and GmTT2B partially overlapped, suggesting conserved but potentially divergent roles for these two TFs in regulating PA accumulation in soybean. The two key enzymes anthocyanidin reductase and leucoanthocyanidin reductase were differentially upregulated, by GmTT2A/GmTT2B and GmMYB5A, respectively. Transgenic soybean plants overexpressing GmTT2B or MtLAP1 (a proven up-regulator of the upstream reactions for production of precursors for PA biosynthesis in legumes) showed increased accumulation of PAs and anthocyanins, respectively, associated with transcriptional reprogramming paralleling the RNA-seq data collected in soybean hairy roots. Collectively, our results show that engineered PA biosynthesis in soybean exhibits qualitative and spatial differences from the better-studied model systems Arabidopsis thaliana and Medicago truncatula, and suggest targets for engineering PAs in soybean plants.

MYB transcription factors GmMYBA2 and GmMYBR function in a feedback loop to control pigmentation of seed coat in soybean

Soybean has undergone extensive selection pressures for seed nutrient composition and seed color during domestication, but the major genetic loci controlling seed coat color have not been completely understood, and the transcriptional regulation relationship among the loci remains elusive. Here, two major regulators, GmMYBA2 and GmMYBR, were functionally characterized as an anthocyanin activator and repressor, respectively. Ectopic expression of GmMYBA2 in soybean hairy roots conferred the enhanced accumulation of delphinidin and cyanidin types of anthocyanins in W1t and w1T backgrounds, respectively, through activating anthocyanin biosynthetic genes in the reported loci. The seed coat pigmentation of GmMYBA2-overexpressing transgenic plants in the W1 background mimicked the imperfect black phenotype (W1/w1, i, R, t), suggesting that GmMYBA2 was responsible for the R locus. Molecular and biochemical analysis showed that GmMYBA2 interacted with GmTT8a to directly activate anthocyanin biosynthetic genes. GmMYBA2 and GmMYBR might form a feedback loop to fine-tune seed coat coloration, which was confirmed in transgenic soybeans. Both GmTT8a and GmMYBR that were activated by GmMYBA2 in turn enhanced and obstructed the formation of the GmMYBA2-GmTT8a module, respectively. The results revealed the sophisticated regulatory network underlying the soybean seed coat pigmentation loci and shed light on the understanding of the seed coat coloration and other seed inclusions.

Antioxidant Contributors in Seed, Seed Coat, and Cotyledon of γ-ray-Induced Soybean Mutant Lines with Different Seed Coat Colors

The development of soybean with high antioxidant activities for use in the food and cosmetics industries is a target of breeding programs. In soybean, antioxidants are associated with seed color, although the metabolic basis for seed coloration remains incompletely understood. We selected six γ-ray-induced mutant lines that exhibited black, partially black, brown, partially brown, or yellowish-white pigmentation in the seed coat. Antioxidant activity and contents of anthocyanins, flavan-3-ols, and isoflavones were evaluated in the seed coat and cotyledons. The lines with black or brown seeds showed the highest antioxidant activities. The cotyledons showed no significant differences in seed coat components or antioxidant activities among lines. Black and brown seed coat components showed the highest antioxidant activities. The black seed coat contained five anthocyanins, whereas seed coats of brown- and yellow-seeded lines entirely lacked anthocyanins. Both black and brown seeds were rich in flavan-3-ols, including catechin and epicatechin, which were the predominant antioxidant contributors in brown seeds. Isoflavone contents showed weaker correlations with antioxidant activity than anthocyanins and flavan-3-ols. These results demonstrated that antioxidant activities were determined by anthocyanins in black seeds and flavan-3-ols in brown and black seeds, whereas relatively low antioxidant activities in yellow seeds reflected their high isoflavone contents.

Overexpression of an ethylene-forming ACC oxidase (ACO) gene precedes the Minute Hilum seed coat phenotype in Glycine max

BACKGROUND

To elucidate features of seed development, we investigated the transcriptome of a soybean isoline from the germplasm collection that contained an introgressed allele known as minute hilum (mi) which confers a smaller hilum region where the seed attaches to the pod and also results in seed coat cracking surrounding the hilum region.

RESULTS

RNAs were extracted from immature seed from an extended hilum region (i.e., the hilum and a small ring of tissue surrounding the hilum in which the cracks form) at three different developmental stages:10-25, 25-50 and 50-100 mg seed fresh weight in two independent replicates for each stage. The transcriptomes of these samples from both the Clark isoline containing the mi allele (PI 547628, UC413, ii R t mi G), and its recurrent Clark 63 parent isoline (PI 548532, UC7, ii R T Mi g), which was used for six generations of backcrossing, were compared for differential expression of 88,648 Glyma models of the soybean genome Wm82.a2. The RNA sequence data obtained from the 12 cDNA libraries were subjected to padj value < 0.05 and at least two-fold expression differences to select with confidence genes differentially expressed in the hilum-containing tissue of the seed coat between the two lines. Glyma.09G008400 annotated as encoding an ethylene forming enzyme, ACC oxidase (ACO), was found to be highly overexpressed in the mi hilum region at 165 RPKMs (reads per kilobase per million mapped reads) compared to the standard line at just 0.03 RPKMs. Evidence of changes in expression of genes downstream of the ethylene pathway included those involved in auxin and gibberellin hormone action and extensive differences in expression of cell wall protein genes. These changes are postulated to determine the restricted hilum size and cracking phenotypes.

CONCLUSIONS

We present transcriptome and phenotypic evidence that substantially higher expression of an ethylene-forming ACO gene likely shifts hormone balance and sets in motion downstream changes resulting in a smaller hilum phenotype and the cracks observed in the minute hilum (mi) isoline as compared to its recurrent parent.

Mapping Quantitative Trait Loci for Soybean Seedling Shoot and Root Architecture Traits in an Inter-Specific Genetic Population

Wild soybean species (Glycine soja Siebold & Zucc.) comprise a unique resource to widen the genetic base of cultivated soybean [Glycine max (L.) Merr.] for various agronomic traits. An inter-specific mapping population derived from a cross of cultivar Williams 82 and PI 483460B, a wild soybean accession, was utilized for genetic characterization of root architecture traits. The objectives of this study were to identify and characterize quantitative trait loci (QTL) for seedling shoot and root architecture traits, as well as to determine additive/epistatic interaction effects of identified QTLs. A total of 16,469 single nucleotide polymorphisms (SNPs) developed for the Illumina beadchip genotyping platform were used to construct a high resolution genetic linkage map. Among the 11 putative QTLs identified, two significant QTLs on chromosome 7 were determined to be associated with total root length (RL) and root surface area (RSA) with favorable alleles from the wild soybean parent. These seedling root traits, RL (BARC_020495_04641 ~ BARC_023101_03769) and RSA (SNP02285 ~ SNP18129_Magellan), could be potential targets for introgression into cultivated soybean background to improve both tap and lateral roots. The RL QTL region harbors four candidate genes with higher expression in root tissues: Phosphofructokinase (Glyma.07g126400), Snf7 protein (Glyma.07g127300), unknown functional gene (Glyma.07g127900), and Leucine Rich-Repeat protein (Glyma.07g127100). The novel alleles inherited from the wild soybean accession could be used as molecular markers to improve root system architecture and productivity in elite soybean lines.

Transcriptomic Insights into Mechanisms of Early Seed Maturation in the Garden Pea (Pisum sativum L.)

The garden pea (Pisum sativum L.) is a legume crop of immense economic value. Extensive breeding has led to the emergence of numerous pea varieties, of which some are distinguished by accelerated development in various stages of ontogenesis. One such trait is rapid seed maturation, which, despite novel insights into the genetic control of seed development in legumes, remains poorly studied. This article presents an attempt to dissect mechanisms of early maturation in the pea line Sprint-2 by means of whole transcriptome RNA sequencing in two developmental stages. By using a de novo assembly approach, we have obtained a reference transcriptome of 25,756 non-redundant entries expressed in pea seeds at either 10 or 20 days after pollination. Differential expression in Sprint-2 seeds has affected 13,056 transcripts. A comparison of the two pea lines with a common maturation rate demonstrates that while at 10 days after pollination, Sprint-2 seeds show development retardation linked to intensive photosynthesis, morphogenesis, and cell division, and those at 20 days show a rapid onset of desiccation marked by the cessation of translation and cell anabolism and accumulation of dehydration-protective and -storage moieties. Further inspection of certain transcript functional categories, including the chromatin constituent, transcription regulation, protein turnover, and hormonal regulation, has revealed transcriptomic trends unique to specific stages and cultivars. Among other remarkable features, Sprint-2 demonstrated an enhanced expression of transposable element-associated open reading frames and an altered expression of major maturation regulators and DNA methyltransferase genes. To the best of our knowledge, this is the first comparative transcriptomic study in which the issue of the seed maturation rate is addressed.

Many Facets of Dynamic Plasticity in Plants

The evolutionary processes that transitioned plants to land-based habitats also incorporated a multiplicity of strategies to enhance resilience to the greater environmental variation encountered on land. The sensing of light, its quality, quantity, and duration, is central to plant survival and, as such, serves as a central network hub. Similarly, plants as sessile organisms that can encounter isolation must continually assess their reproductive options, requiring plasticity in propagation by self- and cross-pollination or asexual strategies. Irregular fluctuations and intermittent extremes in temperature, soil fertility, and moisture conditions have given impetus to genetic specializations for network resiliency, protein neofunctionalization, and internal mechanisms to accelerate their evolution. We review some of the current advancements made in understanding plant resiliency and phenotypic plasticity mechanisms. These mechanisms incorporate unusual nuclear-cytoplasmic interactions, various transposable element (TE) activities, and epigenetic plasticity of central gene networks that are broadly pleiotropic to influence resiliency phenotypes.

A Large-Scale Genome-Wide Association Analyses of Ethiopian Sorghum Landrace Collection Reveal Loci Associated With Important Traits

The eastern Africa region, Ethiopia and its surroundings, is considered as the center of origin and diversity for sorghum, and has contributed to global sorghum genetic improvement. The germplasm from this region harbors enormous genetic variation for various traits but little is known regarding the genetic architecture of most traits. Here, 1425 Ethiopian landrace accessions were phenotyped under field conditions for presence or absence of awns, panicle compactness and shape, panicle exsertion, pericarp color, glume cover, plant height and smut resistance under diverse environmental conditions in Ethiopia. In addition, F1 hybrids obtained from a subset of 1341 accessions crossed to an A1 cytoplasmic male sterile line, ATx623, were scored for fertility/sterility reactions. Subsequently, genotyping-by-sequencing generated a total of 879,407 SNPs from which 72,190 robust SNP markers were selected after stringent quality control (QC). Pairwise distance-based hierarchical clustering identified 11 distinct groups. Of the genotypes assigned to either one of the 11 sub-populations, 65% had high ancestry membership coefficient with the likelihood of more than 0.60 and the remaining 35% represented highly admixed accessions. A genome-wide association study (GWAS) identified loci and SNPs associated with aforementioned traits. GWAS based on compressed mixed linear model (CMLM) identified SNPs with significant association (FDR ≤ 0.05) to the different traits studied. The percentage of total phenotypic variation explained with significant SNPs across traits ranged from 2 to 43%. Candidate genes showing significant association with different traits were identified. The sorghum bHLH transcription factor, ABORTED MICROSPORES was identified as a strong candidate gene conditioning male fertility. Notably, sorghum CLAVATA1 receptor like kinase, known for regulation of plant growth, and the ETHYLENE RESPONSIVE TRANSCRIPTION FACTOR gene RAP2-7, known to suppress transition to flowering, were significantly associated with plant height. In addition, the YELLOW SEED1 like MYB transcription factor and TANNIN1 showed strong association with pericarp color validating previous observations. Overall, the genetic architecture of natural variation representing the complex Ethiopian sorghum germplasm was established. The study contributes to the characterization of genes and alleles controlling agronomic traits, and will serve as a source of markers for molecular breeding.

Contribution of transposable elements in the plant's genome

Plants maintain extensive growth flexibility under different environmental conditions, allowing them to continuously and rapidly adapt to alterations in their environment. A large portion of many plant genomes consists of transposable elements (TEs) that create new genetic variations within plant species. Different types of mutations may be created by TEs in plants. Many TEs can avoid the host's defense mechanisms and survive alterations in transposition activity, internal sequence and target site. Thus, plant genomes are expected to utilize a variety of mechanisms to tolerate TEs that are near or within genes. TEs affect the expression of not only nearby genes but also unlinked inserted genes. TEs can create new promoters, leading to novel expression patterns or alternative coding regions to generate alternate transcripts in plant species. TEs can also provide novel cis-acting regulatory elements that act as enhancers or inserts within original enhancers that are required for transcription. Thus, the regulation of plant gene expression is strongly managed by the insertion of TEs into nearby genes. TEs can also lead to chromatin modifications and thereby affect gene expression in plants. TEs are able to generate new genes and modify existing gene structures by duplicating, mobilizing and recombining gene fragments. They can also facilitate cellular functions by sharing their transposase-coding regions. Hence, TE insertions can not only act as simple mutagens but can also alter the elementary functions of the plant genome. Here, we review recent discoveries concerning the contribution of TEs to gene expression in plant genomes and discuss the different mechanisms by which TEs can affect plant gene expression and reduce host defense mechanisms.

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.

The endogenous transposable element Tgm9 is suitable for generating knockout mutants for functional analyses of soybean genes and genetic improvement in soybean

In soybean, variegated flowers can be caused by somatic excision of the CACTA-type transposable element Tgm9 from Intron 2 of the DFR2 gene encoding dihydroflavonol-4-reductase of the anthocyanin pigment biosynthetic pathway. DFR2 was mapped to the W4 locus, where the allele containing Tgm9 was termed w4-m. In this study we have demonstrated that previously identified morphological mutants (three chlorophyll deficient mutants, one male sterile-female fertile mutant, and three partial female sterile mutants) were caused by insertion of Tgm9 following its excision from DFR2. Analyses of Tgm9 insertion sites among 105 independent mutants demonstrated that Tgm9 hops to all 20 soybean chromosomes from its original location on Chromosome 17. Some genomic regions are prone to increased Tgm9-insertions. Tgm9 transposed over 25% of the time into exon or intron sequences. Tgm9 is therefore suitable for generating an indexed insertional mutant collection for functional analyses of most soybean genes. Furthermore, desirable Tgm9-induced stable knockout mutants can be utilized in generating improved traits for commercial soybean cultivars.

Genome-wide Association Mapping of Qualitatively Inherited Traits in a Germplasm Collection

Genome-wide association (GWA) has been used as a tool for dissecting the genetic architecture of quantitatively inherited traits. We demonstrate here that GWA can also be highly useful for detecting many major genes governing categorically defined phenotype variants that exist for qualitatively inherited traits in a germplasm collection. Genome-wide association mapping was applied to categorical phenotypic data available for 10 descriptive traits in a collection of ∼13,000 soybean [ (L.) Merr.] accessions that had been genotyped with a 50,000 single nucleotide polymorphism (SNP) chip. A GWA on a panel of accessions of this magnitude can offer substantial statistical power and mapping resolution, and we found that GWA mapping resulted in the identification of strong SNP signals for 24 classical genes as well as several heretofore unknown genes controlling the phenotypic variants in those traits. Because some of these genes had been cloned, we were able to show that the narrow GWA mapping SNP signal regions that we detected for the phenotypic variants had chromosomal bp spans that, with just one exception, overlapped the bp region of the cloned genes, despite local variation in SNP number and nonuniform SNP distribution in the chip set.

Relationship among color development, anthocyanin and pigment-related gene expression in 'Crimson Seedless' grapes treated with abscisic acid and sucrose

'Crimson Seedless' is one of the most important table grape varieties in Chile, but under certain environmental conditions, the fruit exhibits inadequate red color development, causing economic losses due to lower product quality. The use of plant growth regulators, such as abscisic acid (ABA) and ethylene, during development increases the anthocyanin content of the skin, improving the color of the berry. Recently, sucrose has been identified as a signaling molecule capable of regulating the expression of genes of the anthocyanin biosynthesis pathway. The aim of this study was to analyze the effect of application of ABA and/or sucrose on color development and their relationship with anthocyanin metabolism. Applications of ABA (400 ppm or 200 ppm) and/or sucrose (90 mM) were performed close to the véraison stage. During development and at harvest, quality attributes such as berry firmness, total soluble solids and titratable acidity were not affected by these treatments. Increased red color development was observed in fruits treated with ABA and/or sucrose, due to accumulation of anthocyanins. Fruits subjected to sucrose treatment showed higher levels of anthocyanins than untreated fruits but lower levels than fruits treated with ABA. Increased expression of genes involved in anthocyanin biosynthesis was observed in ABA- and sucrose-treated fruits compared to untreated fruits. Based on these findings, we demonstrated that sucrose improved fruit color development by increasing synthesis and accumulation of anthocyanins, thus allowing earlier harvests and improving table grape quality.

Mutations in Argonaute5 Illuminate Epistatic Interactions of the K1 and I Loci Leading to Saddle Seed Color Patterns in Glycine max

The soybean (Glycine max) seed coat has distinctive, genetically programmed patterns of pigmentation, and the recessive k1 mutation can epistatically overcome the dominant I and ii alleles, which inhibit seed color by producing small interfering RNAs (siRNAs) targeting chalcone synthase (CHS) mRNAs. Small RNA sequencing of dissected regions of immature seed coats demonstrated that CHS siRNA levels cause the patterns produced by the ii and ik alleles of the I locus, which restrict pigment to the hilum or saddle region of the seed coat, respectively. To identify the K1 locus, we compared RNA-seq data from dissected regions of two Clark isolines having similar saddle phenotypes mediated by CHS siRNAs but different genotypes (homozygous ik K1 versus homozygous ii k1). By examining differentially expressed genes, mapping information, and genome resequencing, we identified a 129-bp deletion in Glyma.11G190900 encoding Argonaute5 (AGO5), a member of the Argonaute family. Amplicon sequencing of several independent saddle pattern mutants from different genetic backgrounds revealed independent lesions affecting AGO5, thus establishing Glyma.11G190900 as the K1 locus. Nonfunctional AGO5 from k1 alleles leads to altered distributions of CHS siRNAs, thus explaining how the k1 mutation reverses the phenotype of the seed coat regions from yellow to pigmented, even in the presence of the normally dominant I or ii alleles.

Transcriptomic Analysis of Seed Coats in Yellow-Seeded Brassica napus Reveals Novel Genes That Influence Proanthocyanidin Biosynthesis

Yellow seeds are a favorable trait for Brassica crops breeding due to better quality than their black-seeded counterparts. Here, we compared the Brassica napus seed coat transcriptomes between yellow- and brown-seeded near-isogenic lines (Y-NIL and B-NIL) that were developed from the resynthesized yellow-seeded line No. 2127-17. A total of 4,974 differentially expressed genes (DEG) were identified during seed development, involving 3,128 up-regulated and 1,835 down-regulated genes in yellow seed coats. Phenylpropanoid and flavonoid biosynthesis pathways were enriched in down-regulated genes, whereas the top two pathways for up-regulated genes were plant-pathogen interaction and plant hormone signal transduction. Twelve biosynthetic genes and three regulatory genes involved in the flavonoid pathway exhibited similar expression patterns in seed coats during seed development, of which the down-regulation mainly contributed to the reduction of proanthocyanidins (PAs) in yellow seed coats, indicating that these genes associated with PA biosynthesis may be regulated by an unreported common regulator, possibly corresponding to the candidate for the dominant black-seeded gene D in the NILs. Three transcription factor (TF) genes, including one bHLH gene and two MYB-related genes that are located within the previous seed coat color quantitative trait locus (QTL) region on chromosome A09, also showed similar developmental expression patterns to the key PA biosynthetic genes and they might thus potentially involved participate in flavonoid biosynthesis regulation. Our study identified novel potential TFs involved in PAs accumulation and will provide pivotal information for identifying the candidate genes for seed coat color in B. napus.

Identification and Validation of Loci Governing Seed Coat Color by Combining Association Mapping and Bulk Segregation Analysis in Soybean

Soybean seed coat exists in a range of colors from yellow, green, brown, black, to bicolor. Classical genetic analysis suggested that soybean seed color was a moderately complex trait controlled by multi-loci. However, only a couple of loci could be detected using a single biparental segregating population. In this study, a combination of association mapping and bulk segregation analysis was employed to identify genes/loci governing this trait in soybean. A total of 14 loci, including nine novel and five previously reported ones, were identified using 176,065 coding SNPs selected from entire SNP dataset among 56 soybean accessions. Four of these loci were confirmed and further mapped using a biparental population developed from the cross between ZP95-5383 (yellow seed color) and NY279 (brown seed color), in which different seed coat colors were further dissected into simple trait pairs (green/yellow, green/black, green/brown, yellow/black, yellow/brown, and black/brown) by continuously developing residual heterozygous lines. By genotyping entire F2 population using flanking markers located in fine-mapping regions, the genetic basis of seed coat color was fully dissected and these four loci could explain all variations of seed colors in this population. These findings will be useful for map-based cloning of genes as well as marker-assisted breeding in soybean. This work also provides an alternative strategy for systematically isolating genes controlling relative complex trait by association analysis followed by biparental mapping.

Transposable element influences on gene expression in plants

Transposable elements (TEs) comprise a major portion of many plant genomes and bursts of TE movements cause novel genomic variation within species. In order to maintain proper gene function, plant genomes have evolved a variety of mechanisms to tolerate the presence of TEs within or near genes. Here, we review our understanding of the interactions between TEs and gene expression in plants by assessing three ways that transposons can influence gene expression. First, there is growing evidence that TE insertions within introns or untranslated regions of genes are often tolerated and have minimal impact on expression level or splicing. However, there are examples in which TE insertions within genes can result in aberrant or novel transcripts. Second, TEs can provide novel alternative promoters, which can lead to new expression patterns or original coding potential of an alternate transcript. Third, TE insertions near genes can influence regulation of gene expression through a variety of mechanisms. For example, TEs may provide novel cis-acting regulatory sites behaving as enhancers or insert within existing enhancers to influence transcript production. Alternatively, TEs may change chromatin modifications in regions near genes, which in turn can influence gene expression levels. Together, the interactions of genes and TEs provide abundant evidence for the role of TEs in changing basic functions within plant genomes beyond acting as latent genomic elements or as simple insertional mutagens. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Identification of Candidate Anthocyanin-Related Genes by Transcriptomic Analysis of 'Furongli' Plum (Prunus salicina Lindl.) during Fruit Ripening Using RNA-Seq

Anthocyanins are important pigments and are responsible for red coloration in plums. However, little is known about the molecular mechanisms underlying anthocyanin accumulation in plum fruits. In this study, the RNA-seq technique was used to analyze the transcriptomic changes during fruit ripening in the red-fleshed plum (Prunus salicina Lindl.) cultivar 'Furongli'. Over 161 million high-quality reads were assembled into 52,093 unigenes and 49.4% of these were annotated using public databases. Of these, 25,681 unigenes had significant hits to the sequences in the NCBI Nr database, 17,203 unigenes showed significant similarity to known proteins in the Swiss-Prot database and 5816 and 8585 unigenes had significant similarity to existing sequences in the Kyoto Encyclopedia of Genes and Genomes and the Cluster of Orthologous Groups databases, respectively. A total of 3548 unigenes were differentially expressed during fruit ripening and 119 of these were annotated as involved in "biosynthesis of other secondary metabolites." Biological pathway analysis and gene ontology term enrichment analysis revealed that 13 differentially expressed genes are involved in anthocyanin biosynthesis. Furthermore, transcription factors such as MYB and bHLH, which may control anthocyanin biosynthesis, were identified through coexpression analysis of transcription factors, and structural genes. Real-time qPCR analysis of candidate genes showed good correlation with the transcriptome data. These results contribute to our understanding of the molecular mechanisms underlying anthocyanin biosynthesis in plum flesh. The transcriptomic data generated in this study provide a basis for further studies of fruit ripening in plum.

CACTA-superfamily transposable element is inserted in MYB transcription factor gene of soybean line producing variegated seeds

The R gene of soybean, presumably encoding a MYB transcription factor, controls seed coat color. The gene consists of multiple alleles, R (black), r-m (black spots and (or) concentric streaks on brown seed), and r (brown seed). This study was conducted to determine the structure of the MYB transcription factor gene in a near-isogenic line (NIL) having r-m allele. PCR amplification of a fragment of the candidate gene Glyma.09G235100 generated a fragment of about 1 kb in the soybean cultivar Clark, whereas a fragment of about 14 kb in addition to fragments of 1 and 1.4 kb were produced in L72-2040, a Clark 63 NIL with the r-m allele. Clark 63 is a NIL of Clark with the rxp and Rps1 alleles. A DNA fragment of 13 060 bp was inserted in the intron of Glyma.09G235100 in L72-2040. The fragment had the CACTA motif at both ends, imperfect terminal inverted repeats (TIR), inverse repetition of short sequence motifs close to the 5' and 3' ends, and a duplication of three nucleotides at the site of integration, indicating that it belongs to a CACTA-superfamily transposable element. We designated the element as Tgm11. Overall nucleotide sequence, motifs of TIR, and subterminal repeats were similar to those of Tgm1 and Tgs1, suggesting that these elements comprise a family.

GREAM: A Web Server to Short-List Potentially Important Genomic Repeat Elements Based on Over-/Under-Representation in Specific Chromosomal Locations, Such as the Gene Neighborhoods, within or across 17 Mammalian Species

Background

Genome-wide repeat sequences, such as LINEs, SINEs and LTRs share a considerable part of the mammalian nuclear genomes. These repeat elements seem to be important for multiple functions including the regulation of transcription initiation, alternative splicing and DNA methylation. But it is not possible to study all repeats and, hence, it would help to short-list before exploring their potential functional significance via experimental studies and/or detailed in silico analyses.

Result

We developed the ‘Genomic Repeat Element Analyzer for Mammals’ (GREAM) for analysis, screening and selection of potentially important mammalian genomic repeats. This web-server offers many novel utilities. For example, this is the only tool that can reveal a categorized list of specific types of transposons, retro-transposons and other genome-wide repetitive elements that are statistically over-/under-represented in regions around a set of genes, such as those expressed differentially in a disease condition. The output displays the position and frequency of identified elements within the specified regions. In addition, GREAM offers two other types of analyses of genomic repeat sequences: a) enrichment within chromosomal region(s) of interest, and b) comparative distribution across the neighborhood of orthologous genes. GREAM successfully short-listed a repeat element (MER20) known to contain functional motifs. In other case studies, we could use GREAM to short-list repetitive elements in the azoospermia factor a (AZFa) region of the human Y chromosome and those around the genes associated with rat liver injury. GREAM could also identify five over-represented repeats around some of the human and mouse transcription factor coding genes that had conserved expression patterns across the two species.

Conclusion

GREAM has been developed to provide an impetus to research on the role of repetitive sequences in mammalian genomes by offering easy selection of more interesting repeats in various contexts/regions. GREAM is freely available at http://resource.ibab.ac.in/GREAM/.