Allele-specific assay reveals functional variation in the chalcone synthase promoter of Arabidopsis thaliana that is compatible with neutral evolution

Transcriptome Analysis and Intraspecific Variation in Spanish Fir (Abies pinsapo Boiss.)

Spanish fir (Abies pinsapo Boiss.) is an endemic, endangered tree that has been scarcely investigated at the molecular level. In this work, the transcriptome of Spanish fir was assembled, providing a large catalog of expressed genes (22,769), within which a high proportion were full-length transcripts (12,545). This resource is valuable for functional genomics studies and genome annotation in this relict conifer species. Two intraspecific variations of A. pinsapo can be found within its largest population at the Sierra de las Nieves National Park: one with standard green needles and another with bluish-green needles. To elucidate the causes of both phenotypes, we studied different physiological and molecular markers and transcriptome profiles in the needles. "Green" trees showed higher electron transport efficiency and enhanced levels of chlorophyll, protein, and total nitrogen in the needles. In contrast, needles from "bluish" trees exhibited higher contents of carotenoids and cellulose. These results agreed with the differential transcriptomic profiles, suggesting an imbalance in the nitrogen status of "bluish" trees. Additionally, gene expression analyses suggested that these differences could be associated with different epigenomic profiles. Taken together, the reported data provide new transcriptome resources and a better understanding of the natural variation in this tree species, which can help improve guidelines for its conservation and the implementation of adaptive management strategies under climatic change.

Compensatory trans-regulatory alleles minimizing variation in TDH3 expression are common within Saccharomyces cerevisiae

Heritable variation in gene expression is common within species. Much of this variation is due to genetic differences outside of the gene with altered expression and is trans-acting. This trans-regulatory variation is often polygenic, with individual variants typically having small effects, making the genetic architecture and evolution of trans-regulatory variation challenging to study. Consequently, key questions about trans-regulatory variation remain, including the variability of trans-regulatory variation within a species, how selection affects trans-regulatory variation, and how trans-regulatory variants are distributed throughout the genome and within a species. To address these questions, we isolated and measured trans-regulatory differences affecting TDH3 promoter activity among 56 strains of Saccharomyces cerevisiae, finding that trans-regulatory backgrounds varied approximately twofold in their effects on TDH3 promoter activity. Comparing this variation to neutral models of trans-regulatory evolution based on empirical measures of mutational effects revealed that despite this variability in the effects of trans-regulatory backgrounds, stabilizing selection has constrained trans-regulatory differences within this species. Using a powerful quantitative trait locus mapping method, we identified ∼100 trans-acting expression quantitative trait locus in each of three crosses to a common reference strain, indicating that regulatory variation is more polygenic than previous studies have suggested. Loci altering expression were located throughout the genome, and many loci were strain specific. This distribution and prevalence of alleles is consistent with recent theories about the genetic architecture of complex traits. In all mapping experiments, the nonreference strain alleles increased and decreased TDH3 promoter activity with similar frequencies, suggesting that stabilizing selection maintained many trans-acting variants with opposing effects. This variation may provide the raw material for compensatory evolution and larger scale regulatory rewiring observed in developmental systems drift among species.

Gene regulatory divergence between locally adapted ecotypes in their native habitats

Local adaptation is a key driver of ecological specialization and the formation of new species. Despite its importance, the evolution of gene regulatory divergence among locally adapted populations is poorly understood, especially how that divergence manifests in nature. Here, we evaluate gene expression divergence and allele-specific gene expression responses for locally adapted coastal perennial and inland annual accessions of the yellow monkeyflower, Mimulus guttatus, in a field reciprocal transplant experiment. Overall, 6765 (73%) of surveyed genes were differentially expressed between coastal and inland habitats, while 7213 (77%) were differentially expressed between the coastal perennial and inland annual accessions. Cis-regulatory variation was pervasive, affecting 79% (5532) of differentially expressed genes. We detected trans effects for 52% (3611) of differentially expressed genes. Expression plasticity of alleles across habitats (G × E interactions) appears to be relatively common (affecting 18% of transcripts) and could minimize fitness trade-offs at loci that contribute to local adaptation. We also found evidence that at least one chromosomal inversion may act as supergene by holding together haplotypes of differentially expressed genes, but this pattern depends on habitat context. Our results highlight multiple key patterns regarding the relationship between gene expression and the evolution of locally adapted populations.

Divergence of annual and perennial species in the Brassicaceae and the contribution of cis-acting variation at FLC orthologues

Variation in life history contributes to reproductive success in different environments. Divergence of annual and perennial angiosperm species is an extreme example that has occurred frequently. Perennials survive for several years and restrict the duration of reproduction by cycling between vegetative growth and flowering, whereas annuals live for 1 year and flower once. We used the tribe Arabideae (Brassicaceae) to study the divergence of seasonal flowering behaviour among annual and perennial species. In perennial Brassicaceae, orthologues of FLOWERING LOCUS C (FLC), a floral inhibitor in Arabidopsis thaliana, are repressed by winter cold and reactivated in spring conferring seasonal flowering patterns, whereas in annuals, they are stably repressed by cold. We isolated FLC orthologues from three annual and two perennial Arabis species and found that the duplicated structure of the A. alpina locus is not required for perenniality. The expression patterns of the genes differed between annuals and perennials, as observed among Arabidopsis species, suggesting a broad relevance of these patterns within the Brassicaceae. Also analysis of plants derived from an interspecies cross of A. alpina and annual A. montbretiana demonstrated that cis-regulatory changes in FLC orthologues contribute to their different transcriptional patterns. Sequence comparisons of FLC orthologues from annuals and perennials in the tribes Arabideae and Camelineae identified two regulatory regions in the first intron whose sequence variation correlates with divergence of the annual and perennial expression patterns. Thus, we propose that related cis-acting changes in FLC orthologues occur independently in different tribes of the Brassicaceae during life history evolution.

Differential fine-tuning of gene expression regulation in coffee leaves by CcDREB1D promoter haplotypes under water deficit

Despite the importance of the DREB1D gene (also known as CBF4) in plant responses to water deficit and cold stress, studies analysing its regulation by transgenic approaches are lacking. In the current work, a functional study of three CcDREB1D promoter haplotypes (named HP15, HP16 and HP17) isolated from drought-tolerant and drought-sensitive clones of Coffea canephora was carried out in plants of C. arabica stably transformed by Agrobacterium tumefaciens by analysing their ability to regulate the expression of the uidA reporter gene in response to water deficit mimicked by polyethylene glycol (-2.0 MPa) and low relative humidity treatments. A deletion analysis of their corresponding 5'-upstream regions revealed increased specificity of β-glucuronidase activity in the polyethylene glycol and low relative humidity treatments, with high expression in leaf mesophyll and guard cells in full-length constructs. RT-qPCR assays also revealed that the HP16 haplotype (specific to clone tolerant to water deficit) had stronger and earlier activity compared with the HP15 and HP17 haplotypes. As most of the cis-regulatory elements involved in ABA-dependent and -independent networks, tissue specificity and light regulation are common to these haplotypes, we propose that their organization, as well as the nucleic acid polymorphisms present outside these boxes, may play a role in modulating activities of DREB1D promoters in guard cells.

Different Alleles of a Gene Encoding Leucoanthocyanidin Reductase (PaLAR3) Influence Resistance against the Fungus Heterobasidion parviporum in Picea abies

Despite the fact that fungal diseases are a growing menace for conifers in modern silviculture, only a very limited number of molecular markers for pathogen resistance have been validated in conifer species. A previous genetic study indicated that the resistance of Norway spruce (Picea abies) to Heterobasidion annosum s.l., a pathogenic basidiomycete species complex, is linked to a quantitative trait loci that associates with differences in fungal growth in sapwood (FGS) that includes a gene, PaLAR3, which encodes a leucoanthocyanidin reductase. In this study, gene sequences showed the presence of two PaLAR3 allelic lineages in P. abies. Higher resistance was associated with the novel allele, which was found in low frequency in the four P. abies populations that we studied. Norway spruce plants carrying at least one copy of the novel allele showed a significant reduction in FGS after inoculation with Heterobasidion parviporum compared to their half-siblings carrying no copies, indicating dominance of this allele. The amount of (+) catechin, the enzymatic product of PaLAR3, was significantly higher in bark of trees homozygous for the novel allele. Although we observed that the in vitro activities of the enzymes encoded by the two alleles were similar, we could show that allele-specific transcript levels were significantly higher for the novel allele, indicating that regulation of gene expression is responsible for the observed effects in resistance, possibly caused by differences in cis-acting elements that we observe in the promoter region of the two alleles.

Differential transcriptional activity of SAD, FAD2 and FAD3 desaturase genes in developing seeds of linseed contributes to varietal variation in α-linolenic acid content

Linseed or flax (Linum usitatissimum L.) varieties differ markedly in their seed α-linolenic acid (ALA) levels. Fatty acid desaturases play a key role in accumulating ALA in seed. We performed fatty acid (FA) profiling of various seed developmental stages of ten Indian linseed varieties including one mutant variety. Depending on their ALA contents, these varieties were grouped under high ALA and low ALA groups. Transcript profiling of six microsomal desaturase genes (SAD1, SAD2, FAD2, FAD2-2, FAD3A and FAD3B), which act sequentially in the fatty acid desaturation pathway, was performed using real-time PCR. We observed gene specific as well as temporal expression pattern for all the desaturases and their differential expression profiles corresponded well with the variation in FA accumulation in the two groups. Our study points to efficient conversion of intermediate FAs [stearic (SA), oleic (OA) and linoleic acids (LA)] to the final product, ALA, due to efficient action of all the desaturases in high ALA group. While in the low ALA group, even though the initial conversion up to OA was efficient, later conversions up to ALA seemed to be inefficient, leading to higher accumulation of OA and LA instead of ALA. We sequenced the six desaturase genes from the ten varieties and observed that variation in the amino acid (AA) sequences of desaturases was not responsible for differential ALA accumulation, except in the mutant variety TL23 with very low (<2%) ALA content. In TL23, a point mutation in the FAD3A gene resulted into a premature stop codon generating a truncated protein with 291 AA.

Social environment influences the relationship between genotype and gene expression in wild baboons

Variation in the social environment can have profound effects on survival and reproduction in wild social mammals. However, we know little about the degree to which these effects are influenced by genetic differences among individuals, and conversely, the degree to which social environmental variation mediates genetic reaction norms. To better understand these relationships, we investigated the potential for dominance rank, social connectedness and group size to modify the effects of genetic variation on gene expression in the wild baboons of the Amboseli basin. We found evidence for a number of gene-environment interactions (GEIs) associated with variation in the social environment, encompassing social environments experienced in adulthood as well as persistent effects of early life social environment. Social connectedness, maternal dominance rank and group size all interacted with genotype to influence gene expression in at least one sex, and either in early life or in adulthood. These results suggest that social and behavioural variation, akin to other factors such as age and sex, can impact the genotype-phenotype relationship. We conclude that GEIs mediated by the social environment are important in the evolution and maintenance of individual differences in wild social mammals, including individual differences in responses to social stressors.

The evolution of phenylpropanoid metabolism in the green lineage

Phenolic secondary metabolites are only produced by plants wherein they play important roles in both biotic and abiotic defense in seed plants as well as being potentially important bioactive compounds with both nutritional and medicinal benefits reported for animals and humans as a consequence of their potent antioxidant activity. During the long evolutionary period in which plants have adapted to the environmental niches in which they exist (and especially during the evolution of land plants from their aquatic algal ancestors), several strategies such as gene duplication and convergent evolution have contributed to the evolution of this pathway. In this respect, diversity and redundancy of several key genes of phenolic secondary metabolism such as polyketide synthases, cytochrome P450s, Fe(2+)/2-oxoglutarate-dependent dioxygenases and UDP-glycosyltransferases have played an essential role. Recent technical developments allowing affordable whole genome sequencing as well as a better inventory of species-by-species chemical diversity have resulted in a dramatic increase in the number of tools we have to assess how these pathways evolved. In parallel, reverse genetics combined with detailed molecular phenotyping is allowing us to elucidate the functional importance of individual genes and metabolites and by this means to provide further mechanistic insight into their biological roles. In this review, phenolic metabolite-related gene sequences (for a total of 65 gene families including shikimate biosynthetic genes) are compared across 23 independent species, and the phenolic metabolic complement of various plant species are compared with one another, in attempt to better understand the evolution of diversity in this crucial pathway.

Flagellin perception varies quantitatively in Arabidopsis thaliana and its relatives

Much is known about the evolution of plant immunity components directed against specific pathogen strains: They show pervasive functional variation and have the potential to coevolve with pathogen populations. However, plants are effectively protected against most microbes by generalist immunity components that detect conserved pathogen-associated molecular patterns (PAMPs) and control the onset of PAMP-triggered immunity. In Arabidopsis thaliana, the receptor kinase flagellin sensing 2 (FLS2) confers recognition of bacterial flagellin (flg22) and activates a manifold defense response. To decipher the evolution of this system, we performed functional assays across a large set of A. thaliana genotypes and Brassicaceae relatives. We reveal extensive variation in flg22 perception, most of which results from changes in protein abundance. The observed variation correlates with both the severity of elicited defense responses and bacterial proliferation. We analyzed nucleotide variation segregating at FLS2 in A. thaliana and detected a pattern of variation suggestive of the rapid fixation of a novel adaptive allele. However, our study also shows that evolution at the receptor locus alone does not explain the evolution of flagellin perception; instead, components common to pathways downstream of PAMP perception likely contribute to the observed quantitative variation. Within and among close relatives, PAMP perception evolves quantitatively, which contrasts with the changes in recognition typically associated with the evolution of R genes.

Evolutionary and functional analyses of cytochrome P4501A promoter polymorphisms in natural populations

The functional importance of variable, transcriptional regulatory sequences within and among natural populations is largely unexplored. We analysed the cytochrome P4501A (CYP1A) promoter in three populations of the minnow, Fundulus heteroclitus, because two SNPs in the promoter and first intron of CYP1A are under selection in populations adapted to pollutants. To define the importance of these SNPs, 1630 bp of the CYP1A promoter and first intron and exon were sequenced in eight individuals from three populations: a population from a polluted environment resistant to some aromatic pollutants and two flanking reference populations. CYP1A is induced by many aromatic pollutants, but in populations adapted to pollutants, CYP1A has been shown to be refractory to induction. We were interested in understanding whether variation in the CYP1A promoter explains mechanism(s) of adaptation to these aromatic pollutants. The CYP1A promoter was extremely variable (an average of 9.3% of the promoter nucleotides varied among all populations) and exhibited no fixed differences between populations. As CYP1A is poorly inducible in adapted fish, we hypothesized that CYP1A promoter regions might vary functionally between populations. Unexpectedly, in vitro analysis showed significantly greater transcription from CYP1A promoters found in the population from the polluted environment relative to promoters found in both reference populations. Thus, despite extensive variation among populations and lack of fixed differences between populations, individuals from a polluted environment have significantly enhanced promoter activity. These data demonstrate that intraspecific variation, which provides the raw material for natural selection to act on, can occur while maintaining promoter function.

Widespread interspecific divergence in cis-regulation of transposable elements in the Arabidopsis genus

Transposable elements (TEs) are so abundant and variable that they count among the most important mutational sources in genomes. Nonetheless, little is known about the genetics of their variation in activity or silencing across closely related species. Here, we demonstrate that regulation of TE genes can differ dramatically between the two closely related Arabidopsis species A. thaliana and A. lyrata. In leaf and floral tissues of F1 interspecific hybrids, about 47% of TEs show allele-specific expression, with the A. lyrata copy being generally expressed at higher level. We confirm that TEs are generally expressed in A. lyrata but not in A. thaliana. Allele-specific differences in TE expression are associated with divergence in epigenetic modifications like DNA and histone methylation between species as well as with sequence divergence. Our data demonstrate that A. thaliana silences TEs much better than A. lyrata. For long terminal repeat retrotransposons, these differences are more pronounced for younger insertions. Interspecific differences in TE silencing may have a great impact on genome size changes.

Comparative analysis of LEA-like 11-24 gene expression and regulation in related plant species within the Linderniaceae that differ in desiccation tolerance

The resurrection plant Craterostigma plantagineum is able to withstand desiccation of its vegetative tissues and is found in areas with variable water availability. The closely related species Lindernia brevidens and Lindernia subracemosa are both endemic to montane rainforests of coastal Africa, but remarkably L. brevidens is tolerant to desiccation. We studied the regulation of the desiccation-related LEA-like 11-24 gene at multiple levels in closely related species in order to investigate the conservation of mechanisms involved in desiccation tolerance. The dehydration-responsive transcription of the LEA-like 11-24 gene is differentially regulated in these plants. Comparison of the LEA-like 11-24 core promoter regions revealed that promoters have different activities, but some functional cis-acting elements are conserved between species. Upon dehydration, LEA-like 11-24 proteins are phosphorylated at different levels and phosphorylation sites are not conserved among the three LEA-like 11-24 proteins. Differences in the regulation of the LEA-like 11-24 gene in the studied plant species appear to be the result of mutations that occurred during evolution. We postulate that L. brevidens will eventually lose the ability to survive vegetative desiccation, given that this trait appears not to be essential for survival.

Promoting the promoter

Recent evolutionary studies clearly indicate that evolution is mainly driven by changes in the complex mechanisms of gene regulation and not solely by polymorphism in protein-encoding genes themselves. After a short description of the cis-regulatory mechanism, we intend in this review to argue that by applying newly available technologies and by merging research areas such as evolutionary and developmental biology, population genetics, ecology and molecular cell biology it is now possible to study evolution in an integrative way. We contend that, by analysing the effects of promoter sequence variation on phenotypic diversity in natural populations, we will soon be able to break the barrier between the study of extant genetic variability and the study of major developmental changes. This will lead to an integrative view of evolution at different scales. Because of their sessile nature and their continuous development, plants must permanently regulate their gene expression to react to their environment, and can, therefore, be considered as a remarkable model for these types of studies.

Allele-specific gene expression in a wild nonhuman primate population

Natural populations hold enormous potential for evolutionary genetic studies, especially when phenotypic, genetic and environmental data are all available on the same individuals. However, untangling the genotype-phenotype relationship in natural populations remains a major challenge. Here, we describe results of an investigation of one class of phenotype, allele-specific gene expression (ASGE), in the well-studied natural population of baboons of the Amboseli basin, Kenya. ASGE measurements identify cases in which one allele of a gene is overexpressed relative to the alternative allele of the same gene, within individuals, thus providing a control for background genetic and environmental effects. Here, we characterize the incidence of ASGE in the Amboseli baboon population, focusing on the genetic and environmental contributions to ASGE in a set of eleven genes involved in immunity and defence. Within this set, we identify evidence for common ASGE in four genes. We also present examples of two relationships between cis-regulatory genetic variants and the ASGE phenotype. Finally, we identify one case in which this relationship is influenced by a novel gene-environment interaction. Specifically, the dominance rank of an individual's mother during its early life (an aspect of that individual's social environment) influences the expression of the gene CCL5 via an interaction with cis-regulatory genetic variation. These results illustrate how environmental and ecological data can be integrated into evolutionary genetic studies of functional variation in natural populations. They also highlight the potential importance of early life environmental variation in shaping the genetic architecture of complex traits in wild mammals.

Cryptic variation between species and the basis of hybrid performance

Studies on natural variation in gene expression and its phenotypic effects provide fresh insights into the origins of vigour and sterility in species hybrids.

Crosses between closely related species give two contrasting results. One result is that species hybrids may be inferior to their parents, for example, being less fertile [1]. The other is that F1 hybrids may display superior performance (heterosis), for example with increased vigour [2]. Although various hypotheses have been proposed to account for these two aspects of hybridisation, their biological basis is still poorly understood [3]. To gain further insights into this issue, we analysed the role that variation in gene expression may play. We took a conserved trait, flower asymmetry in Antirrhinum, and determined the extent to which the underlying regulatory genes varied in expression among closely related species. We show that expression of both genes analysed, CYC and RAD, varies significantly between species because of cis-acting differences. By making a quantitative genotype-phenotype map, using a range of mutant alleles, we demonstrate that the species lie on a plateau in gene expression-morphology space, so that the variation has no detectable phenotypic effect. However, phenotypic differences can be revealed by shifting genotypes off the plateau through genetic crosses. Our results can be readily explained if genomes are free to evolve within an effectively neutral zone in gene expression space. The consequences of this drift will be negligible for individual loci, but when multiple loci across the genome are considered, we show that the variation may have significant effects on phenotype and fitness, causing a significant drift load. By considering these consequences for various gene-expression–fitness landscapes, we conclude that F1 hybrids might be expected to show increased performance with regard to conserved traits, such as basic physiology, but reduced performance with regard to others. Thus, our study provides a new way of explaining how various aspects of hybrid performance may arise through natural variation in gene activity.

A major conundrum in biology is why hybrids between species display two opposing features. On the one hand, hybrids are often more vigorous or productive than their parents, a phenomenon called hybrid vigor or hybrid superiority. On the other hand they often show reduced vigour and fertility, known as hybrid inferiority. Various theories have been proposed to account for these two aspects of hybrid performance, yet we still lack a coherent account of how these conflicting characteristics arise. To address this issue, we looked at the role that variation in gene expression between parental species may play. By measuring this variation and its effect on phenotype, we show that expression for specific genes may be free to vary during evolution within particular bounds. Although such variation may have little phenotypic effect when each locus is considered individually, the collective effect of variation across multiple genes may become highly significant. Using arguments from theoretical population genetics we show how these effects might lead to both hybrid superiority and inferiority, providing fresh insights into the age-old problem of hybrid performance.

Key considerations for measuring allelic expression on a genomic scale using high-throughput sequencing

Differences in gene expression are thought to be an important source of phenotypic diversity, so dissecting the genetic components of natural variation in gene expression is important for understanding the evolutionary mechanisms that lead to adaptation. Gene expression is a complex trait that, in diploid organisms, results from transcription of both maternal and paternal alleles. Directly measuring allelic expression rather than total gene expression offers greater insight into regulatory variation. The recent emergence of high-throughput sequencing offers an unprecedented opportunity to study allelic transcription at a genomic scale for virtually any species. By sequencing transcript pools derived from heterozygous individuals, estimates of allelic expression can be directly obtained. The statistical power of this approach is influenced by the number of transcripts sequenced and the ability to unambiguously assign individual sequence fragments to specific alleles on the basis of transcribed nucleotide polymorphisms. Here, using mathematical modelling and computer simulations, we determine the minimum sequencing depth required to accurately measure relative allelic expression and detect allelic imbalance via high-throughput sequencing under a variety of conditions. We conclude that, within a species, a minimum of 500-1000 sequencing reads per gene are needed to test for allelic imbalance, and consequently, at least five to 10 millions reads are required for studying a genome expressing 10 000 genes. Finally, using 454 sequencing, we illustrate an application of allelic expression by testing for cis-regulatory divergence between closely related Drosophila species.

Assessing the influence of adjacent gene orientation on the evolution of gene upstream regions in Arabidopsis thaliana

The orientation of flanking genes may influence the evolution of intergenic regions in which cis-regulatory elements are likely to be located: divergently transcribed genes share their 5' regions, resulting either in smaller "private" spaces or in overlapping regulatory elements. Thus, upstream sequences of divergently transcribed genes (bi-directional upstream regions, or URs) may be more constrained than those of uni-directional gene pairs. We investigated this effect by analyzing nucleotide variation segregating within and between Arabidopsis species. Compared to uni-directional URs, bi-directional URs indeed display lower population mutation rate, as well as more low-frequency polymorphisms. Furthermore, we find that bi-directional regions undergo selection for the maintenance of intergenic distance. Altogether, however, we observe considerable variation in evolutionary rates, with putative signatures of selection on two uni-directional upstream regions.

Genomic features that predict allelic imbalance in humans suggest patterns of constraint on gene expression variation

Variation in gene expression is an important contributor to phenotypic diversity within and between species. Although this variation often has a genetic component, identification of the genetic variants driving this relationship remains challenging. In particular, measurements of gene expression usually do not reveal whether the genetic basis for any observed variation lies in cis or in trans to the gene, a distinction that has direct relevance to the physical location of the underlying genetic variant, and which may also impact its evolutionary trajectory. Allelic imbalance measurements identify cis-acting genetic effects by assaying the relative contribution of the two alleles of a cis-regulatory region to gene expression within individuals. Identification of patterns that predict commonly imbalanced genes could therefore serve as a useful tool and also shed light on the evolution of cis-regulatory variation itself. Here, we show that sequence motifs, polymorphism levels, and divergence levels around a gene can be used to predict commonly imbalanced genes in a human data set. Reduction of this feature set to four factors revealed that only one factor significantly differentiated between commonly imbalanced and nonimbalanced genes. We demonstrate that these results are consistent between the original data set and a second published data set in humans obtained using different technical and statistical methods. Finally, we show that variation in the single allelic imbalance-associated factor is partially explained by the density of genes in the region of a target gene (allelic imbalance is less probable for genes in gene-dense regions), and, to a lesser extent, the evenness of expression of the gene across tissues and the magnitude of negative selection on putative regulatory regions of the gene. These results suggest that the genomic distribution of functional cis-regulatory variants in the human genome is nonrandom, perhaps due to local differences in evolutionary constraint.

Asymmetric allele-specific expression in relation to developmental variation and drought stress in barley hybrids

In the present study, we analysed allele-specific expression (ASE) in the selfing species barley to assess the frequency of cis-acting regulatory variation and the effects of genetic background, developmental differences and drought stress on allelic expression levels. We measured ASE ratios in 30 genes putatively involved in stress responses in five hybrids and their reciprocals, namely Hordeum spontaneum 41-1/Alexis (HAl), Hordeum spontaneum 41-1/Arta (HAr), Sloop/WI3408 (SW), Tadmor/Sloop (TS) and Tadmor/WI3408 (TW). In order to detect cis-acting variation related to drought and developmental changes, the barley hybrids were grown under control and water-limited conditions, and leaf tissue was harvested at two developmental stages. The analysis demonstrated that more than half of the genes measured (63%) showed allelic differences in expression of up to 19-fold due to cis-regulatory variation in at least one cross by treatment/stage combination. Drought stress induced changes in allelic expression ratios, indicating differences between drought responsive cis-elements. In addition, ASE differences between developmental stages suggested the presence of cis-acting elements interacting with developmental cues. We were also able to demonstrate that the levels and frequency of allelic imbalance and hence differences in cis-regulatory elements are correlated with the genetic divergence between the parental lines, but may also arise as an adaptation to diverse habitats. Our findings suggest that cis-regulatory variation is a common phenomenon in barley, and may provide a molecular basis of transgression. Differential expression of near-isogenic members of the same gene family could potentially result in hybrid lines out performing their parents in terms of expression level, timing and response to developmental and environmental cues. Identification and targeted manipulation of cis-regulatory elements will assist in breeding improved crops with a better adaptation to changing environments.

A sterile-female technique proposed for control of Striga hermonthica and other intractable weeds: advantages, shortcomings and risk management

Weeds have posed intractable challenges to farmers since the dawn of agriculture. This article describes in detail a proposed control strategy based on the introduction of genes conferring female sterility into the genome of an intractable target weed. Spread of these genes through target populations via pollen would be facilitated by their incorporation within active transposable elements. Advantages (e.g. self-dissemination, self-proliferation, target specificity) and shortcomings (e.g. high cost, long project incubation period, limited range of possible targets) of this strategy are discussed in depth, as are assessment and management of its attendant biological and ecological risks, such as the risk of introduced genes spreading to non-target species. The parasitic weed Striga hermonthica (Del.) Benth. is examined as a potential target.

Cis-regulatory variation is typically polyallelic in Drosophila

Gene expression levels vary heritably, with approximately 25-35% of the loci affecting expression acting in cis. We characterized standing cis-regulatory variation among 16 wild-derived strains of Drosophila melanogaster. Our experiment's robust biological and technical replication enabled precise estimates of variation in allelic expression on a high-throughput SNP genotyping platform. We observed concordant, significant differential allelic expression (DAE) in 7/10 genes queried with multiple SNPs, and every member of a set of eight additional, one-assay genes suggest significant DAE. Four of the high-confidence, multiple-assay genes harbor three or more statistically distinguishable allelic classes, often at intermediate frequency. Numerous intermediate-frequency, detectable regulatory polymorphisms cast doubt on a model in which cis-acting variation is a product of deleterious mutations of large effect. Comparing our data to predictions of population genetics theory using coalescent simulations, we estimate that a typical gene harbors 7-15 cis-regulatory sites (nucleotides) at which a selectively neutral mutation would elicit an observable expression phenotype. If standing cis-regulatory variation is actually slightly deleterious, the true mutational target size is larger.

Gametophytically alloplasmic CMS line of rice (Oryza sativa L.) with variant orfH79 haplotype corresponds to specific fertility restorer

For years discovery and identification of the cytoplasmic male sterility (CMS) resource in wild rice is the most intriguing events in breeding field. orfH79, a chimeric gene in mitochondria, has been suggested being the determinant for Honglian CMS in rice. In this report orfH79 gene as molecular marker to screen the wild rice, we found eight accessions with orfH79 gene in the total 42 investigated objects. Sequence analysis revealed that there were a total of nine nucleotide substitutions resulting in the change of nine amino acids in the newly identified orfH79 in wild rice, which further fell into seven haplotypes. In order to investigate the underlying relationship between orfH79 haplotypes and the corresponding fertility restorers, four accessions were selected with different orfH79 haplotype as female parents to hybridize the Honglian maintainer line, Yuetai B. After eight consecutive recurrent backcrosses, four alloplasmic CMS lines with different orfH79 haplotype were developed. Microscopic observation exhibited that their pollen grains were spherical and clear in 1% I(2)-KI solution same as that of Honglian CMS line. Moreover, these four CMS lines displayed various fertility restoring model through test cross, suggesting that each orfH79 haplotye represents a new CMS type and corresponds to their specific Rf allele.

Evolution of primate gene expression: drift and corrective sweeps?

Changes in gene expression play an important role in species' evolution. Earlier studies uncovered evidence that the effect of mutations on expression levels within the primate order is skewed, with many small downregulations balanced by fewer but larger upregulations. In addition, brain-expressed genes appeared to show an increased rate of evolution on the branch leading to human. However, the lack of a mathematical model adequately describing the evolution of gene expression precluded the rigorous establishment of these observations. Here, we develop mathematical tools that allow us to revisit these earlier observations in a model-testing and inference framework. We introduce a model for skewed gene-expression evolution within a phylogenetic tree and use a separate model to account for biological or experimental outliers. A Bayesian Markov chain Monte Carlo inference procedure allows us to infer the phylogeny and other evolutionary parameters, while quantifying the confidence in these inferences. Our results support previous observations; in particular, we find strong evidence for a sustained positive skew in the distribution of gene-expression changes in primate evolution. We propose a "corrective sweep" scenario to explain this phenomenon.

The new mutation theory of phenotypic evolution

Recent studies of developmental biology have shown that the genes controlling phenotypic characters expressed in the early stage of development are highly conserved and that recent evolutionary changes have occurred primarily in the characters expressed in later stages of development. Even the genes controlling the latter characters are generally conserved, but there is a large component of neutral or nearly neutral genetic variation within and between closely related species. Phenotypic evolution occurs primarily by mutation of genes that interact with one another in the developmental process. The enormous amount of phenotypic diversity among different phyla or classes of organisms is a product of accumulation of novel mutations and their conservation that have facilitated adaptation to different environments. Novel mutations may be incorporated into the genome by natural selection (elimination of preexisting genotypes) or by random processes such as genetic and genomic drift. However, once the mutations are incorporated into the genome, they may generate developmental constraints that will affect the future direction of phenotypic evolution. It appears that the driving force of phenotypic evolution is mutation, and natural selection is of secondary importance.

In vitro assays fail to predict in vivo effects of regulatory polymorphisms

A typical paradigm in the investigation of complex human disease is to assess the effects of cis-regulatory polymorphisms implicated in association studies on transcription in cellular expression systems. Evidence from in vitro transfection studies is often assumed to be sufficient evidence for the in vivo functional importance of a polymorphism in the context of human disease, even though many confounding effects (e.g. temporal regulation, tissue specificity, genetic background) are not considered. In this study, we evaluate this assumption directly by examining the translation of in vitro results on allele-specific expression to an in vivo system using four genes that have been well documented through reporter assays to have promoter polymorphisms affecting transcription level: monoamine oxidase A (MAOA), neuropeptide Y (NPY), endothelial nitric oxide synthase (NOS3), and prodynorphin (PDYN). In our study, MAOA was found to have large allelic imbalances, which indicates that there is in vivo variation in the expression of this gene. However, the imbalances observed were not correlated with genotype at the putatively functional polymorphism. PDYN, NOS3 and NPY did not have large allelic imbalances. Overall, there was no statistically significant effect of these polymorphisms on expression level as measured by imbalance ratios in any of these genes. These results suggest that the functional effects of a polymorphism on gene expression may be more complicated and context dependent than is often assumed and also imply that the use of cell-based expression studies to support the role of such polymorphisms in disease etiology should be treated with caution.

Recombinational micro-evolution of functionally different metallothionein promoter alleles from Orchesella cincta

Background

Metallothionein (mt) transcription is elevated in heavy metal tolerant field populations of Orchesella cincta (Collembola). This suggests that natural selection acts on transcriptional regulation of mt in springtails at sites where cadmium (Cd) levels in soil reach toxic values This study investigates the nature and the evolutionary origin of polymorphisms in the metallothionein promoter (pmt) and their functional significance for mt expression.

Results

We sequenced approximately 1600 bp upstream the mt coding region by genome walking. Nine pmt alleles were discovered in NW-European populations. They differ in the number of some indels, consensus transcription factor binding sites and core promoter elements. Extensive recombination events between some of the alleles can be inferred from the alignment. A deviation from neutral expectations was detected in a cadmium tolerant population, pointing towards balancing selection on some promoter stretches. Luciferase constructs were made from the most abundant alleles, and responses to Cd, paraquat (oxidative stress inducer) and moulting hormone were studied in cell lines. By using paraquat we were able to dissect the effect of oxidative stress from the Cd specific effect, and extensive differences in mt induction levels between these two stressors were observed.

Conclusion

The pmt alleles evolved by a number of recombination events, and exhibited differential inducibilities by Cd, paraquat and molting hormone. In a tolerant population from a metal contaminated site, promoter allele frequencies differed significantly from a reference site and nucleotide polymorphisms in some promoter stretches deviated from neutral expectations, revealing a signature of balancing selection. Our results suggest that the structural differences in the Orchesella cincta metallothionein promoter alleles contribute to the metallothionein -over-expresser phenotype in cadmium tolerant populations.

Diverse selective modes among orthologs/paralogs of the chalcone synthase (Chs) gene family of Arabidopsis thaliana and its relative A. halleri ssp. gemmifera

As a model system, Arabidopsis thaliana and its wild relatives have played an important role in the study of genomics and evolution in plants. In this study, we examined the genetic diversity of the chalcone synthase (Chs) gene, which encodes a key enzyme of the flavonoid pathway and is located on chromosome five, as well as two Chs-like genes on the first and fourth chromosomes of Arabidopsis. The objectives of the study are to determine if natural selection operates differentially on the paralogs of the Chs gene family in A. thaliana and Arabidopsis halleri ssp. gemmifera. The mode of selection was inferred from Tajima's D values from noncoding and coding regions, as well as from the ratio of nonsynonymous to synonymous substitutions. Both McDonald-Kreitman and HKA tests revealed the effects of selection on the allelic distribution, except for the chromosome 1 paralog in ssp. gemmifera. The Chs gene on chromosome 5 was under purifying selection in both species. Significant, negative Tajima's D values at synonymous sites and positive Fay and Wu's H values within coding region, plus reduced genetic variability in introns, indicated effects of background selection in shaping the evolution of this gene region in A. thaliana. The Chs paralog on chromosome 1 was under positive selection in A. thaliana, while interspecific introgression and balancing selection determined the fates of the paralog and resulted in high heterogeneity in ssp. gemmifera. Local adaptation differentiated populations of Japan and China at the locus. In contrast, the other Chs-paralog of chromosome 4 was shaped by purifying selection in A. thaliana, while under positive selection in ssp. gemmifera, as indicated by dn/ds>1. Moreover, these contrasting patterns of selection have likely resulted in functional divergence in Arabidopsis, as indicated by radical amino acid substitutions at the chalcone synthase/stilbene synthase motif of the Chs genes. Unlike previous studies of the evolutionary history of A. thaliana, the high levels of genetic diversity in most gene regions of Chs paralogs and nonsignificant Tajima's D in the intron sequences of the Chs gene family in A. thaliana did not reflect the effects of a recent demographic expansion.

Tissue-specific activity of the blind mole rat and the two nucleotide-mutated mouse alphaB-crystallin promoter in transgenic mice

The alphaB-crystallin and HspB2 genes are located approximately 0.9 kb apart in a head-to-head arrangement in mammals. Previous experiments have shown that a truncated -668/+45 alphaB-crystallin enhancer/promoter fragment from blind mole rats (Spalax ehrenbergi), which have nonfunctional lenses, lacks lens activity and has enhanced muscle activity in transgenic mice. Here we show that the full-length mole rat alphaB-crystallin intergenic region behaves similarly in transgenic mice. A two-nucleotide mutation ((-273)CA-->G) in the mouse alphaB-crystallin enhancer/promoter fragment mimicking the wild-type mole rat sequence functionally converted the mouse promoter fragment to that of the wild-type mole rat promoter when tested in transgenic mice. The reciprocal mutation in the mole rat promoter fragment ((-272)G-->CA) did not affect its activity. Oligonucleotides from the wild-type mouse and mole rat alphaB-crystallin promoter region under study formed distinct complexes with nuclear proteins from cultured cells. The mouse mutant sequence lost binding ability, whereas the mutated mole rat sequence gained the ability to form a complex similar in size to that of the wild-type mouse oligonucleotide. Our data support the idea that blind mole rats' alphaB-crystallin promoter activity was modified during the evolution of subterranean life and shows that tissue-specific promoter activity can be modulated by changing as few as two apparently neutral nucleotides in the mouse alphaB-crystallin enhancer region, implying the importance of the context of regulatory sequences for promoter activity.

Cis-regulatory evolution of chalcone-synthase expression in the genus Arabidopsis

The contribution of cis-regulation to adaptive evolutionary change is believed to be essential, yet little is known about the evolutionary rules that govern regulatory sequences. Here, we characterize the short-term evolutionary dynamics of a cis-regulatory region within and among two closely related species, A. lyrata and A. halleri, and compare our findings to A. thaliana. We focused on the cis-regulatory region of chalcone synthase (CHS), a key enzyme involved in the synthesis of plant secondary metabolites. We observed patterns of nucleotide diversity that differ among species but do not depart from neutral expectations. Using intra- and interspecific F1 progeny, we have evaluated functional cis-regulatory variation in response to light and herbivory, environmental cues, which are known to induce CHS expression. We find that substantial cis-regulatory variation segregates within and among populations as well as between species, some of which results from interspecific genetic introgression. We further demonstrate that, in A. thaliana, CHS cis-regulation in response to herbivory is greater than in A. lyrata or A. halleri. Our work indicates that the evolutionary dynamics of a cis-regulatory region is characterized by pervasive functional variation, achieved mostly by modification of response modules to one but not all environmental cues. Our study did not detect the footprint of selection on this variation.

Cloning of DOG1, a quantitative trait locus controlling seed dormancy in Arabidopsis

Genetic variation for seed dormancy in nature is a typical quantitative trait controlled by multiple loci on which environmental factors have a strong effect. Finding the genes underlying dormancy quantitative trait loci is a major scientific challenge, which also has relevance for agriculture and ecology. In this study we describe the identification of the DELAY OF GERMINATION 1 (DOG1) gene previously identified as a quantitative trait locus involved in the control of seed dormancy. This gene was isolated by a combination of positional cloning and mutant analysis and is absolutely required for the induction of seed dormancy. DOG1 is a member of a small gene family of unknown molecular function, with five members in Arabidopsis. The functional natural allelic variation present in Arabidopsis is caused by polymorphisms in the cis-regulatory region of the DOG1 gene and results in considerable expression differences between the DOG1 alleles of the accessions analyzed.

Genetic mechanisms and evolutionary significance of natural variation in Arabidopsis

Genomic studies of natural variation in model organisms provide a bridge between molecular analyses of gene function and evolutionary investigations of adaptation and natural selection. In the model plant species Arabidopsis thaliana, recent studies of natural variation have led to the identification of genes underlying ecologically important complex traits, and provided new insights about the processes of genome evolution, geographic population structure, and the selective mechanisms shaping complex trait variation in natural populations. These advances illustrate the potential for a new synthesis to elucidate mechanisms for the adaptive evolution of complex traits from nucleotide sequences to real-world environments.

Evidence for the monophyletic evolution of benzylisoquinoline alkaloid biosynthesis in angiosperms

Benzylisoquinoline alkaloids (BIAs) consist of more than 2500 diverse structures largely restricted to the order Ranunculales and the eumagnoliids. However, BIAs also occur in the Rutaceae, Lauraceae, Cornaceae and Nelumbonaceae, and sporadically throughout the order Piperales. Several of these alkaloids function in the defense of plants against herbivores and pathogens - thus, the capacity for BIA biosynthesis is expected to play an important role in the reproductive fitness of certain plants. Biochemical and molecular phylogenetic approaches were used to investigate the evolution of BIA biosynthesis in basal angiosperms. The occurrence of (S)-norcoclaurine synthase (NCS; EC 4.2.1.78) activity in 90 diverse plant species was compared to the distribution of BIAs superimposed onto a molecular phylogeny. These results support the monophyletic origin of BIA biosynthesis prior to the emergence of the eudicots. Phylogenetic analyses of NCS, berberine bridge enzyme and several O-methyltransferases suggest a latent molecular fingerprint for BIA biosynthesis in angiosperms not known to accumulate such alkaloids. The limited occurrence of BIAs outside the Ranunculales and eumagnoliids suggests the requirement for a highly specialized, yet evolutionarily unstable cellular platform to accommodate or reactivate the pathway in divergent taxa. The molecular cloning and functional characterization of NCS from opium poppy (Papaver somniferum L.) is also reported. Pathogenesis-related (PR)10 and Bet v 1 major allergen proteins share homology with NCS, but recombinant polypeptides were devoid of NCS activity.