Bioconductor: open software development for computational biology and bioinformatics

The Bioconductor project is an initiative for the collaborative creation of extensible software for computational biology and bioinformatics. The goals of the project include: fostering collaborative development and widespread use of innovative software, reducing barriers to entry into interdisciplinary scientific research, and promoting the achievement of remote reproducibility of research results. We describe details of our aims and methods, identify current challenges, compare Bioconductor to other open bioinformatics projects, and provide working examples.

Hallucinogens, including mescaline, psilocybin, and lysergic acid diethylamide (LSD), profoundly affect perception, cognition, and mood. All known drugs of this class are 5-HT(2A) receptor (2AR) agonists, yet closely related 2AR agonists such as lisuride lack comparable psychoactive properties. Why only certain 2AR agonists are hallucinogens and which neural circuits mediate their effects are poorly understood. By genetically expressing 2AR only in cortex, we show that 2AR-regulated pathways on cortical neurons are sufficient to mediate the signaling pattern and behavioral response to hallucinogens. Hallucinogenic and nonhallucinogenic 2AR agonists both regulate signaling in the same 2AR-expressing cortical neurons. However, the signaling and behavioral responses to the hallucinogens are distinct. While lisuride and LSD both act at 2AR expressed by cortex neurons to regulate phospholipase C, LSD responses also involve pertussis toxin-sensitive heterotrimeric G(i/o) proteins and Src. These studies identify the long-elusive neural and signaling mechanisms responsible for the unique effects of hallucinogens.

SUMMARY The embryonic stem (ES) cell transcriptional and epigenetic networks are critical for the maintenance of ES cell self-renewal. However, it remains unclear whether components of these networks functionally interact and if so, what factors mediate such interactions. Here we show that WD-repeat protein-5 (Wdr5), a core member of the mammalian Trithorax (trxG) complex, positively correlates with the undifferentiated state and is a novel regulator of ES cell self-renewal. We demonstrate that Wdr5, an ‘effector’ of H3K4 methylation, interacts with the pluripotency transcription factor Oct4. Genome-wide protein localization and transcriptome analyses demonstrate overlapping gene regulatory functions between Oct4 and Wdr5. We show that the Oct4-Sox2-Nanog circuitry and trxG cooperate in activating transcription of key self-renewal regulators. Furthermore, Wdr5 expression is required for the efficient formation of induced pluripotent stem (iPS) cells. We propose an integrated model of transcriptional and epigenetic control, mediated by select trxG members, for maintenance of ES cell self-renewal and somatic cell reprogramming.

multiple testing, family-wise error rate, false discovery rate, adjustedp-value, fast algorithm, minP, microarray, 62J15, 62G09, 62P10,

Advanced paternal age (APA) has been shown to be a significant risk factor in the offspring for neurodevelopmental psychiatric disorders, such as schizophrenia and autism spectrum disorders. During aging, de novo mutations accumulate in the male germline and are frequently transmitted to the offspring with deleterious effects. In addition, DNA methylation during spermatogenesis is an active process, which is susceptible to errors that can be propagated to subsequent generations. Here we test the hypothesis that the integrity of germline DNA methylation is compromised during the aging process. A genome-wide DNA methylation screen comparing sperm from young and old mice revealed a significant loss of methylation in the older mice in regions associated with transcriptional regulation. The offspring of older fathers had reduced exploratory and startle behaviors and exhibited similar brain DNA methylation abnormalities as observed in the paternal sperm. Offspring from old fathers also had transcriptional dysregulation of developmental genes implicated in autism and schizophrenia. Our findings demonstrate that DNA methylation abnormalities arising in the sperm of old fathers are a plausible mechanism to explain some of the risks that APA poses to resulting offspring.

Summary Background Chronic injury in kidney transplants remains a major cause of allograft loss. The aim of this study was to identify a gene set capable of predicting renal allografts at risk of progressive injury due to fibrosis. Methods This Genomics of Chronic Allograft Rejection (GoCAR) study is a prospective, multicentre study. We prospectively collected biopsies from renal allograft recipients (n=204) with stable renal function 3 months after transplantation. We used microarray analysis to investigate gene expression in 159 of these tissue samples. We aimed to identify genes that correlated with the Chronic Allograft Damage Index (CADI) score at 12 months, but not fibrosis at the time of the biopsy. We applied a penalised regression model in combination with permutation-based approach to derive an optimal gene set to predict allograft fibrosis. The GoCAR study is registered with ClinicalTrials.gov, number NCT00611702. Findings We identified a set of 13 genes that was independently predictive for the development of fibrosis at 1 year (ie, CADI-12 ≥2). The gene set had high predictive capacity (area under the curve [AUC] 0·967), which was superior to that of baseline clinical variables (AUC 0·706) and clinical and pathological variables (AUC 0·806). Furthermore routine pathological variables were unable to identify which histologically normal allografts would progress to fibrosis (AUC 0·754), whereas the predictive gene set accurately discriminated between transplants at high and low risk of progression (AUC 0·916). The 13 genes also accurately predicted early allograft loss (AUC 0·842 at 2 years and 0·844 at 3 years). We validated the predictive value of this gene set in an independent cohort from the GoCAR study (n=45, AUC 0·866) and two independent, publically available expression datasets (n=282, AUC 0·831 and n=24, AUC 0·972). Interpretation Our results suggest that this set of 13 genes could be used to identify kidney transplant recipients at risk of allograft loss before the development of irreversible damage, thus allowing therapy to be modified to prevent progression to fibrosis. Funding National Institutes of Health.

Heterotrimeric guanine nucleotide–binding protein (G protein)–coupled receptors (GPCRs) can form multiprotein complexes (heteromers), which can alter the pharmacology and functions of the constituent receptors. Previous findings demonstrated that the Gq/11-coupled serotonin 5-HT2A receptor and the Gi/o-coupled metabotropic glutamate 2 (mGlu2) receptor—GPCRs that are involved in signaling alterations associated with psychosis—assemble into a heteromeric complex in the mammalian brain. In single-cell experiments with various mutant versions of the mGlu2 receptor, we showed that stimulation of cells expressing mGlu2–5-HT2A heteromers with an mGlu2 agonist led to activation of Gq/11 proteins by the 5-HT2A receptors. For this crosstalk to occur, one of the mGlu2 subunits had to couple to Gi/o proteins, and we determined the relative location of the Gi/o-contacting subunit within the mGlu2 homodimer of the heteromeric complex. Additionally, mGlu2-dependent activation of Gq/11, but not Gi/o, was reduced in the frontal cortex of 5-HT2A knockout mice and was reduced in the frontal cortex of postmortem brains from schizophrenic patients. These findings offer structural insights into this important target in molecular psychiatry.

Background Whether young adults who are infected with SARS-CoV-2 are at risk of subsequent infection is uncertain. We investigated the risk of subsequent SARS-CoV-2 infection among young adults seropositive for a previous infection. Methods This analysis was performed as part of the prospective COVID-19 Health Action Response for Marines study (CHARM). CHARM included predominantly male US Marine recruits, aged 18–20 years, following a 2-week unsupervised quarantine at home. After the home quarantine period, upon arrival at a Marine-supervised 2-week quarantine facility (college campus or hotel), participants were enrolled and were assessed for baseline SARS-CoV-2 IgG seropositivity, defined as a dilution of 1:150 or more on receptor-binding domain and full-length spike protein ELISA. Participants also completed a questionnaire consisting of demographic information, risk factors, reporting of 14 specific COVID-19-related symptoms or any other unspecified symptom, and brief medical history. SARS-CoV-2 infection was assessed by PCR at weeks 0, 1, and 2 of quarantine and participants completed a follow-up questionnaire, which included questions about the same COVID-19-related symptoms since the last study visit. Participants were excluded at this stage if they had a positive PCR test during quarantine. Participants who had three negative swab PCR results during quarantine and a baseline serum serology test at the beginning of the supervised quarantine that identified them as seronegative or seropositive for SARS-CoV-2 then went on to basic training at Marine Corps Recruit Depot—Parris Island. Three PCR tests were done at weeks 2, 4, and 6 in both seropositive and seronegative groups, along with the follow-up symptom questionnaire and baseline neutralising antibody titres on all subsequently infected seropositive and selected seropositive uninfected participants (prospective study period). Findings Between May 11, 2020, and Nov 2, 2020, we enrolled 3249 participants, of whom 3168 (98%) continued into the 2-week quarantine period. 3076 (95%) participants, 2825 (92%) of whom were men, were then followed up during the prospective study period after quarantine for 6 weeks. Among 189 seropositive participants, 19 (10%) had at least one positive PCR test for SARS-CoV-2 during the 6-week follow-up (1·1 cases per person-year). In contrast, 1079 (48%) of 2247 seronegative participants tested positive (6·2 cases per person-year). The incidence rate ratio was 0·18 (95% CI 0·11–0·28; p<0·001). Among seropositive recruits, infection was more likely with lower baseline full-length spike protein IgG titres than in those with higher baseline full-length spike protein IgG titres (hazard ratio 0·45 [95% CI 0·32–0·65]; p<0·001). Infected seropositive participants had viral loads that were about 10-times lower than those of infected seronegative participants (ORF1ab gene cycle threshold difference 3·95 [95% CI 1·23–6·67]; p=0·004). Among seropositive participant...