The rph1 Gene Is a Common Contributor to the Evolution of Phosphine Resistance in Independent Field Isolates of Rhyzopertha Dominica

Phosphine is the only economically viable fumigant for routine control of insect pests of stored food products, but its continued use is now threatened by the world-wide emergence of high-level resistance in key pest species. Phosphine has a unique mode of action relative to well-characterised contact pesticides. Similarly, the selective pressures that lead to resistance against field sprays differ dramatically from those encountered during fumigation. The consequences of these differences have not been investigated adequately. We determine the genetic basis of phosphine resistance in Rhyzopertha dominica strains collected from New South Wales and South Australia and compare this with resistance in a previously characterised strain from Queensland. The resistance levels range from 225 and 100 times the baseline response of a sensitive reference strain. Moreover, molecular and phenotypic data indicate that high-level resistance was derived independently in each of the three widely separated geographical regions. Despite the independent origins, resistance was due to two interacting genes in each instance. Furthermore, complementation analysis reveals that all three strains contain an incompletely recessive resistance allele of the autosomal rph1 resistance gene. This is particularly noteworthy as a resistance allele at rph1 was previously proposed to be a necessary first step in the evolution of high-level resistance. Despite the capacity of phosphine to disrupt a wide range of enzymes and biological processes, it is remarkable that the initial step in the selection of resistance is so similar in isolated outbreaks.

Genetic diversity and microcystin production by Anabaena in the Gulf of Finland, Baltic Sea

Syanobakteerit (sinilevät) ovat olleet Itämeressä koko nykymuotoisen Itämeren ajan, sillä paleolimnologiset todisteet niiden olemassaolosta Itämeren alueella ovat noin 7000 vuoden takaa. Syanobakteerien massaesiintymät eli kukinnat ovat kuitenkin sekä levinneet laajemmille alueille että tulleet voimakkaimmiksi viimeisten vuosikymmenien aikana. Tähän on osasyynä ihmisten aiheuttama kuormitus, joka rehevöittää Itämerta. Suomenlahti, jota tämä tutkimus käsittelee, on kärsinyt tästä rehevöitymiskehityksestä muita Itämeren altaita enemmän. Syanobakteerit muodostavat jokakesäisiä kukintoja Suomenlahdella - niin sen avomerialueilla kuin rannoillakin. Yleisimmät kukintoja muodostavat syanobakteerisuvut ovat Nodularia, Anabaena ja Aphanizomenon. Kukinnat aiheuttavat paitsi esteettistä haittaa myös terveydellisen riskitekijän. Niiden myrkyllisyys liitetään usein Nodularia-suvun tuottamaan nodulariini-maksamyrkkyyn. Itämeren Aphanizomenon-suvun on todettu olevan myrkytön. Vaikka Itämeren kukintoja aiheuttavista Nodularia- ja Aphanizomenon-syanobakteereista tiedetään varsin paljon, on molekyylimenetelmiin pohjautuva syanobakteeritutkimus ohittanut Itämeren Anabaena-suvun monelta osin. Tämän työn tarkoituksena oli syventää käsitystämme Itämeren Anabaena-syanobakteerista, sen mahdollisesta myrkyllisyydestä, geneettisestä monimuotoisuudesta ja fylogeneettisista sukulaisuussuhteista. Tässä työssä eristettiin 49 planktista Anabaena-kantaa, joista viisi tuottivat mikrokystiinejä. Tämä oli ensimmäinen yksiselitteinen todiste, että Itämeren Anabaena tuottaa maksamyrkyllisiä mikrokystiini-yhdisteitä. Jokainen eristetty myrkyllinen Anabaena-kanta tuotti useita mikrokystiini-variantteja. Lisäksi mikrokystiinejä löydettiin kukintanäytteistä, joissa oli myrkkyä syntetisoivia geenejä sisältäneitä Anabaena-syanobakteereita. Myrkkyjä löydettiin molempina tutkimusvuosina 2003 ja 2004. Myrkkyjen esiintyminen ei siten ollut vain yksittäinen ilmiö. Tässä työssä saimme viitteitä siitä, että maksamyrkyllinen Anabaena-syanobakteeri esiintyisi vähäsuolaisissa vesissä. Tämä riippuvuussuhde jää kuitenkin tulevien tutkimuksien selvitettäväksi. Tässä työssä havaittiin mikrokystiinisyntetaasi-geenien inaktivoituminen Itämeren Anabaena-kannassa ja kukintanäytteissä. Kuvasimme Anabaena-kannan mikrokystiinisyntetaasigeenien sisältä insertioita, jotka hyvin todennäköisesti inaktivoivat myrkyntuoton. Insertion sisältäneeltä kannalta löysimme kuitenkin kaikki mikrokystiinisyntetaasigeenit osoittaen, että geenien olemassaolo ei välttämättä varmista kannan mikrokystiinintuottoa. Mielenkiintoista oli se, että inaktivaation aiheuttavia insertioita löytyi kukintanäytteistä molemmilta tutkimusvuosilta. Vastaavia insertioita ei kuitenkaan löydetty makean veden Anabaena-kannoista tai järvinäytteistä. On yleistä, että syanobakteerikukinnoista löytyy usean syanobakteerisuvun edustajia. Myrkyllisiä sukuja tai lajeja ei voida kuitenkaan erottaa mikroskooppisesti myrkyttömistä. Käsillä olevassa tutkimuksessa kehitettiin molekyylimenetelmä, jolla on mahdollista määrittää kukinnan mahdollisesti maksamyrkylliset syanobakteerisuvut. Tätä menetelmää sovellettiin Itämeren kukintojen tutkimiseen. Itämeren pintavesistä ja ranta-alueiden pohjasta eristetyt Anabaena-kannat osoittautuivat geneettisesti monimuotoisiksi. Tämä Anabaena-syanobakteerien geneettinen monimuotoisuus vahvistettiin monistamalla geenejä suoraan kukintanäytteistä ilman kantojen eristystä. Makeiden vesien ja Itämeren Anabaena-kannat ovat geneettisesti hyvin samankaltaisia. Geneettisissä vertailuissa kävi kuitenkin ilmi, että pohjassa elävien Anabaena-kantojen geneettinen monimuotoisuus oli suurempaa kuin pintavesistä eristettyjen kantojen. Itämeren Anabaena-kantojen sekvenssit muodostivat omia ryhmiä sukupuun sisällä, jolloin on mahdollista, että nämä edustavat Itämeren omia Anabaena-ekotyyppejä. Tämä tutkimus oli ensimmäinen, jossa uusin molekyylimenetelmin systemaattisesti selvitettiin Itämeren Anabaena-syanobakteerin geneettistä populaatiorakennetta, fylogeniaa ja myrkyntuottoa. Tulevaisuudessa monitorointitutkimuksissa on otettava huomioon myös Itämeren Anabaena-syanobakteerin mahdollinen maksamyrkyntuotto – erityisesti vähäsuolaisemmilla rannikkovesillä.

Pleistocene isolation, secondary introgression and restricted contemporary gene flow in the pig-eye shark, Carcharhinus amboinensis across northern Australia

We examine the structure and phylogeography of the pig-eye shark (Carcharhinus amboinensis) common in shallow coastal environments in northern Australia using two types of genetic markers, two mitochondrial (control region and NADH hydrogenase 4) and two nuclear (microsatellite and Rag 1) DNA. Two populations were defined within northern Australia on the basis of mitochondrial DNA evidence, but this result was not supported by nuclear microsatellite or Rag 1 markers. One possibility for this structure might be sex-specific behaviours such as female philopatry, although we argue it is doubtful that sufficient time has elapsed for any potential signatures from this behaviour to be expressed in nuclear markers. It is more likely that the observed pattern represents ancient populations repeatedly isolated and connected during episodic sea level changes during the Pleistocene epoch, until current day with restricted contemporary gene flow maintaining population genetic structure. Our results show the need for an understanding of both the history and ecology of a species in order to interpret patterns in genetic structure.

Associations of gene polymorphisms and nutrition with calcium homeostasis and bone mineral density : Studies on skeletal nutrigenetics

Heredity explains a major part of the variation in calcium homeostasis and bone strength, and the susceptibility to osteoporosis is polygenetically regulated. Bone phenotype results from the interplay between lifestyle and genes, and several nutritional factors modulate bone health throughout life. Thus, nutrigenetics examining the genetic variation in nutrient intake and homeostatic control is an important research area in the etiology of osteoporosis. Despite continuing progress in the search for candidate genes for osteoporosis, the results thus far have been inconclusive. The main objective of this thesis was to investigate the associations of lactase, vitamin D receptor (VDR), calcium sensing receptor (CaSR) and parathyroid hormone (PTH) gene polymorphisms and lifestyle factors and their interactions with bone health in Finns at varying stages of the skeletal life span. Markers of calcium homeostasis and bone remodelling were measured from blood and urine samples. Bone strength was measured at peripheral and central bone sites. Lifestyle factors were assessed with questionnaires and interviews. Genetic lactase non-persistence (the C/C-13910 genotype) was associated with lower consumption of milk from childhood, predisposing females in particular to inadequate calcium intake. Consumption of low-lactose milk and milk products was shown to decrease the risk for inadequate calcium intake. In young adulthood, bone loss was more common in males than in females. Males with the lactase C/C-13910 genotype may be more susceptible to bone loss than males with the other lactase genotypes, although calcium intake predicts changes in bone mass more than the lactase genotype. The BsmI and FokI polymorphisms of the VDR gene were associated with bone mass in growing adolescents, but the associations weakened with age. In young adults, the A986S polymorphism of the calcium sensing receptor gene was associated with serum ionized calcium concentrations, and the BstBI polymorphism of the parathyroid gene was related to bone strength. The FokI polymorphism and sodium intake showed an interaction effect on urinary calcium excretion. A novel gene-gene interaction between the VDR FokI and PTH BstBI gene polymorphisms was found in the regulation of PTH secretion and urinary calcium excretion. Further research should be carried out with more number of Finns at varying stages of the skeletal life span and more detailed measurements of bone strength. Research should concern mechanisms by which genetic variants affect calcium homeostasis and bone strength, and the role of diet-gene and gene-gene interactions in the pathogenesis of osteoporosis.

L-rhamnose-1-dehydrogenase gene and L-rhamnose catabolism in the yeast Pichia stipitis

The purpose of this work was to identify some of the genes of the catabolic route of L-rhamnose in the yeast Pichia stipitis. There are at least two distinctly different pathways for L-rhamnose catabolism. The one described in bacteria has phosphorylated intermediates and the enzymes and the genes of this route have been described. The pathway described in yeast does not have phosphorylated intermediates. The intermediates and the enzymes of this pathway are known but none of the genes have been identified. The work was started by purifying the L-rhamnose dehydrogenase, which oxidates L-rhamnose to rhamnonic acid-gamma-lactone. NAD is used as a cofactor in this reaction. A DEAE ion exchange column was used for purification. The active fraction was further purified using a non-denaturing PAGE and the active protein identified by zymogram staining. In the last step the protein was separated in a SDS-PAGE, the protein band trypsinated and analysed by MALDI-TOF MS. This resulted in the identification of the corresponding gene, RHA1, which was then, after a codon change, expressed in Saccharomyces cerevisiae. Also C- or N-terminal histidine tags were added but as the activity of the enzyme was lost or strongly reduced these were not used. The kinetic properties of the protein were analysed in the cell extract. Substrate specifity was tested with different sugars; L-rhamnose, L-lyxose and L-mannose were oxidated by the enzyme. Vmax values were 180 nkat/mg, 160 nkat/mg and 72 nkat/mg, respectively. The highest affinity was towards L-rhamnose, the Km value being 0.9 mM. Lower affinities were obtained with L-lyxose, Km 4.3 mM, and L-mannose Km 25 mM. Northern analysis was done to study the transcription of RHA1 with different carbon sources. Transcription was observed only on L-rhamnose suggesting that RHA1 expression is L-rhamnose induced. A RHA1 deletion cassette for P. stipitis was constructed but the cassette had integrated randomly and not targeted to delete the RHA1 gene. Enzyme assays for L-lactaldehyde dehydrogenase were done similarly to L-rhamnose dehydrogenase assays. NAD is used as a cofactor also in this reaction where L-lactaldehyde is oxidised to L-lactate. The observed enzyme activities were very low and the activity was lost during the purification procedures.

Sequence variation in the putative effector gene SIX8 facilitates molecular differentiation of Fusarium oxysporum f. sp. cubense

Fusarium oxysporum f. sp. cubense (Foc), causal agent of fusarium wilt of banana, is among the most destructive pathogens of banana and plantain. The development of a molecular diagnostic capable of reliably distinguishing between the various races of the pathogen is of key importance to disease management. However, attempts to distinguish isolates using the standard molecular loci typically used for fungal phylogenetics have been complicated by a poor correlation between phylogeny and pathogenicity. Among the available alternative loci are several putative effector genes, known as SIX genes, which have been successfully used to differentiate the three races of F. oxysporum f. sp. lycopersici. In this study, an international collection of Foc isolates was screened for the presence of the putative effector SIX8. Using a PCR and sequencing approach, variation in Foc-SIX8 was identified which allowed race 4 to be differentiated from race 1 and 2 isolates, and tropical and subtropical race 4 isolates to be distinguished from one another.

Fungal community structure in disease suppressive soils assessed by 28S LSU gene sequencing

Natural biological suppression of soil-borne diseases is a function of the activity and composition of soil microbial communities. Soil microbe and phytopathogen interactions can occur prior to crop sowing and/or in the rhizosphere, subsequently influencing both plant growth and productivity. Research on suppressive microbial communities has concentrated on bacteria although fungi can also influence soil-borne disease. Fungi were analyzed in co-located soils 'suppressive' or 'non-suppressive' for disease caused by Rhizoctonia solani AG 8 at two sites in South Australia using 454 pyrosequencing targeting the fungal 28S LSU rRNA gene. DNA was extracted from a minimum of 125 g of soil per replicate to reduce the micro-scale community variability, and from soil samples taken at sowing and from the rhizosphere at 7 weeks to cover the peak Rhizoctonia infection period. A total of ∼994,000 reads were classified into 917 genera covering 54% of the RDP Fungal Classifier database, a high diversity for an alkaline, low organic matter soil. Statistical analyses and community ordinations revealed significant differences in fungal community composition between suppressive and non-suppressive soil and between soil type/location. The majority of differences associated with suppressive soils were attributed to less than 40 genera including a number of endophytic species with plant pathogen suppression potentials and mycoparasites such as Xylaria spp. Non-suppressive soils were dominated by Alternaria , Gibberella and Penicillum. Pyrosequencing generated a detailed description of fungal community structure and identified candidate taxa that may influence pathogen-plant interactions in stable disease suppression. © 2014 Penton et al.

Resistance gene analogues in sugarcane and sorghum and their association with quantitative trait loci for rust resistance

Fifty-four different sugarcane resistance gene analogue (RGA) sequences were isolated, characterized, and used to identify molecular markers linked to major disease-resistance loci in sugarcane. Ten RGAs were identified from a sugarcane stem expressed sequence tag (EST) library; the remaining 44 were isolated from sugarcane stem, leaf, and root tissue using primers designed to conserved RGA motifs. The map location of 31 of the RGAs was determined in sugarcane and compared with the location of quantitative trait loci (QTL) for brown rust resistance. After 2 years of phenotyping, 3 RGAs were shown to generate markers that were significantly associated with resistance to this disease. To assist in the understanding of the complex genetic structure of sugarcane, 17 of the 31 RGAs were also mapped in sorghum. Comparative mapping between sugarcane and sorghum revealed syntenic localization of several RGA clusters. The 3 brown rust associated RGAs were shown to map to the same linkage group (LG) in sorghum with 2 mapping to one region and the third to a region previously shown to contain a major rust-resistance QTL in sorghum. These results illustrate the value of using RGAs for the identification of markers linked to disease resistance loci and the value of simultaneous mapping in sugarcane and sorghum.

Methods to improve gene signal: Application to cDNA microarrays

Microarrays are high throughput biological assays that allow the screening of thousands of genes for their expression. The main idea behind microarrays is to compute for each gene a unique signal that is directly proportional to the quantity of mRNA that was hybridized on the chip. A large number of steps and errors associated with each step make the generated expression signal noisy. As a result, microarray data need to be carefully pre-processed before their analysis can be assumed to lead to reliable and biologically relevant conclusions. This thesis focuses on developing methods for improving gene signal and further utilizing this improved signal for higher level analysis. To achieve this, first, approaches for designing microarray experiments using various optimality criteria, considering both biological and technical replicates, are described. A carefully designed experiment leads to signal with low noise, as the effect of unwanted variations is minimized and the precision of the estimates of the parameters of interest are maximized. Second, a system for improving the gene signal by using three scans at varying scanner sensitivities is developed. A novel Bayesian latent intensity model is then applied on these three sets of expression values, corresponding to the three scans, to estimate the suitably calibrated true signal of genes. Third, a novel image segmentation approach that segregates the fluorescent signal from the undesired noise is developed using an additional dye, SYBR green RNA II. This technique helped in identifying signal only with respect to the hybridized DNA, and signal corresponding to dust, scratch, spilling of dye, and other noises, are avoided. Fourth, an integrated statistical model is developed, where signal correction, systematic array effects, dye effects, and differential expression, are modelled jointly as opposed to a sequential application of several methods of analysis. The methods described in here have been tested only for cDNA microarrays, but can also, with some modifications, be applied to other high-throughput technologies. Keywords: High-throughput technology, microarray, cDNA, multiple scans, Bayesian hierarchical models, image analysis, experimental design, MCMC, WinBUGS.

Gene flow from transgenic plant populations: models and applications for risk assessment

The future use of genetically modified (GM) plants in food, feed and biomass production requires a careful consideration of possible risks related to the unintended spread of trangenes into new habitats. This may occur via introgression of the transgene to conventional genotypes, due to cross-pollination, and via the invasion of GM plants to new habitats. Assessment of possible environmental impacts of GM plants requires estimation of the level of gene flow from a GM population. Furthermore, management measures for reducing gene flow from GM populations are needed in order to prevent possible unwanted effects of transgenes on ecosystems. This work develops modeling tools for estimating gene flow from GM plant populations in boreal environments and for investigating the mechanisms of the gene flow process. To describe spatial dimensions of the gene flow, dispersal models are developed for the local and regional scale spread of pollen grains and seeds, with special emphasis on wind dispersal. This study provides tools for describing cross-pollination between GM and conventional populations and for estimating the levels of transgenic contamination of the conventional crops. For perennial populations, a modeling framework describing the dynamics of plants and genotypes is developed, in order to estimate the gene flow process over a sequence of years. The dispersal of airborne pollen and seeds cannot be easily controlled, and small amounts of these particles are likely to disperse over long distances. Wind dispersal processes are highly stochastic due to variation in atmospheric conditions, so that there may be considerable variation between individual dispersal patterns. This, in turn, is reflected to the large amount of variation in annual levels of cross-pollination between GM and conventional populations. Even though land-use practices have effects on the average levels of cross-pollination between GM and conventional fields, the level of transgenic contamination of a conventional crop remains highly stochastic. The demographic effects of a transgene have impacts on the establishment of trangenic plants amongst conventional genotypes of the same species. If the transgene gives a plant a considerable fitness advantage in comparison to conventional genotypes, the spread of transgenes to conventional population can be strongly increased. In such cases, dominance of the transgene considerably increases gene flow from GM to conventional populations, due to the enhanced fitness of heterozygous hybrids. The fitness of GM plants in conventional populations can be reduced by linking the selectively favoured primary transgene to a disfavoured mitigation transgene. Recombination between these transgenes is a major risk related to this technique, especially because it tends to take place amongst the conventional genotypes and thus promotes the establishment of invasive transgenic plants in conventional populations.

Inactivation of Aspergillus flavus spores by curcumin-mediated photosensitization

Minimizing fungal infection is essential to the control of mycotoxin contamination of foods and feeds but many potential control methods are not without their own safety concerns for the consumers. Photodynamic inactivation is a novel light-based approach which offers a promising alternative to conventional methods for the control of mycotoxigenic fungi. This study describes the use of curcumin to inactivate spores of Aspergillus flavus, one of the major aflatoxin producing fungi in foods and feeds. Curcumin is a natural polyphenolic compound from the spice turmeric (Curcuma longa). In this study the plant has shown to be an effective photosensitiser when combined with visible light (420 nm). The experiment was conducted in in vitro and in vivo where A. flavus spores were treated with different photosensitiser concentration and light dose both in buffer solution and on maize kernels. Comparison of fungal load from treated and untreated samples was determined, and reductions of fungal spore counts of up to 3 log CFU ml−1 in suspension and 2 log CFU g−1 in maize kernels were obtained using optimal dye concentrations and light dose combinations. The results in this study indicate that curcumin-mediated photosensitization is a potentially effective method to decontaminate A. flavus spores in foods and feeds.

Water stress induced alterations in ornithine aminotransferase of ragi (Eleusine coracana): Protection by proline against heat inactivation and denaturation by urea and guanidinium chloride

Water stress resulted in a specific response leading to a large and significant increase (80-fold) in free proline content of ragi (Eleusine coracana) leaves and seedlings. L-Proline protected ornithine aminotransferase, an enzyme in the pathway for proline biosynthesis, isolated from normal and stressed ragi leaves against heat inactivation and denaturation by urea and guanidinium chloride. The protection of the stressed enzyme by L-proline was much more complete than that of the enzyme isolated from normal leaves. While L-ornithine, one of the substrates, protected the stressed enzyme against inactivation, it enhanced the rate of inactivation of the normal enzyme. α-Ketoglutarate protected both the normal and stressed enzyme against inactivation and denaturation. These results support the suggestion that ornithine aminotransferase has undergone a structural alteration during water stress. In view of the causal relationship between elevated temperature and water stress of plants under natural conditions, the protection afforded by proline against inactivation and denaturation of the enzyme from stressed leaves assumes significance. These results provide an explanation for a possible functional importance of proline accumulation during water stress.

Senataxin controls meiotic silencing through ATR activation and chromatin remodeling

Senataxin, defective in ataxia oculomotor apraxia type 2, protects the genome by facilitating the resolution of RNA–DNA hybrids (R-loops) and other aspects of RNA processing. Disruption of this gene in mice causes failure of meiotic recombination and defective meiotic sex chromosome inactivation, leading to male infertility. Here we provide evidence that the disruption of Setx leads to reduced SUMOylation and disruption of protein localization across the XY body during meiosis. We demonstrate that senataxin and other DNA damage repair proteins, including ataxia telangiectasia and Rad3-related protein-interacting partner, are SUMOylated, and a marked downregulation of both ataxia telangiectasia and Rad3-related protein-interacting partner and TopBP1 leading to defective activation and signaling through ataxia telangiectasia and Rad3-related protein occurs in the absence of senataxin. Furthermore, chromodomain helicase DNA-binding protein 4, a component of the nucleosome remodeling and deacetylase chromatin remodeler that interacts with both ataxia telangiectasia and Rad3-related protein and senataxin was not recruited efficiently to the XY body, triggering altered histone acetylation and chromatin conformation in Setx−/− pachytene-staged spermatocytes. These results demonstrate that senataxin has a critical role in ataxia telangiectasia and Rad3-related protein- and chromodomain helicase DNA-binding protein 4-mediated transcriptional silencing and chromatin remodeling during meiosis providing greater insight into its critical role in gene regulation to protect against neurodegeneration.

Construction of a global map of human gene expression : the process, tools and analysis

This thesis studies human gene expression space using high throughput gene expression data from DNA microarrays. In molecular biology, high throughput techniques allow numerical measurements of expression of tens of thousands of genes simultaneously. In a single study, this data is traditionally obtained from a limited number of sample types with a small number of replicates. For organism-wide analysis, this data has been largely unavailable and the global structure of human transcriptome has remained unknown. This thesis introduces a human transcriptome map of different biological entities and analysis of its general structure. The map is constructed from gene expression data from the two largest public microarray data repositories, GEO and ArrayExpress. The creation of this map contributed to the development of ArrayExpress by identifying and retrofitting the previously unusable and missing data and by improving the access to its data. It also contributed to creation of several new tools for microarray data manipulation and establishment of data exchange between GEO and ArrayExpress. The data integration for the global map required creation of a new large ontology of human cell types, disease states, organism parts and cell lines. The ontology was used in a new text mining and decision tree based method for automatic conversion of human readable free text microarray data annotations into categorised format. The data comparability and minimisation of the systematic measurement errors that are characteristic to each lab- oratory in this large cross-laboratories integrated dataset, was ensured by computation of a range of microarray data quality metrics and exclusion of incomparable data. The structure of a global map of human gene expression was then explored by principal component analysis and hierarchical clustering using heuristics and help from another purpose built sample ontology. A preface and motivation to the construction and analysis of a global map of human gene expression is given by analysis of two microarray datasets of human malignant melanoma. The analysis of these sets incorporate indirect comparison of statistical methods for finding differentially expressed genes and point to the need to study gene expression on a global level.