Multiple loci and genetic interactions involving flowering time genes regulate stem branching among natural variants of Arabidopsis

Flowering time: From physiology, through genetics to mechanism

Plant species have evolved different requirements for environmental/endogenous cues to induce flowering. Originally, these varying requirements were thought to reflect the action of different molecular mechanisms. Thinking changed when genetic and molecular analysis in Arabidopsis thaliana revealed that a network of environmental and endogenous signaling input pathways converge to regulate a common set of "floral pathway integrators." Variation in the predominance of the different input pathways within a network can generate the diversity of requirements observed in different species. Many genes identified by flowering time mutants were found to encode general developmental and gene regulators, with their targets having a specific flowering function. Studies of natural variation in flowering were more successful at identifying genes acting as nodes in the network central to adaptation and domestication. Attention has now turned to mechanistic dissection of flowering time gene function and how that has changed during adaptation. This will inform breeding strategies for climate-proof crops and help define which genes act as critical flowering nodes in many other species.

Involvement of CONSTANS-like Proteins in Plant Flowering and Abiotic Stress Response

The process of flowering in plants is a pivotal stage in their life cycle, and the CONSTANS-like (COL) protein family, known for its photoperiod sensing ability, plays a crucial role in regulating plant flowering. Over the past two decades, homologous genes of COL have been identified in various plant species, leading to significant advancements in comprehending their involvement in the flowering pathway and response to abiotic stress. This article presents novel research progress on the structural aspects of COL proteins and their regulatory patterns within transcription complexes. Additionally, we reviewed recent information about their participation in flowering and abiotic stress response, aiming to provide a more comprehensive understanding of the functions of COL proteins.

Artificial selection reveals complex genetic architecture of shoot branching and its response to nitrate supply in Arabidopsis

Quantitative traits may be controlled by many loci, many alleles at each locus, and subject to genotype-by-environment interactions, making them difficult to map. One example of such a complex trait is shoot branching in the model plant Arabidopsis, and its plasticity in response to nitrate. Here, we use artificial selection under contrasting nitrate supplies to dissect the genetic architecture of this complex trait, where loci identified by association mapping failed to explain heritability estimates. We found a consistent response to selection for high branching, with correlated responses in other traits such as plasticity and flowering time. Genome-wide scans for selection and simulations suggest that at least tens of loci control this trait, with a distinct genetic architecture between low and high nitrate treatments. While signals of selection could be detected in the populations selected for high branching on low nitrate, there was very little overlap in the regions selected in three independent populations. Thus the regulatory network controlling shoot branching can be tuned in different ways to give similar phenotypes.

Genetic Mapping of Genotype-by-Ploidy Effects in Arabidopsis thaliana

Plants can express different phenotypic responses following polyploidization, but ploidy-dependent phenotypic variation has so far not been assigned to specific genetic factors. To map such effects, segregating populations at different ploidy levels are required. The availability of an efficient haploid inducer line in Arabidopsis thaliana allows for the rapid development of large populations of segregating haploid offspring. Because Arabidopsis haploids can be self-fertilised to give rise to homozygous doubled haploids, the same genotypes can be phenotyped at both the haploid and diploid ploidy level. Here, we compared the phenotypes of recombinant haploid and diploid offspring derived from a cross between two late flowering accessions to map genotype × ploidy (G × P) interactions. Ploidy-specific quantitative trait loci (QTLs) were detected at both ploidy levels. This implies that mapping power will increase when phenotypic measurements of monoploids are included in QTL analyses. A multi-trait analysis further revealed pleiotropic effects for a number of the ploidy-specific QTLs as well as opposite effects at different ploidy levels for general QTLs. Taken together, we provide evidence of genetic variation between different Arabidopsis accessions being causal for dissimilarities in phenotypic responses to altered ploidy levels, revealing a G × P effect. Additionally, by investigating a population derived from late flowering accessions, we revealed a major vernalisation-specific QTL for variation in flowering time, countering the historical bias of research in early flowering accessions.

Natural variation of Dt2 determines branching in soybean

Shoot branching is fundamentally important in determining soybean yield. Here, through genome-wide association study, we identify one predominant association locus on chromosome 18 that confers soybean branch number in the natural population. Further analyses determine that Dt2 is the corresponding gene and the natural variations in Dt2 result in significant differential transcriptional levels between the two major haplotypes. Functional characterization reveals that Dt2 interacts with GmAgl22 and GmSoc1a to physically bind to the promoters of GmAp1a and GmAp1d and to activate their transcription. Population genetic investigation show that the genetic differentiation of Dt2 display significant geographic structure. Our study provides a predominant gene for soybean branch number and may facilitate the breeding of high-yield soybean varieties.

A MADS-box gene with similarity to FLC is induced by cold and correlated with epigenetic changes to control budbreak in kiwifruit

Temperate perennials require exposure to chilling temperatures to resume growth in the following spring. Growth and dormancy cycles are controlled by complex genetic regulatory networks and are governed by epigenetic mechanisms, but the specific genes and mechanisms remain poorly understood. To understand how seasonal changes and chilling regulate dormancy and growth in the woody perennial vine kiwifruit (Ac, Actinidia chinensis), a transcriptome study of kiwifruit buds in the field and controlled conditions was performed. A MADS-box gene with homology to Arabidopsis FLOWERING LOCUS C (FLC) was identified and characterized. Elevated expression of AcFLC-like (AcFLCL) was detected during bud dormancy and chilling. A long noncoding (lnc) antisense transcript with an expression pattern opposite to AcFLCL and shorter sense noncoding RNAs were identified. Chilling induced an increase in trimethylation of lysine-4 of histone H3 (H3K4me3) in the 5' end of the gene, indicating multiple layers of epigenetic regulation in response to cold. Overexpression of AcFLCL in kiwifruit gave rise to plants with earlier budbreak, whilst gene editing using CRISPR-Cas9 resulted in transgenic lines with substantially delayed budbreak, suggesting a role in activation of growth. These results have implications for the future management and breeding of perennials for resilience to changing climate.

Natural variation in autumn expression is the major adaptive determinant distinguishing Arabidopsis FLC haplotypes

In Arabidopsis thaliana, winter is registered during vernalization through the temperature-dependent repression and epigenetic silencing of floral repressor FLOWERING LOCUS C (FLC). Natural Arabidopsis accessions show considerable variation in vernalization. However, which aspect of the FLC repression mechanism is most important for adaptation to different environments is unclear. By analysing FLC dynamics in natural variants and mutants throughout winter in three field sites, we find that autumnal FLC expression, rather than epigenetic silencing, is the major variable conferred by the distinct Arabidopsis FLChaplotypes. This variation influences flowering responses of Arabidopsis accessions resulting in an interplay between promotion and delay of flowering in different climates to balance survival and, through a post-vernalization effect, reproductive output. These data reveal how expression variation through non-coding cis variation at FLC has enabled Arabidopsis accessions to adapt to different climatic conditions and year-on-year fluctuations.

Genotype networks of 80 quantitative Arabidopsis thaliana phenotypes reveal phenotypic evolvability despite pervasive epistasis

We study the genotype-phenotype maps of 80 quantitative phenotypes in the model plant Arabidopsis thaliana, by representing the genotypes affecting each phenotype as a genotype network. In such a network, each vertex or node corresponds to an individual's genotype at all those genomic loci that affect a given phenotype. Two vertices are connected by an edge if the associated genotypes differ in exactly one nucleotide. The 80 genotype networks we analyze are based on data from genome-wide association studies of 199 A. thaliana accessions. They form connected graphs whose topography differs substantially among phenotypes. We focus our analysis on the incidence of epistasis (non-additive interactions among mutations) because a high incidence of epistasis can reduce the accessibility of evolutionary paths towards high or low phenotypic values. We find epistatic interactions in 67 phenotypes, and in 51 phenotypes every pairwise mutant interaction is epistatic. Moreover, we find phenotype-specific differences in the fraction of accessible mutational paths to maximum phenotypic values. However, even though epistasis affects the accessibility of maximum phenotypic values, the relationships between genotypic and phenotypic change of our analyzed phenotypes are sufficiently smooth that some evolutionary paths remain accessible for most phenotypes, even where epistasis is pervasive. The genotype network representation we use can complement existing approaches to understand the genetic architecture of polygenic traits in many different organisms.

Epistasis Detection and Modeling for Genomic Selection in Cowpea (Vigna unguiculata L. Walp.)

Genetic architecture reflects the pattern of effects and interaction of genes underlying phenotypic variation. Most mapping and breeding approaches generally consider the additive part of variation but offer limited knowledge on the benefits of epistasis which explains in part the variation observed in traits. In this study, the cowpea multiparent advanced generation inter-cross (MAGIC) population was used to characterize the epistatic genetic architecture of flowering time, maturity, and seed size. In addition, consideration for epistatic genetic architecture in genomic-enabled breeding (GEB) was investigated using parametric, semi-parametric, and non-parametric genomic selection (GS) models. Our results showed that large and moderate effect-sized two-way epistatic interactions underlie the traits examined. Flowering time QTL colocalized with cowpea putative orthologs of Arabidopsis thaliana and Glycine max genes like PHYTOCLOCK1 (PCL1 [Vigun11g157600]) and PHYTOCHROME A (PHY A [Vigun01g205500]). Flowering time adaptation to long and short photoperiod was found to be controlled by distinct and common main and epistatic loci. Parametric and semi-parametric GS models outperformed non-parametric GS model, while using known quantitative trait nucleotide(s) (QTNs) as fixed effects improved prediction accuracy when traits were controlled by large effect loci. In general, our study demonstrated that prior understanding of the genetic architecture of a trait can help make informed decisions in GEB.

Altered Expression of an FT Cluster Underlies a Major Locus Controlling Domestication-Related Changes to Chickpea Phenology and Growth Habit

Flowering time is a key trait in breeding and crop evolution, due to its importance for adaptation to different environments and for yield. In the particular case of chickpea, selection for early phenology was essential for the successful transition of this species from a winter to a summer crop. Here, we used genetic and expression analyses in two different inbred populations to examine the genetic control of domestication-related differences in flowering time and growth habit between domesticated chickpea and its wild progenitor Cicer reticulatum. A single major quantitative trait locus for flowering time under short-day conditions [Days To Flower (DTF)3A] was mapped to a 59-gene interval on chromosome three containing a cluster of three FT genes, which collectively showed upregulated expression in domesticated relative to wild parent lines. An equally strong association with growth habit suggests a pleiotropic effect of the region on both traits. These results indicate the likely molecular explanation for the characteristic early flowering of domesticated chickpea, and the previously described growth habit locus Hg. More generally, they point to de-repression of this specific gene cluster as a conserved mechanism for achieving adaptive early phenology in temperate legumes.

Seed traits are pleiotropically regulated by the flowering time gene PERPETUAL FLOWERING 1 (PEP1) in the perennial Arabis alpina

The life cycles of plants are characterized by two major life history transitions-germination and the initiation of flowering-the timing of which are important determinants of fitness. Unlike annuals, which make the transition from the vegetative to reproductive phase only once, perennials iterate reproduction in successive years. The floral repressor PERPETUAL FLOWERING 1 (PEP1), an ortholog of FLOWERING LOCUS C, in the alpine perennial Arabis alpina ensures the continuation of vegetative growth after flowering and thereby restricts the duration of the flowering episode. We performed greenhouse and garden experiments to compare flowering phenology, fecundity and seed traits between A. alpina accessions that have a functional PEP1 allele and flower seasonally and pep1 mutants and accessions that carry lesions in PEP1 and flower perpetually. In the garden, perpetual genotypes flower asynchronously and show higher winter mortality than seasonal ones. PEP1 also pleiotropically regulates seed dormancy and longevity in a way that is functionally divergent from FLC. Seeds from perpetual genotypes have shallow dormancy and reduced longevity regardless of whether they after-ripened in plants grown in the greenhouse or in the experimental garden. These results suggest that perpetual genotypes have higher mortality during winter but compensate by showing higher seedling establishment. Differences in seed traits between seasonal and perpetual genotypes are also coupled with differences in hormone sensitivity and expression of genes involved in hormonal pathways. Our study highlights the existence of pleiotropic regulation of seed traits by hub developmental regulators such as PEP1, suggesting that seed and flowering traits in perennial plants might be optimized in a coordinated fashion.

Assessment of heterosis in two Arabidopsis thaliana common-reference mapping populations

Hybrid vigour, or heterosis, has been of tremendous importance in agriculture for the improvement of both crops and livestock. Notwithstanding large efforts to study the phenomenon of heterosis in the last decades, the identification of common molecular mechanisms underlying hybrid vigour remain rare. Here, we conducted a systematic survey of the degree of heterosis in Arabidopsis thaliana hybrids. For this purpose, two overlapping Arabidopsis hybrid populations were generated by crossing a large collection of naturally occurring accessions to two common reference lines. In these Arabidopsis hybrid populations the range of heterosis for several developmental and yield related traits was examined, and the relationship between them was studied. The traits under study were projected leaf area at 17 days after sowing, flowering time, height of the main inflorescence, number of side branches from the main stem or from the rosette base, total seed yield, seed weight, seed size and the estimated number of seeds per plant. Predominantly positive heterosis was observed for leaf area and height of the main inflorescence, whereas mainly negative heterosis was observed for rosette branching. For the other traits both positive and negative heterosis was observed in roughly equal amounts. For flowering time and seed size only low levels of heterosis were detected. In general the observed heterosis levels were highly trait specific. Furthermore, no correlation was observed between heterosis levels and the genetic distance between the parental lines. Since all selected lines were a part of the Arabidopsis genome wide association (GWA) mapping panel, a genetic mapping approach was applied to identify possible regions harbouring genetic factors causal for heterosis, with separate calculations for additive and dominance effects. Our study showed that the genetic mechanisms underlying heterosis were highly trait specific in our hybrid populations and greatly depended on the genetic background, confirming the elusive character of heterosis.

Intraspecific variation of residual heterozygosity and its utility for quantitative genetic studies in maize

BACKGROUND

Residual heterozygosity (RH) in advanced inbred lines of plants benefits quantitative trait locus (QTL) mapping studies. However, knowledge of factors affecting the genome-wide distribution of RH remains limited.

RESULTS

A set of 2196 heterogeneous inbred family (HIF) maize lines derived from 12 recombinant inbred line (RIL) populations was genotyped using the Maize50K SNP chip. A total of 18,615 unique RH intervals were identified, ranging from 505 to 2095 intervals per population, with average maize genome coverage of 94.8%. Across all populations, there were 8.6 RH intervals per HIF line on average, ranging from 1.8 to 14 intervals; the average size of an RH interval was approximately 58.7 Mb, ranging from 7.2 to 74.1 Mb. A given RH region was present in an average of 5 different individuals within a population. Seven RH hotspots, where RH segments were enriched in the genome, were found to be subject to selection during population development. The RH patterns varied significantly across populations, presumably reflecting differences in the genetic background of each population, and 8 QTLs were found to affect heterozygosity levels in the RH hotspots. The potential use of this HIF library for the fine mapping of QTLs was assessed based on publicly available QTL information, achieving a ≤ 1 Mb resolution on average.

CONCLUSION

The examined library of HIF lines offers insight into the RH landscape and its intraspecific variation and provides a useful resource for the QTL cloning of important agronomic traits in maize.

The genetics underlying natural variation of plant-plant interactions, a beloved but forgotten member of the family of biotic interactions

Despite the importance of plant-plant interactions on crop yield and plant community dynamics, our understanding of the genetic and molecular bases underlying natural variation of plant-plant interactions is largely limited in comparison with other types of biotic interactions. By listing 63 quantitative trait loci (QTL) mapping and global gene expression studies based on plants directly challenged by other plants, we explored whether the genetic architecture and the function of the candidate genes underlying natural plant-plant interactions depend on the type of interactions between two plants (competition versus commensalism versus reciprocal helping versus asymmetry). The 16 transcriptomic studies are unevenly distributed between competitive interactions (n = 12) and asymmetric interactions (n = 4, all focusing on response to parasitic plants). By contrast, 17 and 30 QTL studies were identified for competitive interactions and asymmetric interactions (either weed suppressive ability or response to parasitic plants), respectively. Surprisingly, no studies have been carried out on the identification of genetic and molecular bases underlying natural variation in positive interactions. The candidate genes underlying natural plant-plant interactions can be classified into seven categories of plant function that have been identified in artificial environments simulating plant-plant interactions either frequently (photosynthesis, hormones), only recently (cell wall modification and degradation, defense pathways against pathogens) or rarely (ABC transporters, histone modification and meristem identity/life history traits). Finally, we introduce several avenues that need to be explored in the future to obtain a thorough understanding of the genetic and molecular bases underlying plant-plant interactions within the context of realistic community complexity.

A genome-wide association study of 23 agronomic traits in Chinese wheat landraces

Uncovering the genetic basis of agronomic traits in wheat landraces is important for ensuring global food security via the development of improved varieties. Here, 723 wheat landraces from 10 Chinese agro-ecological zones were evaluated for 23 agronomic traits in six environments. All accessions could be clustered into five subgroups based on phenotypic data via discriminant function analysis, which was highly consistent with genotypic classification. A genome-wide association study was conducted for these traits using 52 303 DArT-seq markers to identify marker-trait associations and candidate genes. Using both the general linear model and the mixed linear model, 149 significant markers were identified for 21 agronomic traits based on best linear unbiased prediction values. Considering the linkage disequilibrium decay distance in this study, significant markers within 10 cM were combined as a quantitative trait locus (QTL), with a total of 29 QTL identified for 15 traits. Of these, five QTL for heading date, flag leaf width, peduncle length, and thousand kernel weight had been reported previously. Twenty-five candidate genes associated with significant markers were identified. These included the known vernalization genes VRN-B1 and vrn-B3 and the photoperiod response genes Ppd and PRR. Overall, this study should be helpful in elucidating the underlying genetic mechanisms of complex agronomic traits and performing marker-assisted selection in wheat.

HiSeeker: Detecting High-Order SNP Interactions Based on Pairwise SNP Combinations

Detecting single nucleotide polymorphisms’ (SNPs) interaction is one of the most popular approaches for explaining the missing heritability of common complex diseases in genome-wide association studies. Many methods have been proposed for SNP interaction detection, but most of them only focus on pairwise interactions and ignore high-order ones, which may also contribute to complex traits. Existing methods for high-order interaction detection can hardly handle genome-wide data and suffer from low detection power, due to the exponential growth of search space. In this paper, we proposed a flexible two-stage approach (called HiSeeker) to detect high-order interactions. In the screening stage, HiSeeker employs the chi-squared test and logistic regression model to efficiently obtain candidate pairwise combinations, which have intermediate or significant associations with the phenotype for interaction detection. In the search stage, two different strategies (exhaustive search and ant colony optimization-based search) are utilized to detect high-order interactions from candidate combinations. The experimental results on simulated datasets demonstrate that HiSeeker can more efficiently and effectively detect high-order interactions than related representative algorithms. On two real case-control datasets, HiSeeker also detects several significant high-order interactions, whose individual SNPs and pairwise interactions have no strong main effects or pairwise interaction effects, and these high-order interactions can hardly be identified by related algorithms.

Manipulation of flowering time and branching by overexpression of the tomato transcription factor SlZFP2

Flowering of higher plants is orchestrated by complex regulatory networks through integration of various environmental signals such as photoperiod, temperature, light quality and developmental cues. In Arabidopsis, transcription of the flowering integrator gene FLOWERING LOCUS T (FT) that several flowering pathways converge to is directly regulated by more than ten transcription factors. However, very little is known about the transcriptional regulation of the FT homolog SINGLE FLOWER TRUESS (SFT) in the day-neutral plant tomato (Solanum lycopersicum). Previously, we showed that the zinc finger transcription factor SlZFP2 plays important roles in regulation of seed germination and fruit ripening in tomato and also found that overexpression of SlZFP2 impacted flowering and branching. Here, we characterized in detail the early flowering and high branching phenotypes by overexpression of this transcription factor. Our data showed that overexpression of SlZFP2 accelerated flowering in an SFT-dependent manner as demonstrated by elevated SFT expression in the leaves and the transcription factor's binding ability to SFT promoter in vitro and in vivo. Furthermore, overexpression of the SlZFP2 gene in the sft plants failed to rescue the mutant's late flowering. Through analysis of grafting phenotype, growth response of branches to auxin application and transcriptome profiling by RNA sequencing, we also showed that overexpression of SlZFP2 affected shoot apical dominance through multiple regulatory pathways. Our results suggest that the transcription factor SlZFP2 has potential applications in genetic modification of plant architecture and flowering time for tomato production and other crops as well.

1,135 Genomes Reveal the Global Pattern of Polymorphism in Arabidopsis thaliana

Arabidopsis thaliana serves as a model organism for the study of fundamental physiological, cellular, and molecular processes. It has also greatly advanced our understanding of intraspecific genome variation. We present a detailed map of variation in 1,135 high-quality re-sequenced natural inbred lines representing the native Eurasian and North African range and recently colonized North America. We identify relict populations that continue to inhabit ancestral habitats, primarily in the Iberian Peninsula. They have mixed with a lineage that has spread to northern latitudes from an unknown glacial refugium and is now found in a much broader spectrum of habitats. Insights into the history of the species and the fine-scale distribution of genetic diversity provide the basis for full exploitation of A. thaliana natural variation through integration of genomes and epigenomes with molecular and non-molecular phenotypes.

Environmental and genetic interactions reveal FLOWERING LOCUS C as a modulator of the natural variation for the plasticity of flowering in Arabidopsis

The timing of flowering initiation depends strongly on the environment, a property termed as the plasticity of flowering. Such plasticity determines the adaptive potential of plants because it provides phenotypic buffer against environmental changes, and its natural variation contributes to evolutionary adaptation. We addressed the genetic mechanisms of the natural variation for this plasticity in Arabidopsis thaliana by analysing a population of recombinant inbred lines derived from Don-0 and Ler accessions collected from distinct climates. Quantitative trait locus (QTL) mapping in four environmental conditions differing in photoperiod, vernalization treatment and ambient temperature detected the folllowing: (i) FLOWERING LOCUS C (FLC) as a large effect QTL affecting flowering time differentially in all environments; (ii) numerous QTL displaying smaller effects specifically in some conditions; and (iii) significant genetic interactions between FLC and other loci. Hence, the variation for the plasticity of flowering is determined by a combination of environmentally sensitive and specific QTL, and epistasis. Analysis of FLC from Don identified a new and more active allele likely caused by a cis-regulatory deletion covering the non-coding RNA COLDAIR. Further characterization of four FLC natural alleles showed different environmental and genetic interactions. Thus, FLC appears as a major modulator of the natural variation for the plasticity of flowering to multiple environmental factors.

Prevalent Exon-Intron Structural Changes in the APETALA1/FRUITFULL, SEPALLATA, AGAMOUS-LIKE6, and FLOWERING LOCUS C MADS-Box Gene Subfamilies Provide New Insights into Their Evolution

AP1/FUL, SEP, AGL6, and FLC subfamily genes play important roles in flower development. The phylogenetic relationships among them, however, have been controversial, which impedes our understanding of the origin and functional divergence of these genes. One possible reason for the controversy may be the problems caused by changes in the exon-intron structure of genes, which, according to recent studies, may generate non-homologous sites and hamper the homology-based sequence alignment. In this study, we first performed exon-by-exon alignments of these and three outgroup subfamilies (SOC1, AG, and STK). Phylogenetic trees reconstructed based on these matrices show improved resolution and better congruence with species phylogeny. In the context of these phylogenies, we traced evolutionary changes of exon-intron structures in each subfamily. We found that structural changes have occurred frequently following gene duplication and speciation events. Notably, exons 7 and 8 (if present) suffered more structural changes than others. With the knowledge of exon-intron structural changes, we generated more reasonable alignments containing all the focal subfamilies. The resulting trees showed that the SEP subfamily is sister to the monophyletic group formed by AP1/FUL and FLC subfamily genes and that the AGL6 subfamily forms a sister group to the three abovementioned subfamilies. Based on this topology, we inferred the evolutionary history of exon-intron structural changes among different subfamilies. Particularly, we found that the eighth exon originated before the divergence of AP1/FUL, FLC, SEP, and AGL6 subfamilies and degenerated in the ancestral FLC-like gene. These results provide new insights into the origin and evolution of the AP1/FUL, FLC, SEP, and AGL6 subfamilies.

A genome-wide association study of plant height and primary branch number in rapeseed (Brassica napus)

Crop plant architecture plays a highly important role in its agronomic performance. Plant height (PH) and primary branch number (PB) are two major factors that affect the plant architecture of rapeseed (Brassica napus). Previous studies have shown that these two traits are controlled by multiple quantitative trait loci (QTL); however, QTLs have not been delimited to regions less than 10cM. Genome-wide association study (GWAS) is a highly efficient approach for identifying genetic loci controlling traits at relatively high resolution. In this study, variations in PH and PB of a panel of 472 rapeseed accessions that had previously been analyzed by a 60k SNP array were investigated for three consecutive years and studied by GWAS. Eight QTLs on chromosome A03, A05, A07 and C07 were identified for PH, and five QTLs on A01, A03, A07 and C07 were identified for PB. Although most QTLs have been detected in previous studies based on linkage analyses, the two QTLs of PH on A05 and the QTL of PB on C07 were novel. In the genomic regions close to the GWAS peaks, orthologs of the genes involved in flower development, phytohormone biosynthesis, metabolism and signaling in Arabidopsis were identified.

A Bayesian model for detection of high-order interactions among genetic variants in genome-wide association studies

BACKGROUND

A central question for disease studies and crop improvements is how genetics variants drive phenotypes. Genome Wide Association Study (GWAS) provides a powerful tool for characterizing the genotype-phenotype relationships in complex traits and diseases. Epistasis (gene-gene interaction), including high-order interaction among more than two genes, often plays important roles in complex traits and diseases, but current GWAS analysis usually just focuses on additive effects of single nucleotide polymorphisms (SNPs). The lack of effective computational modelling of high-order functional interactions often leads to significant under-utilization of GWAS data.

RESULTS

We have developed a novel Bayesian computational method with a Markov Chain Monte Carlo (MCMC) search, and implemented the method as a Bayesian High-order Interaction Toolkit (BHIT) for detecting epistatic interactions among SNPs. BHIT first builds a Bayesian model on both continuous data and discrete data, which is capable of detecting high-order interactions in SNPs related to case--control or quantitative phenotypes. We also developed a pipeline that enables users to apply BHIT on different species in different use cases.

CONCLUSIONS

Using both simulation data and soybean nutritional seed composition studies on oil content and protein content, BHIT effectively detected some high-order interactions associated with phenotypes, and it outperformed a number of other available tools. BHIT is freely available for academic users at http://digbio.missouri.edu/BHIT/.

Multiple pathways regulate shoot branching

Shoot branching patterns result from the spatio-temporal regulation of axillary bud outgrowth. Numerous endogenous, developmental and environmental factors are integrated at the bud and plant levels to determine numbers of growing shoots. Multiple pathways that converge to common integrators are most probably involved. We propose several pathways involving not only the classical hormones auxin, cytokinins and strigolactones, but also other signals with a strong influence on shoot branching such as gibberellins, sugars or molecular actors of plant phase transition. We also deal with recent findings about the molecular mechanisms and the pathway involved in the response to shade as an example of an environmental signal controlling branching. We propose the TEOSINTE BRANCHED1, CYCLOIDEA, PCF transcription factor TB1/BRC1 and the polar auxin transport stream in the stem as possible integrators of these pathways. We finally discuss how modeling can help to represent this highly dynamic system by articulating knowledges and hypothesis and calculating the phenotype properties they imply.

QTL analysis of root morphology, flowering time, and yield reveals trade-offs in response to drought in Brassica napus

Drought escape and dehydration avoidance represent alternative strategies for drought adaptation in annual crops. The mechanisms underlying these two strategies are reported to have a negative correlation, suggesting a trade-off. We conducted a quantitative trait locus (QTL) analysis of flowering time and root mass, traits representing each strategy, in Brassica napus to understand if a trade-off exists and what the genetic basis might be. Our field experiment used a genotyped population of doubled haploid lines and included both irrigated and rainfed treatments, allowing analysis of plasticity in each trait. We found strong genetic correlations among all traits, suggesting a trade-off among traits may exist. Summing across traits and treatments we found 20 QTLs, but many of these co-localized to two major QTLs, providing evidence that the trade-off is genetically constrained. To understand the mechanistic relationship between root mass, flowering time, and QTLs, we analysed the data by conditioning upon correlated traits. Our results suggest a causal model where such QTLs affect root mass directly as well as through their impacts on flowering time. Additionally, we used draft Brassica genomes to identify orthologues of well characterized Arabidopsis thaliana flowering time genes as candidate genes. This research provides valuable clues to breeding for drought adaptation as it is the first to analyse the inheritance of the root system in B. napus in relation to drought.

Induced and natural variation of promoter length modulates the photoperiodic response of FLOWERING LOCUS T

FLOWERING LOCUS T (FT) regulates the floral transition in many plant species by integrating environmental seasonal signals and internal cues. Here we show that two interdependent regulatory regions are necessary and sufficient to convey photoperiod responsiveness to FT. While a minimal distance between the regulatory regions is required to fully suppress FT expression under short days, increased distance reduces promoter response to long days. Natural variation at FT creating promoter length differences is widespread, correlates with longitudinal and latitudinal clines and affects a promoter region physically interacting with both photoperiod control regions. Three major FT promoter variants correlate with differences in FT allele usage in F1 hybrids. We propose that FT variation in cis could be adaptive by conferring differences in FT transcriptional control ultimately translating to increased fitness.

Gene expression at FLOWERING LOCUS T controls floral transition in many plants and is regulated by both environmental signals and internal cues. Liu et al. show that the distance between two regulatory sequences in the FT promoter varies with geographical location and determines responsiveness to photoperiod.

QTL meta-analysis in Arabidopsis reveals an interaction between leaf senescence and resource allocation to seeds

Sequential and monocarpic senescence are observed at vegetative and reproductive stages, respectively. Both facilitate nitrogen (N) remobilization and control the duration of carbon (C) fixation. Genetic and environmental factors control N and C resource allocation to seeds. Studies of natural variation in Arabidopsis thaliana revealed differences between accessions for leaf senescence phenotypes, seed N and C contents, and N remobilization efficiency-related traits. Here, a quantitative genetics approach was used to gain a better understanding of seed filling regulation in relation to leaf senescence and resource allocation. For that purpose, three Arabidopsis recombinant inbred line populations (Ct-1×Col-0, Cvi-0×Col-0, Bur-0×Col-0) were used to map QTL (quantitative trait loci) for ten traits related to senescence, resource allocation, and seed filling. The use of common markers across the three different maps allowed direct comparisons of the positions of the detected QTL in a single consensus map. QTL meta-analysis was then used to identify interesting regions (metaQTL) where QTL for several traits co-localized. MetaQTL were compared with positions of candidate genes known to be involved in senescence processes and flowering time. Finally, investigation of the correlation between yield and seed N concentration in the three populations suggests that leaf senescence disrupts the negative correlation generally observed between these two traits.

Genetic architecture of naturally occurring quantitative traits in plants: an updated synthesis

Deciphering the genetic and molecular bases of quantitative variation is a long-standing challenge in plant biology because it is essential for understanding evolution and for accelerating plant breeding. Recent multi-trait analyses at different phenotypic levels are uncovering the pleiotropy and the genetic regulation underlying high-level complex traits. Thus, the number of known causal loci, genes and nucleotide polymorphisms is expanding. Current plant causal catalogs contain ∼400 genes and natural polymorphisms revealing several dysfunctional allelic series that involve multiple mutations. In addition, repeated evolution of quantitative traits mediated by large effect alleles is found across plant phylogeny. Finally, systematic analyses of genetic and environmental interactions are beginning to elucidate the molecular mechanisms of relevant interactions.

Post-fertilization expression of FLOWERING LOCUS T suppresses reproductive reversion

FLOWERING LOCUS T (FT) encodes a systemic signal communicating the perception of long day photoperiod from leaves to the shoot apex to induce the floral transition. Transient expression of FT in the phloem companion cells of rosette leaves for one to several days was previously shown to be sufficient to commit plants to flowering. Here we show that partial commitment results in pleiotropic inflorescence meristem reversion phenotypes. FT expression is much stronger in organs formed after the floral transition such as cauline leaves, sepals, and developing siliques. We show that expression of FT and its paralog TWIN SISTER OF FT (TSF) after the floral transition plays a role in inflorescence meristem stabilization even if plants flower very late in development. CONSTANS (CO), the major activator of FT, is not required to prevent late reproductive reversion. The requirement for FT is temporal since reproductive reversion to a vegetative state occurs only in recently formed inflorescence meristems. Unlike for the expression of FT in leaves, neither the distal putative FT enhancer nor long-day photoperiod is required for FT expression in developing siliques. Expression of FT in developing siliques and their supporting stems is sufficient to stabilize flowering during the sensitive developmental window indicating that fruit generated FT participates in inflorescence stabilization.