Illumina mate-paired DNA sequencing-library preparation using Cre-Lox recombination

Standard Illumina mate-paired libraries are constructed from 3- to 5-kb DNA fragments by a blunt-end circularization. Sequencing reads that pass through the junction of the two joined ends of a 3-5-kb DNA fragment are not easy to identify and pose problems during mapping and de novo assembly. Longer read lengths increase the possibility that a read will cross the junction. To solve this problem, we developed a mate-paired protocol for use with Illumina sequencing technology that uses Cre-Lox recombination instead of blunt end circularization. In this method, a LoxP sequence is incorporated at the junction site. This sequence allows screening reads for junctions without using a reference genome. Junction reads can be trimmed or split at the junction. Moreover, the location of the LoxP sequence in the reads distinguishes mate-paired reads from spurious paired-end reads. We tested this new method by preparing and sequencing a mate-paired library with an insert size of 3 kb from Saccharomyces cerevisiae. We present an analysis of the library quality statistics and a new bio-informatics tool called DeLoxer that can be used to analyze an IlluminaCre-Lox mate-paired data set. We also demonstrate how the resulting data significantly improves a de novo assembly of the S. cerevisiae genome

The role of integrated databases in microbial genome sequence analysis and metabolic reconstruction

This paper provides an overview of the PUMA system which provides access to data about metabolic pathways, enzymes, compounds, organisms, encoded activity, and assay condition information for enzymes in particular organisms and multiple sequence alignments

Estrogen pathway polymorphisms in relation to primary open angle glaucoma: An analysis accounting for gender from the United States

Purpose Circulating estrogen levels are relevant in glaucoma phenotypic traits. We assessed the association between an estrogen metabolism single nucleotide polymorphism (SNP) panel in relation to primary open angle glaucoma (POAG), accounting for gender. Methods: We included 3,108 POAG cases and 3,430 controls of both genders from the Glaucoma Genes and Environment (GLAUGEN) study and the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) consortium genotyped on the Illumina 660W-Quad platform. We assessed the relation between the SNP panels representative of estrogen metabolism and POAG using pathway- and gene-based approaches with the Pathway Analysis by Randomization Incorporating Structure (PARIS) software. PARIS executes a permutation algorithm to assess statistical significance relative to the pathways and genes of comparable genetic architecture. These analyses were performed using the meta-analyzed results from the GLAUGEN and NEIGHBOR data sets. We evaluated POAG overall as well as two subtypes of POAG defined as intraocular pressure (IOP) ≥22 mmHg (high-pressure glaucoma [HPG]) or IOP 0.99). Among women, gene-based analyses revealed that the catechol-O-methyltransferase gene showed strong associations with HTG (permuted gene p≤0.001) and NPG (permuted gene p=0.01). Conclusions: The estrogen SNP pathway was associated with POAG among women

GALC Deletions Increase the Risk of Primary Open-Angle Glaucoma: The Role of Mendelian Variants in Complex Disease

DNA copy number variants (CNVs) have been reported in many human diseases including autism and schizophrenia. Primary Open Angle Glaucoma (POAG) is a complex adult-onset disorder characterized by progressive optic neuropathy and vision loss. Previous studies have identified rare CNVs in POAG; however, their low frequencies prevented formal association testing. We present here the association between POAG risk and a heterozygous deletion in the galactosylceramidase gene (GALC). This CNV was initially identified in a dataset containing 71 Caucasian POAG cases and 478 ethnically matched controls obtained from dbGAP (study accession phs000126.v1.p1.) (p = 0.017, fisher's exact test). It was validated with array comparative genomic hybridization (arrayCGH) and realtime PCR, and replicated in an independent POAG dataset containing 959 cases and 1852 controls (p = 0.021, OR (odds ratio) = 3.5, 95% CI −1.1–12.0). Evidence for association was strengthened when the discovery and replication datasets were combined (p = 0.002; OR = 5.0, 95% CI 1.6–16.4). Several deletions with different endpoints were identified by array CGH of POAG patients. Homozygous deletions that eliminate GALC enzymatic activity cause Krabbe disease, a recessive Mendelian disorder of childhood displaying bilateral optic neuropathy and vision loss. Our findings suggest that heterozygous deletions that reduce GALC activity are a novel mechanism increasing risk of POAG. This is the first report of a statistically-significant association of a CNV with POAG risk, contributing to a growing body of evidence that CNVs play an important role in complex, inherited disorders. Our findings suggest an attractive biomarker and potential therapeutic target for patients with this form of POAG

A proteogenomic update to Yersinia: enhancing genome annotation

<p>Abstract</p> <p>Background</p> <p>Modern biomedical research depends on a complete and accurate proteome. With the widespread adoption of new sequencing technologies, genome sequences are generated at a near exponential rate, diminishing the time and effort that can be invested in genome annotation. The resulting gene set contains numerous errors in even the most basic form of annotation: the primary structure of the proteins.</p> <p>Results</p> <p>The application of experimental proteomics data to genome annotation, called proteogenomics, can quickly and efficiently discover misannotations, yielding a more accurate and complete genome annotation. We present a comprehensive proteogenomic analysis of the plague bacterium, <it>Yersinia pestis KIM</it>. We discover non-annotated genes, correct protein boundaries, remove spuriously annotated ORFs, and make major advances towards accurate identification of signal peptides. Finally, we apply our data to 21 other <it>Yersinia </it>genomes, correcting and enhancing their annotations.</p> <p>Conclusions</p> <p>In total, 141 gene models were altered and have been updated in RefSeq and Genbank, which can be accessed seamlessly through any NCBI tool (e.g. blast) or downloaded directly. Along with the improved gene models we discover new, more accurate means of identifying signal peptides in proteomics data.</p

A Common Variant in MIR182 Is Associated With Primary Open-Angle Glaucoma in the NEIGHBORHOOD Consortium

PURPOSE. Noncoding microRNAs (miRNAs) have been implicated in the pathogenesis of glaucoma. We aimed to identify common variants in miRNA coding genes (MIR) associated with primary open-angle glaucoma (POAG). METHODS. Using the NEIGHBORHOOD data set (3853 cases/33,480 controls with European ancestry), we first assessed the relation between 85 variants in 76 MIR genes and overall POAG. Subtype-specific analyses were performed in high-tension glaucoma (HTG) and normal-tension glaucoma subsets. Second, we examined the expression of miR-182, which was associated with POAG, in postmortem human ocular tissues (ciliary body, cornea, retina, and trabecular meshwork [TM]), using miRNA sequencing (miRNA-Seq) and droplet digital PCR (ddPCR). Third, miR-182 expression was also examined in human aqueous humor (AH) by using miRNA-Seq. Fourth, exosomes secreted from primary human TM cells were examined for miR-182 expression by using miRNA-Seq. Fifth, using ddPCR we compared miR182 expression in AH between five HTG cases and five controls. RESULTS. Only rs76481776 in MIR182 gene was associated with POAG after adjustment for multiple comparisons (odds ratio [OR] ¼ 1.23, 95% confidence interval [CI]: 1.11–1.42, P ¼ 0.0002). Subtype analysis indicated that the association was primarily in the HTG subset (OR ¼ 1.26, 95% CI: 1.08–1.47, P ¼ 0.004). The risk allele T has been associated with elevated miR-182 expression in vitro. Data from ddPCR and miRNA-Seq confirmed miR-182 expression in all examined ocular tissues and TM-derived exosomes. Interestingly, miR-182 expression in AH was 2-fold higher in HTG patients than nonglaucoma controls (P ¼ 0.03) without controlling for medication treatment. CONCLUSIONS. Our integrative study is the first to associate rs76481776 with POAG via elevated miR-182 expression

A Common Variant in MIR182 Is Associated With Primary Open-Angle Glaucoma in the NEIGHBORHOOD Consortium

PURPOSE. Noncoding microRNAs (miRNAs) have been implicated in the pathogenesis of glaucoma. We aimed to identify common variants in miRNA coding genes (MIR) associated with primary open-angle glaucoma (POAG). METHODS. Using the NEIGHBORHOOD data set (3853 cases/33,480 controls with European ancestry), we first assessed the relation between 85 variants in 76 MIR genes and overall POAG. Subtype-specific analyses were performed in high-tension glaucoma (HTG) and normal-tension glaucoma subsets. Second, we examined the expression of miR-182, which was associated with POAG, in postmortem human ocular tissues (ciliary body, cornea, retina, and trabecular meshwork [TM]), using miRNA sequencing (miRNA-Seq) and droplet digital PCR (ddPCR). Third, miR-182 expression was also examined in human aqueous humor (AH) by using miRNA-Seq. Fourth, exosomes secreted from primary human TM cells were examined for miR-182 expression by using miRNA-Seq. Fifth, using ddPCR we compared miR182 expression in AH between five HTG cases and five controls. RESULTS. Only rs76481776 in MIR182 gene was associated with POAG after adjustment for multiple comparisons (odds ratio [OR] ¼ 1.23, 95% confidence interval [CI]: 1.11–1.42, P ¼ 0.0002). Subtype analysis indicated that the association was primarily in the HTG subset (OR ¼ 1.26, 95% CI: 1.08–1.47, P ¼ 0.004). The risk allele T has been associated with elevated miR-182 expression in vitro. Data from ddPCR and miRNA-Seq confirmed miR-182 expression in all examined ocular tissues and TM-derived exosomes. Interestingly, miR-182 expression in AH was 2-fold higher in HTG patients than nonglaucoma controls (P ¼ 0.03) without controlling for medication treatment. CONCLUSIONS. Our integrative study is the first to associate rs76481776 with POAG via elevated miR-182 expression

Association of Genetic Variants With Primary Open-Angle Glaucoma Among Individuals With African Ancestry.

Importance:Primary open-angle glaucoma presents with increased prevalence and a higher degree of clinical severity in populations of African ancestry compared with European or Asian ancestry. Despite this, individuals of African ancestry remain understudied in genomic research for blinding disorders. Objectives:To perform a genome-wide association study (GWAS) of African ancestry populations and evaluate potential mechanisms of pathogenesis for loci associated with primary open-angle glaucoma. Design, Settings, and Participants:A 2-stage GWAS with a discovery data set of 2320 individuals with primary open-angle glaucoma and 2121 control individuals without primary open-angle glaucoma. The validation stage included an additional 6937 affected individuals and 14 917 unaffected individuals using multicenter clinic- and population-based participant recruitment approaches. Study participants were recruited from Ghana, Nigeria, South Africa, the United States, Tanzania, Britain, Cameroon, Saudi Arabia, Brazil, the Democratic Republic of the Congo, Morocco, Peru, and Mali from 2003 to 2018. Individuals with primary open-angle glaucoma had open iridocorneal angles and displayed glaucomatous optic neuropathy with visual field defects. Elevated intraocular pressure was not included in the case definition. Control individuals had no elevated intraocular pressure and no signs of glaucoma. Exposures:Genetic variants associated with primary open-angle glaucoma. Main Outcomes and Measures:Presence of primary open-angle glaucoma. Genome-wide significance was defined as P &lt; 5 × 10-8 in the discovery stage and in the meta-analysis of combined discovery and validation data. Results:A total of 2320 individuals with primary open-angle glaucoma (mean [interquartile range] age, 64.6 [56-74] years; 1055 [45.5%] women) and 2121 individuals without primary open-angle glaucoma (mean [interquartile range] age, 63.4 [55-71] years; 1025 [48.3%] women) were included in the discovery GWAS. The GWAS discovery meta-analysis demonstrated association of variants at amyloid-β A4 precursor protein-binding family B member 2 (APBB2; chromosome 4, rs59892895T&gt;C) with primary open-angle glaucoma (odds ratio [OR], 1.32 [95% CI, 1.20-1.46]; P = 2 × 10-8). The association was validated in an analysis of an additional 6937 affected individuals and 14 917 unaffected individuals (OR, 1.15 [95% CI, 1.09-1.21]; P &lt; .001). Each copy of the rs59892895*C risk allele was associated with increased risk of primary open-angle glaucoma when all data were included in a meta-analysis (OR, 1.19 [95% CI, 1.14-1.25]; P = 4 × 10-13). The rs59892895*C risk allele was present at appreciable frequency only in African ancestry populations. In contrast, the rs59892895*C risk allele had a frequency of less than 0.1% in individuals of European or Asian ancestry. Conclusions and Relevance:In this genome-wide association study, variants at the APBB2 locus demonstrated differential association with primary open-angle glaucoma by ancestry. If validated in additional populations this finding may have implications for risk assessment and therapeutic strategies

Lens intracellular hydrostatic pressure is generated by the circulation of sodium and modulated by gap junction coupling

We recently modeled fluid flow through gap junction channels coupling the pigmented and nonpigmented layers of the ciliary body. The model suggested the channels could transport the secretion of aqueous humor, but flow would be driven by hydrostatic pressure rather than osmosis. The pressure required to drive fluid through a single layer of gap junctions might be just a few mmHg and difficult to measure. In the lens, however, there is a circulation of Na+ that may be coupled to intracellular fluid flow. Based on this hypothesis, the fluid would cross hundreds of layers of gap junctions, and this might require a large hydrostatic gradient. Therefore, we measured hydrostatic pressure as a function of distance from the center of the lens using an intracellular microelectrode-based pressure-sensing system. In wild-type mouse lenses, intracellular pressure varied from ∼330 mmHg at the center to zero at the surface. We have several knockout/knock-in mouse models with differing levels of expression of gap junction channels coupling lens fiber cells. Intracellular hydrostatic pressure in lenses from these mouse models varied inversely with the number of channels. When the lens’ circulation of Na+ was either blocked or reduced, intracellular hydrostatic pressure in central fiber cells was either eliminated or reduced proportionally. These data are consistent with our hypotheses: fluid circulates through the lens; the intracellular leg of fluid circulation is through gap junction channels and is driven by hydrostatic pressure; and the fluid flow is generated by membrane transport of sodium