Search

Background There is an incompletely understood increased risk for cardiovascular disease (CVD) among people with HIV (PWH). We investigated if a collection of biomarkers were associated with CVD among PWH. Mendelian randomization (MR) was used to identify potentially causal associations. Methods Data from follow-up in 4 large trials among PWH were used to identify 131 incident CVD cases and they were matched to 259 participants without incident CVD (controls). Tests of associations between 460 baseline protein levels and case status were conducted. Results Univariate analysis found CLEC6A, HGF, IL-6, IL-10RB, and IGFBP7 as being associated with case status and a multivariate model identified 3 of these: CLEC6A (odds ratio [OR] = 1.48, P = .037), HGF (OR = 1.83, P = .012), and IL-6 (OR = 1.45, P = .016). MR methods identified 5 significantly associated proteins: AXL, CHI3L1, GAS6, IL-6RA, and SCGB3A2. Conclusions These results implicate inflammatory and fibrotic processes as contributing to CVD. While some of these biomarkers are well established in the general population and in PWH (IL-6 and its receptor), some are novel to PWH (HGF, AXL, and GAS6) and some are novel overall (CLEC6A). Further investigation into the uniqueness of these biomarkers in PWH and the role of these biomarkers as targets among PWH is warranted.

Background: Levels of plasma SARS-CoV-2 nucleocapsid (N) antigen may be an important biomarker in patients with COVID-19 and enhance our understanding of the pathogenesis of COVID-19. Objective: To evaluate whether levels of plasma antigen can predict short-term clinical outcomes and identify clinical and viral factors associated with plasma antigen levels in hospitalized patients with SARS-CoV-2. Design: Cross-sectional study of baseline plasma antigen level from 2540 participants enrolled in the TICO (Therapeutics for Inpatients With COVID-19) platform trial from August 2020 to November 2021, with additional data on day 5 outcome and time to discharge. Setting: 114 centers in 10 countries. Participants: Adults hospitalized for acute SARS-CoV-2 infection with 12 days or less of symptoms. Measurements: Baseline plasma viral N antigen level was measured at a central laboratory. Delta variant status was determined from baseline nasal swabs using reverse transcriptase polymerase chain reaction. Associations between baseline patient characteristics and viral factors and baseline plasma antigen levels were assessed using both unadjusted and multivariable modeling. Association between elevated baseline antigen level of 1000 ng/L or greater and outcomes, including worsening of ordinal pulmonary scale at day 5 and time to hospital discharge, were evaluated using logistic regression and Fine-Gray regression models, respectively. Results: Plasma antigen was below the level of quantification in 5% of participants at enrollment, and 1000 ng/L or greater in 57%. Baseline pulmonary severity of illness was strongly associated with plasma antigen level, with mean plasma antigen level 3.10-fold higher among those requiring noninvasive ventilation or high-flow nasal cannula compared with room air (95% CI, 2.22 to 4.34). Plasma antigen level was higher in those who lacked antispike antibodies (6.42 fold; CI, 5.37 to 7.66) and in those with the Delta variant (1.73 fold; CI, 1.41 to 2.13). Additional factors associated with higher baseline antigen level included male sex, shorter time since hospital admission, decreased days of remdesivir, and renal impairment. In contrast, race, ethnicity, body mass index, and immunocompromising conditions were not associated with plasma antigen levels. Plasma antigen level of 1000 ng/L or greater was associated with a markedly higher odds of worsened pulmonary status at day 5 (odds ratio, 5.06 [CI, 3.41 to 7.50]) and longer time to hospital discharge (median, 7 vs. 4 days; subhazard ratio, 0.51 [CI, 0.45 to 0.57]), with subhazard ratios similar across all levels of baseline pulmonary severity. Limitations: Plasma samples were drawn at enrollment, not hospital presentation. No point-of-care test to measure plasma antigen is currently available. Conclusion: Elevated plasma antigen is highly associated with both severity of pulmonary illness and clinically important patient outcomes. Multiple clinical and viral factors are associated with plasma antigen level at presentation. These data support a potential role of ongoing viral replication in the pathogenesis of SARS-CoV-2 in hospitalized patients.

Objectives Elevated levels of interleukin-6 (IL-6), D-dimer, and C-reactive protein (hsCRP) are associated with increased incidence of comorbid disease and mortality among people living with HIV (PLWH). Prior studies suggest a genetic basis for these biomarker elevations in the general population. The study objectives are to identify the genetic basis for these biomarkers among PLWH. Methods Baseline levels of hsCRP, D-dimer and IL-6, and single nucleotide polymorphisms (SNPs) were determined for 7,768 participants in three HIV treatment trials. Single variant analysis was performed for each biomarker on samples from each of three ethnic groups (African [AFR], Admixed American [AMR], European [EUR]) within each trial including covariates relevant to biomarker levels. For each ethnic group, the results were pooled across trials, then further pooled across ethnicities. Results The transethnic analysis identified three, two and one known loci associated with hsCRP, D-dimer and IL-6 levels, respectively, and two novel loci, FGB and GCNT1, associated with D-dimer levels. Lead SNPs exhibited similar effects across ethnicities. Additionally, three novel, ethnic-specific loci were identified: CATSPERG associated with D-dimer in AFR and PROX1-AS1 and TRAPPC9 associated with IL-6 in AFR and AMR, respectively. Conclusions Eleven loci associated with three biomarker levels were identified in PLWH from the three studies including six loci known in the general population and five novel loci associated with D-dimer and IL-6 levels. These findings support the hypothesis that host genetics may partially contribute to chronic inflammation in PLWH and help to identify potential targets for intervention of serious non-AIDS complications.

Supplemental Digital Content is Available in the Text., Objectives: Interleukin-27 (IL-27) is known as an anti-HIV cytokine. We have recently demonstrated that IL-27-pretreatment promotes phytohemagglutinin-stimulated CD4(+) T cells into HIV-1-resistant cells by inhibiting an uncoating step. Purpose: To further characterize the function of the HIV resistant T cells, we investigated profiles of microRNA in the cells using microRNA sequencing (miRNA-seq) and assessed anti-HIV effect of the microRNAs. Methods: Phytohemagglutinin-stimulated CD4(+) T cells were treated with or without IL-27 for 3 days. MicroRNA profiles were analyzed using miRNA-seq. To assess anti-HIV effect, T cells or macrophages were transfected with synthesized microRNA mimics and then infected with HIVNL4.3 or HIVAD8. Anti-HIV effect was monitored by a p24 antigen enzyme-linked immunosorbent assay kit. interferon (IFN)-α, IFN-β, or IFN-λ production was quantified using each subtype-specific enzyme-linked immunosorbent assay kit. Results: A comparative analysis of microRNA profiles indicated that expression of known miRNAs was not significantly changed in IL-27-treated cells compared with untreated T cells; however, a total of 15 novel microRNAs (miRTC1 ∼ miRTC15) were identified. Anti-HIV assay using overexpression of each novel microRNA revealed that 10 nM miRTC14 (GenBank accession number: MF281439) remarkably suppressed HIV infection by (99.3 ± 0.27%, n = 9) in macrophages but not in T cells. The inhibition was associated through induction of >1000 pg/mL of IFN-αs and IFN-λ1. Conclusion: We discovered a total of 15 novel microRNAs in T cells and characterized that miRTC14, one of the novel microRNAs, was a potent IFN-inducing anti-HIV miRNA, implicating that regulation of the expression of miRTC14 may be a potent therapeutic tool for not only HIV but also other virus infection.

Significance In HIV-infected patients on combination antiretroviral therapy (cART), greater than 95% of proviruses in the peripheral blood are “defective.” Historically, these defective proviruses have been thought to be dead-end products with no real pathophysiological significance, as they do not encode replication-competent viruses. Contrary to this view, we have identified cells in tissue culture and from cART-treated patients that harbor defective proviruses and produce viral proteins. Features found in these translationally competent yet defective proviruses suggest that HIV-1 infection results in modification of the CD4+ T cell genome analogous to human endogenous retroviruses. We propose that these defective HIV-1 proviruses are biologically significant, despite being “replication incompetent,” have the potential to elicit immune activation, and may serve as a barrier to HIV-1 cure., HIV-1 proviruses persist in the CD4+ T cells of HIV-infected individuals despite years of combination antiretroviral therapy (cART) with suppression of HIV-1 RNA levels

BACKGROUND: In a randomized, placebo-controlled, clinical trial, bamlanivimab, a SARS-CoV-2-neutralizing monoclonal antibody, given in combination with remdesivir, did not improve outcomes among hospitalized patients with COVID-19 based on an early futility assessment.OBJECTIVE: To evaluate the a priori hypothesis that bamlanivimab has greater benefit in patients without detectable levels of endogenous neutralizing antibody (nAb) at study entry than in those with antibodies, especially if viral levels are high.DESIGN: Randomized, placebo-controlled trial. (ClinicalTrials.gov: NCT04501978).SETTING: Multicenter trial.PATIENTS: Hospitalized patients with COVID-19 without end-organ failure.INTERVENTION: Bamlanivimab (7000 mg) or placebo.MEASUREMENTS: Antibody, antigen, and viral RNA levels were centrally measured on stored specimens collected at baseline. Patients were followed for 90 days for sustained recovery (defined as discharge to home and remaining home for 14 consecutive days) and a composite safety outcome (death, serious adverse events, organ failure, or serious infections).RESULTS: Among 314 participants (163 receiving bamlanivimab and 151 placebo), the median time to sustained recovery was 19 days and did not differ between the bamlanivimab and placebo groups (subhazard ratio [sHR], 0.99 [95% CI, 0.79 to 1.22]; sHR > 1 favors bamlanivimab). At entry, 50% evidenced production of anti-spike nAbs; 50% had SARS-CoV-2 nucleocapsid plasma antigen levels of at least 1000 ng/L. Among those without and with nAbs at study entry, the sHRs were 1.24 (CI, 0.90 to 1.70) and 0.74 (CI, 0.54 to 1.00), respectively (nominal P for interaction = 0.018). The sHR (bamlanivimab vs. placebo) was also more than 1 for those with plasma antigen or nasal viral RNA levels above median level at entry and was greatest for those without antibodies and with elevated levels of antigen (sHR, 1.48 [CI, 0.99 to 2.23]) or viral RNA (sHR, 1.89 [CI, 1.23 to 2.91]). Hazard ratios for the composite safety outcome (LIMITATION: Subgroup analysis of a trial prematurely stopped because of futility; small sample size; multiple subgroups analyzed.CONCLUSION: Efficacy and safety of bamlanivimab may differ depending on whether an endogenous nAb response has been mounted. The limited sample size of the study does not allow firm conclusions based on these findings, and further independent trials are required that assess other types of passive immune therapies in the same patient setting.PRIMARY FUNDING SOURCE: U.S. government Operation Warp Speed and National Institute of Allergy and Infectious Diseases.

Macrophages play an essential role in the immune system. Recent studies have shown that long non-coding RNAs (lncRNAs) can regulate genes encoding products involved in the immune response. Interleukin (IL)-27 is a member of the IL-6/IL-12 family of cytokines with broad anti-viral effects that inhibits human immunodeficiency virus (HIV) type-1 and herpes simplex virus (HSV). However, little is known about the role of lncRNAs in macrophages affected by IL-27. Therefore, we investigated the expression profiles of mRNA and lncRNA in human monocyte-derived macrophages (MDMs) regulated by IL-27. Monocytes were differentiated in the presence of macrophage-colony stimulatory factor (M-CSF)- or human AB serum with or without IL-27, and these cells were the subject for the profile analysis using RNA-Seq. We identified 146 lncRNAs (including 88 novel ones) and 434 coding genes were differentially regulated by IL-27 in both M-CSF- and AB serum-induced macrophages. Using weighted gene co-expression network analysis, we obtained four modules. The immune system, cell cycle, and regulation of complement cascade pathways were enriched in different modules. The network of mRNAs and lncRNAs in the pathways suggest that lncRNAs might regulate immune activity in macrophages. This study provides potential insight into the roles of lncRNA in macrophages regulated by IL-27.

Interleukin-27 (IL-27) is a pleiotropic cytokine that influences the innate and adaptive immune systems. It inhibits viral infection and regulates the expression of microRNAs (miRNAs). We recently reported that macrophages differentiated from human primary monocytes in the presence of IL-27 and human AB serum resisted human immunodeficiency virus (HIV) infection and showed significant autophagy induction. In the current study, the miRNA profiles in these cells were investigated, especially focusing on the identification of novel miRNAs regulated by IL-27-treatment. The miRNA sequencing analysis detected 38 novel miRNAs. Real-time reverse transcription polymerase chain reaction (RT-PCR) analysis confirmed that IL-27 differentially regulated the expression of 16 of the 38 miRNAs. Overexpression of the synthesized miRNA mimics by transfection revealed that miRAB40 had potent HIV-inhibiting and autophagy-inducing properties. B18R, an interferon (IFN)-neutralization protein, partially suppressed both activities, indicating that the two functions were induced via IFN-dependent and -independent pathways. Although the target mRNA(s) of miRAB40 involving in the induction of both functions was unable to identify in this study, the discovery of miRAB40, a potential HIV-inhibiting and autophagy inducing miRNA, may provide novel insights into the miRNA (small none-coding RNA)-mediated regulation of HIV inhibition and autophagy induction as an innate immune response.

Macrophages play an essential role in the immune system. Recent studies have shown that long non-coding RNAs (lncRNAs) can regulate genes encoding products involved in the immune response. Interleukin (IL)-27 is a member of the IL-6/IL-12 family of cytokines with broad anti-viral effects that inhibits human immunodeficiency virus (HIV) type-1 and herpes simplex virus (HSV). However, little is known about the role of lncRNAs in macrophages affected by IL-27. Therefore, we investigated the expression profiles of mRNA and lncRNA in human monocyte-derived macrophages (MDMs) regulated by IL-27. Monocytes were differentiated in the presence of macrophage-colony stimulatory factor (M-CSF)- or human AB serum with or without IL-27, and these cells were the subject for the profile analysis using RNA-Seq. We identified 146 lncRNAs (including 88 novel ones) and 434 coding genes were differentially regulated by IL-27 in both M-CSF- and AB serum-induced macrophages. Using weighted gene co-expression network analysis, we obtained four modules. The immune system, cell cycle, and regulation of complement cascade pathways were enriched in different modules. The network of mRNAs and lncRNAs in the pathways suggest that lncRNAs might regulate immune activity in macrophages. This study provides potential insight into the roles of lncRNA in macrophages regulated by IL-27.

Herpes simplex virus 1 (HSV-1) strain McKrae is highly virulent and relatively neuroinvasive in animal models compared with other wild-type HSV-1 strains. To identify the genetic determinants that lead to the unique phenotypes of the McKrae strain, we sequenced its genome with PacBio single-molecule real-time (SMRT) technology and resolved the complete sequence.

BACKGROUND: HIV-infected patients with poor virologic control and multidrug-resistant virus have limited therapeutic options. The current study was undertaken to evaluate the safety, immunologic effects, and antiviral activity of peripheral lymphocytes transferred from an elite controller, whose immune system is able to control viral replication without antiretroviral medications, to an HLA-B*2705–matched progressor. METHODS: Approximately 22 billion cells were collected from an elite controller by lymphapheresis and infused within 6 hours into a recipient with a preinfusion CD4(+) T cell count of 10 cells/μL (1%) and HIV plasma viral load of 114,993 copies/mL. RESULTS: Donor cells were cleared from the recipient’s peripheral blood by day 8. A transient decrease in viral load to 58,421 (day 3) was followed by a rebound to 702,972 (day 6) before returning to baseline values by day 8. The decreased viral load was temporally associated with peak levels of donor T cells, including CD8(+) T cells that had high levels of expression of Ki67, perforin, and granzyme B. Notably, recipient CD8(+) T cells also showed increased expression of these markers, especially in HIV-specific tetramer–positive cells. CONCLUSION: These results suggest that the adoptive transfer of lymphocytes from an HIV-infected elite controller to an HIV-infected patient with progressive disease may be able to perturb the immune system of the recipient in both positive and negative ways. TRIAL REGISTRATION: ClinicalTrials.gov NCT00559416. FUNDING: Intramural Research Programs of the US NIH Clinical Center and the National Institute of Allergy and Infectious Diseases (NIAID); the National Cancer Institute.

The impact of variation in host genetics on replication of human immunodeficiency virus type 1 (HIV-1) in demographically diverse populations remains uncertain. In the current study, we performed a genome-wide screen for associations of single-nucleotide polymorphisms (SNPs) to viral load (VL) in antiretroviral therapy–naive participants (n = 2440) with varying demographics from the Strategic Timing of AntiRetroviral Treatment (START) trial. Associations were assessed using genotypic data generated by a customized SNP array, imputed HLA alleles, and multiple linear regression. Genome-wide significant associations between SNPs and VL were observed in the major histocompatibility complex class I region (MHC I), with effect sizes ranging between 0.14 and 0.39 log10 VL (copies/mL). Supporting the SNP findings, we identified several HLA alleles significantly associated with VL, extending prior observations that the (MHC I) is a major host determinant of HIV-1 control with shared genetic variants across diverse populations and underscoring the limitations of genome-wide association studies as being merely a screening tool.

After initiating antiretroviral therapy (ART), a rapid decline in plasma viremia is followed by reservoir stabilization. Viral outgrowth assay suggests the reservoir continues to decline slowly, but variation over time and among individuals complicates our understanding of selective pressures during ART. We used full-length sequencing to study more than 800 HIV proviruses of two subjects on ART at four time points over nine years to investigate the selection pressures influencing the dynamics of the reservoir. We found that intact as well as defective proviruses capable of significant protein expression decrease over time. Moreover, proviruses lacking genetic elements to promote viral protein expression, yet containing strong splice donor sequences increase relative to other defectives over time, especially among clones. Our work suggests that HIV expression occurs to a significant extent during ART and results in HIV clearance, but this is obscured by clones generated by donor splice site-enhanced clonal expansion.

We have previously demonstrated that Interleukin-27 differentially regulates the expression of seven novel microRNAs. Here we elucidate the functional significance of these novel microRNAs. Of the seven microRNAs, over expression of miRNA-6852 (miR-SX4) mimic induces cell cycle arrest at G2/M phase and induces necrosis in HEK293 and panel of cervical cancer cells (Human Papilloma Virus (HPV) infected cell lines; HeLa, CaSki and SiHa cells). To define the mechanism of the miR-SX4-mediated G2/M arrest, a microarray gene chip array and western blot analysis were performed. FoxM1, a transcription factor is identified as a key protein down-regulated by miR-SX4, even though the miR-SX4 does not target 3’UTR of FoxM1. Knock down of FoxM1 using si-RNA demonstrate that FoxM1 silenced cell induces G2/M cell cycle arrest and necrosis. Our data demonstrated for the first time that miR-SX4 could be a potent anti-cancer microRNA.

Many of the infectious diseases which have jeopardized and still are a threat to public health are caused by RNA viruses, including HIV, HCV, Influenza virus, Ebola virus and Zika virus. Because of a high rate of mutations and recombination events, rapidly evolving RNA viruses prevail within a host as a collection of closely related variants, referred to as viral quasi-species. Uncovering the genetic diversity (i.e., inferring viral haplotypes and their proportions in the population) of an RNA virus can significantly benefit the study of disease progression, antiviral drug design, vaccine design and viral pathogenesis. Recent advances in PacBio single-molecule sequencing offers sufficient throughput and contiguous long reads (>10kb) covering the full length of most genes and RNA viral genomes, providing a potential to reliably profile the viral populations. However, the relatively high error rate (2~15%) in the long-read data requires novel analysis methods to deconvolute sequences derived from complex viral mixtures. We examined samples containing complex mixtures of near-full-length HIV-1 genomes, single molecules sequenced as near-full-length (9kb) amplicons directly from PCR products, and developed a novel signature-based self-tuning spectral clustering method called SigClust to accurately determine the identity (above 99.5%) and relative abundances of viral genomes in the mixtures. Results on real HIV-1 and influenza benchmark data sets demonstrate efficacy and superior performance of SigClust.

We previously identified Ku70, a subunit of a DNA repair protein complex, as a cytosolic DNA sensor that induces the production of interferon-λ1 (IFN-λ1) by human primary cells and cell lines. IFN-λ1 is a type III IFN and has similar antiviral activity to that of the type I IFNs (IFN-α and IFN-β). We observed that human embryonic kidney (HEK) 293T cells, which are deficient in the innate immune adaptor protein STING (stimulator of IFN genes), did not produce IFN-λ1 in response to DNA unless they were reconstituted with STING. Conversely, parental HEK 293 cells produced IFN-λ1 after they were exposed to exogenous DNA; however, when STING was knocked out in the HEK 293 cells through the CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 genome editing system, they lost this response. Through confocal microscopy, we demonstrated that endogenous Ku70 was located in the nucleus and then translocated to the cytoplasm upon DNA exposure to form a complex with STING. Additionally, the DNA binding domain of Ku70 was essential for formation of the Ku70-STING complex. Knocking down STING in primary human macrophages inhibited their ability to produce IFN-λ1 in response to transfection with DNA or infection with the DNA virus HSV-2 (herpes simplex virus–2). Together, these data suggest that STING mediates the Ku70-mediated IFN-λ1 innate immune response to exogenous DNA or DNA virus infection.

MicroRNAs (miRNAs) regulate gene expression and thereby influence cell fate and function. Recent studies suggest that an abundant class of miRNAs play important roles in immune cells, such as T cells, natural killer (NK) cells, B cells, and dendritic cells (DCs). Interleukin (IL)-27 is a member of the IL-12 family of cytokines with broad anti-viral effects. It is a potent inhibitor of HIV-1 infection in CD4+ T cells and macrophages, as well as monocyte-derived immature dendritic cells (iDCs). This pilot study compared miRNA profiles between iDCs and IL-27-treated iDCs (27DCs) using deep sequencing methods and identified 46 known miRNAs that were significantly differentially expressed in 27DCs: 36 were upregulated and 10 downregulated by IL-27. Many of the potential target genes of these miRNAs are involved in IL-27 associated pathways, such as JAK/STAT, MAPKs, and PI3K and several were also previously reported to be involved in the regulation of human DC function. This study found that these miRNAs also potentially target several viral genomes and therefore may have antiviral effects. Four of these differential miRNAs (miR-99a-5p, miR-222-3p, miR-138-5p, and miR-125b-5p) were validated using quantitative reverse transcription PCR (RT-qPCR). Twenty-two novel miRNAs were discovered from deep sequencing and confirmed using RT-qPCR. This study furthers the understanding of the role of IL-27 in immunity and lays a foundation for future characterization of the role of specific miRNAs in DCs.

Pneumocystis jirovecii is an important opportunistic pathogen associated with AIDS and other immunodeficient conditions. Currently, very little is known about its nuclear and mitochondrial genomes. In this study, we sequenced the complete mitochondrial genome (mtDNA) of this organism and its closely related species Pneumocystis carinii and Pneumocystis murina by a combination of sequencing technologies. Our study shows that P. carinii and P. murina mtDNA share a nearly identical number and order of genes in a linear configuration, whereas P. jirovecii has a circular mtDNA containing nearly the same set of genes but in a different order. Detailed studies of the mtDNA terminal structures of P. murina and P. carinii suggest a unique replication mechanism for linear mtDNA. Phylogenetic analysis supports a close association of Pneumocystis species with Taphrina, Saitoella, and Schizosaccharomyces, and divergence within Pneumocystis species, with P. murina and P. carinii being more closely related to each other than either is to P. jirovecii. Comparative analysis of four complete P. jirovecii mtDNA sequences in this study and previously reported mtDNA sequences for diagnosing and genotyping suggests that the current diagnostic and typing methods can be improved using the complete mtDNA data. The availability of the complete P. jirovecii mtDNA also opens the possibility of identifying new therapeutic targets.—Ma, L., Huang, D. W., Cuomo, C. A., Sykes, S., Fantoni, G., Das, B., Sherman, B. T., Yang, J., Huber, C., Xia, Y., Davey, E., Kutty, G., Bishop, L., Sassi, M., Lempicki, R. A., Kovacs, J. A. Sequencing and characterization of the complete mitochondrial genomes of three Pneumocystis species provide new insights into divergence between human and rodent Pneumocystis.

β-glucans, which can activate innate immune responses, are a major component in the cell wall of the cyst form of Pneumocystis. In the current study, we examined whether β-1,3-glucans are masked by surface proteins in Pneumocystis and what role β-glucans play in Pneumocystis-associated inflammation. For 3 species, including Pneumocystis jirovecii, which causes Pneumocystis pneumonia in humans, Pneumocystis carinii, and Pneumocystis murina, β-1,3-glucans were masked in most organisms, as demonstrated by increased exposure following trypsin treatment. Using quantitative polymerase chain reaction and microarray techniques, we demonstrated in a mouse model of Pneumocystis pneumonia that treatment with caspofungin, an inhibitor of β-1,3-glucan synthesis, for 21 days decreased expression of a broad panel of inflammatory markers, including interferon γ, tumor necrosis factor α, interleukin 1β, interleukin 6, and multiple chemokines/chemokine ligands. Thus, β-glucans in Pneumocystis cysts are largely masked, which likely decreases innate immune activation; this mechanism presumably was developed for interactions with immunocompetent hosts, in whom organism loads are substantially lower. In immunosuppressed hosts with a high organism burden, organism death and release of glucans appears to be an important contributor to deleterious host inflammatory responses.

Pneumocystis jirovecii is a major cause of life-threatening pneumonia in immunosuppressed patients including transplant recipients and those with HIV/AIDS, yet surprisingly little is known about the biology of this fungal pathogen. Here we report near complete genome assemblies for three Pneumocystis species that infect humans, rats and mice. Pneumocystis genomes are highly compact relative to other fungi, with substantial reductions of ribosomal RNA genes, transporters, transcription factors and many metabolic pathways, but contain expansions of surface proteins, especially a unique and complex surface glycoprotein superfamily, as well as proteases and RNA processing proteins. Unexpectedly, the key fungal cell wall components chitin and outer chain N-mannans are absent, based on genome content and experimental validation. Our findings suggest that Pneumocystis has developed unique mechanisms of adaptation to life exclusively in mammalian hosts, including dependence on the lungs for gas and nutrients and highly efficient strategies to escape both host innate and acquired immune defenses., Pneumocystis jirovecii is a fungus that can cause life-threatening pneumonia in immunocompromised patients. Here, the authors sequence the genomes of P. jirovecii and two other Pneumocystis species, and show the unexpected absence of chitin (a near universal fungal cell wall component).

Background HIV-1 set-point viral load (SPVL) is predictive of disease progression and shows variability across HIV-1-positive (HIV+) persons. Various factors may influence SPVL including viral features, environmental exposure and host genetics. To identify single nucleotide polymorphisms (SNPs) associated with SPVL, we performed a genome-wide association study (GWAS) on a subset of participants from the Strategic Timing of AntiRetroviral Treatment (START) study covering a demographically diverse population. Methods. Consenting participants were antiretroviral therapy (ART)-naïve and SPVL was taken as log10(HIV RNA) at study entry. Genotypic data were generated on a custom content Affymetrix Axiom SNP array covering 770,558 probes. The Ensembl Gene database, assembly GRCh37.p13, was used for annotation. Principal component analysis (PCA) was used to identify population structures, and analysis of variance (ANOVA) was performed to detect associations between SNPs and SPVL. SNPs with zero variance or minor allele frequency (MAF) ≤0.05 were removed. Results. Among the 2,544 participants, PCA showed distinct population structures with strong separation between black (n = 578) and nonblack (n = 1,966) participants, Figure 1. ANOVA was performed independently on both subsets. Two SNPs located in the Major Histocompatibility Complex (MHC) class I region of chromosome six reached genome-wide significance (P < 5 × 10–8) in the non-black population: rs4418214 (P = 1.74 × 10-10), and rs57989216 (P = 3.96 × 10–8), Figure 2. Two additional SNPs, rs9264942 (P = 5.99 × 10–8) and rs7356880 (P = 9.69 × 10–8), in the same region approached significance. The minor alleles of all four SNPs were associated with lower SPVL, Figure 3. While no SNPs reached genome-wide significance in the black group, we observed similar trends toward lower SPVL for both rs4418214 and rs57989216. Conclusion. In this study we confirm the association of a previously reported SNP (rs4418214) and identify a novel candidate SNP (rs57989216) associated with lower SPVL in a population of nonblack, ART-naïve HIV+ persons. Current findings suggest that the effects of these SNPs are consistent across race groups, but further studies are required to confirm this. Our results support previous findings that variation in the MHC class I region is a major host determinant of HIV-1 control. Disclosures D. D. Murray, Centre of Excellence for Health, Immunity and Infectious diseases (CHIP), Department of Infectious Diseases, Rigshospitalet, University of Copenhagen, Copenhagen, Denmark: Employee, Salary. J. M. Molina, Gilead: Scientific Advisor, Consulting fee. Merck: Scientific Advisor, Consulting fee. ViiV: Scientific Advisor, Consulting fee. Teva: Scientific Advisor, Consulting fee.

DAVID bioinformatics resources consists of an integrated biological knowledgebase and analytic tools aimed at systematically extracting biological meaning from large gene/protein lists. This protocol explains how to use DAVID, a high-throughput and integrated data-mining environment, to analyze gene lists derived from high-throughput genomic experiments. The procedure first requires uploading a gene list containing any number of common gene identifiers followed by analysis using one or more text and pathway-mining tools such as gene functional classification, functional annotation chart or clustering and functional annotation table. By following this protocol, investigators are able to gain an in-depth understanding of the biological themes in lists of genes that are enriched in genome-scale studies.

Functional analysis of large gene lists, derived in most cases from emerging high-throughput genomic, proteomic and bioinformatics scanning approaches, is still a challenging and daunting task. The gene-annotation enrichment analysis is a promising high-throughput strategy that increases the likelihood for investigators to identify biological processes most pertinent to their study. Approximately 68 bioinformatics enrichment tools that are currently available in the community are collected in this survey. Tools are uniquely categorized into three major classes, according to their underlying enrichment algorithms. The comprehensive collections, unique tool classifications and associated questions/issues will provide a more comprehensive and up-to-date view regarding the advantages, pitfalls and recent trends in a simpler tool-class level rather than by a tool-by-tool approach. Thus, the survey will help tool designers/developers and experienced end users understand the underlying algorithms and pertinent details of particular tool categories/tools, enabling them to make the best choices for their particular research interests.

Development of HIV-1 drug resistance mutations (HDRMs) is one of the major reasons for the clinical failure of antiretroviral therapy. Treatment success rates can be improved by applying personalized anti-HIV regimens based on a patient’s HDRM profile. However, the sensitivity and specificity of the HDRM profile is limited by the methods used for detection. Sanger-based sequencing technology has traditionally been used for determining HDRM profiles at the single nucleotide variant (SNV) level, but with a sensitivity of only ≥ 20% in the HIV population of a patient. Next Generation Sequencing (NGS) technologies offer greater detection sensitivity (~ 1%) and larger scope (hundreds of samples per run). However, NGS technologies produce reads that are too short to enable the detection of the physical linkages of individual SNVs across the haplotype of each HIV strain present. In this article, we demonstrate that the single-molecule long reads generated using the Third Generation Sequencer (TGS), PacBio RS II, along with the appropriate bioinformatics analysis method, can resolve the HDRM profile at a more advanced quasispecies level. The case studies on patients’ HIV samples showed that the quasispecies view produced using the PacBio method offered greater detection sensitivity and was more comprehensive for understanding HDRM situations, which is complement to both Sanger and NGS technologies. In conclusion, the PacBio method, providing a promising new quasispecies level of HDRM profiling, may effect an important change in the field of HIV drug resistance research.

UNLABELLED Our current biological knowledge is spread over many independent bioinformatics databases where many different types of gene and protein identifiers are used. The heterogeneous and redundant nature of these identifiers limits data analysis across different bioinformatics resources. It is an even more serious bottleneck of data analysis for larger datasets, such as gene lists derived from microarray and proteomic experiments. The DAVID Gene ID Conversion Tool (DICT), a web-based application, is able to convert user's input gene or gene product identifiers from one type to another in a more comprehensive and high-throughput manner with a uniquely enhanced ID-ID mapping database. AVAILABILITY http://david.abcc.ncifcrf.gov/conversion.jsp.

PacBio RS, a newly emerging third-generation DNA sequencing platform, is based on a real-time, single-molecule, nano-nitch sequencing technology that can generate very long reads (up to 20-kb) in contrast to the shorter reads produced by the first and second generation sequencing technologies. As a new platform, it is important to assess the sequencing error rate, as well as the quality control (QC) parameters associated with the PacBio sequence data. In this study, a mixture of 10 prior known, closely related DNA amplicons were sequenced using the PacBio RS sequencing platform. After aligning Circular Consensus Sequence (CCS) reads derived from the above sequencing experiment to the known reference sequences, we found that the median error rate was 2.5% without read QC, and improved to 1.3% with an SVM based multi-parameter QC method. In addition, a De Novo assembly was used as a downstream application to evaluate the effects of different QC approaches. This benchmark study indicates that even though CCS reads are post error-corrected it is still necessary to perform appropriate QC on CCS reads in order to produce successful downstream bioinformatics analytical results.

Background Gene duplication is a normal evolutionary process. If there is no selective advantage in keeping the duplicated gene, it is usually reduced to a pseudogene and disappears from the genome. However, some paralogs are retained. These gene products are likely to be beneficial to the organism, e.g. in adaptation to new environmental conditions. The aim of our analysis is to investigate the properties of paralog-forming genes in prokaryotes, and to analyse the role of these retained paralogs by relating gene properties to life style of the corresponding prokaryotes. Results Paralogs were identified in a number of prokaryotes, and these paralogs were compared to singletons of persistent orthologs based on functional classification. This showed that the paralogs were associated with for example energy production, cell motility, ion transport, and defence mechanisms. A statistical overrepresentation analysis of gene and protein annotations was based on paralogs of the 200 prokaryotes with the highest fraction of paralog-forming genes. Biclustering of overrepresented gene ontology terms versus species was used to identify clusters of properties associated with clusters of species. The clusters were classified using similarity scores on properties and species to identify interesting clusters, and a subset of clusters were analysed by comparison to literature data. This analysis showed that paralogs often are associated with properties that are important for survival and proliferation of the specific organisms. This includes processes like ion transport, locomotion, chemotaxis and photosynthesis. However, the analysis also showed that the gene ontology terms sometimes were too general, imprecise or even misleading for automatic analysis. Conclusions Properties described by gene ontology terms identified in the overrepresentation analysis are often consistent with individual prokaryote lifestyles and are likely to give a competitive advantage to the organism. Paralogs and singletons dominate different categories of functional classification, where paralogs in particular seem to be associated with processes involving interaction with the environment.

Background Many species form extraembryonic tissues during embryogenesis, such as the placenta of humans and other viviparous mammals. Extraembryonic tissues have various roles in protecting, nourishing and patterning embryos. Prior to gastrulation in zebrafish, the yolk syncytial layer - an extraembryonic nuclear syncytium - produces signals that induce mesoderm and endoderm formation. Mesoderm and endoderm precursor cells are situated in the embryonic margin, an external ring of cells along the embryo-yolk interface. The yolk syncytial layer initially forms below the margin, in a domain called the external yolk syncytial layer (E-YSL). Results We hypothesize that key components of the yolk syncytial layer's mesoderm and endoderm inducing activity are expressed as mRNAs in the E-YSL. To identify genes expressed in the E-YSL, we used microarrays to compare the transcription profiles of intact pre-gastrula embryos with pre-gastrula embryonic cells that we had separated from the yolk and yolk syncytial layer. This identified a cohort of genes with enriched expression in intact embryos. Here we describe our whole mount in situ hybridization analysis of sixty-eight of them. This includes ten genes with E-YSL expression (camsap1l1, gata3, znf503, hnf1ba, slc26a1, slc40a1, gata6, gpr137bb, otop1 and cebpa), four genes with expression in the enveloping layer (EVL), a superficial epithelium that protects the embryo (zgc:136817, zgc:152778, slc14a2 and elovl6l), three EVL genes whose expression is transiently confined to the animal pole (elovl6l, zgc:136359 and clica), and six genes with transient maternal expression (mtf1, wu:fj59f04, mospd2, rftn2, arrdc1a and pho). We also assessed the requirement of Nodal signaling for the expression of selected genes in the E-YSL, EVL and margin. Margin expression was Nodal dependent for all genes we tested, including the concentrated margin expression of an EVL gene: zgc:110712. All other instances of EVL and E-YSL expression that we tested were Nodal independent. Conclusion We have devised an effective strategy for enriching and identifying genes expressed in the E-YSL of pre-gastrula embryos. To our surprise, maternal genes and genes expressed in the EVL were also enriched by this strategy. A number of these genes are promising candidates for future functional studies on early embryonic patterning.

High-throughput genomics screening studies, such as microarray, proteomics, etc., often result in large, "interesting" gene lists, ranging in size from hundreds to thousands of genes. Given the challenges of functionally interpreting such large gene lists, it is necessary to incorporate bioinformatics tools in the analysis. DAVID is a Web-based application that provides a high-throughput and integrative gene functional annotation environment to systematically extract biological themes behind large gene lists. High-throughput gene functional analysis with DAVID will provide important insights that allow investigators to understand the biological themes within their given genomic study. This unit will describe step-by-step procedures to use DAVID tools, as well as a brief rationale and key parameters in the DAVID analysis.

Background Due to the complex and distributed nature of biological research, our current biological knowledge is spread over many redundant annotation databases maintained by many independent groups. Analysts usually need to visit many of these bioinformatics databases in order to integrate comprehensive annotation information for their genes, which becomes one of the bottlenecks, particularly for the analytic task associated with a large gene list. Thus, a highly centralized and ready-to-use gene-annotation knowledgebase is in demand for high throughput gene functional analysis. Description The DAVID Knowledgebase is built around the DAVID Gene Concept, a single-linkage method to agglomerate tens of millions of gene/protein identifiers from a variety of public genomic resources into DAVID gene clusters. The grouping of such identifiers improves the cross-reference capability, particularly across NCBI and UniProt systems, enabling more than 40 publicly available functional annotation sources to be comprehensively integrated and centralized by the DAVID gene clusters. The simple, pair-wise, text format files which make up the DAVID Knowledgebase are freely downloadable for various data analysis uses. In addition, a well organized web interface allows users to query different types of heterogeneous annotations in a high-throughput manner. Conclusion The DAVID Knowledgebase is designed to facilitate high throughput gene functional analysis. For a given gene list, it not only provides the quick accessibility to a wide range of heterogeneous annotation data in a centralized location, but also enriches the level of biological information for an individual gene. Moreover, the entire DAVID Knowledgebase is freely downloadable or searchable at http://david.abcc.ncifcrf.gov/knowledgebase/.

DAVID, a web-accessible program that integrates functional genomic annotations with intuitive graphical summaries, has been described and will assist in the interpretation of genome-scale datasets by facilitating the transition from data collection to biological meaning., The distributed nature of biological knowledge poses a major challenge to the interpretation of genome-scale datasets, including those derived from microarray and proteomic studies. This report describes DAVID, a web-accessible program that integrates functional genomic annotations with intuitive graphical summaries. Lists of gene or protein identifiers are rapidly annotated and summarized according to shared categorical data for Gene Ontology, protein domain, and biochemical pathway membership. DAVID assists in the interpretation of genome-scale datasets by facilitating the transition from data collection to biological meaning.

Summary: The database for annotation, visualization and integrated discovery (DAVID), which can be freely accessed at http://david.abcc.ncifcrf.gov/, is a web-based online bioinformatics resource that aims to provide tools for the functional interpretation of large lists of genes/proteins. It has been used by researchers from more than 5000 institutes worldwide, with a daily submission rate of ∼1200 gene lists from ∼400 unique researchers, and has been cited by more than 6000 scientific publications. However, the current web interface does not support programmatic access to DAVID, and the uniform resource locator (URL)-based application programming interface (API) has a limit on URL size and is stateless in nature as it uses URL request and response messages to communicate with the server, without keeping any state-related details. DAVID-WS (web service) has been developed to automate user tasks by providing stateful web services to access DAVID programmatically without the need for human interactions. Availability: The web service and sample clients (written in Java, Perl, Python and Matlab) are made freely available under the DAVID License at http://david.abcc.ncifcrf.gov/content.jsp?file=WS.html. Contact: xiaoli.jiao@nih.gov; rlempicki@nih.gov

Background Microarray experimentation requires the application of complex analysis methods as well as the use of non-trivial computer technologies to manage the resultant large data sets. This, together with the proliferation of tools and techniques for microarray data analysis, makes it very challenging for a laboratory scientist to keep up-to-date with the latest developments in this field. Our aim was to develop a distributed e-support system for microarray data analysis and management. Results EMAAS (Extensible MicroArray Analysis System) is a multi-user rich internet application (RIA) providing simple, robust access to up-to-date resources for microarray data storage and analysis, combined with integrated tools to optimise real time user support and training. The system leverages the power of distributed computing to perform microarray analyses, and provides seamless access to resources located at various remote facilities. The EMAAS framework allows users to import microarray data from several sources to an underlying database, to pre-process, quality assess and analyse the data, to perform functional analyses, and to track data analysis steps, all through a single easy to use web portal. This interface offers distance support to users both in the form of video tutorials and via live screen feeds using the web conferencing tool EVO. A number of analysis packages, including R-Bioconductor and Affymetrix Power Tools have been integrated on the server side and are available programmatically through the Postgres-PLR library or on grid compute clusters. Integrated distributed resources include the functional annotation tool DAVID, GeneCards and the microarray data repositories GEO, CELSIUS and MiMiR. EMAAS currently supports analysis of Affymetrix 3' and Exon expression arrays, and the system is extensible to cater for other microarray and transcriptomic platforms. Conclusion EMAAS enables users to track and perform microarray data management and analysis tasks through a single easy-to-use web application. The system architecture is flexible and scalable to allow new array types, analysis algorithms and tools to be added with relative ease and to cope with large increases in data volume.

EASE is a customizable software application for rapid biological interpretation of gene lists that result from the analysis of microarray, proteomics, SAGE and other high-throughput genomic data., EASE is a customizable software application for rapid biological interpretation of gene lists that result from the analysis of microarray, proteomics, SAGE and other high-throughput genomic data. The biological themes returned by EASE recapitulate manually determined themes in previously published gene lists and are robust to varying methods of normalization, intensity calculation and statistical selection of genes. EASE is a powerful tool for rapidly converting the results of functional genomics studies from 'genes' to 'themes'.

Human pygmy populations inhabit different regions of the world, from Africa to Melanesia. In Asia, short-statured populations are often referred to as "negritos." Their short stature has been interpreted as a consequence of thermoregulatory, nutritional, and/or locomotory adaptations to life in tropical forests. A more recent hypothesis proposes that their stature is the outcome of a life history trade-off in high-mortality environments, where early reproduction is favored and, consequently, early sexual maturation and early growth cessation have coevolved. Some serological evidence of deficiencies in the growth hormone/insulin-like growth factor axis have been previously associated with pygmies' short stature. Using genome-wide single-nucleotide polymorphism genotype data, we first tested whether different negrito groups living in the Philippines and Papua New Guinea are closely related and then investigated genomic signals of recent positive selection in African, Asian, and Papuan pygmy populations. We found that negritos in the Philippines and Papua New Guinea are genetically more similar to their nonpygmy neighbors than to one another and have experienced positive selection at different genes. These results indicate that geographically distant pygmy groups are likely to have evolved their short stature independently. We also found that selection on common height variants is unlikely to explain their short stature and that different genes associated with growth, thyroid function, and sexual development are under selection in different pygmy groups.

Background Gene isoforms are commonly found in both prokaryotes and eukaryotes. Since each isoform may perform a specific function in response to changing environmental conditions, studying the dynamics of gene isoforms is important in understanding biological processes and disease conditions. However, genome-wide identification of gene isoforms is technically challenging due to the high degree of sequence identity among isoforms. Traditional targeted sequencing approach, involving Sanger sequencing of plasmid-cloned PCR products, has low throughput and is very tedious and time-consuming. Next-generation sequencing technologies such as Illumina and 454 achieve high throughput but their short read lengths are a critical barrier to accurate assembly of highly similar gene isoforms, and may result in ambiguities and false joining during sequence assembly. More recently, the third generation sequencer represented by the PacBio platform offers sufficient throughput and long reads covering the full length of typical genes, thus providing a potential to reliably profile gene isoforms. However, the PacBio long reads are error-prone and cannot be effectively analyzed by traditional assembly programs. Results We present a clustering-based analysis pipeline integrated with PacBio sequencing data for profiling highly similar gene isoforms. This approach was first evaluated in comparison to de novo assembly of 454 reads using a benchmark admixture containing 10 known, cloned msg genes encoding the major surface glycoprotein of Pneumocystis jirovecii. All 10 msg isoforms were successfully reconstructed with the expected length (~1.5 kb) and correct sequence by the new approach, while 454 reads could not be correctly assembled using various assembly programs. When using an additional benchmark admixture containing 22 known P. jirovecii msg isoforms, this approach accurately reconstructed all but 4 these isoforms in their full-length (~3 kb); these 4 isoforms were present in low concentrations in the admixture. Finally, when applied to the original clinical sample from which the 22 known msg isoforms were cloned, this approach successfully identified not only all known isoforms accurately (~3 kb each) but also 48 novel isoforms. Conclusions PacBio sequencing integrated with the clustering-based analysis pipeline achieves high-throughput and high-resolution discrimination of highly similar sequences, and can serve as a new approach for genome-wide characterization of gene isoforms and other highly repetitive sequences. Electronic supplementary material The online version of this article (doi:10.1186/s13040-016-0090-8) contains supplementary material, which is available to authorized users.

The DAVID gene functional classification tool uses a novel fuzzy clustering algorithm to condense a list of genes or associated biological terms into organized classes of related genes or biology, called biological modules., The DAVID Gene Functional Classification Tool uses a novel agglomeration algorithm to condense a list of genes or associated biological terms into organized classes of related genes or biology, called biological modules. This organization is accomplished by mining the complex biological co-occurrences found in multiple sources of functional annotation. It is a powerful method to group functionally related genes and terms into a manageable number of biological modules for efficient interpretation of gene lists in a network context.