A first de novo transcriptome assembly of feijoa (Acca sellowiana [Berg] Burret) reveals key genes involved in flavonoid biosynthesis.

Acca sellowiana [Berg] Burret, a cultivated fruit tree originating from South America, is gaining the attention of the nutraceutical and pharmaceutical industries due to their high content of flavonoids and other phenolic compounds in fruits, leaves, and flowers. Flavonoids are a diverse group of secondary metabolites with antioxidant, anti-inflammatory, and antimicrobial properties. They also play a crucial role in plant immune response. Despite their importance, the lack of research on A. sellowiana genomics and transcriptomics hinders a deeper understanding of the molecular mechanisms behind flavonoid biosynthesis and its regulation. Here, we de novo assembled and benchmarked 11 A. sellowiana transcriptomes from leaves and floral tissues at three developmental stages using high-throughput sequencing. We selected and annotated the best assembly according to commonly used metrics and databases. This reference transcriptome consisted of 221,649 nonredundant transcripts, of which 107,612 were functionally annotated. We then used this reference transcriptome to explore the expression profiling of key secondary metabolite genes. Transcripts from genes involved in the flavonoid and anthocyanin biosynthesis pathways were identified. We also identified 4068 putative transcription factors, with the most abundant families being bHLH, C2H2, NAC, MYB, and MYB-related. Transcript expression profiling revealed distinct patterns of gene expression during flower development. Particularly, we found 71 differentially expressed transcripts representing 14 enzymes of the flavonoid pathway, suggesting major changes in flavonoid accumulation across floral stages. Our findings will contribute to understanding the genetic basis of flavonoids and provide a foundation for further research and exploitation of the economic potential of this species.

Alfalfa (Medicago sativa L.) pho2 mutant plants hyperaccumulate phosphate.

In this article, we describe a set of novel alfalfa (Medicago sativa L.) plants that hyper-accumulate Phosphate ion (Pi) at levels 3- to 6-fold higher than wild-type. This alfalfa germplasm will have practical applications reclaiming Pi from contaminated or enriched soil or be used in conservation buffer strips to protect waterways from Pi run-off. Hyper-accumulating alfalfa plants were generated by targeted mutagenesis of PHOSPHATE2 (PHO2) using newly created CRISPR/Cas9 reagents and an improved mutant screening strategy. PHO2 encodes a ubiquitin conjugating E2 enzyme (UBC24) previously characterized in Arabidopsis thaliana, Medicago truncatula, and Oryza sativa. Mutations of PHO2 disrupt Pi homeostasis resulting in Pi hyper-accumulation. Successful CRISPR/Cas9 editing of PHO2 demonstrates that this is an efficient mutagenesis tool in alfalfa despite its complex autotetraploid genome structure. Arabidopsis and M. truncatula ortholog genes were used to identify PHO2 haplotypes in outcrossing tetraploid M. sativa with the aim of generating heritable mutations in both PHO2-like genes (PHO2-B and PHO2-C). After delivery of the reagent and regeneration from transformed leaf explants, plants with mutations in all haplotypes of PHO2-B and PHO2-C were identified. These plants were evaluated for morphology, Pi accumulation, heritable transmission of targeted mutations, segregation of mutant haplotypes and removal of T-DNA(s). The Agrobacterium-mediated transformation assay and gene editing reagents reported here were also evaluated for further optimization for future alfalfa functional genomic studies.

Potato improvement through genetic engineering.

Potato (Solanum tuberosum L.) is the third most important crop worldwide and a staple food for many people worldwide. Genetically, it poses many challenges for traditional breeding due to its autotetraploid nature and its tendency toward inbreeding depression. Breeding programs have focused on productivity, nutritional quality, and disease resistance. Some of these traits exist in wild potato relatives but their introgression into elite cultivars can take many years and, for traits such as pest resistance, their effect is often short-lasting. These problems can be addressed by genetic modification (GM) or gene editing (GE) and open a wide horizon for potato crop improvement. Current genetically modified and gene edited varieties include those with Colorado potato beetle and late blight resistance, reduction in acrylamide, and modified starch content. RNAi hairpin technology can be used to silence the haplo-alleles of multiple genes simultaneously, whereas optimization of newer gene editing technologies such as base and prime editing will facilitate the routine generation of advanced edits across the genome. These technologies will likely gain further relevance as increased target specificity and decreased off-target effects are demonstrated. In this Review, we discuss recent work related to these technologies in potato improvement.

De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes.

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

Multiple wheat genomes reveal global variation in modern breeding.

Advances in genomics have expedited the improvement of several agriculturally important crops but similar efforts in wheat (Triticum spp.) have been more challenging. This is largely owing to the size and complexity of the wheat genome1, and the lack of genome-assembly data for multiple wheat lines2,3. Here we generated ten chromosome pseudomolecule and five scaffold assemblies of hexaploid wheat to explore the genomic diversity among wheat lines from global breeding programs. Comparative analysis revealed extensive structural rearrangements, introgressions from wild relatives and differences in gene content resulting from complex breeding histories aimed at improving adaptation to diverse environments, grain yield and quality, and resistance to stresses4,5. We provide examples outlining the utility of these genomes, including a detailed multi-genome-derived nucleotide-binding leucine-rich repeat protein repertoire involved in disease resistance and the characterization of Sm16, a gene associated with insect resistance. These genome assemblies will provide a basis for functional gene discovery and breeding to deliver the next generation of modern wheat cultivars.

Wildfire effects on diversity and composition in soil bacterial communities.

In recent years, the Mediterranean area has witnessed an increase of both the frequency and severity of large fires, which appears to be intimately associated with climate and land use changes. To measure the impact of wildfires on living organisms, diverse indicators have been proposed. These indicators of fire severity traditionally rely on quantifying the damage caused to the vegetal component of ecosystems. However, the use of bacterial communities as severity indicators has received less attention. Here, we studied the differences between bacterial communities of three different Mediterranean ecosystems, two shrubby and one arboreal, two months after a large wildfire. Two levels of severity were compared to a control unburnt soil. The results showed that greater fire severity triggers a reduction in the diversity of soil bacterial communities. In high-severity fires, this reduction reached 40.6 and 58.6% of the control values for richness and Shannon's diversity, respectively. We also found that the greatest differences between communities could be attributed first to the severity of the fire, and second to the ecosystem from which they originated. Importantly, species of just five families of bacteria: Oxalobacteraceae, Micrococcaceae, Paenibacillaceae, Bacillaceae and Planococcaceae, became dominant in all three ecosystems. The average frequency increase for particular species was 100 times. However, due to random uncontrolled factors, the species that became dominant in each community were not always the same.

Heritable temporal gene expression patterns correlate with metabolomic seed content in developing hexaploid oat seed.

Oat ranks sixth in world cereal production and has a higher content of health-promoting compounds compared with other cereals. However, there is neither a robust oat reference genome nor transcriptome. Using deeply sequenced full-length mRNA libraries of oat cultivar Ogle-C, a de novo high-quality and comprehensive oat seed transcriptome was assembled. With this reference transcriptome and QuantSeq 3' mRNA sequencing, gene expression was quantified during seed development from 22 diverse lines across six time points. Transcript expression showed higher correlations between adjacent time points. Based on differentially expressed genes, we identified 22 major temporal co-expression (TCoE) patterns of gene expression and revealed enriched gene ontology biological processes. Within each TCoE set, highly correlated transcripts, putatively commonly affected by genetic background, were clustered and termed genetic co-expression (GCoE) sets. Seventeen of the 22 TCoE sets had GCoE sets with median heritabilities higher than 0.50, and these heritability estimates were much higher than that estimated from permutation analysis, with no divergence observed in cluster sizes between permutation and non-permutation analyses. Linear regression between 634 metabolites from mature seeds and the PC1 score of each of the GCoE sets showed significantly lower p-values than permutation analysis. Temporal expression patterns of oat avenanthramides and lipid biosynthetic genes were concordant with previous studies of avenanthramide biosynthetic enzyme activity and lipid accumulation. This study expands our understanding of physiological processes that occur during oat seed maturation and provides plant breeders the means to change oat seed composition through targeted manipulation of key pathways.

Dense genotyping-by-sequencing linkage maps of two Synthetic W7984×Opata reference populations provide insights into wheat structural diversity.

Wheat (Triticum aestivum) genetic maps are a key enabling tool for genetic studies. We used genotyping-by-sequencing-(GBS) derived markers to map recombinant inbred line (RIL) and doubled haploid (DH) populations from crosses of W7984 by Opata, and used the maps to explore features of recombination control. The RIL and DH populations, SynOpRIL and SynOpDH, were composed of 906 and 92 individuals, respectively. Two high-density genetic linkage framework maps were constructed of 2,842 and 2,961 cM, harboring 3,634 and 6,580 markers, respectively. Using imputation, we added 43,013 and 86,042 markers to the SynOpRIL and SynOpDH maps. We observed preferential recombination in telomeric regions and reduced recombination in pericentromeric regions. Recombination rates varied between subgenomes, with the D genomes of the two populations exhibiting the highest recombination rates of 0.26-0.27 cM/Mb. QTL mapping identified two additive and three epistatic loci associated with crossover number. Additionally, we used published POPSEQ data from SynOpDH to explore the structural variation in W7984 and Opata. We found that chromosome 5AS is missing from W7984. We also found 2,332 variations larger than 100 kb. Structural variants were more abundant in distal regions, and overlapped 9,196 genes. The two maps provide a resource for trait mapping and genomic-assisted breeding.

De Novo Transcriptome Assembly in Polyploid Species.

In the absence of a reference genome, the ultimate goal of a de novo transcriptome assembly is to accurately and comprehensively reconstruct the set of messenger RNA transcripts represented in the sample. Non-reference assembly of the transcriptome of polyploid species poses a particular challenge because of the presence of homeologs that are difficult to disentangle at the sequence level. This is especially true for hexaploid oats, which have three highly similar subgenomes, two of which are thought to be nearly identical. Under these circumstances, most software packages and established pipelines encounter difficulties in rendering an accurate transcriptome because they are typically developed, refined, and tested for diploid organisms. We present a protocol for transcriptome assembly in oats that can be extended both to other polyploids and species with highly duplicated genomes.

Genome Assembly of the Fungus Cochliobolus miyabeanus, and Transcriptome Analysis during Early Stages of Infection on American Wildrice (Zizania palustris L.).

The fungus Cochliobolus miyabeanus causes severe leaf spot disease on rice (Oryza sativa) and two North American specialty crops, American wildrice (Zizania palustris) and switchgrass (Panicum virgatum). Despite the importance of C. miyabeanus as a disease-causing agent in wildrice, little is known about either the mechanisms of pathogenicity or host defense responses. To start bridging these gaps, the genome of C. miyabeanus strain TG12bL2 was shotgun sequenced using Illumina technology. The genome assembly consists of 31.79 Mbp in 2,378 scaffolds with an N50 = 74,921. It contains 11,000 predicted genes of which 94.5% were annotated. Approximately 10% of total gene number is expected to be secreted. The C. miyabeanus genome is rich in carbohydrate active enzymes, and harbors 187 small secreted peptides (SSPs) and some fungal effector homologs. Detoxification systems were represented by a variety of enzymes that could offer protection against plant defense compounds. The non-ribosomal peptide synthetases and polyketide synthases (PKS) present were common to other Cochliobolus species. Additionally, the fungal transcriptome was analyzed at 48 hours after inoculation in planta. A total of 10,674 genes were found to be expressed, some of which are known to be involved in pathogenicity or response to host defenses including hydrophobins, cutinase, cell wall degrading enzymes, enzymes related to reactive oxygen species scavenging, PKS, detoxification systems, SSPs, and a known fungal effector. This work will facilitate future research on C. miyabeanus pathogen-associated molecular patterns and effectors, and in the identification of their corresponding wildrice defense mechanisms.

Subgenome-specific assembly of vitamin E biosynthesis genes and expression patterns during seed development provide insight into the evolution of oat genome.

Vitamin E is essential for humans and thus must be a component of a healthy diet. Among the cereal grains, hexaploid oats (Avena sativa L.) have high vitamin E content. To date, no gene sequences in the vitamin E biosynthesis pathway have been reported for oats. Using deep sequencing and orthology-guided assembly, coding sequences of genes for each step in vitamin E synthesis in oats were reconstructed, including resolution of the sequences of homeologs. Three homeologs, presumably representing each of the three oat subgenomes, were identified for the main steps of the pathway. Partial sequences, likely representing pseudogenes, were recovered in some instances as well. Pairwise comparisons among homeologs revealed that two of the three putative subgenome-specific homeologs are almost identical for each gene. Synonymous substitution rates indicate the time of divergence of the two more similar subgenomes from the distinct one at 7.9-8.7 MYA, and a divergence between the similar subgenomes from a common ancestor 1.1 MYA. A new proposed evolutionary model for hexaploid oat formation is discussed. Homeolog-specific gene expression was quantified during oat seed development and compared with vitamin E accumulation. Homeolog expression largely appears to be similar for most of genes; however, for some genes, homoeolog-specific transcriptional bias was observed. The expression of HPPD, as well as certain homoeologs of VTE2 and VTE4, is highly correlated with seed vitamin E accumulation. Our findings expand our understanding of oat genome evolution and will assist efforts to modify vitamin E content and composition in oats.

MicroRNA Maturation and MicroRNA Target Gene Expression Regulation Are Severely Disrupted in Soybean dicer-like1 Double Mutants.

Small nonprotein-coding microRNAs (miRNAs) are present in most eukaryotes and are central effectors of RNA silencing-mediated mechanisms for gene expression regulation. In plants, DICER-LIKE1 (DCL1) is the founding member of a highly conserved family of RNase III-like endonucleases that function as core machinery proteins to process hairpin-like precursor transcripts into mature miRNAs, small regulatory RNAs, 21-22 nucleotides in length. Zinc finger nucleases (ZFNs) were used to generate single and double-mutants of putative soybean DCL1 homologs, DCL1a and DCL1b, to confirm their functional role(s) in the soybean miRNA pathway. Neither DCL1 single mutant, dcl1a or dcl1b plants, exhibited a pronounced morphological or molecular phenotype. However, the dcl1a/dcl1b double mutant expressed a strong morphological phenotype, characterized by reduced seed size and aborted seedling development, in addition to defective miRNA precursor transcript processing efficiency and deregulated miRNA target gene expression. Together, these findings indicate that the two soybean DCL1 paralogs, DCL1a and DCL1b, largely play functionally redundant roles in the miRNA pathway and are essential for normal plant development.

Overexpression of AtDREB1D transcription factor improves drought tolerance in soybean.

Drought is one of the major abiotic stresses that affect productivity in soybean (Glycine max L.) Several genes induced by drought stress include functional genes and regulatory transcription factors. The Arabidopsis thaliana DREB1D transcription factor driven by the constitutive and ABA-inducible promoters was introduced into soybean through Agrobacterium tumefaciens-mediated gene transfer. Several transgenic lines were generated and molecular analysis was performed to confirm transgene integration. Transgenic plants with an ABA-inducible promoter showed a 1.5- to two-fold increase of transgene expression under severe stress conditions. Under well-watered conditions, transgenic plants with constitutive and ABA-inducible promoters showed reduced total leaf area and shoot biomass compared to non-transgenic plants. No significant differences in root length or root biomass were observed between transgenic and non-transgenic plants under non-stress conditions. When subjected to gradual water deficit, transgenic plants maintained higher relative water content because the transgenic lines used water more slowly as a result of reduced total leaf area. This caused them to wilt slower than non-transgenic plants. Transgenic plants showed differential drought tolerance responses with a significantly higher survival rate compared to non-transgenic plants when subjected to comparable severe water-deficit conditions. Moreover, the transgenic plants also showed improved drought tolerance by maintaining 17-24 % greater leaf cell membrane stability compared to non-transgenic plants. The results demonstrate the feasibility of engineering soybean for enhanced drought tolerance by expressing stress-responsive genes.

Analysis and annotation of the hexaploid oat seed transcriptome.

BACKGROUND: Next generation sequencing provides new opportunities to explore transcriptomes. However, challenges remain for accurate differentiation of homoeoalleles and paralogs, particularly in polyploid organisms with no supporting genome sequence. In this study, RNA-Seq was employed to generate and characterize the first gene expression atlas for hexaploid oat. RESULTS: The software packages Trinity and Oases were used to produce a transcript assembly from nearly 134 million 100-bp paired-end reads from developing oat seeds. Based on the quality-parameters employed, Oases assemblies were superior. The Oases 67-kmer assembly, denoted dnOST (de novo Oat Seed Transcriptome), is over 55 million nucleotides in length and the average transcript length is 1,043 nucleotides. The 74.8× sequencing depth was adequate to differentiate a large proportion of putative homoeoalleles and paralogs. To assess the robustness of dnOST, we successfully identified gene transcripts associated with the biosynthetic pathways of three compounds with health-promoting properties (avenanthramides, tocols, ß-glucans), and quantified their expression. CONCLUSIONS: To our knowledge, this study provides the first direct performance comparison between two major assemblers in a polyploid organism. The workflow we developed provides a useful guide for comparable analyses in other organisms. The transcript assembly developed here is a major advance. It expands the number of oat ESTs 3-fold, and constitutes the first comprehensive transcriptome study in oat. This resource will be a useful new tool both for analysis of genes relevant to nutritional enhancement of oat, and for improvement of this crop in general.

Evaluation of candidate reference genes for normalization of quantitative RT-PCR in soybean tissues under various abiotic stress conditions.

Quantitative RT-PCR can be a very sensitive and powerful technique for measuring differential gene expression. Changes in gene expression induced by abiotic stresses are complex and multifaceted, which make determining stably expressed genes for data normalization difficult. To identify the most suitable reference genes for abiotic stress studies in soybean, 13 candidate genes collected from literature were evaluated for stability of expression under dehydration, high salinity, cold and ABA (abscisic acid) treatments using delta CT and geNorm approaches. Validation of reference genes indicated that the best reference genes are tissue- and stress-dependent. With respect to dehydration treatment, the Fbox/ABC, Fbox/60s gene pairs were found to have the highest expression stability in the root and shoot tissues of soybean seedlings, respectively. Fbox and 60s genes are the most suitable reference genes across dehydrated root and shoot tissues. Under salt stress the ELF1b/IDE and Fbox/ELF1b are the most stably expressed gene pairs in roots and shoots, respectively, while 60s/Fbox is the best gene pair in both tissues. For studying cold stress in roots or shoots, IDE/60s and Fbox/Act27 are good reference gene pairs, respectively. With regard to gene expression analysis under ABA treatment in either roots, shoots or across these tissues, 60s/ELF1b, ELF1b/Fbox and 60s/ELF1b are the most suitable reference genes, respectively. The expression of ELF1b/60s, 60s/Fbox and 60s/Fbox genes was most stable in roots, shoots and both tissues, respectively, under various stresses studied. Among the genes tested, 60s was found to be the best reference gene in different tissues and under various stress conditions. The highly ranked reference genes identified from this study were proved to be capable of detecting subtle differences in expression rates that otherwise would be missed if a less stable reference gene was used.

Major locus and other novel additive and epistatic loci involved in modulation of isoflavone concentration in soybean seeds.

Seeds of soybean [Glycine max (L.) Merr.] accumulate more isoflavones than any tissue of any plant species. In other plant parts, isoflavones are usually released to counteract the effects of various biotic and abiotic stresses. Because of the benefits to the plant and positive implications that consumption may have on human health, increasing isoflavones is a goal of many soybean breeding programs. However, altering isoflavone levels through marker-assisted selection (MAS) has been impractical due to the small and often environmentally variable contributions that each individual quantitative trait locus (QTL) has on total isoflavones. In this study, we developed a Magellan × PI 437654 F(7)-RIL population to construct a highly saturated non-redundant linkage map that encompassed 451 SNP and SSR molecular markers and used it to locate genomic regions that govern accumulation of isoflavones in the seeds of soybean. Five QTLs were found that contribute to the concentration of isoflavones, having single or multiple additive effects on isoflavone component traits. We also validated a major locus which alone accounted for up to 10% of the phenotypic variance for glycitein, and 35-37% for genistein, daidzein and the sum of all three soybean isoflavones. This QTL was consistently associated with increased concentration of isoflavones across different locations, years and crosses. It was the most important QTL in terms of net increased amounts of all isoflavone forms. Our results suggest that this locus would be an excellent candidate to target for MAS. Also, several minor QTLs were identified that interacted in an additive-by-additive epistatic manner, to increase isoflavone concentration.

Intricate environment-modulated genetic networks control isoflavone accumulation in soybean seeds.

BACKGROUND: Soybean (Glycine max [L] Merr.) seed isoflavones have long been considered a desirable trait to target in selection programs for their contribution to human health and plant defense systems. However, attempts to modify seed isoflavone contents have not always produced the expected results because their genetic basis is polygenic and complex. Undoubtedly, the extreme variability that seed isoflavones display over environments has obscured our understanding of the genetics involved. RESULTS: In this study, a mapping population of RILs with three replicates was analyzed in four different environments (two locations over two years). We found a total of thirty-five main-effect genomic regions and many epistatic interactions controlling genistein, daidzein, glycitein and total isoflavone accumulation in seeds. The use of distinct environments permitted detection of a great number of environment-modulated and minor-effect QTL. Our findings suggest that isoflavone seed concentration is controlled by a complex network of multiple minor-effect loci interconnected by a dense epistatic map of interactions. The magnitude and significance of the effects of many of the nodes and connections in the network varied depending on the environmental conditions. In an attempt to unravel the genetic architecture underlying the traits studied, we searched on a genome-wide scale for genomic regions homologous to the most important identified isoflavone biosynthetic genes. We identified putative candidate genes for several of the main-effect and epistatic QTL and for QTL reported by other groups. CONCLUSIONS: To better understand the underlying genetics of isoflavone accumulation, we performed a large scale analysis to identify genomic regions associated with isoflavone concentrations. We not only identified a number of such regions, but also found that they can interact with one another and with the environment to form a complex adaptable network controlling seed isoflavone levels. We also found putative candidate genes in several regions and overall we advanced the knowledge of the genetics underlying isoflavone synthesis.

Differential expression of isoflavone biosynthetic genes in soybean during water deficits.

Numerous environmental factors influence isoflavone accumulation and have long hampered their genetic dissection. Temperature and water regimes are two of the most significant abiotic factors. However, while the effects of temperature have been widely studied, little is known about how water scarcity might affect isoflavone concentration in seeds. Studies have shown that accumulation of isoflavones is promoted by well-watered conditions, but the molecular basis remains elusive. The length and severity of the water stress required to induce changes are also still unknown. In the present work, several intensities of water stress were evaluated at various critical stages for soybean [Glycine max (L.) Merr.] seed development, in both field and controlled environments. The results suggested that only long-term progressive drought, spanning most of the seed developmental stages, significantly decreased isoflavone content in seeds. The reduction is proportional to the intensity of the stress and appears to occur in a genotype-dependent manner. However, regardless of water regime, isoflavone compounds were mainly accumulated in the later seed developmental stages. Transcripts of the most important genes for isoflavone biosynthesis were also quantified from samples collected at key seed developmental stages under well-watered and long-term water deficit conditions. Expression of CHS7, CHS8 and IFS2 correlated with isoflavone accumulation under well-watered conditions. Interestingly, we found that the two isoflavone synthase genes in soybean (IFS1 and IFS2) showed different patterns of expression. The abundance of IFS1 transcripts was maintained at a constant rate, whereas IFS2 was down-regulated and highly correlated with isoflavone accumulation under both water deficit and well-watered conditions, suggesting IFS2 as a main contributor to isoflavone diminution under drought.

Genetic control of soybean seed isoflavone content: importance of statistical model and epistasis in complex traits.

A major objective for geneticists is to decipher genetic architecture of traits associated with agronomic importance. However, a majority of such traits are complex, and their genetic dissection has been traditionally hampered not only by the number of minor-effect quantitative trait loci (QTL) but also by genome-wide interacting loci with little or no individual effect. Soybean (Glycine max [L.] Merr.) seed isoflavonoids display a broad range of variation, even in genetically stabilized lines that grow in a fixed environment, because their synthesis and accumulation are affected by many biotic and abiotic factors. Due to this complexity, isoflavone QTL mapping has often produced conflicting results especially with variable growing conditions. Herein, we comparatively mapped soybean seed isoflavones genistein, daidzein, and glycitein by using several of the most commonly used mapping approaches: interval mapping, composite interval mapping, multiple interval mapping and a mixed-model based composite interval mapping. In total, 26 QTLs, including many novel regions, were found bearing additive main effects in a population of RILs derived from the cross between Essex and PI 437654. Our comparative approach demonstrates that statistical mapping methodologies are crucial for QTL discovery in complex traits. Despite a previous understanding of the influence of additive QTL on isoflavone production, the role of epistasis is not well established. Results indicate that epistasis, although largely dependent on the environment, is a very important genetic component underlying seed isoflavone content, and suggest epistasis as a key factor causing the observed phenotypic variability of these traits in diverse environments.