Jason L. Weirather

Serial Quantitative PCR Assay for Detection, Species Discrimination, and Quantification of Leishmania spp. in Human Samples

ABSTRACT The Leishmania species cause a variety of human disease syndromes. Methods for diagnosis and species differentiation are insensitive and many require invasive sampling. Although quantitative PCR (qPCR) methods are reported for leishmania detection, no systematic method to quantify parasites and determine the species in clinical specimens is established. We developed a serial qPCR strategy to identify and rapidly differentiate Leishmania species and quantify parasites in clinical or environmental specimens. SYBR green qPCR is mainly employed, with corresponding TaqMan assays for validation. The screening primers recognize kinetoplast minicircle DNA of all Leishmania species. Species identification employs further qPCR set(s) individualized for geographic regions, combining species-discriminating probes with melt curve analysis. The assay was sufficient to detect Leishmania parasites, make species determinations, and quantify Leishmania spp. in sera, cutaneous biopsy specimens, or cultured isolates from subjects from Bangladesh or Brazil with different forms of leishmaniasis. The multicopy kinetoplast DNA (kDNA) probes were the most sensitive and useful for quantification based on promastigote standard curves. To test their validity for quantification, kDNA copy numbers were compared between Leishmania species, isolates, and life stages using qPCR. Maxicircle and minicircle copy numbers differed up to 6-fold between Leishmania species, but the differences were smaller between strains of the same species. Amastigote and promastigote leishmania life stages retained similar numbers of kDNA maxi- or minicircles. Thus, serial qPCR is useful for leishmania detection and species determination and for absolute quantification when compared to a standard curve from the same Leishmania species.

Interferon Regulatory Factor 6 promotes differentiation of the periderm by activating expression of Grainyhead-like 3

Interferon Regulatory Factor 6 (IRF6) is a transcription factor that, in mammals, is required for the differentiation of skin, breast epithelium, and oral epithelium. However, the transcriptional targets that mediate these effects are currently unknown. In zebrafish and frog embryos Irf6 is necessary for differentiation of the embryonic superficial epithelium, or periderm. Here we use microarrays to identify genes that are expressed in the zebrafish periderm and whose expression is inhibited by a dominant-negative variant of Irf6 (dnIrf6). These methods identify Grhl3, an ancient regulator of the epidermal permeability barrier, as acting downstream of Irf6. In human keratinocytes, IRF6 binds conserved elements near the GHRL3 promoter. We show that one of these elements has enhancer activity in human keratinocytes and zebrafish periderm, suggesting that Irf6 directly stimulates Grhl3 expression in these tissues. Simultaneous inhibition of grhl1 and grhl3 disrupts periderm differentiation in zebrafish, and, intriguingly, forced grhl3 expression restores periderm markers in both zebrafish injected with dnIrf6 and frog embryos depleted of Irf6. Finally, in Irf6 deficient mouse embryos, Grhl3 expression in the periderm and oral epithelium is virtually absent. These results indicate that Grhl3 is a key effector of Irf6 in periderm differentiation.

Study of parasite kinetics with antileishmanial drugs using real-time quantitative PCR in Indian visceral leishmaniasis

OBJECTIVES This study describes parasite kinetics in the blood of visceral leishmaniasis patients treated with liposomal amphotericin B (L-AmB) or a preformed fat emulsion of amphotericin B (ApL) using real-time quantitative PCR (qPCR). METHODS Forty-six patients were treated with a single dose (15 mg/kg of body weight) of either L-AmB (n = 13) or ApL (n = 33). qPCR was used to estimate parasite kinetics by detection of Leishmania donovani DNA using kinetoplast DNA-specific primers in peripheral blood samples using an absolute quantification method. RESULTS The mean parasite load decreased from baseline (day 0) values of 894.07 and 980.48 to 71.72 and 211.52 parasite genomes/mL at day 7 in L-AmB and ApL groups, respectively, and at day 30 these further declined to 8.30 and 133.98 parasite genomes/mL, respectively. At day 30 post-treatment evaluation, the decline in parasite load was significantly greater (P = 0.024) with L-AmB compared with ApL. Four of 33 patients in the ApL group failed treatment (1 primary failure and 3 relapses) with the presence of parasites, whereas all patients in the L-AmB group were cured at 6 month follow-up. CONCLUSIONS qPCR can be a tool to measure parasite dynamics accurately and provide a marker to measure the efficacy of various drugs. It can be used as a test of cure, allowing us to do away with invasive and risky methods such as splenic or bone marrow aspiration.

Characterization of fusion genes and the significantly expressed fusion isoforms in breast cancer by hybrid sequencing

We developed an innovative hybrid sequencing approach, IDP-fusion, to detect fusion genes, determine fusion sites and identify and quantify fusion isoforms. IDP-fusion is the first method to study gene fusion events by integrating Third Generation Sequencing long reads and Second Generation Sequencing short reads. We applied IDP-fusion to PacBio data and Illumina data from the MCF-7 breast cancer cells. Compared with the existing tools, IDP-fusion detects fusion genes at higher precision and a very low false positive rate. The results show that IDP-fusion will be useful for unraveling the complexity of multiple fusion splices and fusion isoforms within tumorigenesis-relevant fusion genes.

Leishmania infantum chagasi in Northeastern Brazil: Asymptomatic Infection at the Urban Perimeter

Visceral leishmaniasis (VL) is endemic in large cities in Brazil, including Natal. We determined the prevalence of asymptomatic human infection with Leishmania infantum chagasi and associated environmental risks around Natal. Infection was detected by Leishmania skin test (LST) and anti-leishmanial antibodies in humans and anti-leishmanial antibodies in dogs. Amongst 345 humans, 24.6% were seropositive, and 38.6% were LST-positive. Prevalence of positive serology was similar in both sexes and across all ages. However, positive LST responses increased with age, suggesting that LST is long-lasting and cumulative. Multinomial logistic analysis showed that LST response varied with location (P = 0.007) and that males were more frequently LST-positive (P = 0.027). Indicators of lower socioeconomic status associated significantly with human infection. Furthermore, there was geographic coincidence of seropositive humans and dogs (r = 0.7926, P = 0.011). These data suggest that dog and human L. i. chagasi infection are intimately interrelated in environmental conditions associated with low income.

Comprehensive comparison of Pacific Biosciences and Oxford Nanopore Technologies and their applications to transcriptome analysis

Background: Given the demonstrated utility of Third Generation Sequencing [Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT)] long reads in many studies, a comprehensive analysis and comparison of their data quality and applications is in high demand. Methods: Based on the transcriptome sequencing data from human embryonic stem cells, we analyzed multiple data features of PacBio and ONT, including error pattern, length, mappability and technical improvements over previous platforms. We also evaluated their application to transcriptome analyses, such as isoform identification and quantification and characterization of transcriptome complexity, by comparing the performance of PacBio, ONT and their corresponding Hybrid-Seq strategies (PacBio+Illumina and ONT+Illumina). Results: PacBio shows overall better data quality, while ONT provides a higher yield. As with data quality, PacBio performs marginally better than ONT in most aspects for both long reads only and Hybrid-Seq strategies in transcriptome analysis. In addition, Hybrid-Seq shows superior performance over long reads only in most transcriptome analyses. Conclusions: Both PacBio and ONT sequencing are suitable for full-length single-molecule transcriptome analysis. As this first use of ONT reads in a Hybrid-Seq analysis has shown, both PacBio and ONT can benefit from a combined Illumina strategy. The tools and analytical methods developed here provide a resource for future applications and evaluations of these rapidly-changing technologies.

Genomic Strategy Identifies a Missense Mutation in WD-Repeat Domain 65 (WDR65) in an Individual With Van der Woude Syndrome

Genetic variation in the transcription factor interferon regulatory factor 6 (IRF6) causes and contributes risk for oral clefting disorders. We hypothesized that genes regulated by IRF6 are also involved in oral clefting disorders. We used five criteria to identify potential IRF6 target genes; differential gene expression in skin taken from wild‐type and Irf6‐deficient murine embryos, localization to the Van der Woude syndrome 2 (VWS2) locus at 1p36‐1p32, overlapping expression with Irf6, presence of a conserved predicted‐binding site in the promoter region, and a mutant murine phenotype that was similar to the Irf6 mutant mouse. Previously, we observed altered expression for 573 genes; 13 were located in the murine region syntenic to the VWS2 locus. Two of these genes, Wdr65 and Stratifin, met 4 of 5 criteria. Wdr65 was a novel gene that encoded a predicted protein of 1,250 amino acids with two WD domains. As potential targets for Irf6 regulation, we hypothesized that disease‐causing mutations will be found in WDR65 and Stratifin in individuals with VWS or VWS‐like syndromes. We identified a potentially etiologic missense mutation in WDR65 in a person with VWS who does not have an exonic mutation in IRF6. The expression and mutation data were consistent with the hypothesis that WDR65 was a novel gene involved in oral clefting.

Non-coding yet non-trivial: a review on the computational genomics of lincRNAs.

Long intergenic non-coding RNAs (lincRNAs) represent one of the most mysterious RNA species encoded by the human genome. Thanks to next generation sequencing (NGS) technology and its applications, we have recently witnessed a surge in non-coding RNA research, including lincRNA research. Here, we summarize the recent advancement in genomics studies of lincRNAs. We review the emerging characteristics of lincRNAs, the experimental and computational approaches to identify lincRNAs, their known mechanisms of regulation, the computational methods and resources for lincRNA functional predictions, and discuss the challenges to understanding lincRNA comprehensively.

IDP-ASE: haplotyping and quantifying allele-specific expression at the gene and gene isoform level by hybrid sequencing

Abstract Allele-specific expression (ASE) is a fundamental problem in studying gene regulation and diploid transcriptome profiles, with two key challenges: (i) haplotyping and (ii) estimation of ASE at the gene isoform level. Existing ASE analysis methods are limited by a dependence on haplotyping from laborious experiments or extra genome/family trio data. In addition, there is a lack of methods for gene isoform level ASE analysis. We developed a tool, IDP-ASE, for full ASE analysis. By innovative integration of Third Generation Sequencing (TGS) long reads with Second Generation Sequencing (SGS) short reads, the accuracy of haplotyping and ASE quantification at the gene and gene isoform level was greatly improved as demonstrated by the gold standard data GM12878 data and semi-simulation data. In addition to methodology development, applications of IDP-ASE to human embryonic stem cells and breast cancer cells indicate that the imbalance of ASE and non-uniformity of gene isoform ASE is widespread, including tumorigenesis relevant genes and pluripotency markers. These results show that gene isoform expression and allele-specific expression cooperate to provide high diversity and complexity of gene regulation and expression, highlighting the importance of studying ASE at the gene isoform level. Our study provides a robust bioinformatics solution to understand ASE using RNA sequencing data only.

Comprehensive candidate gene analysis for symptomatic or asymptomatic outcomes of Leishmania infantum infection in Brazil

Genetic risk factors contribute to asymptomatic versus symptomatic visceral leishmaniasis (VL) outcomes following infection with Leishmania infantum. We therefore carried out a family‐based (n = 918 post–quality control fully genotyped and phenotyped individuals) candidate gene study for symptomatic VL or asymptomatic delayed‐type hypersensitivity (DTH) skin test phenotypes in highly endemic neighborhoods of northeast Brazil. A total of 248 SNPs were genotyped in 42 genes selected as candidates on the basis of prior genetic, immunological, and transcriptional profiling studies. The most significant association with the VL phenotype was with SNP rs6785358 (P = 5.7e‐04; pcorrected = 0.026) 3.8 kb upstream of TGFBR2, the gene encoding the type 2 receptor for transforming growth factor beta (TGFβ). A second inhibitory member of the TGBβ superfamily signaling pathway, SMAD7, was associated with the DTH phenotype (SNP rs7238442: P = 0.001; pcorrected = 0.051). The most significant association for the DTH phenotype was with SNP rs10800309 (P = –8.4e‐06; pcorrected = 3.9e‐04) situated 3.1 kb upstream of FCGR2A, the gene encoding the low‐affinity IIa receptor for the Fc fragment of IgG. Overall, our results imply a role for IgG‐mediated inflammation in determining DTH associated with asymptomatic infection and contribute to growing evidence that the TGFβ pathway is important in the immunopathogenesis of VL.