Quantitative trait loci mapping of floral and leaf morphology traits in Arabidopsis thaliana: evidence for modular genetic architecture

Summary Morphological variation within organisms is integrated and often modular in nature. That is to say, the size and shape of traits tend to vary in a coordinated and structured manner across sets of organs or parts of an organism. The genetic basis of this morphological integration is largely unknown. Here, we report on quantitative trait loci (QTL) analysis of leaf and floral organ size in Arabidopsis thaliana. We evaluate patterns of genetic correlations among traits and perform whole-genome scans using QTL mapping methods. We detected significant genetic variation for the size and shape of each floral and leaf trait in our study. Moreover, we found large positive genetic correlations among sets of either flower or leaf traits, but low and generally nonsignificant genetic correlations between flower and leaf traits. These results support the hypothesis of independent floral and vegetative modules. We consider co-localization of QTL for different traits as support for a pleiotropic basis of morphological integration and modularity. A total of eight QTL affecting flower and three QTL affecting leaf traits were identified. Most QTL affected either floral or leaf traits, providing a general explanation for high correlations within and low correlations between modules. Only two genomic locations affected both flower and leaf growth. These results are discussed in the context of the evolution of modules, pleiotropy, and the putative homologous relationship between leaves and flowers.     

Quantitative trait loci for inflorescence development in Arabidopsis thaliana

Variation in inflorescence development patterns is a central factor in the evolutionary ecology of plants. The genetic architectures of 13 traits associated with inflorescence developmental timing, architecture, rosette morphology, and fitness were investigated in Arabidopsis thaliana, a model plant system. There is substantial naturally occurring genetic variation for inflorescence development traits, with broad sense heritabilities computed from 21 Arabidopsis ecotypes ranging from 0.134 to 0.772. Genetic correlations are significant for most (64/78) pairs of traits, suggesting either pleiotropy or tight linkage among loci. Quantitative trait locus (QTL) mapping indicates 47 and 63 QTL for inflorescence developmental traits in Ler x Col and Cvi x Ler recombinant inbred mapping populations, respectively. Several QTL associated with different developmental traits map to the same Arabidopsis chromosomal regions, in agreement with the strong genetic correlations observed. Epistasis among QTL was observed only in the Cvi x Ler population, and only between regions on chromosomes 1 and 5. Examination of the completed Arabidopsis genome sequence in three QTL regions revealed between 375 and 783 genes per region. Previously identified flowering time, inflorescence architecture, floral meristem identity, and hormone signaling genes represent some of the many candidate genes in these regions.     

Quantitative trait loci analysis of leaf and plant longevity in Arabidopsis thaliana

The natural variation in leaf and plant longevity in Arabidopsis thaliana was analysed in a set of 45 ecotypes and 155 recombinant inbred lines derived from a Cape Verde Islands (Cvi) x Landsberg erecta (Ler) cross. Post-bolting longevity was inversely related to time to flowering and rosette leaf number in the set of 45 ecotypes, with Cvi having the longest and Ler the shortest post-bolting longevity. The recombinant inbred line population was tested under low or high soil nutrient levels (LN or HN, respectively). Three quantitative trait loci (QTL), one in chromosome 3 and two in chromosomes 1 and 5, were associated with longevity of the 6th rosette leaf under LN and HN, respectively. Four QTL for post-bolting longevity were found in chromosomes 1, 3, 4, and 5, and two in chromosomes 1 and 5 under LN and HN, respectively. An epistatic interaction affecting post-bolting longevity under LN, but not HN, was detected. Ler and Cvi carry a mix of increasing and decreasing alleles for the QTL affecting longevity of the 6th leaf and post-bolting longevity. Longevity of the 6th rosette leaf was associated with different QTL than post-bolting longevity, and it was affected by different QTL depending on nutrient availability. By contrast, the major QTL affecting post-bolting longevity exerted significant effects irrespective of soil nutrient availability.     

Quantitative trait loci for partial resistance to Pseudomonas syringae pv. maculicola in Arabidopsis thaliana

Segregation of partial resistance to Pseudomonas syringae pv. maculicola (Psm) ES4326 was studied in the recombinant inbred population created from accessions (ecotypes) Columbia (Col‐4), the more susceptible parent, and Landsberg (Ler‐0). Plants were spray inoculated with lux‐transformed bacteria in experiments to measure susceptibility. The amount of disease produced on a range of Col × Ler lines by spray inoculation was highly correlated with that produced by pressure infiltration of bacteria into the apoplast. Quantitative trait locus (QTL) analysis identified four loci that contributed to partial resistance: QRps.JIC‐1.1, QRps.JIC‐2.1, QRps.JIC‐3.1 and QRps.JIC‐5.1 on chromosomes 1, 2, 3 and 5, respectively. QRps.JIC‐3.1, located 8.45 cM from the top of the consensus genetic map of chromosome 3, had a large, approximately additive effect on partial resistance, explaining 50% of the genetic variation in this population. Fine mapping narrowed the region within which this QTL was located to 62 genes. A list of candidate genes included several major classes of resistance gene.     

Quantitative trait loci controlling light and hormone response in two accessions of Arabidopsis thaliana

We have mapped quantitative trait loci (QTL) responsible for natural variation in light and hormone response between the Cape Verde Islands (Cvi) and Landsberg erecta (Ler) accessions of Arabidopsis thaliana using recombinant inbred lines (RILs). Hypocotyl length was measured in four light environments: white, blue, red, and far-red light and in the dark. In addition, white light plus gibberellin (GA) and dark plus the brassinosteroid biosynthesis inhibitor brassinazole (BRZ) were used to detect hormone effects. Twelve QTL were identified that map to loci not previously known to affect light response, as well as loci where candidate genes have been identified from known mutations. Some QTL act in all environments while others show genotype-by-environment interaction. A global threshold was established to identify a significant epistatic interaction between two loci that have few main effects of their own. LIGHT1, a major QTL, has been confirmed in a near isogenic line (NIL) and maps to a new locus with effects in all light environments. The erecta mutation can explain the effect of the HYP2 QTL in the blue, BRZ, and dark environments, but not in far-red. LIGHT2, also confirmed in an NIL, has effects in white and red light and shows interaction with GA. The phenotype and map position of LIGHT2 suggest the photoreceptor PHYB as a candidate gene. Natural variation in light and hormone response thus defines both new genes and known genes that control light response in wild accessions.     

Quantitative trait loci controlling leaf venation in Arabidopsis

Leaf veins provide the mechanical support and are responsible for the transport of nutrients and water to the plant. High vein density is a prerequisite for plants to have C4 photosynthesis. We investigated the genetic variation and genetic architecture of leaf venation traits within the species Arabidopsis thaliana using natural variation. Leaf venation traits, including leaf vein density (LVD) were analysed in 66 worldwide accessions and 399 lines of the multi‐parent advanced generation intercross population. It was shown that there is no correlation between LVD and photosynthesis parameters within A. thaliana. Association mapping was performed for LVD and identified 16 and 17 putative quantitative trait loci (QTLs) in the multi‐parent advanced generation intercross and worldwide sets, respectively. There was no overlap between the identified QTLs suggesting that many genes can affect the traits. In addition, linkage mapping was performed using two biparental recombinant inbred line populations. Combining linkage and association mapping revealed seven candidate genes. For one of the candidate genes, RCI2c, we demonstrated its function in leaf venation patterning.     

Quantitative trait loci affecting differences in floral morphology between two species of monkeyflower (Mimulus).

Conspicuous differences in floral morphology are partly responsible for reproductive isolation between two sympatric species of monkeyflower because of their effect on visitation of the flowers by different pollinators. Mimulus lewisii flowers are visited primarily by bumblebees, whereas M. cardinalis flowers are visited mostly by hummingbirds. The genetic control of 12 morphological differences between the flowers of M. lewisii and M. cardinalis was explored in a large linkage mapping population of F2 plants n = 465 to provide an accurate estimate of the number and magnitude of effect of quantitative trait loci (QTLs) governing each character. Between one and six QTLs were identified for each trait. Most (9/12) traits appear to be controlled in part by at least one major QTL explaining >/=25% of the total phenotypic variance. This implies that either single genes of individually large effect or linked clusters of genes with a large cumulative effect can play a role in the evolution of reproductive isolation and speciation.     

Quantitative trait loci controlling root growth and architecture in Arabidopsis thaliana confirmed by heterogeneous inbred family

Arabidopsis thaliana provides a scientifically attractive and simple model for studying root growth and architecture and, subsequently, for discovering new genes involved in the control of these characters in plants. We have used the natural variation available in Arabidopsis accessions and mapped quantitative trait loci (QTLs) for primary root length (PRL), lateral root number (LRN) and density (LRD) and for total length of the lateral root system (LRL) in the Bay-0 × Shahdara population. Total phenotypic variation was very large, and despite the importance of the environmental component we were able to map 13 QTLs and one epistatic interaction between QTLs. Our results highlight the biological relevance and genetic control of lateral root density in this material. We were also able to show that variation in the extent of the lateral root system depends mainly on the growth of the existing lateral roots rather than in a change in LRN. Factors controlling lateral root growth seemed to have no major effect on primary root growth. Moreover, Shahdara QTL alleles always increased the length of the lateral roots, which may be taken as an adaptation to its very dry natural environment in Tadjikistan. A QTL for PRL was confirmed using a type of near-isogenic line called a heterogeneous inbred family (HIF), and this QTL is a candidate for further fine-mapping and cloning.     

Quantitative trait loci affecting delta13C and response to differential water availibility in Arabidopsis thaliana

Phenotypic plasticity is an important response mechanism of plants to environmental heterogeneity. Here, we explored the genetic basis of plastic responses of Arabidopsis thaliana to water deficit by experimentally mapping quantitative trait loci (QTL) in two recombinant inbred populations (Cvi x Ler and Ler x Col). We detected genetic variation and significant genotype-by-environment interactions for many traits related to water use. We also mapped 26 QTL, including six for carbon isotope composition (delta13C). Negative genetic correlations between fruit length and fruit production as well as between flowering time and branch production were corroborated by QTL colocalization, suggesting these correlations are due to pleiotropy or physical linkage. Water-limited plants were more apically dominant with greater root:shoot ratios and higher delta13C (higher water-use efficiency) when compared to well-watered plants. Many of the QTL effects for these traits interacted significantly with the irrigation treatment, suggesting that the observed phenotypic plasticity is genetically based. We specifically searched for epistatic (QTL-QTL) interactions using a two-dimensional genome scan, which allowed us to detect epistasis regardless of additive genetic effects. We found several significant QTL-QTL interactions including three that exhibited environmental dependence. These results provide preliminary evidence for proposed genetic mechanisms underlying phenotypic plasticity.     

Quantitative trait loci analysis of water and anion contents in interaction with nitrogen availability in Arabidopsis thaliana

In plants, water and anion parameters are linked, for example through the integration of nutritional signaling and the response to diverse stress. In this work, Arabidopsis thaliana is used as a model system to dissect the genetic variation of these parameters by quantitative trait loci (QTL) mapping in the 415 recombinant inbred lines of the Bay-0 x Shahdara population. Water, nitrate, chloride, and phosphate contents were measured at the vegetative stage in the shoots of plants grown in controlled conditions. Two contrasting nitrogen (N) conditions were studied, one leading to the complete depletion of the nitrate pool in the plants. Most of the observed genetic variation was identified as QTL, with medium but also large phenotypic contributions. QTL colocalization provides a genetic basis for the correlation between water and nitrate contents in nonlimiting N conditions and water and chloride contents in limiting N conditions. The 34 new QTL described here represent at least 19 loci polymorphic between Bay-0 and Shahdara; some may correspond to known genes from water/anion transport systems, while others clearly identify new genes controlling or interacting with water/anion absorption and accumulation. Interestingly, flowering-time genes probably play a role in the regulation of water content in our conditions.     

Allele-specific effects of PDF2 on floral morphology in Arabidopsis thaliana

The class IV Homeodomain-leucine zipper (HD-ZIP IV) gene family includes several genes that are functionally significant in epidermal development. Our recent study revealed that double mutants of the epidermis-expressed HD-ZIP IV members, PROTODERMAL FACTOR2 (PDF2) in combination with some HOMEODOMAIN GLABROUS (HDG, pronounced “hedge”) genes, affect stamen development and specification of petal and stamen identity, possibly in a non cell-autonomous manner. However, the effect of the pdf2 mutations on the floral development was largely different depending on T-DNA insertion locations: pdf2–1 hdg flowers exhibited homeotic conversion of petals and stamens, while pdf2–2 hdg flowers had only a reduced number of stamens. Here, we used 2 additional pdf2 alleles to make double mutants and found that their floral phenotypes were rather similar to those of pdf2–2 hdg. The allele-specific effect caused by pdf2–1, which carries a T-DNA in a steroidogenic acute regulatory protein-related lipid transfer (START) domain-encoding region, suggests the importance of the START domain in proper function of HD-ZIP IV proteins.     

Mapping quantitative trait loci associated with selenate tolerance in Arabidopsis thaliana

Selenium is essential for many organisms, but is toxic at higher levels. To investigate the genetic basis of selenate tolerance in Arabidopsis thaliana, quantitative trait loci (QTL) associated with selenate tolerance in accessions Landsberg erecta and Columbia were mapped using recombinant inbred lines (RILs). The selenate tolerance index (TI(D10) = root growth + 30 microm selenate/root growth control x 100%) was fourfold higher for parental line Col-4 (59%) than for parent Ler-0 (15%). Among the 96 F8 RILs, TI(D10) ranged from 11 to 75% (mean 37%). Using composite interval mapping, three QTL were found on chromosomes 1, 3 and 5, which together explained 24% of variation in TI(D10) and 32% of the phenotypic variation for the difference in root length +/- Se (RL(D10)). Highly significant epistatic interactions between the QTL and markers on chromosome 2 explained additional variation for both traits. Potential candidate genes for Se tolerance in each of the QTL regions are discussed. These results offer insight into the genetic basis of selenate tolerance, and may be useful for identification of selenate-tolerance genes.     

Quantitative trait loci analysis of powdery mildew disease resistance in the Arabidopsis thaliana accession kashmir-1.

Powdery mildew diseases are economically important diseases, caused by obligate biotrophic fungi of the Erysiphales. To understand the complex inheritance of resistance to the powdery mildew disease in the model plant Arabidopsis thaliana, quantitative trait loci analysis was performed using a set of recombinant inbred lines derived from a cross between the resistant accession Kashmir-1 and the susceptible accession Columbia glabrous1. We identified and mapped three independent powdery mildew quantitative disease resistance loci, which act additively to confer disease resistance. The locus with the strongest effect on resistance was mapped to a 500-kbp interval on chromosome III.     

Natural variation in GA1 associates with floral morphology in Arabidopsis thaliana

? The genetic architecture of floral traits is evolutionarily important due to the fitness consequences of quantitative variation in floral morphology. Yet, little is known about the genes underlying these traits in natural populations. Using Arabidopsis thaliana, we examine molecular variation at GIBBERELLIC ACID REQUIRING 1 (GA1) and test for associations with floral morphology. ? We examined full-length sequence in 32 accessions and describe two haplotypes (comprising four nonsynonymous polymorphisms) in GA1 that segregate at intermediate frequencies. In 133 A. thaliana accessions, we test for genotype-phenotype associations and corroborate these findings in segregating progenies. ? The two common GA1 haplotypes were associated with the length of petals, stamens, and to a lesser extent style-stigma length. Associations were confirmed in a segregating progeny developed from 19 accessions. We find analogous results in recombinant inbred lines of the Bayreuth × Shahdara cross, which differ only at one of 4 SNPs, suggesting that this SNP may contribute to the observed association. ? Assuming GA1 causally affects floral organ size, it is interesting that adjacent petal and stamen whorls are most strongly affected. This pattern suggests that GA1 could contribute to the greater strength of petal-stamen correlations relative to other floral-length correlations observed in some Brassicaceous species.     

Quantitative trait loci analysis of nitrogen use efficiency in Arabidopsis

Improving plant nitrogen (N) use efficiency or controlling soil N requires a better knowledge of the regulation of plant N metabolism. This could be achieved using Arabidopsis as a model genetic system, taking advantage of the natural variation available among ecotypes. Here, we describe an extensive study of N metabolism variation in the Bay-0 x Shahdara recombinant inbred line population, using quantitative trait locus (QTL) mapping. We mapped QTL for traits such as shoot growth, total N, nitrate, and free-amino acid contents, measured in two contrasting N environments (contrasting nitrate availability in the soil), in controlled conditions. Genetic variation and transgression were observed for all traits, and most of the genetic variation was identified through QTL and QTL x QTL epistatic interactions. The 48 significant QTL represent at least 18 loci that are polymorphic between parents; some may correspond to known genes from the N metabolic pathway, but others represent new genes controlling or interacting with N physiology. The correlations between traits are dissected through QTL colocalizations: The identification of the individual factors contributing to the regulation of different traits sheds new light on the relations among these characters. We also point out that the regulation of our traits is mostly specific to the N environment (N availability). Finally, we describe four interesting loci at which positional cloning is feasible.     

Genotype-environment interactions at quantitative trait loci affecting inflorescence development in Arabidopsis thaliana

Phenotypic plasticity and genotype-environment interactions (GEI) play a prominent role in plant morphological diversity and in the potential functional capacities of plant life-history traits. The genetic basis of plasticity and GEI, however, is poorly understood in most organisms. In this report, inflorescence development patterns in Arabidopsis thaliana were examined under different, ecologically relevant photoperiod environments for two recombinant inbred mapping populations (Ler x Col and Cvi x Ler) using a combination of quantitative genetics and quantitative trait locus (QTL) mapping. Plasticity and GEI were regularly observed for the majority of 13 inflorescence traits. These observations can be attributable (at least partly) to variable effects of specific QTL. Pooled across traits, 12/44 (27.3%) and 32/62 (51.6%) of QTL exhibited significant QTL x environment interactions in the Ler x Col and Cvi x Ler lines, respectively. These interactions were attributable to changes in magnitude of effect of QTL more often than to changes in rank order (sign) of effect. Multiple QTL x environment interactions (in Cvi x Ler) clustered in two genomic regions on chromosomes 1 and 5, indicating a disproportionate contribution of these regions to the phenotypic patterns observed. High-resolution mapping will be necessary to distinguish between the alternative explanations of pleiotropy and tight linkage among multiple genes.     

Epistasis and genotype-environment interaction for quantitative trait loci affecting flowering time in Arabidopsis thaliana

A major goal of evolutionary biology is to understand the genetic architecture of the complex quantitative traits that may lead to adaptations in natural populations. Of particular relevance is the evaluation of the frequency and magnitude of epistasis (gene–gene and gene–environment interaction) as it plays a controversial role in models of adaptation within and among populations. Here, we explore the genetic basis of flowering time in Arabidopsis thaliana using a series of quantitative trait loci (QTL) mapping experiments with two recombinant inbred line (RIL) mapping populations [Columbia (Col) x Landsberg erecta (Ler), Ler x Cape Verde Islands (Cvi)]. We focus on the response of RILs to a series of environmental conditions including drought stress, leaf damage, and apical damage. These data were explicitly evaluated for the presence of epistasis using Bayesian based multiple-QTL genome scans. Overall, we mapped fourteen QTL affecting flowering time. We detected two significant QTL–QTL interactions and several QTL–environment interactions for flowering time in the Ler x Cvi population. QTL–environment interactions were due to environmentally induced changes in the magnitude of QTL effects and their interactions across environments – we did not detect antagonistic pleiotropy. We found no evidence for QTL interactions in the Ler x Col population. We evaluate these results in the context of several other studies of flowering time in Arabidopsis thaliana and adaptive evolution in natural populations.     

Quantitative trait loci analysis of primary cell wall composition in Arabidopsis

Quantitative trait loci (QTL) analysis was used to identify genes underlying natural variation in primary cell wall composition in Arabidopsis (Arabidopsis thaliana). The cell walls of dark-grown seedlings of a Bay-0 x Shahdara recombinant inbred line population were analyzed using three miniaturized global cell wall fingerprinting techniques: monosaccharide composition analysis by gas chromatography, xyloglucan oligosaccharide mass profiling, and whole-wall Fourier-transform infrared microspectroscopy. Heritable variation and transgression were observed for the arabinose-rhamnose ratio, xyloglucan side-chain composition (including O-acetylation levels), and absorbance for a subset of Fourier-transform infrared wavenumbers. In total, 33 QTL, corresponding to at least 11 different loci controlling dark-grown hypocotyl length, pectin composition, and levels of xyloglucan fucosylation and O-acetylation, were identified. One major QTL, accounting for 51% of the variation in the arabinose-rhamnose ratio, affected the number of arabinan side chains presumably attached to the pectic polysaccharide rhamnogalacturonan I, paving the way to positional cloning of the first gene underlying natural variation in pectin structure. Several QTL were found to be colocalized, which may have implications for the regulation of xyloglucan metabolism. These results demonstrate the feasibility of combining fingerprinting techniques, natural variation, and quantitative genetics to gain original insight into the molecular mechanisms underlying the structure and metabolism of cell wall polysaccharides.     

Quantitative trait loci and candidate genes underlying genotype by environment interaction in the response of Arabidopsis thaliana to drought

Drought stress was imposed on two sets of Arabidopsis thaliana genotypes grown in sand under short‐day conditions and analysed for several shoot and root growth traits. The response to drought was assessed for quantitative trait locus (QTL) mapping in a genetically diverse set of Arabidopsis accessions using genome‐wide association (GWA) mapping, and conventional linkage analysis of a recombinant inbred line (RIL) population. Results showed significant genotype by environment interaction (G×E) for all traits in response to different watering regimes. For the RIL population, the observed G×E was reflected in 17 QTL by environment interactions (Q×E), while 17 additional QTLs were mapped not showing Q×E. GWA mapping identified 58 single nucleotide polymorphism (SNPs) associated with loci displaying Q×E and an additional 16 SNPs associated with loci not showing Q×E. Many candidate genes potentially underlying these loci were suggested. The genes for RPS3C and YLS7 were found to contain conserved amino acid differences when comparing Arabidopsis accessions with strongly contrasting drought response phenotypes, further supporting their candidacy. One of these candidate genes co‐located with a QTL mapped in the RIL population.