Sequencing of TGF-beta pathway genes in familial cases of intracranial aneurysm.

BACKGROUND AND PURPOSE: Familial aggregation of intracranial aneurysms (IA) strongly suggests a genetic contribution to pathogenesis. However, genetic risk factors have yet to be defined. For families affected by aortic aneurysms, specific gene variants have been identified, many affecting the receptors to transforming growth factor-beta (TGF-beta). In recent work, we found that aortic and intracranial aneurysms may share a common genetic basis in some families. We hypothesized, therefore, that mutations in TGF-beta receptors might also play a role in IA pathogenesis. METHODS: To identify genetic variants in TGF-beta and its receptors, TGFB1, TGFBR1, TGFBR2, ACVR1, TGFBR3, and ENG were directly sequenced in 44 unrelated patients with familial IA. Novel variants were confirmed by restriction digestion analyses, and allele frequencies were analyzed in cases versus individuals without known intracranial disease. Similarly, allele frequencies of a subset of known SNPs in each gene were also analyzed for association with IA. RESULTS: No mutations were found in TGFB1, TGFBR1, TGFBR2, or ACVR1. Novel variants identified in ENG (p.A60E) and TGFBR3 (p.W112R) were not detected in at least 892 reference chromosomes. ENG p.A60E showed significant association with familial IA in case-control studies (P=0.0080). No association with IA could be found for any of the known polymorphisms tested. CONCLUSIONS: Mutations in TGF-beta receptor genes are not a major cause of IA. However, we identified rare variants in ENG and TGFBR3 that may be important for IA pathogenesis in a subset of families.

A trilocus sequence typing scheme for hospital epidemiology and subspecies differentiation of an important nosocomial pathogen, Enterococcus faecalis.

In this study, we present a trilocus sequence typing (TLST) scheme based on intragenic regions of two antigenic genes, ace and salA (encoding a collagen/laminin adhesin and a cell wall-associated antigen, respectively), and a gene associated with antibiotic resistance, lsa (encoding a putative ABC transporter), for subspecies differentiation of Enterococcus faecalis. Each of the alleles was analyzed using 50 E. faecalis isolates representing 42 diverse multilocus sequence types (ST(M); based on seven housekeeping genes) and four groups of clonally linked (by pulsed-field gel electrophoresis [PFGE]) isolates. The allelic profiles and/or concatenated sequences of the three genes agreed with multilocus sequence typing (MLST) results for typing of 49 of the 50 isolates; in addition to the one exception, two isolates were found to have identical TLST types but were single-locus variants (differing by a single nucleotide) by MLST and were therefore also classified as clonally related by MLST. TLST was also comparable to PFGE for establishing short-term epidemiological relationships, typing all isolates classified as clonally related by PFGE with the same type. TLST was then applied to representative isolates (of each PFGE subtype and isolation year) of a collection of 48 hospital isolates and demonstrated the same relationships between isolates of an outbreak strain as those found by MLST and PFGE. In conclusion, the TLST scheme described here was shown to be successful for investigating short-term epidemiology in a hospital setting and may provide an alternative to MLST for discriminating isolates.

Computational identification and functional validation of regulatory motifs in cartilage-expressed genes.

Chondrocyte gene regulation is important for the generation and maintenance of cartilage tissues. Several regulatory factors have been identified that play a role in chondrogenesis, including the positive transacting factors of the SOX family such as SOX9, SOX5, and SOX6, as well as negative transacting factors such as C/EBP and delta EF1. However, a complete understanding of the intricate regulatory network that governs the tissue-specific expression of cartilage genes is not yet available. We have taken a computational approach to identify cis-regulatory, transcription factor (TF) binding motifs in a set of cartilage characteristic genes to better define the transcriptional regulatory networks that regulate chondrogenesis. Our computational methods have identified several TFs, whose binding profiles are available in the TRANSFAC database, as important to chondrogenesis. In addition, a cartilage-specific SOX-binding profile was constructed and used to identify both known, and novel, functional paired SOX-binding motifs in chondrocyte genes. Using DNA pattern-recognition algorithms, we have also identified cis-regulatory elements for unknown TFs. We have validated our computational predictions through mutational analyses in cell transfection experiments. One novel regulatory motif, N1, found at high frequency in the COL2A1 promoter, was found to bind to chondrocyte nuclear proteins. Mutational analyses suggest that this motif binds a repressive factor that regulates basal levels of the COL2A1 promoter.

Analysis of mutations in $\beta$-tubulin genes affecting tubulin stability and assembly

Cmd4 is a colcemid-sensitive CHO cell line that is temperature sensitive for growth and expresses an altered $\beta$-tubulin, $\beta\sb1$. One revertant of this cell line, D2, exhibits a further alteration in $\beta\sb1$ resulting in an acidic shift in its isoelectric point and a decrease in its molecular weight to 40 kD, as measured by two dimensional gel electrophoresis. This $\beta$-tubulin variant has been shown to be assembly-defective and unstable. Characterization of the mutant $\beta\sb1$ in D2 by high pressure liquid chromatography (HPLC) revealed the loss of methionine containing tryptic peptides 7,8,9, and 10. Southern analysis of the genomic DNA digested with several different restriction enzymes resulted in the appearance of new restriction fragments 250 base pairs shorter than the corresponding fragments from the wild-type $\beta\sb1$-tubulin gene. Northern analysis on mRNA from D2 revealed two new message products that also differed by 250 bases from the corresponding wild type $\beta$-tubulin transcripts. To precisely define the region of the alteration, cloning and sequencing of the mutant and wild type genomic $\beta$-tubulin genes were conducted. A size-selected EcoRI genomic library was prepared using the Stratagene lambda Zap II phage cloning system. Using subclones of CHO $\beta$-tubulin cDNA as probes, a 2.5 kb wild type clone and a 2.3 kb mutant clone were identified from this library. Each of these was shown to contain a portion of the gene extending from intron 3 through the end of the coding sequence in exon 4 and into the 3$\sp\prime$ untranslated region on the basis of alignment with the published human $\beta$-tubulin sequence. Sequencing of the mutant 2.3 kb clone revealed that the mutation is due to a 246 base pair internal deletion in exon 4 (base pair 756-1001) that encodes amino acids 253-334. This deletion results in the loss of a putative binding site for GTP which could potentially explain the phenotype of this mutant $\beta$-tubulin. Also sequence comparison of the 3$\sp\prime$ untranslated region between different species revealed the conservation of 200 base pairs with 78% homology. It is proposed that this region could play an important role in the regulation of $\beta$-tubulin gene expression. ^

Mapping and cloning of genes from portions of the human genome using somatic cell hybrids

Human x rodent somatic cell hybrids have played an important role in human genetics research. They have been especially useful for assigning genes to chromosomes and isolating DNA markers from specific regions of the human genome.^ By employing a combination of somatic cell genetic, recombinant DNA, and cytogenetic techniques, human DNA excision repair gene ERCC4 was mapped regionally to human 16p13.13-13.2, even though the gene has not been cloned. Human x Chinese hamster ovary (CHO) cell hybrids selected for human ERCC4 activity and containing 16p13.1-p13.3 as the only human genetic material were identified. These hybrids were used to order DNA markers located in 16p13.1-p13.3. New DNA markers physically close to ERCC4 were isolated from such hybrids. Using amplified human DNA from the hybrids as probe in fluorescent in situ hybridization, the short arm breakpoint in the chromosome 16 inversion associated with acute myelomonocytic leukemia (AMML) was found to be physically close to the ERCC4 gene. The physical mapping and eventually, the cloning of the ERCC4 gene, will benefit the understanding of the DNA repair system and the study of other important biomedical problems such as tumorigenesis.^ To facilitate the cloning of ERCC4 gene and, in general, the cloning of genes from any defined regions of the human genome, a method was developed for the direct isolation of human transcribed genes ffom somatic cell hybrids. cDNA was prepared from human x rodent hybrid by using consensus 5$\sp\prime$ splice site sequences as primers. These primers were designed to select immature, unspliced messenger RNA (still retaining species specific repeat sequences) as templates. Screening of a derived cDNA library for human repeat sequences resulted in the isolation of human clones at the anticipated frequency with characteristics expected of exons of transcribed human genes. The usefulness of the splice site specific primers was analyzed and the cDNA synthesis conditions with these primers were optimized. The procedure was shown to be sensitive enough to clone weakly expressed genes. Studying the expression of the represented genes with the isolated clones was shown to be feasible. Such regional specific human gene fragments will be very valuable for many human genetic studies such as the search of inherited disease genes and the construction of a cDNA map of the human genome. ^

Aboral ectoderm-specific expression and regulation of the LpS1 genes of {\it Lytechinus pictus\/}

The Spec genes serve as molecular markers for examining the ontogeny of the aboral ectoderm lineage of sea urchin embryos. These genes are activated at late-cleavage stage only in cells contributing to the aboral ectoderm of Strongylocentrotus purpuratus and encode 14,000-17,000 Da calcium-binding proteins. A comparative analysis was undertaken to better understand the mechanisms underlying the activation and function of the Spec genes by investigating Spec homologues from Lytechinus pictus, a distantly related sea urchin. Spec antibodies cross-reacted with 34,000 Da proteins in L. pictus embryos that displayed a similar ontogenetic pattern to that of Spec proteins. One cDNA clone, LpS1, was isolated by hybridization to a synthetic oligonucleotide corresponding to a calcium-binding domain or EF-hand. The LpS1 mRNA has developmental properties similar to those of the Spec mRNAs. LpS1 encodes a 34,000 Da protein containing eight EF-hand domains, which share structural homology with the Spec EF-hands; however, little else in the protein sequence is conserved, implying that calcium-binding is important for Spec protein function. Genomic DNA blot analysis showed two LpS1 genes, LpS1$\alpha$ and LpS1$\beta$, in L. pictus. Partial gene structures for both LpS1$\alpha$ and $\beta$ were constructed based on genomic clones isolated from an L. pictus genomic library. These revealed internal duplications of the LpS1 genes that accounted for the eight EF-hand domains in the LpS1 proteins. Sequencing analysis showed there was little in common among the 5$\sp\prime$-flanking regions of the LpS1 and Spec genes except for the presence of a binding site for the transcription factor USF.^ A sea urchin gene-transfer expression system showed that 762 base pairs (bp) of 5$\sp\prime$-flanking DNA from the LpS1$\beta$ gene were sufficient for correct temporal and spatial expression of reporter genes in sea urchin embryos. Deletions at the 5$\sp\prime$ end to 511, 368, or 108bp resulted in a 3-4 fold decrease in chloramphenicol acetyltransferase (CAT) activity and disrupted the restricted activation of the lac Z gene in aboral ectoderm cells.^ A full-length Spec1 protein and a truncated LpS1 protein were induced and partially purified from an in vitro expression system. (Abstract shortened with permission of author.) ^

Transcriptional analysis of liver-specific expression and differential interleukin-6 response of rat kininogen genes

Analyses of rat T1 kininogen gene/chloramphenicol acetyltransferase (T1K/CAT) constructs revealed two regions important for tissue-specific and induced regulation of T1 kininogen.^ Although the T1 kininogen gene is inducible by inflammatory cytokines, a highly homologous K kininogen gene is minimally responsive. Moreover, the basal expression of a KK/CAT construct was 5- to 7-fold higher than that of the analogous T1K/CAT construct. To examine the molecular basis of this differential regulation, a series of promoter swapping experiments was carried out. Our transfection results showed that at least two regions in the K kininogen gene are important for its high basal expression: a distal 19-bp region (C box) constituted a binding site for CCAAT/enhancer binding protein (C/EBP) family proteins and a proximal 66-bp region contained two adjacent binding sites for hepatocyte nuclear factor-3 (HNF-3). The distal HNF-3 binding site from the K kininogen promoter demonstrated a stronger affinity than that from the T1 kininogen promoter. Since C/EBP and HNF-3 are highly enriched in the liver and known to enhance transcription of liver-specific genes, differential binding affinities of these factors accounted for the higher basal expression of the K kininogen gene.^ In contrast to the K kininogen C box, the T1 kininogen C box does not bind C/EBP presumably due to their two-nucleotide divergence. This sequence divergence, however, converts it to a consensus binding sequence for two IL-6-inducible transcription factors--IL-6 response element binding protein and acute-phase response factor. To functionally determine whether C box sequences are important for their differential acute-phase response, T1 and K kininogen C boxes were swapped and analyzed after transfection into Hep3B cells. Our results showed that the T1 kininogen C box is indeed one of the IL-6 response elements in T1 kininogen promoter. Furthermore, its function can be modulated by a 5$\sp\prime$-adjacent C/EBP-binding site (B box) whose mutation significantly reduced the overall induced activity. Moreover, this B box is the target site for binding and transactivation of another IL-6 inducible transcription factor C/EBP$\delta.$ Evolutionary divergence of a few critical nucleotides can either lead to subtle changes in the binding affinities of a given transcription factor or convert a binding sequence for a constitutive factor to a site recognized by an inducible factor. (Abstract shortened by UMI.) ^

Mitotic mapping of a novel tumor suppressor locus on human chromosome 3q important in osteosarcoma tumorigenesis

Tumor-specific loss of constitutional heterozygosity by deletion, mitotic recombination or nondisjunction is a common mechanism for tumor suppressor allele inactivation. When loss of heterozygosity is the result of mitotic recombination, or a segmental deletion event, only a portion of the chromosome is lost. This information can be used to map the location of new tumor suppressor genes. In osteosarcoma, the highest frequencies of loss of heterozygosity have been reported for chromosomes 3q, 13q, 17p. On chromosomes 13q and 17p, allelic losses are associated with loss of function at the retinoblastoma susceptibility locus (RB1) and the p53 locus, respectively. Chromosome 3q is also of particular interest because the high percent of loss of heterozygosity (62%-75%) suggests the presence of another tumor suppressor important for osteosarcoma tumorigenesis. To localize this putative tumor suppressor gene, we used polymorphic markers on chromosome 3q to find the smallest common region of allele loss. This putative tumor suppressor was localized to a 700 kb region on chromosome 3q26.2 between the polymorphic loci D3S1282 and D3S1246. ^

{\it Hox\/} genes and specification of regional identity in the axial skeleton: Analysis of the individual and combined mutant phenotypes of the paralogous group 4 murine {\it Hox\/} genes; {\it hoxa\/}-4, {\it hoxb\/}-4 and {\it hoxd\/}-4

The Hox gene products are transcription factors involved in specifying regional identity along the anteroposterior body axis. In Drosophila, where these genes are known as HOM-C (Homeotic-complex) genes and where they have been most extensively studied, they are expressed in restricted domains along the anteroposterior axis with different anterior limits. Genetic analysis of a large number of gain- and loss-of-function alleles of these genes has revealed that these genes are important in specifying segmental identity at their anterior limits of expression. Furthermore, there is a functional dominance of posterior genes over anterior genes, such that posterior genes can dominantly specify their developmental programs in spite of the expression of more anterior genes in the same segment. In the mouse, there are four clusters of HOM-C genes, called Hox genes. Thus, there may be up to four genes, called paralogs, that are more highly homologous to each other and to their Drosophila homolog than they are to the other mouse Hox genes. The single mutants for two paralogous genes, hoxa-4 and hoxd-4, presented in this dissertation, are similar to several other mouse Hox mutants in that they show partial, incompletely penetrant homeotic transformations of vertebrae at their anterior limit of expression. These mutants were then bred with hoxb-4 mutants (Ramirez-Solis, et al. 1993) to generate the three possible double mutant combinations as well as the triple mutant. The skeletal phenotypes of these group 4 Hox compound mutants displayed clear alterations in regional identity, such that a nearly complete transformation towards the morphology of the first cervical vertebra occurs. These results suggest a certain degree of functional redundancy among paralogous genes in specifying regional identity. Furthermore, there was a remarkable dose-dependent increase in the number of vertebrae transformed to a first cervical vertebra identity, including the second through the fifth cervical vertebrae in the triple mutant. Thus, these genes are required in a larger anteroposterior domain than is revealed by the single mutant phenotypes alone, such that multiple mutations in these genes result in transformations of vertebrae that are not at their anterior limit of expression. ^

COMPARISON OF DNA TRANSFECTION AND MICROCELL-MEDIATED CHROMOSOME TRANSFER FOR DETECTING TRANSFORMING GENES IN HUMAN TUMOR CELL LINES

DNA mediated gene transfection is an important tool for moving and isolating genes from one cell type and putting them into a foreign genetic background. DNA transfection studies have been done routinely in many laboratories to identify and isolate transforming sequences in human tumors and tumor cell lines. A second technique, microcell-mediated chromosome transfer, allows the transfer of small numbers of intact human chromosome from one cell to another. This work was done to compare the efficiency of these two techniques in the transformation of NIH 3T3 mouse fibroblast cells.^ My intent in comparing these two techniques was to see if there was a difference in the transforming capability of DNA which has been purified of all associated protein and RNAs, and that of DNA which is introduced into a cell in its native form, the chromosome. If chromosomal sequences were capable of transforming the 3T3 cells in culture, the method could then be used as a way to isolate the relevant tumorigenic chromosomes from human tumors.^ The study shows, however, that even for those cell lines that contain transforming sequences identified by DNA-mediated gene transfer, those same sequences were unable to transform 3T3 cells when introduced to the cells by somatic fusion of human tumor microcells. I believe that the human transforming sequences in their original genetic conformation are not recognized by the mouse cell as genes which should be expressed; therefore, no noticeable transformation event was selected by this technique. ^

Additive effects of HLA alleles and innate immune genes determine viral outcome in HCV infection

BACKGROUND Chronic HCV infection is a leading cause of liver-related morbidity globally. The innate and adaptive immune responses are thought to be important in determining viral outcomes. Polymorphisms associated with the IFNL3 (IL28B) gene are strongly associated with spontaneous clearance and treatment outcomes. OBJECTIVE This study investigates the importance of HLA genes in the context of genetic variation associated with the innate immune genes IFNL3 and KIR2DS3. DESIGN We assess the collective influence of HLA and innate immune genes on viral outcomes in an Irish cohort of women (n=319) who had been infected from a single source as well as a more heterogeneous cohort (Swiss Cohort, n=461). In the Irish cohort, a number of HLA alleles are associated with different outcomes, and the impact of IFNL3-linked polymorphisms is profound. RESULTS Logistic regression was performed on data from the Irish cohort, and indicates that the HLA-A*03 (OR 0.36 (0.15 to 0.89), p=0.027) -B*27 (OR 0.12 (0.03 to 0.45), p=<0.001), -DRB1*01:01 (OR 0.2 (0.07 to 0.61), p=0.005), -DRB1*04:01 (OR 0.31 (0.12 to 0.85, p=0.02) and the CC IFNL3 rs12979860 genotypes (OR 0.1 (0.04 to 0.23), p<0.001) are significantly associated with viral clearance. Furthermore, DQB1*02:01 (OR 4.2 (2.04 to 8.66), p=0.008), KIR2DS3 (OR 4.36 (1.62 to 11.74), p=0.004) and the rs12979860 IFNL3 'T' allele are associated with chronic infection. This study finds no interactive effect between IFNL3 and these Class I and II alleles in relation to viral clearance. There is a clear additive effect, however. Data from the Swiss cohort also confirms independent and additive effects of HLA Class I, II and IFNL3 genes in their prediction of viral outcome. CONCLUSIONS This data supports a critical role for the adaptive immune response in the control of HCV in concert with the innate immune response.

A novel universal DNA labeling and amplification system for rapid microarray-based detection of 117 antibiotic resistance genes in Gram-positive bacteria.

A rapid and simple DNA labeling system has been developed for disposable microarrays and has been validated for the detection of 117 antibiotic resistance genes abundant in Gram-positive bacteria. The DNA was fragmented and amplified using phi-29 polymerase and random primers with linkers. Labeling and further amplification were then performed by classic PCR amplification using biotinylated primers specific for the linkers. The microarray developed by Perreten et al. (Perreten, V., Vorlet-Fawer, L., Slickers, P., Ehricht, R., Kuhnert, P., Frey, J., 2005. Microarray-based detection of 90 antibiotic resistance genes of gram-positive bacteria. J.Clin.Microbiol. 43, 2291-2302.) was improved by additional oligonucleotides. A total of 244 oligonucleotides (26 to 37 nucleotide length and with similar melting temperatures) were spotted on the microarray, including genes conferring resistance to clinically important antibiotic classes like β-lactams, macrolides, aminoglycosides, glycopeptides and tetracyclines. Each antibiotic resistance gene is represented by at least 2 oligonucleotides designed from consensus sequences of gene families. The specificity of the oligonucleotides and the quality of the amplification and labeling were verified by analysis of a collection of 65 strains belonging to 24 species. Association between genotype and phenotype was verified for 6 antibiotics using 77 Staphylococcus strains belonging to different species and revealed 95% test specificity and a 93% predictive value of a positive test. The DNA labeling and amplification is independent of the species and of the target genes and could be used for different types of microarrays. This system has also the advantage to detect several genes within one bacterium at once, like in Staphylococcus aureus strain BM3318, in which up to 15 genes were detected. This new microarray-based detection system offers a large potential for applications in clinical diagnostic, basic research, food safety and surveillance programs for antimicrobial resistance.

Human NR5A1/SF-1 mutations show decreased activity on BDNF (brain-derived neurotrophic factor), an important regulator of energy balance: testing impact of novel SF-1 mutations beyond steroidogenesis

CONTEXT Human NR5A1/SF-1 mutations cause 46,XY disorder of sex development (DSD) with broad phenotypic variability, and rarely cause adrenal insufficiency although SF-1 is an important transcription factor for many genes involved in steroidogenesis. In addition, the Sf-1 knockout mouse develops obesity with age. Obesity might be mediated through Sf-1 regulating activity of brain-derived neurotrophic factor (BDNF), an important regulator of energy balance in the ventromedial hypothalamus. OBJECTIVE To characterize novel SF-1 gene variants in 4 families, clinical, genetic and functional studies were performed with respect to steroidogenesis and energy balance. PATIENTS 5 patients with 46,XY DSD were found to harbor NR5A1/SF-1 mutations including 2 novel variations. One patient harboring a novel mutation also suffered from adrenal insufficiency. METHODS SF-1 mutations were studied in cell systems (HEK293, JEG3) for impact on transcription of genes involved in steroidogenesis (CYP11A1, CYP17A1, HSD3B2) and in energy balance (BDNF). BDNF regulation by SF-1 was studied by promoter assays (JEG3). RESULTS Two novel NR5A1/SF-1 mutations (Glu7Stop, His408Profs*159) were confirmed. Glu7Stop is the 4th reported SF-1 mutation causing DSD and adrenal insufficiency. In vitro studies revealed that transcription of the BDNF gene is regulated by SF-1, and that mutant SF-1 decreased BDNF promoter activation (similar to steroid enzyme promoters). However, clinical data from 16 subjects carrying SF-1 mutations showed normal birth weight and BMI. CONCLUSIONS Glu7Stop and His408Profs*159 are novel SF-1 mutations identified in patients with 46,XY DSD and adrenal insufficiency (Glu7Stop). In vitro, SF-1 mutations affect not only steroidogenesis but also transcription of BDNF which is involved in energy balance. However, in contrast to mice, consequences on weight were not found in humans with SF-1 mutations.

Extensive Copy-Number Variation of Young Genes across Stickleback Populations

Duplicate genes emerge as copy-number variations (CNVs) at the population level, and remain copy-number polymorphic until they are fixed or lost. The successful establishment of such structural polymorphisms in the genome plays an important role in evolution by promoting genetic diversity, complexity and innovation. To characterize the early evolutionary stages of duplicate genes and their potential adaptive benefits, we combine comparative genomics with population genomics analyses to evaluate the distribution and impact of CNVs across natural populations of an eco-genomic model, the three-spined stickleback. With whole genome sequences of 66 individuals from populations inhabiting three distinct habitats, we find that CNVs generally occur at low frequencies and are often only found in one of the 11 populations surveyed. A subset of CNVs, however, displays copy-number differentiation between populations, showing elevated within-population frequencies consistent with local adaptation. By comparing teleost genomes to identify lineage-specific genes and duplications in sticklebacks, we highlight rampant gene content differences among individuals in which over 30% of young duplicate genes are CNVs. These CNV genes are evolving rapidly at the molecular level and are enriched with functional categories associated with environmental interactions, depicting the dynamic early copy-number polymorphic stage of genes during population differentiation.

The transporter GAT1 plays an important role in GABA-mediated carbon-nitrogen interactions in Arabidopsis

Glutamate derived γ-aminobutyric acid (GABA) is synthetized in the cytosol prior to delivery to the mitochondria where it is catabolized via the TCA cycle. GABA accumulates under various environmental conditions, but an increasing number of studies show its involvement at the crossroad between C and N metabolism. To assess the role of GABA in modulating cellular metabolism, we exposed seedlings of A. thaliana GABA transporter gat1 mutant to full nutrition medium and media deficient in C and N combined with feeding of different concentrations (0.5 and 1 mM) of exogenous GABA. GC-MS based metabolite profiling showed an expected effect of medium composition on the seedlings metabolism of mutant and wild type alike. That being said, a significant interaction between GAT1 deficiency and medium composition was determined with respect to magnitude of change in relative amino acid levels. The effect of exogenous GABA treatment on metabolism was contingent on both the medium and the genotype, leading for instance to a drop in asparagine under full nutrition and low C conditions and glucose under all tested media, but not to changes in GABA content. We additionally assessed the effect of GAT1 deficiency on the expression of glutamate metabolism related genes and genes involved in abiotic stress responses. These results suggest a role for GAT1 in GABA-mediated metabolic alterations in the context of the C-N equilibrium of plant cells.