An atlas of combinatorial transcriptional regulation in mouse and man.

Combinatorial interactions among transcription factors are critical to directing tissue-specific gene expression. To build a global atlas of these combinations, we have screened for physical interactions among the majority of human and mouse DNA-binding transcription factors (TFs). The complete networks contain 762 human and 877 mouse interactions. Analysis of the networks reveals that highly connected TFs are broadly expressed across tissues, and that roughly half of the measured interactions are conserved between mouse and human. The data highlight the importance of TF combinations for determining cell fate, and they lead to the identification of a SMAD3/FLI1 complex expressed during development of immunity. The availability of large TF combinatorial networks in both human and mouse will provide many opportunities to study gene regulation, tissue differentiation, and mammalian evolution.

The association between single-nucleotide polymorphisms within type 1 interferon pathway genes and human immunodeficiency virus type 1 viral load in antiretroviral-naïve participants.

BACKGROUND: Human genetic contribution to HIV progression remains inadequately explained. The type 1 interferon (IFN) pathway is important for host control of HIV and variation in type 1 IFN genes may contribute to disease progression. This study assessed the impact of variations at the gene and pathway level of type 1 IFN on HIV-1 viral load (VL). METHODS: Two cohorts of antiretroviral (ART) naïve participants living with HIV (PLWH) with either early (START) or advanced infection (FIRST) were analysed separately. Type 1 IFN genes (n = 17) and receptor subunits (IFNAR1, IFNAR2) were examined for both cumulated type 1 IFN pathway analysis and individual gene analysis. SKAT-O was applied to detect associations between the genotype and HIV-1 study entry viral load (log10 transformed) as a proxy for set point VL; P-values were corrected using Bonferroni (P < 0.0025). RESULTS: The analyses among those with early infection included 2429 individuals from five continents. The median study entry HIV VL was 14,623 (IQR 3460-45100) copies/mL. Across 673 SNPs within 19 type 1 IFN genes, no significant association with study entry VL was detected. Conversely, examining individual genes in START showed a borderline significant association between IFNW1, and study entry VL (P = 0.0025). This significance remained after separate adjustments for age, CD4+ T-cell count, CD4+/CD8+ T-cell ratio and recent infection. When controlling for population structure using linear mixed effects models (LME), in addition to principal components used in the main model, this was no longer significant (p = 0.0244). In subgroup analyses stratified by geographical region, the association between IFNW1 and study entry VL was only observed among African participants, although, the association was not significant when controlling for population structure using LME. Of the 17 SNPs within the IFNW1 region, only rs79876898 (A > G) was associated with study entry VL (p = 0.0020, beta = 0.32; G associated with higher study entry VL than A) in single SNP association analyses. The findings were not reproduced in FIRST participants. CONCLUSION: Across 19 type 1 IFN genes, only IFNW1 was associated with HIV-1 study entry VL in a cohort of ART-naïve individuals in early stages of their infection, however, this was no longer significant in sensitivity analyses that controlled for population structures using LME.

Impact of Antibiotic Treatment on the Gut Microbiome and its Resistome in Hematopoietic Stem Cell Transplant Recipients.

Antibiotic-resistant bacterial infections are increasingly an issue in allogenic hematopoietic stem cell transplant patients. How antibiotic treatment impacts antibiotic resistance in the human gut microbiome remains poorly understood in vivo. Here, a total of 577 fecal samples from 233 heavily antibiotic-treated transplant patients were examined using high-resolution prescription data and shotgun metagenomics. The 13 most frequently used antibiotics were significantly associated with 154 (40% of tested associations) microbiome features. Use of broad-spectrum ß-lactam antibiotics was most markedly associated with microbial disruption and increase in resistome features. The enterococcal vanA gene was positively associated with 8 of the 13 antibiotics, and in particular piperacillin/tazobactam and vancomycin. Here, we highlight the need for a high-resolution approach in understanding the development of antibiotic resistance in the gut microbiome. Our findings can be used to inform antibiotic stewardship and combat the increasing threat of antibiotic resistance.

Readmissions, Postdischarge Mortality, and Sustained Recovery Among Patients Admitted to Hospital With Coronavirus Disease 2019 (COVID-19).

BACKGROUND: Many interventional in-patient coronavirus disease 2019 (COVID-19) trials assess primary outcomes through day 28 post-randomization. Since a proportion of patients experience protracted disease or relapse, such follow-up period may not fully capture the course of the disease, even when randomization occurs a few days after hospitalization. METHODS: Among adults hospitalized with COVID-19 in eastern Denmark from 18 March 2020-12 January 2021 we assessed all-cause mortality, recovery, and sustained recovery 90 days after admission, and readmission and all-cause mortality 90 days after discharge. Recovery was defined as hospital discharge and sustained recovery as recovery and alive without readmissions for 14 consecutive days. RESULTS: Among 3386 patients included in the study, 2796 (82.6%) reached recovery and 2600 (77.0%) achieved sustained recovery. Of those discharged from hospital, 556 (19.9%) were readmitted and 289 (10.3%) died. Overall, the median time to recovery was 6 days (interquartile range [IQR]: 3-10), and 19 days (IQR: 11-33) among patients in intensive care in the first 2 days of admission. CONCLUSIONS: Postdischarge readmission and mortality rates were substantial. Therefore, sustained recovery should be favored to recovery outcomes in clinical COVID-19 trials. A 28-day follow-up period may be too short for the critically ill.

Nutrient sensing modulates malaria parasite virulence.

The lifestyle of intracellular pathogens, such as malaria parasites, is intimately connected to that of their host, primarily for nutrient supply. Nutrients act not only as primary sources of energy but also as regulators of gene expression, metabolism and growth, through various signalling networks that enable cells to sense and adapt to varying environmental conditions. Canonical nutrient-sensing pathways are presumed to be absent from the causative agent of malaria, Plasmodium, thus raising the question of whether these parasites can sense and cope with fluctuations in host nutrient levels. Here we show that Plasmodium blood-stage parasites actively respond to host dietary calorie alterations through rearrangement of their transcriptome accompanied by substantial adjustment of their multiplication rate. A kinome analysis combined with chemical and genetic approaches identified KIN as a critical regulator that mediates sensing of nutrients and controls a transcriptional response to the host nutritional status. KIN shares homology with SNF1/AMPKα, and yeast complementation studies suggest that it is part of a functionally conserved cellular energy-sensing pathway. Overall, these findings reveal a key parasite nutrient-sensing mechanism that is critical for modulating parasite replication and virulence.

Metabolic Potential of the Gut Microbiome Is Significantly Impacted by Conditioning Regimen in Allogeneic Hematopoietic Stem Cell Transplantation Recipients.

Allogeneic hematopoietic stem cell transplantation (aHSCT) is a putative curative treatment for malignant hematologic disorders. During transplantation, the immune system is suppressed/eradicated through a conditioning regimen (non-myeloablative or myeloablative) and replaced with a donor immune system. In our previous study, we showed changes in gut taxonomic profiles and a decrease in bacterial diversity post-transplant. In this study, we expand the cohort with 114 patients and focus on the impact of the conditioning regimens on taxonomic features and the metabolic functions of the gut bacteria. This is, to our knowledge, the first study to examine the metabolic potential of the gut microbiome in this patient group. Adult aHSCT recipients with shotgun sequenced stool samples collected day -30 to +28 relative to aHSCT were included. One sample was selected per patient per period: pre-aHSCT (day -30-0) and post-aHSCT (day 1-28). In total, 254 patients and 365 samples were included. Species richness, alpha diversity, gene richness and metabolic richness were all lower post-aHSCT than pre-aHSCT and the decline was more pronounced for the myeloablative group. The myeloablative group showed a decline in 36 genera and an increase in 15 genera. For the non-myeloablative group, 30 genera decreased and 16 increased with lower fold changes than observed in the myeloablative group. For the myeloablative group, 32 bacterial metabolic functions decreased, and one function increased. For the non-myeloablative group, three functions decreased, and two functions increased. Hence, the changes in taxonomy post-aHSCT caused a profound decline in bacterial metabolic functions especially in the myeloablative group, thus providing new evidence for associations of myeloablative conditioning and gut dysbiosis from a functional perspective.

Clipped histone H3 is integrated into nucleosomes of DNA replication genes in the human malaria parasite Plasmodium falciparum.

Post-translational modifications of histone H3 N-terminal tails are key epigenetic regulators of virulence gene expression and sexual commitment in the human malaria parasite Plasmodium falciparum Here, we identify proteolytic clipping of the N-terminal tail of nucleosome-associated histone H3 at amino acid position 21 as a new chromatin modification. A cathepsin C-like proteolytic clipping activity is observed in nuclear parasite extracts. Notably, an ectopically expressed version of clipped histone H3, PfH3p-HA, is targeted to the nucleus and integrates into mononucleosomes. Furthermore, chromatin immunoprecipitation and next-generation sequencing analysis identified PfH3p-HA as being highly enriched in the upstream region of six genes that play a key role in DNA replication and repair: In these genes, PfH3p-HA demarcates a specific 1.5 kb chromatin island adjacent to the open reading frame. Our results indicate that, in P. falciparum, the process of histone clipping may precede chromatin integration hinting at preferential targeting of pre-assembled PfH3p-containing nucleosomes to specific genomic regions. The discovery of a protease-directed mode of chromatin organization in P. falciparum opens up new avenues to develop new anti-malarials.

Distinctive roles of age, sex, and genetics in shaping transcriptional variation of human immune responses to microbial challenges.

The contribution of host genetic and nongenetic factors to immunological differences in humans remains largely undefined. Here, we generated bacterial-, fungal-, and viral-induced immune transcriptional profiles in an age- and sex-balanced cohort of 1,000 healthy individuals and searched for the determinants of immune response variation. We found that age and sex affected the transcriptional response of most immune-related genes, with age effects being more stimulus-specific relative to sex effects, which were largely shared across conditions. Although specific cell populations mediated the effects of age and sex on gene expression, including CD8+ T cells for age and CD4+ T cells and monocytes for sex, we detected a direct effect of these intrinsic factors for the majority of immune genes. The mapping of expression quantitative trait loci (eQTLs) revealed that genetic factors had a stronger effect on immune gene regulation than age and sex, yet they affected a smaller number of genes. Importantly, we identified numerous genetic variants that manifested their regulatory effects exclusively on immune stimulation, including a Candida albicans-specific master regulator at the CR1 locus. These response eQTLs were enriched in disease-associated variants, particularly for autoimmune and inflammatory disorders, indicating that differences in disease risk may result from regulatory variants exerting their effects only in the presence of immune stress. Together, this study quantifies the respective effects of age, sex, genetics, and cellular heterogeneity on the interindividual variability of immune responses and constitutes a valuable resource for further exploration in the context of different infection risks or disease outcomes.

Association Between Single-Nucleotide Polymorphisms in HLA Alleles and Human Immunodeficiency Virus Type 1 Viral Load in Demographically Diverse, Antiretroviral Therapy-Naive Participants From the Strategic Timing of AntiRetroviral Treatment Trial.

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.

Identification of the centromeres of Leishmania major: revealing the hidden pieces.

Leishmania affects millions of people worldwide. Its genome undergoes constitutive mosaic aneuploidy, a type of genomic plasticity that may serve as an adaptive strategy to survive distinct host environments. We previously found high rates of asymmetric chromosome allotments during mitosis that lead to the generation of such ploidy. However, the underlying molecular events remain elusive. Centromeres and kinetochores most likely play a key role in this process, yet their identification has failed using classical methods. Our analysis of the unconventional kinetochore complex recently discovered in Trypanosoma brucei (KKTs) leads to the identification of a Leishmania KKT gene candidate (LmKKT1). The GFP-tagged LmKKT1 displays "kinetochore-like" dynamics of intranuclear localization throughout the cell cycle. By ChIP-Seq assay, one major peak per chromosome is revealed, covering a region of 4 ±2 kb. We find two largely conserved motifs mapping to 14 of 36 chromosomes while a higher density of retroposons are observed in 27 of 36 centromeres. The identification of centromeres and of a kinetochore component of Leishmania chromosomes opens avenues to explore their role in mosaic aneuploidy.

First efficient CRISPR-Cas9-mediated genome editing in Leishmania parasites.

Protozoan pathogens that cause leishmaniasis in humans are relatively refractory to genetic manipulation. In this work, we implemented the CRISPR-Cas9 system in Leishmania parasites and demonstrated its efficient use for genome editing. The Cas9 endonuclease was expressed under the control of the Dihydrofolate Reductase-Thymidylate Synthase (DHFR-TS) promoter and the single guide RNA was produced under the control of the U6snRNA promoter and terminator. As a proof of concept, we chose to knockout a tandemly repeated gene family, the paraflagellar rod-2 locus. We were able to obtain null mutants in a single round of transfection. In addition, we confirmed the absence of off-target editions by whole genome sequencing of two independent clones. Our work demonstrates that CRISPR-Cas9-mediated gene knockout represents a major improvement in comparison with existing methods. Beyond gene knockout, this genome editing tool opens avenues for a multitude of functional studies to speed up research on leishmaniasis.

IAA-Ala Resistant3, an evolutionarily conserved target of miR167, mediates Arabidopsis root architecture changes during high osmotic stress.

The functions of microRNAs and their target mRNAs in Arabidopsis thaliana development have been widely documented; however, roles of stress-responsive microRNAs and their targets are not as well understood. Using small RNA deep sequencing and ATH1 microarrays to profile mRNAs, we identified IAA-Ala Resistant3 (IAR3) as a new target of miR167a. As expected, IAR3 mRNA was cleaved at the miR167a complementary site and under high osmotic stress miR167a levels decreased, whereas IAR3 mRNA levels increased. IAR3 hydrolyzes an inactive form of auxin (indole-3-acetic acid [IAA]-alanine) and releases bioactive auxin (IAA), a central phytohormone for root development. In contrast with the wild type, iar3 mutants accumulated reduced IAA levels and did not display high osmotic stress-induced root architecture changes. Transgenic plants expressing a cleavage-resistant form of IAR3 mRNA accumulated high levels of IAR3 mRNAs and showed increased lateral root development compared with transgenic plants expressing wild-type IAR3. Expression of an inducible noncoding RNA to sequester miR167a by target mimicry led to an increase in IAR3 mRNA levels, further confirming the inverse relationship between the two partners. Sequence comparison revealed the miR167 target site on IAR3 mRNA is conserved in evolutionarily distant plant species. Finally, we showed that IAR3 is required for drought tolerance.

Development of a machine learning model for early prediction of plasma leakage in suspected dengue patients.

BACKGROUND: At least a third of dengue patients develop plasma leakage with increased risk of life-threatening complications. Predicting plasma leakage using laboratory parameters obtained in early infection as means of triaging patients for hospital admission is important for resource-limited settings. METHODS: A Sri Lankan cohort including 4,768 instances of clinical data from N = 877 patients (60.3% patients with confirmed dengue infection) recorded in the first 96 hours of fever was considered. After excluding incomplete instances, the dataset was randomly split into a development and a test set with 374 (70%) and 172 (30%) patients, respectively. From the development set, five most informative features were selected using the minimum description length (MDL) algorithm. Random forest and light gradient boosting machine (LightGBM) were used to develop a classification model using the development set based on nested cross validation. An ensemble of the learners via average stacking was used as the final model to predict plasma leakage. RESULTS: Lymphocyte count, haemoglobin, haematocrit, age, and aspartate aminotransferase were the most informative features to predict plasma leakage. The final model achieved the area under the receiver operating characteristics curve, AUC = 0.80 with positive predictive value, PPV = 76.9%, negative predictive value, NPV = 72.5%, specificity = 87.9%, and sensitivity = 54.8% on the test set. CONCLUSION: The early predictors of plasma leakage identified in this study are similar to those identified in several prior studies that used non-machine learning based methods. However, our observations strengthen the evidence base for these predictors by showing their relevance even when individual data points, missing data and non-linear associations were considered. Testing the model on different populations using these low-cost observations would identify further strengths and limitations of the presented model.

Association between ten-eleven methylcytosine dioxygenase 2 genetic variation and viral load in people with HIV.

INTRODUCTION: Identifying genetic factors that influence HIV-pathogenesis is critical for understanding disease pathways. Previous studies have suggested a role for the human gene ten-eleven methylcytosine dioxygenase 2 (TET2) in modulating HIV-pathogenesis. METHODS: We assessed whether genetic variation in TET2 was associated with markers of HIV-pathogenesis using both gene level and single nucleotide polymorphism (SNP) level association in 8512 HIV-positive persons across five clinical trial cohorts. RESULTS: Variation at both the gene and SNP-level of TET2 was found to be associated with levels of HIV viral load (HIV-VL) consistently in the two cohorts that recruited antiretroviral-naïve participants. The SNPs occurred in two clusters of high linkage disequilibrium (LD), one associated with high HIV-VL and the other low HIV-VL, and were predominantly found in Black participants. CONCLUSION: Genetic variation in TET2 was associated with HIV-VL in two large antiretroviral therapy (ART)-naive clinical trial cohorts. The role of TET2 in HIV-pathogenesis warrants further investigation.

Metabolic Profiling Early Post-Allogeneic Haematopoietic Cell Transplantation in the Context of CMV Infection.

Immune dysfunction resulting from allogeneic haematopoietic stem cell transplantation (aHSCT) predisposes one to an elevated risk of cytomegalovirus (CMV) infection. Changes in metabolism have been associated with adverse outcomes, and in this study, we explored the associations between metabolic profiles and post-transplantation CMV infection using plasma samples collected 7-33 days after aHSCT. We included 68 aHSCT recipients from Rigshospitalet, Denmark, 50% of whom experienced CMV infection between days 34-100 post-transplantation. First, we investigated whether 12 metabolites selected based on the literature were associated with an increased risk of post-transplantation CMV infection. Second, we conducted an exploratory network-based analysis of the complete metabolic and lipidomic profiles in relation to clinical phenotypes and biological pathways. Lower levels of trimethylamine N-oxide were associated with subsequent CMV infection (multivariable logistic regression: OR = 0.63; 95% CI = [0.41; 0.87]; p = 0.01). Explorative analysis revealed 12 clusters of metabolites or lipids, among which one was predictive of CMV infection, and the others were associated with conditioning regimens, age upon aHSCT, CMV serostatus, and/or sex. Our results provide evidence for an association between the metabolome and CMV infection post-aHSCT that is independent of known risk factors.

Associations of functional human leucocyte antigen class I groups with HIV viral load in a heterogeneous cohort.

OBJECTIVE: Human leucocyte antigen (HLA) class I alleles are the main host genetic factors involved in controlling HIV-1 viral load (VL). Nevertheless, HLA diversity has proven a significant challenge in association studies. We assessed how accounting for binding affinities of HLA class I alleles to HIV-1 peptides facilitate association testing of HLA with HIV-1 VL in a heterogeneous cohort. DESIGN: Cohort from the Strategic Timing of AntiRetroviral Treatment (START) study. METHODS: We imputed HLA class I alleles from host genetic data (2546 HIV+ participants) and sampled immunopeptidomes from 2079 host-paired viral genomes (targeted amplicon sequencing). We predicted HLA class I binding affinities to HIV-1 and unspecific peptides, grouping alleles into functional clusters through consensus clustering. These functional HLA class I clusters were used to test associations with HIV VL. RESULTS: We identified four clades totaling 30 HLA alleles accounting for 11.4% variability in VL. We highlight HLA-B∗57:01 and B∗57:03 as functionally similar but yet overrepresented in distinct ethnic groups, showing when combined a protective association with HIV+ VL (log, ß -0.25; adj. P-value < 0.05). We further demonstrate only a slight power reduction when using unspecific immunopeptidomes, facilitating the use of the inferred functional HLA groups in other studies. CONCLUSION: The outlined computational approach provides a robust and efficient way to incorporate HLA function and peptide diversity, aiding clinical association studies in heterogeneous cohorts. To facilitate access to the proposed methods and results we provide an interactive application for exploring data.

Predicting type 2 diabetes risk before and after solid organ transplantation using polygenic scores in a Danish cohort.

Type 2 diabetes mellitus (T2DM) can be multifactorial where both genetics and environmental factors play a role. We aimed to investigate the use of polygenic risk scores (PRS) in the prediction of pre-transplant T2DM and post-transplant diabetes mellitus (PTDM) among solid organ transplant (SOT) patients. Using non-genetic risk scores alone; and the combination with PRS, separate logistic regression models were built and compared using receiver operator curves. Patients were assessed pre-transplant and in three post-transplant periods: 0-45, 46-365 and >365 days. A higher PRS was significantly associated with increased odds of pre-transplant T2DM. However, no improvement was observed for pre-transplant T2DM prediction when comparing PRS combined with non-genetic risk scores to using non-genetic risk scores alone. This was also true for predictions of PTDM in all three post-transplant periods. This study demonstrated that polygenic risk was only associated with the risk of T2DM among SOT recipients prior to transplant and not for PTDM. Combining PRS with a clinical model of non-genetic risk scores did not significantly improve the predictive ability, indicating its limited clinical utility in identifying patients at high risk for T2DM before transplantation, suggesting that non-genetic or different genetic factors may contribute to PTDM.

Erratum: Author Correction: Early stimulated immune responses predict clinical disease severity in hospitalized COVID-19 patients.

[This corrects the article DOI: 10.1038/s43856-022-00178-5.].

The response and recovery of the Arabidopsis thaliana transcriptome to phosphate starvation.

BACKGROUND: Over application of phosphate fertilizers in modern agriculture contaminates waterways and disrupts natural ecosystems. Nevertheless, this is a common practice among farmers, especially in developing countries as abundant fertilizers are believed to boost crop yields. The study of plant phosphate metabolism and its underlying genetic pathways is key to discovering methods of efficient fertilizer usage. The work presented here describes a genome-wide resource on the molecular dynamics underpinning the response and recovery in roots and shoots of Arabidopsis thaliana to phosphate-starvation. RESULTS: Genome-wide profiling by micro- and tiling-arrays (accessible from GEO: GSE34004) revealed minimal overlap between root and shoot transcriptomes suggesting two independent phosphate-starvation regulons. Novel gene expression patterns were detected for over 1000 candidates and were classified as either initial, persistent, or latent responders. Comparative analysis to AtGenExpress identified cohorts of genes co-regulated across multiple stimuli. The hormone ABA displayed a dominant role in regulating many phosphate-responsive candidates. Analysis of co-regulation enabled the determination of specific versus generic members of closely related gene families with respect to phosphate-starvation. Thus, among others, we showed that PHR1-regulated members of closely related phosphate-responsive families (PHT1;1, PHT1;7-9, SPX1-3, and PHO1;H1) display greater specificity to phosphate-starvation than their more generic counterparts. CONCLUSION: Our results uncover much larger, staged responses to phosphate-starvation than previously described. To our knowledge, this work describes the most complete genome-wide data on plant nutrient stress to-date.