Allelic variation at a single gene increases food value in a drought-tolerant staple cereal

The production of adequate agricultural outputs to support the growing human population places great demands on agriculture, especially in light of ever-greater restrictions on input resources. Sorghum is a drought-adapted cereal capable of reliable production where other cereals fail, and thus represents a good candidate to address food security as agricultural inputs of water and arable land grow scarce. A long-standing issue with sorghum grain is that it has an inherently lower digestibility. Here we show that a low-frequency allele type in the starch metabolic gene, pullulanase, is associated with increased digestibility, regardless of genotypic background. We also provide evidence that the beneficial allele type is not associated with deleterious pleiotropic effects in the modern field environment. We argue that increasing the digestibility of an adapted crop is a viable way forward towards addressing food security while maximizing water and land-use efficiency.

Extending genetic linkage analysis to diffusion tensor images to map single gene effects on brain fiber architecture

We extended genetic linkage analysis - an analysis widely used in quantitative genetics - to 3D images to analyze single gene effects on brain fiber architecture. We collected 4 Tesla diffusion tensor images (DTI) and genotype data from 258 healthy adult twins and their non-twin siblings. After high-dimensional fluid registration, at each voxel we estimated the genetic linkage between the single nucleotide polymorphism (SNP), Val66Met (dbSNP number rs6265), of the BDNF gene (brain-derived neurotrophic factor) with fractional anisotropy (FA) derived from each subject's DTI scan, by fitting structural equation models (SEM) from quantitative genetics. We also examined how image filtering affects the effect sizes for genetic linkage by examining how the overall significance of voxelwise effects varied with respect to full width at half maximum (FWHM) of the Gaussian smoothing applied to the FA images. Raw FA maps with no smoothing yielded the greatest sensitivity to detect gene effects, when corrected for multiple comparisons using the false discovery rate (FDR) procedure. The BDNF polymorphism significantly contributed to the variation in FA in the posterior cingulate gyrus, where it accounted for around 90-95% of the total variance in FA. Our study generated the first maps to visualize the effect of the BDNF gene on brain fiber integrity, suggesting that common genetic variants may strongly determine white matter integrity.

Allelic variation at a single gene increases food value in a drought-tolerant staple cereal

The production of adequate agricultural outputs to support the growing human population places great demands on agriculture, especially in light of ever-greater restrictions on input resources. Sorghum is a drought-adapted cereal capable of reliable production where other cereals fail, and thus represents a good candidate to address food security as agricultural inputs of water and arable land grow scarce. A long-standing issue with sorghum grain is that it has an inherently lower digestibility. Here we show that a low-frequency allele type in the starch metabolic gene, pullulanase, is associated with increased digestibility, regardless of genotypic background. We also provide evidence that the beneficial allele type is not associated with deleterious pleiotropic effects in the modern field environment. We argue that increasing the digestibility of an adapted crop is a viable way forward towards addressing food security while maximizing water and land-use efficiency.

Unravelling the genetic basis of renal diseases; from single gene to multifactorial disorders

Chronic kidney disease is common with up to 5% of the adult population reported to have an estimated glomerular filtration rate of

Harnessing the complexity of gene expression data from cancer: from single gene to structural pathway methods

High-dimensional gene expression data provide a rich source of information because they capture the expression level of genes in dynamic states that reflect the biological functioning of a cell. For this reason, such data are suitable to reveal systems related properties inside a cell, e.g., in order to elucidate molecular mechanisms of complex diseases like breast or prostate cancer. However, this is not only strongly dependent on the sample size and the correlation structure of a data set, but also on the statistical hypotheses tested. Many different approaches have been developed over the years to analyze gene expression data to (I) identify changes in single genes, (II) identify changes in gene sets or pathways, and (III) identify changes in the correlation structure in pathways. In this paper, we review statistical methods for all three types of approaches, including subtypes, in the context of cancer data and provide links to software implementations and tools and address also the general problem of multiple hypotheses testing. Further, we provide recommendations for the selection of such analysis methods.

Transcriptional functions of DNA Topoisomerases at a genome-wide scale and a single gene levels

The DNA topology is an important modifier of DNA functions. Torsional stress is generated when right handed DNA is either over- or underwound, producing structural deformations which drive or are driven by processes such as replication, transcription, recombination and repair. DNA topoisomerases are molecular machines that regulate the topological state of the DNA in the cell. These enzymes accomplish this task by either passing one strand of the DNA through a break in the opposing strand or by passing a region of the duplex from the same or a different molecule through a double-stranded cut generated in the DNA. Because of their ability to cut one or two strands of DNA they are also target for some of the most successful anticancer drugs used in standard combination therapies of human cancers. An effective anticancer drug is Camptothecin (CPT) that specifically targets DNA topoisomerase 1 (TOP 1). The research project of the present thesis has been focused on the role of human TOP 1 during transcription and on the transcriptional consequences associated with TOP 1 inhibition by CPT in human cell lines. Previous findings demonstrate that TOP 1 inhibition by CPT perturbs RNA polymerase (RNAP II) density at promoters and along transcribed genes suggesting an involvement of TOP 1 in RNAP II promoter proximal pausing site. Within the transcription cycle, promoter pausing is a fundamental step the importance of which has been well established as a means of coupling elongation to RNA maturation. By measuring nascent RNA transcripts bound to chromatin, we demonstrated that TOP 1 inhibition by CPT can enhance RNAP II escape from promoter proximal pausing site of the human Hypoxia Inducible Factor 1 (HIF-1) and c-MYC genes in a dose dependent manner. This effect is dependent from Cdk7/Cdk9 activities since it can be reversed by the kinases inhibitor DRB. Since CPT affects RNAP II by promoting the hyperphosphorylation of its Rpb1 subunit the findings suggest that TOP 1inhibition by CPT may increase the activity of Cdks which in turn phosphorylate the Rpb1 subunit of RNAP II enhancing its escape from pausing. Interestingly, the transcriptional consequences of CPT induced topological stress are wider than expected. CPT increased co-transcriptional splicing of exon1 and 2 and markedly affected alternative splicing at exon 11. Surprisingly despite its well-established transcription inhibitory activity, CPT can trigger the production of a novel long RNA (5’aHIF-1) antisense to the human HIF-1 mRNA and a known antisense RNA at the 3’ end of the gene, while decreasing mRNA levels. The effects require TOP 1 and are independent from CPT induced DNA damage. Thus, when the supercoiling imbalance promoted by CPT occurs at promoter, it may trigger deregulation of the RNAP II pausing, increased chromatin accessibility and activation/derepression of antisense transcripts in a Cdks dependent manner. A changed balance of antisense transcripts and mRNAs may regulate the activity of HIF-1 and contribute to the control of tumor progression After focusing our TOP 1 investigations at a single gene level, we have extended the study to the whole genome by developing the “Topo-Seq” approach which generates a map of genome-wide distribution of sites of TOP 1 activity sites in human cells. The preliminary data revealed that TOP 1 preferentially localizes at intragenic regions and in particular at 5’ and 3’ ends of genes. Surprisingly upon TOP 1 downregulation, which impairs protein expression by 80%, TOP 1 molecules are mostly localized around 3’ ends of genes, thus suggesting that its activity is essential at these regions and can be compensate at 5’ ends. The developed procedure is a pioneer tool for the detection of TOP 1 cleavage sites across the genome and can open the way to further investigations of the enzyme roles in different nuclear processes.

Deoxycytidine kinase and deoxyguanosine kinase of Lactobacillus acidophilus R-26 are colinear products of a single gene

Three of the four deoxynucleoside kinases required for growth of Lactobacillus acidophilus R-26 exist as heterodimeric pairs specific for deoxyadenosine (dAK) and deoxycytidine (dCK) or dAK and deoxyguanosine (dGK). However, only two tandem genes, dak/dgk, are found, and are expressed only as dAK/dGK in transformed Escherichia coli. Sequencing peptides spanning 63% of the native dCK subunit revealed a sequence identical to that deduced from dgk (beginning MTVIVL···), except that dCK lacks residues 2 and 3 (dCK is M··IVL; dGK is ·TVIVL). Also, mass spectrometry indicates that native dCK and dGK subunits are identical in mass adjusted for the first three residues. Furthermore, the native enzymes have identical isoelectric pH values, indicating an equal number of charged residues. To enable E. coli to express peptide having the native dCK sequence, codons 2 and 3 were deleted from the dgk portion of the tandem genes, resulting in expression of protein having the specificities and regulatory properties of native dAK/dCK, including heterotropic stimulation of dAK activity by deoxycytidine or dCTP (not deoxyguanosine or dGTP) and end-product inhibition of the respective activities by dATP and dCTP. Subcloning normal and mutant dgk yielded homodimeric dGK and dCK, respectively. The dCK homodimer strongly resembles human dCK, with a low Km for deoxycytidine, the ability to phosphorylate deoxyadenosine and deoxyguanosine at much higher Km values, and end-product inhibition by dCTP. Thus two distinct and specific enzymes evidently are derived from a single Lactobacillus gene. The mechanism by which this occurs in vivo has yet to be elucidated.

Quantifying single gene copy number by measuring fluorescent probe lengths on combed genomic DNA

An approach was developed for the quantification of subtle gains and losses of genomic DNA. The approach relies on a process called molecular combing. Molecular combing consists of the extension and alignment of purified molecules of genomic DNA on a glass coverslip. It has the advantage that a large number of genomes can be combed per coverslip, which allows for a statistically adequate number of measurements to be made on the combed DNA. Consequently, a high-resolution approach to mapping and quantifying genomic alterations is possible. The approach consists of applying fluorescence hybridization to the combed DNA by using probes to identify the amplified region. Measurements then are made on the linear hybridization signals to ascertain the region's exact size. The reliability of the approach first was tested for low copy number amplifications by determining the copy number of chromosome 21 in a normal and trisomy 21 cell line. It then was tested for high copy number amplifications by quantifying the copy number of an oncogene amplified in the tumor cell line GTL-16. These results demonstrate that a wide range of amplifications can be accurately and reliably quantified. The sensitivity and resolution of the approach likewise was assessed by determining the copy number of a single allele (160 kb) alteration.

A single gene of chloroplast origin codes for mitochondrial and chloroplastic methionyl–tRNA synthetase in Arabidopsis thaliana

One-fifth of the tRNAs used in plant mitochondrial translation is coded for by chloroplast-derived tRNA genes. To understand how aminoacyl–tRNA synthetases have adapted to the presence of these tRNAs in mitochondria, we have cloned an Arabidopsis thaliana cDNA coding for a methionyl–tRNA synthetase. This enzyme was chosen because chloroplast-like elongator tRNAMet genes have been described in several plant species, including A. thaliana. We demonstrate here that the isolated cDNA codes for both the chloroplastic and the mitochondrial methionyl–tRNA synthetase (MetRS). The protein is transported into isolated chloroplasts and mitochondria and is processed to its mature form in both organelles. Transient expression assays using the green fluorescent protein demonstrated that the N-terminal region of the MetRS is sufficient to address the protein to both chloroplasts and mitochondria. Moreover, characterization of MetRS activities from mitochondria and chloroplasts of pea showed that only one MetRS activity exists in each organelle and that both are indistinguishable by their behavior on ion exchange and hydrophobic chromatographies. The high degree of sequence similarity between A. thaliana and Synechocystis MetRS strongly suggests that the A. thaliana MetRS gene described here is of chloroplast origin.

Bidirectional imprinting of a single gene: GNAS1 encodes maternally, paternally, and biallelically derived proteins

The GNAS1 gene encodes the α subunit of the guanine nucleotide-binding protein Gs, which couples signaling through peptide hormone receptors to cAMP generation. GNAS1 mutations underlie the hormone resistance syndrome pseudohypoparathyroidism type Ia (PHP-Ia), so the maternal inheritance displayed by PHP-Ia has raised suspicions that GNAS1 is imprinted. Despite this suggestion, in most tissues Gsα is biallelically encoded. In contrast, the large G protein XLαs, also encoded by GNAS1, is paternally derived. Because the inheritance of PHP-Ia predicts the existence of maternally, rather than paternally, expressed transcripts, we have investigated the allelic origin of other mRNAs derived from GNAS1. We find this gene to be remarkable in the complexity of its allele-specific regulation. Two upstream promoters, each associated with a large coding exon, lie only 11 kb apart, yet show opposite patterns of allele-specific methylation and monoallelic transcription. The more 5′ of these exons encodes the neuroendocrine secretory protein NESP55, which is expressed exclusively from the maternal allele. The NESP55 exon is 11 kb 5′ to the paternally expressed XLαs exon. The transcripts from these two promoters both splice onto GNAS1 exon 2, yet share no coding sequences. Despite their structural unrelatedness, the encoded proteins, of opposite allelic origin, both have been implicated in regulated secretion in neuroendocrine tissues. Remarkably, maternally (NESP55), paternally (XLαs), and biallelically (Gsα) derived proteins all are produced by different patterns of promoter use and alternative splicing of GNAS1, a gene showing simultaneous imprinting in both the paternal and maternal directions.

The glyoxysomal and plastid molecular chaperones (70-kDa heat shock protein) of watermelon cotyledons are encoded by a single gene

The monoclonal a-70-kDa heat shock protein (hsp70) antibody recognizes in crude extracts from watermelon (Citrullus vulgaris) cotyledons two hsps with molecular masses of 70 and 72 kDa. Immunocytochemistry on watermelon cotyledon tissue and on isolated glyoxysomes identified hsp70s in the matrix of glyoxysomes and plastids. Affinity purification and partial amino acid determination revealed the 70-kDa protein to share high sequence identity with cytosolic hsp70s from a number of plant species, while the 72 kDa protein was very similar to plastid hsp70s from pea and cucumber. A full-length cDNA clone encoding the 72-kDa hsp70 was isolated and identified two start methionines in frame within the N-terminal presequence leading either to an N-terminal extension of 67 amino acids or to a shorter one of 47 amino acids. The longer presequence was necessary and sufficient to target a reporter protein into watermelon proplastids in vitro. The shorter extension starting from the second methionine within the long version harbored a consensus peroxisomal targeting signal (RT-X5-KL) that directed in vivo a reporter protein into peroxisomes of the yeast Hansenula polymorpha. Peroxisomal targeting was however prevented, when the 67-residue presequence was fused to the reporter protein, indicating that the peroxisomal targeting signal 2 information is hidden in this context. We propose that the 72-kDa hsp70 is encoded by a single gene, but targeted alternatively into two organelles by the modulated use of its presequence.