Enriching targeted sequencing experiments for rare disease alleles

MOTIVATION: Next-generation targeted resequencing of genome-wide association study (GWAS)-associated genomic regions is a common approach for follow-up of indirect association of common alleles. However, it is prohibitively expensive to sequence all the samples from a well-powered GWAS study with sufficient depth of coverage to accurately call rare genotypes. As a result, many studies may use next-generation sequencing for single nucleotide polymorphism (SNP) discovery in a smaller number of samples, with the intent to genotype candidate SNPs with rare alleles captured by resequencing. This approach is reasonable, but may be inefficient for rare alleles if samples are not carefully selected for the resequencing experiment. RESULTS: We have developed a probability-based approach, SampleSeq, to select samples for a targeted resequencing experiment that increases the yield of rare disease alleles substantially over random sampling of cases or controls or sampling based on genotypes at associated SNPs from GWAS data. This technique allows for smaller sample sizes for resequencing experiments, or allows the capture of rarer risk alleles. When following up multiple regions, SampleSeq selects subjects with an even representation of all the regions. SampleSeq also can be used to calculate the sample size needed for the resequencing to increase the chance of successful capture of rare alleles of desired frequencies. SOFTWARE: http://biostat.mc.vanderbilt.edu/SampleSeq     

Molecular identification and physiological characterization of a novel monosaccharide transporter from Arabidopsis involved in vacuolar sugar transport

The tonoplast monosaccharide transporter (TMT) family comprises three isoforms in Arabidopsis thaliana, and TMT-green fluorescent protein fusion proteins are targeted to the vacuolar membrane. TMT promoter-beta-glucuronidase plants revealed that the TONOPLAST MONOSACCHARIDE TRANSPORTER1 (TMT1) and TMT2 genes exhibit a tissue- and cell type-specific expression pattern, whereas TMT3 is only weakly expressed. TMT1 and TMT2 expression is induced by drought, salt, and cold treatments and by sugar. During cold adaptation, tmt knockout lines accumulated less glucose and fructose compared with wild-type plants, whereas no differences were observed for sucrose. Cold adaptation of wild-type plants substantially promoted glucose uptake into isolated leaf mesophyll vacuoles. Glucose uptake into isolated vacuoles was inhibited by NH(4)(+), fructose, and phlorizin, indicating that transport is energy-dependent and that both glucose and fructose were taken up by the same carrier. Glucose import into vacuoles from two cold-induced tmt1 knockout lines or from triple knockout plants was substantially lower than into corresponding wild-type vacuoles. Monosaccharide feeding into leaf discs revealed the strongest response to sugar in tmt1 knockout lines compared with wild-type plants, suggesting that TMT1 is required for cytosolic glucose homeostasis. Our results indicate that TMT1 is involved in vacuolar monosaccharide transport and plays a major role during stress responses.     

Molecular characterization and comparative analysis of two waxy alleles in barley

Two alleles of the barley waxy locus were characterized from non-waxy cultivar Bowman and waxy cultivar CDC Candle, respectively. Their nucleotide and protein sequences were compared with other known waxy genes. The comparison results indicated that there were 100 polymorphic sites, among which 69 were in the non-coding region and 31 were in the coding region. Out of 100 polymorphic sites, 45 were trans-version, 35 were transition and 20 were indels. A 397 bp deletion and a 193 bp insertion in the promoter region and a 15 bp insertion in the coding region were found in CDC Candle, but not in Bowman. A deletion (11 bp) was detected in Bowman, which exhibited no effects on normal waxy expression. In summary, the 397 bp deletion was supposed to account for the reduction of GBSS I, resulting in the low amylose in CDC Candle; whereas other polymorphic sites might be not correlated with amylose synthesis.     

Molecular characterization of organelle-type Nudix hydrolases in Arabidopsis

Nudix (for nucleoside diphosphates linked to some moiety X) hydrolases act to hydrolyze ribonucleoside and deoxyribonucleoside triphosphates, nucleotide sugars, coenzymes, or dinucleoside polyphosphates. Arabidopsis (Arabidopsis thaliana) contains 27 genes encoding Nudix hydrolase homologues (AtNUDX1 to -27) with a predicted distribution in the cytosol, mitochondria, and chloroplasts. Previously, cytosolic Nudix hydrolases (AtNUDX1 to -11 and -25) were characterized. Here, we conducted a characterization of organelle-type AtNUDX proteins (AtNUDX12 to -24, -26, and -27). AtNUDX14 showed pyrophosphohydrolase activity toward both ADP-ribose and ADP-glucose, although its K(m) value was approximately 100-fold lower for ADP-ribose (13.0+/-0.7 microm) than for ADP-glucose (1,235+/-65 microm). AtNUDX15 hydrolyzed not only reduced coenzyme A (118.7+/-3.4 microm) but also a wide range of its derivatives. AtNUDX19 showed pyrophosphohydrolase activity toward both NADH (335.3+/-5.4 microm) and NADPH (36.9+/-3.5 microm). AtNUDX23 had flavin adenine dinucleotide pyrophosphohydrolase activity (9.1+/-0.9 microm). Both AtNUDX26 and AtNUDX27 hydrolyzed diadenosine polyphosphates (n=4-5). A confocal microscopic analysis using a green fluorescent protein fusion protein showed that AtNUDX15 is distributed in mitochondria and AtNUDX14 -19, -23, -26, and -27 are distributed in chloroplasts. These AtNUDX mRNAs were detected ubiquitously in various Arabidopsis tissues. The T-DNA insertion mutants of AtNUDX13, -14, -15, -19, -20, -21, -25, -26, and -27 did not exhibit any phenotypical differences under normal growth conditions. These results suggest that Nudix hydrolases in Arabidopsis control a variety of metabolites and are pertinent to a wide range of physiological processes.     

Molecular characterization and targeted quantitative profiling of the sphingolipidome in rice

Recent advances in comprehensive metabolite profiling techniques, the foundation of metabolomics, is facilitating our understanding of the functions, regulation and complex networks of various metabolites in organisms. Here, we report a quantitative metabolomics technique for complex plant sphingolipids, composed of various polar head groups as well as structural isomers of hydrophobic ceramide moieties. Rice (Oryza sativa L.) was used as an experimental model of monocotyledonous plants and has been demonstrated to possess a highly complex sphingolipidome including hundreds of molecular species with a wide range of abundance. We established a high‐throughput scheme for lipid preparation and mass spectrometry‐based characterization of complex sphingolipid structures, which provided basic information to create a comprehensive theoretical library for targeted quantitative profiling of complex sphingolipids in rice. The established sphingolipidomic approach combined with multivariate analyses of the large dataset obtained clearly showed that different classes of rice sphingolipids, particularly including subclasses of glycosylinositol phosphoceramide with various sugar‐chain head groups, are distributed with distinct quantitative profiles in various rice tissues, indicating tissue‐dependent metabolism and biological functions of the lipid classes and subclasses. The sphingolipidomic analysis also highlighted that disruption of a lipid‐associated gene causes a typical sphingolipidomic change in a gene‐dependent manner. These results clearly support the utility of the sphingolipidomic approach in application to wide screening of sphingolipid‐metabolic phenotypes as well as deeper investigation of metabolism and biological functions of complex sphingolipid species in plants.     

Molecular and functional characterization of Arabidopsis Cullin 3A

Cullin proteins, which belong to multigenic families in all eukaryotes, associate with other proteins to form ubiquitin protein ligases (E3s) that target substrates for proteolysis by the 26S proteasome. Here, we present the molecular and genetic characterization of a plant Cullin3. In contrast to fungi and animals, the genome of the model plant Arabidopsis thaliana contains two related CUL3 genes, called CUL3A and CUL3B. We found that CUL3A is ubiquitously expressed in plants and is able to interact with the ring-finger protein RBX1. A genomic search revealed the existence of at least 76 BTB-domain proteins in Arabidopsis belonging to 11 major families. Yeast two-hybrid experiments indicate that representative members of certain families are able to physically interact with both CUL3A and CUL3B, suggesting that Arabidopsis CUL3 forms E3 protein complexes with certain BTB domain proteins. In order to determine the function of CUL3A, we used a reverse genetic approach. The cul3a null mutant flowers slightly later than the control plants. Furthermore, this mutant exhibits a reduced sensitivity of the inhibition of hypocotyl growth in far-red light and miss-expresses COP1. The viability of the mutant plants suggests functional redundancy between the two CUL3 genes in Arabidopsis.     

Molecular and physiological analysis of Arabidopsis mutants exhibiting altered sensitivities to aluminium

t Arabidopsis mutants with altered Al sensitivities have been isolated with the goal of identifying genes important for Al resistance and toxicity. By screening in a high-Al environment, seven Al-resistant mutants (alr) were isolated. The alr mutants were semi-dominant and constitute at least two unique loci in Arabidopsis. Nine Al-sensitive mutants (als) were also isolated. Complementation analysis revealed that of the nine als mutants, eight represent unique loci, indicating that Al sensitivity is genetically complex in Arabidopsis. The characterization of mutants with altered Al sensitivity will provide greater insights into the mechanisms of Al toxicity and resistance on a molecular and biochemical level.     

Synergistic derepression of gibberellin signaling by removing RGA and GAI function in Arabidopsis thaliana

RGA and GAI are negative regulators of the gibberellin (GA) signal transduction pathway in Arabidopsis thaliana. These genes may have partially redundant functions because they are highly homologous, and plants containing single null mutations at these loci are phenotypically similar to wild type. Previously, rga loss-of-function mutations were shown to partially suppress defects of the GA-deficient ga1-3 mutant. Phenotypes rescued include abaxial trichome initiation, rosette radius, flowering time, stem elongation, and apical dominance. Here we present work showing that the rga-24 and gai-t6 null mutations have a synergistic effect on plant growth. Although gai-t6 alone has little effect, when combined with rga-24, they completely rescued the above defects of ga1-3 to wild-type or GA-overdose phenotype. However, seed germination and flower development defects were not restored. Additionally, rga-24 and rga-24/gai-t6 but not gai-t6 alone caused increased feedback inhibition of expression of a GA biosynthetic gene in both the ga1-3 and wild-type backgrounds. These results demonstrate that RGA and GAI have partially redundant functions in maintaining the repressive state of the GA-signaling pathway, but RGA plays a more dominant role than GAI. Removing both RGA and GAI function allows for complete derepression of many aspects of GA signaling.     

Molecular and physiological analysis of Arabidopsis mutants defective in cytosolic or chloroplastic aspartate aminotransferase

A single, recessive mutation in soybean (Glycine max L. Merr.), which confers a seed phenotype of increased inorganic phosphate, decreased phytic acid, and a decrease in total raffinosaccharides, has been previously disclosed (S.A. Sebastian, P.S. Kerr, R.W. Pearlstein, W.D. Hitz [2000] Soy in Animal Nutrition, pp 56-74). The genetic lesion causing the multiple changes in seed phenotype is a single base change in the third base of the codon for what is amino acid residue 396 of the mature peptide encoding a seed-expressed myo-inositol 1-phospate synthase gene. The base change causes residue 396 to change from lysine to asparagine. That amino acid change decreases the specific activity of the seed-expressed myo-inositol 1-phosphate synthase by about 90%. Radio tracer experiments indicate that the supply of myo-inositol to the reaction, which converts UDP-galactose and myo-inositol to galactinol is a controlling factor in the conversion of total carbohydrate into the raffinosaccharides in both wild-type and mutant lines. That same decrease in myo-inositol 1-phosphate synthetic capacity leads to a decreased capacity for the synthesis of myo-inositol hexaphosphate (phytic acid) and a concomitant increase in inorganic phosphate.     

Molecular and physiological characterization of the pyrokinin insect neuropeptide family

Peptides from the pyrokinin (PK) family are a large, structurally and functionally diverse group of the insect neuropeptides produced by neurosecretory cells of the insect nervous system. This family contains short and long peptides which share C-terminal -FXPRLa amino acid sequence. Pyrokinins regulate the visceral muscle contractions, pheromone biosynthesis, pupariation and diapause duration in insects. They are encoded by two genes PBAN and capa, which are mainly expressed in the suboesophageal ganglion. Peptides are then transported to the retrocerebral complex and released into haemolymph. Recent studies are focused on application of pyrokinins as biopesticides in the regulation of insect pests growth and development.     

Comparative proteomics and physiological characterization of Arabidopsis thaliana seedlings in responses to Ochratoxin A

Ochratoxin A (OTA) is a mycotoxin that is primarily produced by Aspergillus ochraceus and Penicillium verrucosum. This mycotoxin is a contaminant of food and feedstock worldwide and may induce cell death in plants. To investigate the dynamic growth process of Arabidopsis seedlings in response to OTA stress and to obtain a better understanding of the mechanism of OTA toxicity towards Arabidopsis, a comparative proteomics study using 2-DE and MALDI-TOF/TOF MS/MS was performed. Mass spectrometry analysis identified 59 and 51 differentially expressed proteins in seedlings exposed to 25 and 45 μM OTA for 7 days, respectively. OTA treatment decreased root elongation and leaf area, increased anthocyanin accumulation, damaged the photosynthetic apparatus and inhibited photosynthesis. Treatment of the seedlings with 25 μM OTA enhanced energy metabolism, whereas higher concentration of OTA (45 μM) inhibited energy metabolism in the seedlings. OTA treatment caused an increase of ROS, an enhancement of antioxidant enzyme defense responses, disturbance of redox homeostasis and activation of lipid oxidation. Glutamine and S-adenosylmethionine metabolism may also play important roles in the response to OTA. In conclusion, our study provided novel insights regarding the response of Arabidopsis to OTA at the level of the proteome. These results are expected to be highly useful for understanding the physiological responses and dissecting the OTA response pathways in higher plants.     

The number of alleles in electrophoretic experiments

We examine the rate of divergence of the number of distinguishable alleles in a large sample. Our calculations are based on the Ohta-Kimura model for electrophoretic experiments, in which types can only take discrete values, and mutations are selectively neutral.     

Molecular characterization of plastidic phosphoserine aminotransferase in serine biosynthesis from Arabidopsis

Serine biosynthesis in plants proceeds by two pathways; a photorespiratory pathway which is associated with photorespiration and a pathway from phosphoglycerate. A cDNA encoding plastidic phosphoserine aminotransferase (PSAT) which catalyzes the formation of phosphoserine from phosphohydroxypyruvate has been isolated from Arabidopsis thaliana . Genomic DNA blot analysis indicated that this enzyme is most probably encoded by a single gene and is mapped on the lower arm of chromosome 4. The deduced protein contains an N-terminal extension exhibiting the general features of a plastidic transit peptide, which was confirmed by subcellular organelle localization using GFP (green flourescence protein). Northern analysis indicated preferential expression of PSAT in roots of light-grown plants, supporting the idea that the phosphorylated pathway may play an important role in supplying the serine requirement of plants in non-green tissues. In situ hybridization analysis of PSAT revealed that the gene is generally expressed in all types of cells with a significantly higher amount in the meristem tissue of root tips.     

Molecular and biochemical characterization of a serine racemase from Arabidopsis thaliana

A cDNA encoding a homolog of mammalian serine racemase, a unique enzyme in eukaryotes, was isolated from Arabidopsis thaliana and expressed in Escherichia coli cells. The gene product, of which the amino acid residues for binding pyridoxal 5'-phosphate (PLP) are conserved in this as well as mammalian serine racemases, catalyzes not only serine racemization but also dehydration of serine to pyruvate. The enzyme is a homodimer and requires PLP and divalent cations, Ca2+, Mg2+, Mn2+, Fe2+, or Ni2+, at alkaline pH for both activities. The racemization process is highly specific toward L-serine, whereas L-alanine, L-arginine, and L-glutamine were poor substrates. The Vmax/Km values for racemase activity of L- and D-serine are 2.0 and 1.4 nmol/mg/min/mM, respectively, and those values for L- and D-serine on dehydratase activity are 13 and 5.3 nmol/mg/min/mM, i.e. consistent with the theory of racemization reaction and the specificity of dehydration toward L-serine. Hybridization analysis showed that the serine racemase gene was expressed in various organs of A. thaliana.     

Molecular characterization of four chitinase cDNAs obtained fromCladosporium fulvum-infected tomato

Complementary DNA clones encoding acidic and basic isoforms of tomato chitinases were isolated fromCladosporium fulvum-infected leaves. The clones were sequenced and found to encode the 30 kDa basic intracellular and the 26 and 27 kDa acidic extracellular tomato chitinases previously purified (M.H.A.J. Joostenet al., in preparation). A fourth truncated cDNA which appears to encode an extracellular chitinase with 82% amino acid similarity to the 30 kDa intracellular chitinase was also isolated. Characterization of the clones revealed that the 30 kDa basic intracellular protein is a class I chitinase and that the 26 and 27 kDa acidic extracellular proteins which have 85% peptide sequence similarity are class II chitinases. The characterized cDNA clones represent four from a family of at least six tomato chitinases. Southern blot analysis indicated that, with the exception of the 30 kDa basic intracellular chitinase, the tomato chitinases are encoded by one or two genes. Northern blot analysis showed that the mRNA encoding the 26 kDa acidic extracellular chitinase is induced more rapidly during an incompatibleC. fulvum-tomato interaction than during a compatible interaction. This difference in timing of mRNA induction was not observed for the 30 kDa basic intracellular chitinase.