Comparison of commercially available target enrichment methods for next-generation sequencing.

Isolating high-priority segments of genomes greatly enhances the efficiency of next-generation sequencing (NGS) by allowing researchers to focus on their regions of interest. For the 2010-11 DNA Sequencing Research Group (DSRG) study, we compared outcomes from two leading companies, Agilent Technologies (Santa Clara, CA, USA) and Roche NimbleGen (Madison, WI, USA), which offer custom-targeted genomic enrichment methods. Both companies were provided with the same genomic sample and challenged to capture identical genomic locations for DNA NGS. The target region totaled 3.5 Mb and included 31 individual genes and a 2-Mb contiguous interval. Each company was asked to design its best assay, perform the capture in replicates, and return the captured material to the DSRG-participating laboratories. Sequencing was performed in two different laboratories on Genome Analyzer IIx systems (Illumina, San Diego, CA, USA). Sequencing data were analyzed for sensitivity, specificity, and coverage of the desired regions. The success of the enrichment was highly dependent on the design of the capture probes. Overall, coverage variability was higher for the Agilent samples. As variant discovery is the ultimate goal for a typical targeted sequencing project, we compared samples for their ability to sequence single-nucleotide polymorphisms (SNPs) as a test of the ability to capture both chromosomes from the sample. In the targeted regions, we detected 2546 SNPs with the NimbleGen samples and 2071 with Agilent's. When limited to the regions that both companies included as baits, the number of SNPs was ∼1000 for each, with Agilent and NimbleGen finding a small number of unique SNPs not found by the other.

Molecular mechanisms of pituitary organogenesis: In search of novel regulatory genes.

Defects in pituitary gland organogenesis are sometimes associated with congenital anomalies that affect head development. Lesions in transcription factors and signaling pathways explain some of these developmental syndromes. Basic research studies, including the characterization of genetically engineered mice, provide a mechanistic framework for understanding how mutations create the clinical characteristics observed in patients. Defects in BMP, WNT, Notch, and FGF signaling pathways affect induction and growth of the pituitary primordium and other organ systems partly by altering the balance between signaling pathways. The PITX and LHX transcription factor families influence pituitary and head development and are clinically relevant. A few later-acting transcription factors have pituitary-specific effects, including PROP1, POU1F1 (PIT1), and TPIT (TBX19), while others, such as NeuroD1 and NR5A1 (SF1), are syndromic, influencing development of other endocrine organs. We conducted a survey of genes transcribed in developing mouse pituitary to find candidates for cases of pituitary hormone deficiency of unknown etiology. We identified numerous transcription factors that are members of gene families with roles in syndromic or non-syndromic pituitary hormone deficiency. This collection is a rich source for future basic and clinical studies.

Genomic characterization of human SEC14L1 splice variants within a 17q25 candidate tumor suppressor gene region and identification of an unrelated embedded expressed sequence tag.

Human SEC14L1 shows partial sequence homology to the budding yeast SEC14 protein and the Japanese flying squid retinal-binding protein and was previously generally localized to 17q25. We more precisely mapped SEC14L1 within a discrete region of 17q25 that likely harbors at least one putative breast and ovarian tumor suppressor gene. We determined that this gene consists of 18 exons ranging in size from 70 bp (exon 11) to 3088 bp (exon 17) and spanning at least 58 kb of DNA. Exon 17 contained a highly polymorphic variable number of tandem repeats (VNTR) and was present only in the larger ubiquitously expressed 5.5-kb transcript. The 3.0-kb ubiquitously expressed transcript included sequences at the beginning of exon 17 (designated exon 17a) and the end of exon 17 (designated exon 18), but lacked the internal 2439 bp of exon 17, including the VNTR. This alternative splicing resulted in a predicted protein of 719 residues from the smaller transcript with four more terminal amino acids than the 715 residue protein predicted from the larger transcript. EST H49244 spanned exon 11 of SEC14L1 and was specifically expressed in human peripheral blood leukocytes. One intragenic single nucleotide polymorphism (SNP) was confirmed. SEC14L1 contained the CRAL/TRIO domain also found in alpha-tocopherol transfer protein (TTPA) and cellular retinaldehyde-binding protein (CRALBP). As retinoids have been shown to inhibit the growth of breast cancer cells, loss of the proposed SEC14L1 retinal-binding function may contribute to breast tumorigenesis. As TTPA and CRALBP have been implicated in retinitis pigmentosa (RP), altered SEC14L1 expression may contribute to RP in previously unlinked families. Coding exon-specific PCR primers were designed to aid in future expression and mutational analyses.

Identification of members of the Wnt signaling pathway in the embryonic pituitary gland.

Prop1 is one of several transcription factors important for the development of the pituitary gland. Downstream targets of PROP1 and other critical pituitary transcription factors remain largely unknown. We have generated a partial expression profile of the developing pituitary gland containing over 350 transcripts, using cDNA subtractive hybridization between Prop1(df/df) and wild-type embryonic pituitary gland primordia. Numerous classes of genes including transcription factors, membrane associated molecules, and cell cycle regulators were identified in this study. Of the transcripts, 34% do not have sequence similarity to known genes, but are similar to ESTs, and 4% represent novel sequences. Pituitary gland expression of a number of clones was verified using in situ hybridization. Several members of the Wnt signaling pathway were identified in the developing pituitary gland. The frizzled2 receptor, Apc, beta-catenin, groucho, and a novel isoform of TCF4 (officially named Tcf7l2) were identified in developing pituitary libraries. Three N-terminal alternatively spliced Tcf7l2 isoforms are reported here, each of which lacks a DNA-binding domain. Functional studies indicate that these isoforms can act as endogenous inhibitors of Wnt signaling in some contexts. This is the first report of Tcf7l2 and Fzd2 expression in the developing pituitary. These molecules may be important in mediating Wnt signaling during pituitary ontogeny. We expect other transcripts from these libraries to be involved in pituitary gland development.

Correction of deafness in shaker-2 mice by an unconventional myosin in a BAC transgene.

The shaker-2 mouse mutation, the homolog of human DFNB3, causes deafness and circling behavior. A bacterial artificial chromosome (BAC) transgene from the shaker-2 critical region corrected the vestibular defects, deafness, and inner ear morphology of shaker-2 mice. An unconventional myosin gene, Myo15, was discovered by DNA sequencing of this BAC. Shaker-2 mice were found to have an amino acid substitution at a highly conserved position within the motor domain of this myosin. Auditory hair cells of shaker-2 mice have very short stereocilia and a long actin-containing protrusion extending from their basal end. This histopathology suggests that Myo15 is necessary for actin organization in the hair cells of the cochlea.

Signal transduction by DR3, a death domain-containing receptor related to TNFR-1 and CD95.

Tumor necrosis factor receptor-1 (TNFR-1) and CD95 (also called Fas or APO-1) are cytokine receptors that engage the apoptosis pathway through a region of intracellular homology, designated the "death domain." Another death domain-containing member of the TNFR family, death receptor 3 (DR3), was identified and was shown to induce both apoptosis and activation of nuclear factor kappaB. Expression of DR3 appears to be restricted to tissues enriched in lymphocytes. DR3 signal transduction is mediated by a complex of intracellular signaling molecules including TRADD, TRAF2, FADD, and FLICE. Thus, DR3 likely plays a role in regulating lymphocyte homeostasis.

Mouse Brca1: localization sequence analysis and identification of evolutionarily conserved domains.

The human genes BRCA1, conferring susceptibility to early-onset breast and ovarian cancer, has recently been isolated. Here we describe isolation of cDNAs, sequence analysis, and genomic localization of the murine homolog, Brac1. The mouse cDNA sequence predicts a protein of 1812 amino acids; a number of small gaps account for the 51 fewer residues in the mouse protein relative to human BRCA1. While the predicted mouse and human proteins display on the whole a high level of homology (58% identity, 73% similarity), the regions of greatest homology are at the respective amino and carboxyl termini. Most reported disease-associated missense mutations in human BCRA1 occurred within these more highly conserved terminal regions. A predicted zinc-building RING finger domain near the amino terminus lies within a 50 amino acid stretch that is perfectly conserved in both species. The strong conservation during mammalian evolution argues for the importance of this domain, perhaps mediating a role for BRCA1 in DNA and/or protein binding. We have also identified a conserved highly acidic domain in the carboxyl terminal half of the BCRA1 protein resembling acidic transactivation domains of certain transcription factors. Using an interspecific backcross panel, Brca1 was mapped to a region of mouse chromosome 11 that exhibits conserved linkage with 17q21. The sequence and isolated cDNAs will provide useful reagents for studying the expression of Brca1 in the mouse, and for testing the importance of the evolutionarily conserved domains.

Evidence that the FK506-binding immunophilin heat shock protein 56 is required for trafficking of the glucocorticoid receptor from the cytoplasm to the nucleus.

The FK506-binding immunophilin hsp56 (FKBP52) is one of several chaperone proteins associated with untrasformed steroid receptors in a multiprotein heterocomplex. The function of heat shock protein 56 (hsp56) with respect to receptor action is unknown. hsp56 is not required for glucocorticoid receptor heterocomplex assembly or for proper folding of the receptor hormone-binding domain into a high affinity steroid-binding conformation. In intact cells, the majority of the hsp56 is located in the nucleus, with a minority colocalizing with microtubules in the cytoplasm. hsp56 contains a conserved negatively charged domain that we speculate might serve as a nuclear localization signal recognition sequence. Here we show that injection of an antibody raised against this negative sequence into intact L cells impedes subsequent dexamethasone-mediated shift of the glucocorticoid receptor into the nucleus. Nonimmune rabbit serum and an antibody raised against another site on hsp56 do not affect receptor movement. Inhibition of receptor movement by the 419 antibody against the negative sequence is blocked by preincubation with purified hsp56, but not by preincubation with purified hsp90, hsp70, or BSA. These observations are consistent with the possibility that hsp56 is involved in receptor trafficking to the nucleus, possibly functioning as the nuclear localization signal recognition protein. Receptor trafficking to the nucleus is not affected by FK506, indicating that the peptidylprolyl isomerase activity of hsp56 is not involved.

Understanding the dynamics of cellular responsiveness to modifications of metabolic substrates in perifusion.

A novel microperifusion system with capabilities for continuous, real-time, potentiometric monitoring of extracellular hydrogen ion concentration has been used to define the response of HeLa cells to abrupt changes in extracellular energy sources or introduction of an inhibitor of glycolysis. Glycolytic inhibition, induced by removal of glucose or introduction of iodoacetate, each led to a rapid, continuous decrease in acid release. The response to iodoacetate took longer than removal of glucose, perhaps due to the time required for binding and activation. Once inhibition began, however, the rate of change was greater than following glucose removal. Conversely, recovery time following iodoacetate inhibition was much slower than with glucose removal. Unlike the response to short-term glucose depletion, a second pulse of iodoacetate resulted in a faster response followed by an even longer recovery time. The response to switching between glucose and glutamine began almost without evident delay. The response patterns revealed that HeLa cells prefer glutamine to glucose, but, in the presence of both energy sources, some glucose continues to be used. In summary, these results indicate that continuous, real-time monitoring of the kinetics of hydrogen-ion release can be used to gain new insights into the dynamics of cellular response to perturbations of extracellular energy sources.

Serum-regulated nuclear localization is signal specific.

The activity of nuclear factors can be regulated by blocking their ability to enter the nucleus, but how the cell achieves this is not yet understood. We demonstrate herein the serum-responsive nuclear localization of adenovirus E1a protein and show that this serum dependence is a property of the nuclear localization signal itself. When E1a protein is microinjected into the cytoplasm of cultured cells, it is found in the nucleus 30 min later only if the cells are serum fed; in serum-starved (growth-arrested) cells, the E1a is still cytoplasmic. Substituting the simian virus-40 T-antigen nuclear localization signal in place of the normal E1a signal abolishes this serum effect, and transferring the E1a signal to a heterologous protein also transfers the serum dependence. The serum effect on signal function is first exerted within 40 min after serum addition, suggesting that this is one of the earliest cellular responses to serum feeding. We conclude that the nuclear accumulation (and probably the function) of a protein is influenced not only by the presence of a nuclear localization signal, but also by the nature of that signal.

Pentapeptide nuclear localization signal in adenovirus E1a.

The adenovirus E1a gene products are nuclear proteins important in transcriptional control of viral functions during infection. By producing normal E1a proteins and derivatives of E1a in bacteria and microinjecting these proteins into cultured cells, we were able to examine their ability to localize to the nucleus. We showed that a short peptide sequence at the carboxyl terminus of E1a is necessary for the rapid (30-min) nuclear localization of that protein. Additionally, we showed that just the last five amino acids of E1a are sufficient to direct nuclear accumulation of a heterologous protein, Escherichia coli galactokinase, with the same kinetics as native E1a. The mechanism by which this pentamer mediates rapid nuclear localization was examined by testing the ability of a galactokinase derivative which has no signal pentamer to exit the nucleus, as well as to enter it. Because neither free entry nor exit was detected, the effect of the signal is unlikely to be through increased nuclear retention of freely diffusible proteins but rather by enhancement of entry into the nucleus.

Requirement for the 3' flanking region of the bovine growth hormone gene for accurate polyadenylylation.

We examined whether the sequence extending 3' to the polyadenylylation site of the bovine growth hormone gene contains any signal that affects the polyadenylylation of the growth hormone mRNA. For this purpose, cloned copies of this gene, each containing a different length of growth hormone-specific sequence 3' to the wild-type polyadenylylation site, were used to transfect COS-1 cells. The polyadenylylation site on the mRNAs produced from the exogenously added growth hormone genes were analyzed with an S1 nuclease mapping procedure. We found that a gene containing 84 base pairs of its own 3' flanking sequence is capable of producing an accurately polyadenylylated mRNA. On the other hand, genes containing only 1, 10, or 13 base pairs of 3' flanking sequence were principally polyadenylylated at discrete sites either upstream or downstream from the wild-type position. Using a computer program, we examined whether secondary structures on the primary growth hormone transcript correlated with the site where the mRNA is polyadenylylated.

Characterization of the bovine prolactin gene.

The gene coding for bovine prolactin was shown to exist as a single copy per haploid genome. Three restriction fragment polymorphisms were detected in the prolactin gene by Southern blot analysis of DNA obtained from the semen of pedigreed bulls representing eight breeds. The organization of the bovine prolactin gene was determined by restriction mapping of a clone isolated from a genomic library and by genomic blots. The 5'-flanking region and two exons were sequenced and the transcription start site mapped by primer extension. Comparison of the bovine prolactin sequence reported here with the published sequence of the rat prolactin gene revealed extensive homology (79%), extending 360 nucleotides upstream from the cap site, after which the sequences diverge. The homology exceeds that of the coding regions. A possible alternate intron-exon splice site was noted within the sequence coding for the signal peptide.

Cloning and nucleotide sequencing of the bovine growth hormone gene.

A gene coding for bovine growth hormone was isolated from a bovine genomic library. The nucleotide sequence of the coding regions of the gene was found to be identical with that of a nearly full-length growth hormone cDNA clone. The gene sequence is approximately 1800 bp in length and contains four intervening sequences. The second intervening sequence of 227 nucleotides does not contain a repetitive element similar to that observed in the rat growth hormone gene. A comparison of the 5' and 3' flanking and untranslated regions of the bovine, human and rat growth hormone genes revealed many areas of highly conserved sequence. Especially noteworthy was the observation that all three genes had a 38 nucleotide homologous sequence within their 5' flanking regions located about 100 bp upstream from their transcription initiation sites.

Immunocytochemical localization of the prostaglandin-forming cyclooxygenase in the cerebellar cortex.

An indirect immunocytofluorescence technique was used to examine the distribution of the prostaglandin-forming cyclooxygenase in the cerebellar cortex of the pig, guinea, rat, mouse, cow, rabbit and sheep. Cyclooxygenase antigenicity was detected (a) in the cell bodies of Bergman glial cells in the Purkinje cell layer of the porcine, ovine and bovine cerebellar cortex; (b) in small arterioles throughout the cerebellar cortex in the sheep and cow; and (c) in the endothelial cells of large arteries in all the species examined. No cyclooxygenase-positive staining was apparent in neuronal cell bodies of granule, basket, stellate or Purkinje cells. Our results establish that prostaglandin endoperoxides can be synthesized by the arterial vasculature and at least certain glial cells in the central nervous system.