Variation In Transcript Abundance Research Articles

Gene Expression Reaction Norms Unravel the Molecular and Cellular Processes Underpinning the Plastic Phenotypes of Alternanthera Philoxeroides in Contrasting Hydrological Conditions.

Alternanthera philoxeroides is an amphibious invasive weed that can colonize both aquatic and terrestrial habitats. Individuals growing in different habitats exhibit extensive phenotypic variation but little genetic differentiation. Little is known about the molecular basis underlying environment-induced phenotypic changes. Variation in transcript abundance in A. philoxeroides was characterized throughout the time-courses of pond and upland treatments using RNA-Sequencing. Seven thousand eight hundred and five genes demonstrated variable expression in response to different treatments, forming 11 transcriptionally coordinated gene groups. Functional enrichment analysis of plastically expressed genes revealed pathway changes in hormone-mediated signaling, osmotic adjustment, cell wall remodeling, and programmed cell death, providing a mechanistic understanding of the biological processes underlying the phenotypic changes in A. philoxeroides. Both transcriptional modulation of environmentally sensitive loci and environmentally dependent control of regulatory loci influenced the plastic responses to the environment. Phenotypic responses and gene expression patterns to contrasting hydrological conditions were compared between A. philoxeroides and its alien congener Alternanthera pungens. The terricolous A. pungens displayed limited phenotypic plasticity to different treatments. It was postulated based on gene expression comparison that the interspecific variation in plasticity between A. philoxeroides and A. pungens was not due to environmentally-mediated changes in hormone levels but to variations in the type and relative abundance of different signal transducers and receptors expressed in the target tissue.

Isolation and expression analysis of 18 CsbZIP genes implicated in abiotic stress responses in the tea plant (Camellia sinensis)

Basic leucine zipper (bZIP) transcription factors (TFs) play essential roles in regulating stress processes in plants. Despite the economic importance of this woody crop, there is little information about bZIP TFs in tea plants. In this study, 18 bZIP genes were isolated from the tea plant (Camellia sinensis) and named sequentially from CsbZIP1 to CsbZIP18. According to the phylogenetic classification as in Arabidopsis, the CsbZIP genes spanned ten subgroups (Group A, B, C, D, E, F, H, I, S and K) of bZIP TFs. When analyzed for organ specific expression, all CsbZIP genes were found to be ubiquitously expressed in roots, stems, leaves and flowers. Expression analysis of CsbZIP genes in response to four abiotic stresses showed that in leaves, 9, 9, 15 and 11 CsbZIPs have 2-fold greater variation in transcript abundance under cold, exogenous ABA, high salinity and dehydration conditions, respectively. In roots, 5, 12, 14 and 11 CsbZIPs were differentially expressed under conditions of cold, exogenous ABA, high salinity and dehydration stresses. Moreover, CsbZIP genes in Groups F, H, S and K exhibited several folds up-and/or down-regulation against the above four stresses. Notably, CsbZIP18 of group K showed significant up-regulation in response to these same stresses, suggesting a vital functional role in stress response. Together, these findings increase our knowledge of bZIP TFs in the tea plant and suggest the significance of CsbZIP genes in plant abiotic responses.

StructureFold: genome-wide RNA secondary structure mapping and reconstruction in vivo.

RNAs fold into complex structures that are integral to the diverse mechanisms underlying RNA regulation of gene expression. Recent development of transcriptome-wide RNA structure profiling through the application of structure-probing enzymes or chemicals combined with high-throughput sequencing has opened a new field that greatly expands the amount of in vitro and in vivo RNA structural information available. The resultant datasets provide the opportunity to investigate RNA structural information on a global scale. However, the analysis of high-throughput RNA structure profiling data requires considerable computational effort and expertise. We present a new platform, StructureFold, that provides an integrated computational solution designed specifically for large-scale RNA structure mapping and reconstruction across any transcriptome. StructureFold automates the processing and analysis of raw high-throughput RNA structure profiling data, allowing the seamless incorporation of wet-bench structural information from chemical probes and/or ribonucleases to restrain RNA secondary structure prediction via the RNAstructure and ViennaRNA package algorithms. StructureFold performs reads mapping and alignment, normalization and reactivity derivation, and RNA structure prediction in a single user-friendly web interface or via local installation. The variation in transcript abundance and length that prevails in living cells and consequently causes variation in the counts of structure-probing events between transcripts is accounted for. Accordingly, StructureFold is applicable to RNA structural profiling data obtained in vivo as well as to in vitro or in silico datasets. StructureFold is deployed via the Galaxy platform. StructureFold is freely available as a component of Galaxy available at: https://usegalaxy.org/. yxt148@psu.edu or sma3@psu.edu Supplementary data are available at Bioinformatics online.

Digital multiplexed mRNA analysis of functionally important genes in single human oocytes and correlation of changes in transcript levels with oocyte protein expression

ObjectiveTo investigate functionally important transcripts in single human oocytes with the use of NanoString technology and determine whether observed differences are biologically meaningful.DesignAnalysis of human oocytes with the use of NanoString and immunoblotting.SettingUniversity-affiliated reproductive medicine unit.PatientsWomen undergoing in vitro fertilization.InterventionHuman oocytes were analyzed with the use of NanoString or immunoblotting.Main Outcome MeasuresThe abundance of transcripts for ten functionally important genes—AURKA, AURKC, BUB1, BUB1B (encoding BubR1), CDK1, CHEK1, FYN, MOS, MAP2K1, and WEE2—and six functionally dispensable genes were analyzed with the use of NanoString. BubR1 protein levels in oocytes from younger and older women were compared with the use of immunoblotting.Result(s)All ten functional genes but none of the six dispensable genes were detectable with the use of NanoString in single oocytes. There was 3- to 5-fold variation in BUB1, BUB1B, and CDK1 transcript abundance among individual oocytes from a single patient. Transcripts for these three genes—all players within the spindle assembly checkpoint surveillance mechanism for preventing aneuploidy—were reduced in the same oocyte from an older patient. Mean BUB1B transcripts were reduced by 1.5-fold with aging and associated with marked reductions in BubR1 protein levels.Conclusion(s)The abundance of functionally important transcripts exhibit marked oocyte-to-oocyte heterogeneity to a degree that is sufficient to affect protein expression. Observed variations in transcript abundance are therefore likely to be biologically meaningful, especially if multiple genes within the same pathway are simultaneously affected.

Differential Expression of a Sodium-Phosphate Cotransporter Among Vibrio vulnificus Strains

Vibrio vulnificus is an estuarine bacterium with pathogenic potential. Its three known biotypes differ in host distribution. We have found the nptA gene for a sodium-phosphate cotransporter, which is rare in bacteria, in each biotype. nptA transcript abundance differed significantly among biotypes, leading to the hypothesis that transcript levels differ under environmental conditions associated with estuarine and host environments. nptA transcript abundance was assessed in V. vulnificus biotypes 1 (C and E genotypes), 2 and 3 strains under varied salinity, phosphate concentration, and pH. Differences in transcript abundance separated strains into two groups. Type C and biotype 3 strains formed Group 1, while type E and biotype 2 strains formed Group 2. Group 2 strains had significantly greater nptA RNA transcript abundance than Group 1. Transcript abundance in the two groups also responded differently to pH and salinity, suggesting differential regulation of nptA in response to environmental conditions. Comparison of the deduced amino acid sequences of NptA among strains resulted in strain grouping similar to that based on transcript abundance. Variation in transcript abundance between groups may affect the ability of V. vulnificus strains to colonize hosts and/or to compete as free-living bacteria in various habitats.

The flavonoid pathway in tomato seedlings: transcript abundance and the modeling of metabolite dynamics.

Flavonoids are secondary metabolites present in all terrestrial plants. The flavonoid pathway has been extensively studied, and many of the involved genes and metabolites have been described in the literature. Despite this extensive knowledge, the functioning of the pathway in vivo is still poorly understood. Here, we study the flavonoid pathway using both experiments and mathematical models. We measured flavonoid metabolite dynamics in two tissues, hypocotyls and cotyledons, during tomato seedling development. Interestingly, the same backbone of interactions leads to very different accumulation patterns in the different tissues. Initially, we developed a mathematical model with constant enzyme concentrations that described the metabolic networks separately in both tissues. This model was unable to fit the measured flavonoid dynamics in the hypocotyls, even if we allowed unrealistic parameter values. This suggested us to investigate the effect of transcript abundance on flavonoid accumulation. We found that the expression of candidate flavonoid genes varies considerably with time. Variation in transcript abundance results in enzymatic variation, which could have a large effect on metabolite accumulation. Candidate transcript abundance was included in the mathematical model as representative for enzyme concentration. We fitted the resulting model to the flavonoid dynamics in the cotyledons, and tested it by applying it to the data from hypocotyls. When transcript abundance is included, we are indeed able to explain flavonoid dynamics in both tissues. Importantly, this is possible under the biologically relevant restriction that the enzymatic properties estimated by the model are conserved between the tissues.

Mendelian and Non-Mendelian Regulation of Gene Expression in Maize

Transcriptome variation plays an important role in affecting the phenotype of an organism. However, an understanding of the underlying mechanisms regulating transcriptome variation in segregating populations is still largely unknown. We sought to assess and map variation in transcript abundance in maize shoot apices in the intermated B73×Mo17 recombinant inbred line population. RNA–based sequencing (RNA–seq) allowed for the detection and quantification of the transcript abundance derived from 28,603 genes. For a majority of these genes, the population mean, coefficient of variation, and segregation patterns could be predicted by the parental expression levels. Expression quantitative trait loci (eQTL) mapping identified 30,774 eQTL including 96 trans-eQTL “hotspots,” each of which regulates the expression of a large number of genes. Interestingly, genes regulated by a trans-eQTL hotspot tend to be enriched for a specific function or act in the same genetic pathway. Also, genomic structural variation appeared to contribute to cis-regulation of gene expression. Besides genes showing Mendelian inheritance in the RIL population, we also found genes whose expression level and variation in the progeny could not be predicted based on parental difference, indicating that non-Mendelian factors also contribute to expression variation. Specifically, we found 145 genes that show patterns of expression reminiscent of paramutation such that all the progeny had expression levels similar to one of the two parents. Furthermore, we identified another 210 genes that exhibited unexpected patterns of transcript presence/absence. Many of these genes are likely to be gene fragments resulting from transposition, and the presence/absence of their transcripts could influence expression levels of their ancestral syntenic genes. Overall, our results contribute to the identification of novel expression patterns and broaden the understanding of transcriptional variation in plants.

Assessing the Genome-Wide Effect of Promoter Region Tandem Repeat Natural Variation on Gene Expression

Copy number polymorphisms of nucleotide tandem repeat (TR) regions, such as microsatellites and minisatellites, are mutationally reversible and highly abundant in eukaryotic genomes. Studies linking TR polymorphism to phenotypic variation have led some to suggest that TR variation modulates and majorly contributes to phenotypic variation; however, studies in which the authors assess the genome-wide impact of TR variation on phenotype are lacking. To address this question, we quantified relationships between polymorphism levels in 143 genome-wide promoter region TRs across 16 isolates of the filamentous fungus Aspergillus flavus and its ecotype Aspergillus oryzae with expression levels of their downstream genes. We found that only 4.3% of relationships tested were significant; these findings were consistent with models in which TRs act as “tuning,” “volume,” or “optimality” “knobs” of phenotype but not with “switch” models. Furthermore, the promoter regions of differentially expressed genes between A. oryzae and A. flavus did not show TR enrichment, suggesting that genome-wide differences in molecular phenotype between the two species are not significantly associated with TRs. Although in some cases TR polymorphisms do contribute to transcript abundance variation, these results argue that at least in this case, TRs might not be major modulators of variation in phenotype.

The Rapid Evolution of X-linked Male-Biased Gene Expression and the Large-X Effect in Drosophila yakuba, D. santomea, and Their Hybrids

The X chromosome has a large effect on hybrid dysfunction, particularly on hybrid male sterility. Although the evidence for this so-called large-X effect is clear, its molecular causes are not yet fully understood. One possibility is that, under certain conditions, evolution proceeds faster in X-linked than in autosomal loci (i.e., faster-X effect) due to both natural selection and their hemizygosity in males, an effect that is expected to be greatest in genes with male-biased expression. Here, I study genome-wide variation in transcript abundance between Drosophila yakuba and D. santomea, within these species and in their hybrid males to evaluate both the faster-X and large-X effects at the level of expression. I find that in X-linked male-biased genes (MBGs) expression evolves faster than in their autosomal counterparts, an effect that is accompanied by a unique reduction in expression polymorphism. This suggests that Darwinian selection is driving expression differences between species, likely enhanced by the hemizygosity of the X chromosome in males. Despite the recent split of the two sister species under study, abundant changes in both cis- and trans-regulatory elements underlie expression divergence in the majority of the genes analyzed, with significant differences in allelic ratios of transcript abundance between the two reciprocal F(1) hybrid males. Cis-trans coevolution at molecular level, evolved shortly after populations become isolated, may therefore contribute to explain the breakdown of the regulation of gene expression in hybrid males. Additionally, the X chromosome plays a large role in this hybrid male misexpression, which affects not only MBG but also, to a lesser degree, nonsex-biased genes. Interestingly, hybrid male misexpression is concentrated mostly in autosomal genes, likely facilitated by the rapid evolution of sex-linked trans-acting factors. I suggest that the faster evolution of X-linked MBGs, at both protein and expression levels, contributes to explain the large effect of the X chromosome on hybrid male sterility, likely mediating widespread autosomal misexpression through the preferential recognition of cis-regulatory elements by conspecific trans-acting factors (i.e., cis-trans conspecific recognition).

Drought induces alterations in the stomatal development program in Populus.

Much is known about the physiological control of stomatal aperture as a means by which plants adjust to water availability. By contrast, the role played by the modulation of stomatal development to limit water loss has received much less attention. The control of stomatal development in response to water deprivation in the genus Populus is explored here. Drought induced declines in stomatal conductance as well as an alteration in stomatal development in two genotypes of Populus balsamifera. Leaves that developed under water-deficit conditions had lower stomatal indices than leaves that developed under well-watered conditions. Transcript abundance of genes that could hypothetically underpin drought-responsive changes in stomatal development was examined, in two genotypes, across six time points, under two conditions, well-watered and with water deficit. Populus homologues of STOMAGEN, ERECTA (ER), STOMATA DENSITY AND DISTRIBUTION 1 (SDD1), and FAMA had variable transcript abundance patterns congruent with their role in the modulation of stomatal development in response to drought. Conversely, there was no significant variation in transcript abundance between genotypes or treatments for the Populus homologues of YODA (YDA) and TOO MANY MOUTHS (TMM). The findings highlight the role that could be played by stomatal development during leaf expansion as a longer term means by which to limit water loss from leaves. Moreover, the results point to the key roles played by the regulation of the homologues of STOMAGEN, ER, SDD1, and FAMA in the control of this response in poplar.

Nuclear RNase P of Trypanosoma brucei: A Single Protein in Place of the Multicomponent RNA-Protein Complex

SummaryRNase P is the endonuclease that removes 5′ extensions from tRNA precursors. In its best-known form, the enzyme is composed of a catalytic RNA and a protein moiety variable in number and mass. This ribonucleoprotein enzyme is widely considered ubiquitous and apparently reached its highest complexity in the eukaryal nucleus, where it is typically composed of at least ten subunits. Here, we show that in the protist Trypanosoma brucei, two proteins are the sole forms of RNase P. They localize to the nucleus and the mitochondrion, respectively, and have RNase P activity each on their own. The protein-RNase P is, moreover, capable of replacing nuclear RNase P in yeast cells. This shows that complex ribonucleoprotein structures and RNA catalysis are not necessarily required to support tRNA 5′ end formation in eukaryal cells.

Systems genetics: challenges and developing strategies

Systems genetics is a new discipline based on the transcription mapping, which is also called “genetical genomics”. In recent years, systems genetics has become more practical because of advances in science and technology. Analysis of expression quantitative trait loci (eQTLs) is an emerging technique in which individuals are genotyped across a panel of genetic markers and, simultaneously, phenotyped using DNA microarrays. Depending on eQTL mapping, one can infer the underlying regulatory network responsible for complex diseases or quantitative trait phenotypes. Systems genetics approaches integrate DNA sequence variation, variation in transcript abundance and other molecular phenotypes and variation in organismal phenotypes in a linkage or association mapping population, and allow us to interpret quantitative genetic variation in terms of biologically meaningful causal networks of correlated transcripts. These approaches have been made possible due to the development of massively parallel technologies for quantifying genome-wide levels of transcript abundance. The predictive power of the networks could be enhanced by more systematically integrating protein-protein interactions, protein-DNA interactions, protein-RNA interactions, RNA-RNA interactions, protein state information, methylation state, and interactions with metabolites. Systems genetics research will change the traditional approaches based on reductionism, and allows us to reconsider the living phenomenon and complex disease mechanism. Systems genetics benefits from varied “omics” researches (such as transcriptomics, metabolomics, and phenomics) and the development of bioinformatics tools and mathematical modeling, and will become mature in the near future like many other branches of genetics. Systems genetics is leading researchers to understand genetics systems from holism’s viewpoint, and will open a wide field of vision for genetics researchers in systems biology era.

A Guideline to Family-Wide Comparative State-of-the-Art Quantitative RT-PCR Analysis Exemplified with a Brassicaceae Cross-Species Seed Germination Case Study

Comparative biology includes the comparison of transcriptome and quantitative real-time RT-PCR (qRT-PCR) data sets in a range of species to detect evolutionarily conserved and divergent processes. Transcript abundance analysis of target genes by qRT-PCR requires a highly accurate and robust workflow. This includes reference genes with high expression stability (i.e., low intersample transcript abundance variation) for correct target gene normalization. Cross-species qRT-PCR for proper comparative transcript quantification requires reference genes suitable for different species. We addressed this issue using tissue-specific transcriptome data sets of germinating Lepidium sativum seeds to identify new candidate reference genes. We investigated their expression stability in germinating seeds of L. sativum and Arabidopsis thaliana by qRT-PCR, combined with in silico analysis of Arabidopsis and Brassica napus microarray data sets. This revealed that reference gene expression stability is higher for a given developmental process between distinct species than for distinct developmental processes within a given single species. The identified superior cross-species reference genes may be used for family-wide comparative qRT-PCR analysis of Brassicaceae seed germination. Furthermore, using germinating seeds, we exemplify optimization of the qRT-PCR workflow for challenging tissues regarding RNA quality, transcript stability, and tissue abundance. Our work therefore can serve as a guideline for moving beyond Arabidopsis by establishing high-quality cross-species qRT-PCR.

Directed partial correlation: inferring large-scale gene regulatory network through induced topology disruptions.

Inferring regulatory relationships among many genes based on their temporal variation in transcript abundance has been a popular research topic. Due to the nature of microarray experiments, classical tools for time series analysis lose power since the number of variables far exceeds the number of the samples. In this paper, we describe some of the existing multivariate inference techniques that are applicable to hundreds of variables and show the potential challenges for small-sample, large-scale data. We propose a directed partial correlation (DPC) method as an efficient and effective solution to regulatory network inference using these data. Specifically for genomic data, the proposed method is designed to deal with large-scale datasets. It combines the efficiency of partial correlation for setting up network topology by testing conditional independence, and the concept of Granger causality to assess topology change with induced interruptions. The idea is that when a transcription factor is induced artificially within a gene network, the disruption of the network by the induction signifies a genes role in transcriptional regulation. The benchmarking results using GeneNetWeaver, the simulator for the DREAM challenges, provide strong evidence of the outstanding performance of the proposed DPC method. When applied to real biological data, the inferred starch metabolism network in Arabidopsis reveals many biologically meaningful network modules worthy of further investigation. These results collectively suggest DPC is a versatile tool for genomics research. The R package DPC is available for download (http://code.google.com/p/dpcnet/).

Identification of candidate genes involved in responses of Arabidopsis to polychlorinated biphenyls based on microarray analysis

Polychlorinated biphenyls (PCBs) are considered as persistent organic pollutants. Their long term accumulation in soils and progressive bioaccumulation in the food chain makes PCBs a risk for human health. Phytoremediation, or the use of plants to enhance pollutant dissipation, is an efficient technique for the removal of organic compounds, including PCBs. However, the molecular basis of plant responses to PCBs has not been clearly elucidated. Here we studied the effects on growth and the transcriptional profile of Arabidopsis thaliana after exposure to 2,2′,3,3′-tetrachlorobiphenyl, a representative PCB. A comprehensive survey of global gene expression response to this PCB was done by using Arabidopsis (V4) gene expression microarray (4 × 44 K) to quantify the spatio-temporal variations in transcript abundance of 42,000 genes. The results revealed a coordinated induction or suppression of 146 or 148 genes, respectively. Of these, expression of 40 genes was validated by RT-PCR analysis. The functional classification of these PCB-responsive genes indicated their involvement in various metabolic pathways, such as ion transport, signal transduction, transcriptional regulation, and other processes related to growth and development.

Complementary genetic and genomic approaches help characterize the linkage group I seed protein QTL in soybean

BackgroundThe nutritional and economic value of many crops is effectively a function of seed protein and oil content. Insight into the genetic and molecular control mechanisms involved in the deposition of these constituents in the developing seed is needed to guide crop improvement. A quantitative trait locus (QTL) on Linkage Group I (LG I) of soybean (Glycine max (L.) Merrill) has a striking effect on seed protein content.ResultsA soybean near-isogenic line (NIL) pair contrasting in seed protein and differing in an introgressed genomic segment containing the LG I protein QTL was used as a resource to demarcate the QTL region and to study variation in transcript abundance in developing seed. The LG I QTL region was delineated to less than 8.4 Mbp of genomic sequence on chromosome 20. Using Affymetrix® Soy GeneChip and high-throughput Illumina® whole transcriptome sequencing platforms, 13 genes displaying significant seed transcript accumulation differences between NILs were identified that mapped to the 8.4 Mbp LG I protein QTL region.ConclusionsThis study identifies gene candidates at the LG I protein QTL for potential involvement in the regulation of protein content in the soybean seed. The results demonstrate the power of complementary approaches to characterize contrasting NILs and provide genome-wide transcriptome insight towards understanding seed biology and the soybean genome.

Response to zinc deficiency of two rice lines with contrasting tolerance is determined by root growth maintenance and organic acid exudation rates, and not by zinc‐transporter activity

*Zinc (Zn)-deficient soils constrain rice (Oryza sativa) production and cause Zn malnutrition. The identification of Zn-deficiency-tolerant rice lines indicates that breeding might overcome these constraints. Here, we seek to identify processes underlying Zn-deficiency tolerance in rice at the physiological and transcriptional levels. *A Zn-deficiency-tolerant line RIL46 acquires Zn more efficiently and produces more biomass than its nontolerant maternal line (IR74) at low [Zn](ext) under field conditions. We tested if this was the result of increased expression of Zn(2+) transporters; increased root exudation of deoxymugineic acid (DMA) or low-molecular-weight organic acids (LMWOAs); and/or increased root production. Experiments were performed in field and controlled environment conditions. *There was little genotypic variation in transcript abundance of Zn-responsive root Zn(2+)-transporters between the RIL46 and IR74. However, root exudation of DMA and LMWOA was greater in RIL46, coinciding with increased root expression of putative ligand-efflux genes. Adventitious root production was maintained in RIL46 at low [Zn](ext), correlating with altered expression of root-specific auxin-responsive genes. *Zinc-deficiency tolerance in RIL46 is most likely the result of maintenance of root growth, increased efflux of Zn ligands, and increased uptake of Zn-ligand complexes at low [Zn](ext); these traits are potential breeding targets.

An eQTL Analysis of Partial Resistance to Puccinia hordei in Barley

BackgroundGenetic resistance to barley leaf rust caused by Puccinia hordei involves both R genes and quantitative trait loci. The R genes provide higher but less durable resistance than the quantitative trait loci. Consequently, exploring quantitative or partial resistance has become a favorable alternative for controlling disease. Four quantitative trait loci for partial resistance to leaf rust have been identified in the doubled haploid Steptoe (St)/Morex (Mx) mapping population. Further investigations are required to study the molecular mechanisms underpinning partial resistance and ultimately identify the causal genes.Methodology/Principal FindingsWe explored partial resistance to barley leaf rust using a genetical genomics approach. We recorded RNA transcript abundance corresponding to each probe on a 15K Agilent custom barley microarray in seedlings from St and Mx and 144 doubled haploid lines of the St/Mx population. A total of 1154 and 1037 genes were, respectively, identified as being P. hordei-responsive among the St and Mx and differentially expressed between P. hordei-infected St and Mx. Normalized ratios from 72 distant-pair hybridisations were used to map the genetic determinants of variation in transcript abundance by expression quantitative trait locus (eQTL) mapping generating 15685 eQTL from 9557 genes. Correlation analysis identified 128 genes that were correlated with resistance, of which 89 had eQTL co-locating with the phenotypic quantitative trait loci (pQTL). Transcript abundance in the parents and conservation of synteny with rice allowed us to prioritise six genes as candidates for Rphq11, the pQTL of largest effect, and highlight one, a phospholipid hydroperoxide glutathione peroxidase (HvPHGPx) for detailed analysis.Conclusions/SignificanceThe eQTL approach yielded information that led to the identification of strong candidate genes underlying pQTL for resistance to leaf rust in barley and on the general pathogen response pathway. The dataset will facilitate a systems appraisal of this host-pathogen interaction and, potentially, for other traits measured in this population.

Geographical genomics of human leukocyte gene expression variation in southern Morocco

Studies of the genetics of gene expression reveal expression SNPs that explain variation in transcript abundance. Here we address the robustness of eSNP associations to environmental geography and population structure in a comparison of 194 Arab and Amazigh individuals from a city and two villages in southern Morocco. Gene expression differed between pairs of locations for up to a third of all transcripts, with notable enrichment for ribosomal biosynthesis and oxidative phosphorylation. Robust associations were observed in the leukocyte samples with cis-eSNPs (P < 10−08) for 346 genes, and trans-eSNPs (P < 10−11) with 10 genes. All of these were consistent across the three sample locations and after controlling for ethnicity and relatedness. No evidence for large-effect trans-acting mediators of the pervasive environmental influence was found and instead genetic and environmental factors acted in a largely additive manner.

Widespread variation in transcript abundance within and across developmental stages of Trypanosoma brucei

BackgroundTrypanosoma brucei, the causative agent of African sleeping sickness, undergoes a complex developmental cycle that takes place in mammalian and insect hosts and is accompanied by changes in metabolism and cellular morphology. While differences in mRNA expression have been described for many genes, genome-wide expression analyses have been largely lacking. Trypanosomatids represent a unique case in eukaryotes in that they transcribe protein-coding genes as large polycistronic units, and rarely regulate gene expression at the level of transcription initiation.ResultsHere we present a comprehensive analysis of mRNA expression in several stages of parasite development. Utilizing microarrays that have multiple copies of multiple probes for each gene, we were able to demonstrate with a high degree of statistical confidence that approximately one-fourth of genes show differences in mRNA expression levels in the stages examined. These include complex patterns of gene expression within gene families, including the large family of variant surface glycoproteins (VSGs) and their relatives, where we have identified a number of constitutively expressed family members. Furthermore, we were able to assess the relative abundance of all transcripts in each stage, identifying the genes that are either weakly or highly expressed. Very few genes show no evidence of expression.ConclusionDespite the lack of gene regulation at the level of transcription initiation, our results reveal extensive regulation of mRNA abundance associated with different life cycle and growth stages. In addition, analysis of variant surface glycoprotein gene expression reveals a more complex picture than previously thought. These data provide a valuable resource to the community of researchers studying this lethal agent.