Polymorphic Genes of Major Effect: Consequences for Variation,Selection and Evolution in Arabidopsis thaliana

The importance of genes of major effect for evolutionary trajectories within and among natural populations has long been the subject of intense debate. For example, if allelic variation at a major-effect locus fundamentally alters the structure of quantitative trait variation, then fixation of a single locus can have rapid and profound effects on the rate or direction of subsequent evolutionary change. Using an Arabidopsis thaliana RIL mapping population, we compare G-matrix structure between lines possessing different alleles at ERECTA, a locus known to affect ecologically relevant variation in plant architecture. We find that the allele present at ERECTA significantly alters G-matrix structure—in particular the genetic correlations between branch number and flowering time traits—and may also modulate the strength of natural selection on these traits. Despite these differences, however, when we extend our analysis to determine how evolution might differ depending on the ERECTA allele, we find that predicted responses to selection are similar. To compare responses to selection between allele classes, we developed a resampling strategy that incorporates uncertainty in estimates of selection that can also be used for statistical comparisons of G matrices.THE structure of the genetic variation that underlies phenotypic traits has important consequences for understanding the evolution of quantitative traits (Fisher 1930; Lande 1979; Bulmer 1980; Kimura 1983; Orr 1998; Agrawal et al. 2001). Despite the infinitesimal model''s allure and theoretical tractability (see Orr and Coyne 1992; Orr 1998, 2005a,b for reviews of its influence), evidence has accumulated from several sources (artificial selection experiments, experimental evolution, and QTL mapping) to suggest that genes of major effect often contribute to quantitative traits. Thus, the frequency and role of genes of major effect in evolutionary quantitative genetics have been a subject of intense debate and investigation for close to 80 years (Fisher 1930; Kimura 1983; Orr 1998, 2005a,b). Beyond the conceptual implications, the prevalence of major-effect loci also affects our ability to determine the genetic basis of adaptations and species differences (e.g., Bradshaw et al. 1995, 1998).Although the existence of genes of major effect is no longer in doubt, we still lack basic empirical data on how segregating variation at such genes affects key components of evolutionary process (but see Carrière and Roff 1995). In other words, How does polymorphism at genes of major effect alter patterns of genetic variation and covariation, natural selection, and the likely response to selection? The lack of data stems, in part, from the methods used to detect genes of major effect: experimental evolution (e.g., Bull et al. 1997; Zeyl 2005) and QTL analysis (see Erickson et al. 2004 for a review) often detect such genes retrospectively after they have become fixed in experimental populations or the species pairs used to generate the mapping population. The consequences of polymorphism at these genes on patterns of variation, covariation, selection, and the response to selection—which can be transient (Agrawal et al. 2001)—are thus often unobserved.A partial exception to the absence of data on the effects of major genes comes from artificial selection experiments, in which a substantial evolutionary response to selection in the phenotype after a plateau is often interpreted as evidence for the fixation of a major-effect locus (Frankham et al. 1968; Yoo 1980a,b; Frankham 1980; Shrimpton and Robertson 1988a,b; Caballero et al. 1991; Keightley 1998; see Mackay 1990 and Hill and Caballero 1992 for reviews). However, many of these experiments report only data on the selected phenotype (e.g., bristle number) or, alternatively, the selected phenotype and some measure of fitness (e.g., Frankham et al. 1968, Yoo 1980b; Caballero et al. 1991; Mackay et al. 1994; Fry et al. 1995; Nuzhdin et al. 1995; Zur Lage et al. 1997), making it difficult to infer how a mutation will affect variation, covariation, selection, and evolutionary responses for a suite of traits that might affect fitness themselves. One approach is to document how variation at individual genes of major effect affects the genetic variance–covariance matrix (“G matrix”; Lande 1979), which represents the additive genetic variance and covariance between traits.Although direct evidence for variation at major-effect genes altering patterns of genetic variation, covariation, and selection is rare, there is abundant evidence for the genetic mechanisms that could produce these dynamics. A gene of major effect could have these consequences due to any of at least three genetic mechanisms: (1) pleiotropy, where a gene of major effect influences several traits, including potentially fitness, simultaneously, (2) physical linkage or linkage disequilibrium (LD), in which a gene of major effect is either physically linked or in LD with other genes that influence other traits under selection, and (3) epistasis, in which the allele present at a major-effect gene alters the phenotypic effect of other loci and potentially phenotypes under selection. Evidence for these three evolutionary genetic mechanisms leading to changes in suites of traits comes from a variety of sources, including mutation accumulation experiments (Clark et al. 1995; Fernandez and Lopez-Fanjul 1996), mutation induction experiments (Keightley and Ohnishi 1998), artificial selection experiments (Long et al. 1995), and transposable element insertions (Rollmann et al. 2006). For pleiotropy in particular, major-effect genes that have consequences on several phenotypic traits are well known from the domestication and livestock breeding literature [e.g., myostatin mutations in Belgian blue cattle and whippets (Arthur 1995; Grobet et al. 1997; Mosher et al. 2007), halothane genes in pigs (Christian and Rothschild 1991; Fujii et al. 1991), and Booroola and Inverdale genes in sheep (Amer et al. 1999; Visscher et al. 2000)]. While these data suggest that variation at major-effect genes could—and probably does—influence variation, covariation, and selection on quantitative traits, data on the magnitude of these consequences remain lacking.Recombinant inbred line (RIL) populations are a promising tool for investigating the influence of major-effect loci. During advancement of the lines from F₂''s to RILs, alternate alleles at major-effect genes (and most of the rest of the genome) will be made homozygous, simplifying comparisons among genotypic classes. Because of the high homozygosity, individuals within RILs are nearly genetically identical, facilitating phenotyping of many genotypes under a range of environments. In addition, because of recombination, alternative alleles are randomized across genetic backgrounds—facilitating robust comparisons between sets of lines differing at a major-effect locus.Here we investigate how polymorphism at an artificially induced mutation, the erecta locus in Arabidopsis thaliana, affects the magnitude of these important evolutionary genetic parameters under ecologically realistic field conditions. We use the Landsberg erecta (Ler) × Columbia (Col) RIL population of A. thaliana to examine how variation at a gene of major effect influences genetic variation, covariation, and selection on quantitative traits in a field setting. The Ler × Col RIL population is particularly suitable, because it segregates for an artificially induced mutation at the erecta locus, which has been shown to influence a wide variety of plant traits. The Ler × Col population thus allows a powerful test of the effects of segregating variation at a gene—chosen a priori—with numerous pleiotropic effects. The ERECTA gene is a leucine-rich receptor-like kinase (LRR-RLK) (Torii et al. 1996) and has been shown to affect plant growth rates (El-Lithy et al. 2004), stomatal patterning and transpiration efficiency (Masle et al. 2005; Shpak et al. 2005), bacterial pathogen resistance (Godiard et al. 2003), inflorescence and floral organ size and shape (Douglas et al. 2002; Shpak et al. 2003, 2004), and leaf polarity (Xu et al. 2003; Qi et al. 2004).Specifically, we sought to answer the following questions: (1) Is variation at erecta significantly associated with changes to the G matrix? (2) Is variation at erecta associated with changes in natural selection on genetically variable traits? And (3) is variation at erecta associated with significantly different projected evolutionary responses to selection?     

Identification of QTLs Conferring Resistance to Deltamethrin in Culex pipiens pallens

Culex pipiens pallens is the most abundant Culex mosquito species in northern China and is an important vector of bancroftian filariasis and West Nile virus. Deltamethrin is an insecticide that is widely used for mosquito control, however resistance to this and other insecticides has become a major challenge in the control of vector-borne diseases that appear to be inherited quantitatively. Furthermore, the genetic basis of insecticide resistance remains poorly understood. In this study, quantitative trait loci (QTL) mapping of resistance to deltamethrin was conducted in F₂ intercross segregation populations using bulked segregation analysis (BSA) and amplified fragment length polymorphism markers (AFLP) in Culex pipiens pallens. A genetic linkage map covering 381 cM was constructed and a total of seven QTL responsible for resistance to deltamethrin were detected by composite interval mapping (CIM), which explained 95% of the phenotypic variance. The major QTL in linkage group 2 accounted for 62% of the variance and is worthy of further study. 12 AFLP markers in the map were cloned and the genomic locations of these marker sequences were determined by applying the Basic Local Alignment Search Tool (BLAST) tool to the genome sequence of the closely related Culex quinquefasciatus. Our results suggest that resistance to deltamethrin is a quantitative trait under the control of a major QTL in Culex pipiens pallens. Cloning of related AFLP markers confirm the potential utility for anchoring the genetic map to the physical map. The results provide insight into the genetic architecture of the trait.     

Identification of Mutations in Three Genes That Interact with Zeste in the Control of White Gene Expression in Drosophila Melanogaster

Three previously described genes, enhancer of yellow, 1, 2 and 3, are shown to cooperate with the zeste gene in the control of white gene expression. The mutations e(y)1(u1), e(y)3(u1), and to a lesser extent e(y)2(u1), enhance the effect of the zeste null allele z(v77h). Different combinations of e(y)1(u1), e(y)2(u1) and e(y)3(u1) mutations with several other z alleles also enhance the white mutant phenotype, but only to levels characteristic of white alleles containing a deletion of the upstream eye enhancer. Loss of zeste protein binding sites from the white locus does not eliminate the effect of e(y)1(u1) and e(y)3(u1) mutations, suggesting that the products of these genes interact with some other nucleotide sequences. Combinations of either e(y)1(u1) or e(y)2(u1) mutations with e(y)3(u1) are lethal. The products of these three genes may represent, together with zeste, a group of proteins involved in the organization of long-distance interactions between DNA sequences.     

The Effect of the Knockout of Major Transsulfuration Genes on the Pattern of Protein Synthesis in D. melanogaster

Molecular Biology - The enzymes involved in the transsulfuration pathway and hydrogen sulfide production—cystathionine-β-synthase (CBS), cystathionine-γ-lyase (CSE) and...     

Chromosome Copy Number Variation and Control in the Ciliate Chilodonella uncinata

Copy number variations are widespread in eukaryotes. The unusual genome architecture of ciliates, in particular, with its process of amitosis in macronuclear division, provides a valuable model in which to study copy number variation. The current model of amitosis envisions stochastic distribution of macronuclear chromosomes during asexual reproduction. This suggests that amitosis is likely to result in high levels of copy number variation in ciliates, as dividing daughter cells can have variable copy numbers of chromosomes if chromosomal distribution during amitosis is a stochastic process. We examined chromosomal distribution during amitosis in Chilodonella uncinata, a ciliate with gene-size macronuclear chromosomes. We quantified 4 chromosomes in evolving populations of C. uncinata and modeled the amitotic distribution process. We found that macronuclear chromosomes differ in copy number from one another but that copy number does not change as expected under a stochastic process. The chromosome carrying SSU increased in copy number, which is consistent with selection to increase abundance; however, two other studied chromosomes displayed much lower than expected among-line variance. Our models suggest that balancing selection is sufficient to explain the observed maintenance of chromosome copy during asexual reproduction.     

Variation of Expression Levels of Seven Housekeeping Genes at Different Life-History Stages in Porphyra yezoensis

In order to identify the optimal internal control for relative real-time PCR when studying target gene expression in the red alga Porphyra yezoensis, we quantified the expression of seven housekeeping genes (18S ribosomal RNA, 30S ribosomal protein S8, Polyubiquitin-2, Glyceraldehyde-3-phosphate dehydrogenase, Elongation factor 1-alpha, Beta-tubulin and Actin 3) at different life-history stages. Absolute quantification was done by normalization to total RNA quantity and by normalization to genomic DNA quantity. We used these two normalization approaches, comparing the differences of expression levels of all candidate housekeeping genes between any two generations and across three life-history stages (filamentous sporophytes, leafy gametophytes and conchospores). We found GAPDH had the best stability in all cases and we recommend that GAPDH be considered as a potential internal control for gene expression studies at different life-history stages in P. yezoensis.     

Effect of Within-field Variation in Soil Texture on Heterodera glycines and Soybean Yield

The influence of soil texture on Soybean yield in the presence of Heterodera glycines was investigated by comparing yields of susceptible cultivars with a resistant cultivar for 2 years. Soybean yield was negatively correlated with increasing sand content (P = 0.05). Yields of susceptible cultivars were suppressed with increasing sand content. Final nematode population densities were lowest in plots with greatest sand content. Soybean infection by SCN, as determined by the number of cysts 30 days after planting, was not consistently related to soil texture over 2 years. Initial nematode population density was positively related to soybean yield the first year and negatively related to soybean yield the second, probably a result of greater yield suppression by H. glycines in plots with greater sand content.     

Polymorphisms in Anopheles gambiae Immune Genes Associated with Natural Resistance to Plasmodium falciparum

Many genes involved in the immune response of Anopheles gambiae, the main malaria vector in Africa, have been identified, but whether naturally occurring polymorphisms in these genes underlie variation in resistance to the human malaria parasite, Plasmodium falciparum, is currently unknown. Here we carried out a candidate gene association study to identify single nucleotide polymorphisms (SNPs) associated with natural resistance to P. falciparum. A. gambiae M form mosquitoes from Cameroon were experimentally challenged with three local wild P. falciparum isolates. Statistical associations were assessed between 157 SNPs selected from a set of 67 A. gambiae immune-related genes and the level of infection. Isolate-specific associations were accounted for by including the effect of the isolate in the analysis. Five SNPs were significantly associated to the infection phenotype, located within or upstream of AgMDL1, CEC1, Sp PPO activate, Sp SNAKElike, and TOLL6. Low overall and local linkage disequilibrium indicated high specificity in the loci found. Association between infection phenotype and two SNPs was isolate-specific, providing the first evidence of vector genotype by parasite isolate interactions at the molecular level. Four SNPs were associated to either oocyst presence or load, indicating that the genetic basis of infection prevalence and intensity may differ. The validity of the approach was verified by confirming the functional role of Sp SNAKElike in gene silencing assays. These results strongly support the role of genetic variation within or near these five A. gambiae immune genes, in concert with other genes, in natural resistance to P. falciparum. They emphasize the need to distinguish between infection prevalence and intensity and to account for the genetic specificity of vector-parasite interactions in dissecting the genetic basis of Anopheles resistance to human malaria.     

Genes as Early Responders Regulate Quorum-Sensing and Control Bacterial Cooperation in Pseudomonas aeruginosa

Quorum-sensing (QS) allows bacterial communication to coordinate the production of extracellular products essential for population fitness at higher cell densities. It has been generally accepted that a significant time duration is required to reach appropriate cell density to activate the relevant quiescent genes encoding these costly but beneficial public goods. Which regulatory genes are involved and how these genes control bacterial communication at the early phases are largely un-explored. By determining time-dependent expression of QS-related genes of the opportunistic pathogen Pseudomonas aerugionsa, we show that the induction of social cooperation could be critically influenced by environmental factors to optimize the density of population. In particular, small regulatory RNAs (RsmY and RsmZ) serving as early responders, can promote the expression of dependent genes (e.g. lasR) to boost the synthesis of intracellular enzymes and coordinate instant cooperative behavior in bacterial cells. These early responders, acting as a rheostat to finely modulate bacterial cooperation, which may be quickly activated under environment threats, but peter off when critical QS dependent genes are fully functional for cooperation. Our findings suggest that RsmY and RsmZ critically control the timing and levels of public goods production, which may have implications in sociomicrobiology and infection control.     

Identification of Novel QTLs for Isolate-Specific Partial Resistance to Plasmodiophora brassicae in Brassica rapa

Plasmodiophora brassicae, the causal agent of clubroot disease of the Brassica crops, is widespread in the world. Quantitative trait loci (QTLs) for partial resistance to 4 different isolates of P. brassicae (Pb2, Pb4, Pb7, and Pb10) were investigated using a BC₁F₁ population from a cross between two subspecies of Brassica rapa, i.e. Chinese cabbage inbred line C59-1 as a susceptible recurrent parent and turnip inbred line ECD04 as a resistant donor parent. The BC₁F₂ families were assessed for resistance under controlled conditions. A linkage map constructed with simple sequence repeats (SSR), unigene-derived microsatellite (UGMS) markers, and specific markers linked to published clubroot resistance (CR) genes of B. rapa was used to perform QTL mapping. A total of 6 QTLs residing in 5 CR QTL regions of the B. rapa chromosomes A01, A03, and A08 were identified to account for 12.2 to 35.2% of the phenotypic variance. Two QTL regions were found to be novel except for 3 QTLs in the respective regions of previously identified Crr1, Crr2, and Crr3. QTL mapping results indicated that 1 QTL region was common for partial resistance to the 2 isolates of Pb2 and Pb7, whereas the others were specific for each isolate. Additionally, synteny analysis between B. rapa and Arabidopsis thaliana revealed that all CR QTL regions were aligned to a single conserved crucifer blocks (U, F, and R) on 3 Arabidopsis chromosomes where 2 CR QTLs were detected in A. thaliana. These results suggest that some common ancestral genomic regions were involved in the evolution of CR genes in B. rapa.     

BTG Interacts with Retinoblastoma to Control Cell Fate in Dictyostelium

Background

In the genesis of many tissues, a phase of cell proliferation is followed by cell cycle exit and terminal differentiation. The latter two processes overlap: genes involved in the cessation of growth may also be important in triggering differentiation. Though conceptually distinct, they are often causally related and functional interactions between the cell cycle machinery and cell fate control networks are fundamental to coordinate growth and differentiation. A switch from proliferation to differentiation may also be important in the life cycle of single-celled organisms, and genes which arose as regulators of microbial differentiation may be conserved in higher organisms. Studies in microorganisms may thus contribute to understanding the molecular links between cell cycle machinery and the determination of cell fate choice networks.

Methodology/Principal Findings

Here we show that in the amoebozoan D. discoideum, an ortholog of the metazoan antiproliferative gene btg controls cell fate, and that this function is dependent on the presence of a second tumor suppressor ortholog, the retinoblastoma-like gene product. Specifically, we find that btg-overexpressing cells preferentially adopt a stalk cell (and, more particularly, an Anterior-Like Cell) fate. No btg-dependent preference for ALC fate is observed in cells in which the retinoblastoma-like gene has been genetically inactivated. Dictyostelium btg is the only example of non-metazoan member of the BTG family characterized so far, suggesting that a genetic interaction between btg and Rb predated the divergence between dictyostelids and metazoa.

Conclusions/Significance

While the requirement for retinoblastoma function for BTG antiproliferative activity in metazoans is known, an interaction of these genes in the control of cell fate has not been previously documented. Involvement of a single pathway in the control of mutually exclusive processes may have relevant implication in the evolution of multicellularity.