NF-κB c-Rel Is Crucial for the Regulatory T Cell Immune Checkpoint in Cancer.

Regulatory T cells (Tregs) play a pivotal role in the inhibition of anti-tumor immune responses. Understanding the mechanisms governing Treg homeostasis may therefore be important for development of effective tumor immunotherapy. We have recently demonstrated a key role for the canonical nuclear factor κB (NF-κB) subunits, p65 and c-Rel, in Treg identity and function. In this report, we show that NF-κB c-Rel ablation specifically impairs the generation and maintenance of the activated Treg (aTreg) subset, which is known to be enriched at sites of tumors. Using mouse models, we demonstrate that melanoma growth is drastically reduced in mice lacking c-Rel, but not p65, in Tregs. Moreover, chemical inhibition of c-Rel function delayed melanoma growth by impairing aTreg-mediated immunosuppression and potentiated the effects of anti-PD-1 immunotherapy. Our studies therefore establish inhibition of NF-κB c-Rel as a viable therapeutic approach for enhancing checkpoint-targeting immunotherapy protocols.

A mouse model of Salmonella typhi infection.

Salmonella spp. are gram-negative flagellated bacteria that can cause food- and waterborne gastroenteritis and typhoid fever in humans. We now report that flagellin from Salmonella spp. is recognized in mouse intestine by Toll-like receptor 11 (TLR11). Absence of TLR11 renders mice more susceptible to infection by S. Typhimurium, with increased dissemination of the bacteria and enhanced lethality. Unlike S. Typhimurium, S. Typhi, a human obligatory pathogen that causes typhoid fever, is normally unable to infect mice. TLR11 is expressed in mice, but not in humans, and remarkably, we find that tlr11(-/-) mice are efficiently infected with orally administered S. Typhi. We also find that tlr11(-/-) mice can be immunized against S. Typhi. Therefore, tlr11(-/-) mice represent a small-animal model for the study of the immune response to S. Typhi and for the development of vaccines against this important human pathogen.

An NF-κB Transcription-Factor-Dependent Lineage-Specific Transcriptional Program Promotes Regulatory T Cell Identity and Function.

Both conventional T (Tconv) cells and regulatory T (Treg) cells are activated through ligation of the T cell receptor (TCR) complex, leading to the induction of the transcription factor NF-κB. In Tconv cells, NF-κB regulates expression of genes essential for T cell activation, proliferation, and function. However the role of NF-κB in Treg function remains unclear. We conditionally deleted canonical NF-κB members p65 and c-Rel in developing and mature Treg cells and found they have unique but partially redundant roles. c-Rel was critical for thymic Treg development while p65 was essential for mature Treg identity and maintenance of immune tolerance. Transcriptome and NF-κB p65 binding analyses demonstrated a lineage specific, NF-κB-dependent transcriptional program, enabled by enhanced chromatin accessibility. These dual roles of canonical NF-κB in Tconv and Treg cells highlight the functional plasticity of the NF-κB signaling pathway and underscores the need for more selective strategies to therapeutically target NF-κB.

A platform for curated products from novel open reading frames prompts reinterpretation of disease variants.

Recent evidence from proteomics and deep massively parallel sequencing studies have revealed that eukaryotic genomes contain substantial numbers of as-yet-uncharacterized open reading frames (ORFs). We define these uncharacterized ORFs as novel ORFs (nORFs). nORFs in humans are mostly under 100 codons and are found in diverse regions of the genome, including in long noncoding RNAs, pseudogenes, 3' UTRs, 5' UTRs, and alternative reading frames of canonical protein coding exons. There is therefore a pressing need to evaluate the potential functional importance of these unannotated transcripts and proteins in biological pathways and human disease on a larger scale, rather than one at a time. In this study, we outline the creation of a valuable nORFs data set with experimental evidence of translation for the community, use measures of heritability and selection that reveal signals for functional importance, and show the potential implications for functional interpretation of genetic variants in nORFs. Our results indicate that some variants that were previously classified as being benign or of uncertain significance may have to be reinterpreted.

Statistical sampling of missing environmental variables improves biophysical genomic prediction in wheat.

KEY MESSAGE: The integration of genomic prediction with crop growth models enabled the estimation of missing environmental variables which improved the prediction accuracy of grain yield. Since the invention of whole-genome prediction (WGP) more than two decades ago, breeding programmes have established extensive reference populations that are cultivated under diverse environmental conditions. The introduction of the CGM-WGP model, which integrates crop growth models (CGM) with WGP, has expanded the applications of WGP to the prediction of unphenotyped traits in untested environments, including future climates. However, CGMs require multiple seasonal environmental records, unlike WGP, which makes CGM-WGP less accurate when applied to historical reference populations that lack crucial environmental inputs. Here, we investigated the ability of CGM-WGP to approximate missing environmental variables to improve prediction accuracy. Two environmental variables in a wheat CGM, initial soil water content (InitlSoilWCont) and initial nitrate profile, were sampled from different normal distributions separately or jointly in each iteration within the CGM-WGP algorithm. Our results showed that sampling InitlSoilWCont alone gave the best results and improved the prediction accuracy of grain number by 0.07, yield by 0.06 and protein content by 0.03. When using the sampled InitlSoilWCont values as an input for the traditional CGM, the average narrow-sense heritability of the genotype-specific parameters (GSPs) improved by 0.05, with GNSlope, PreAnthRes, and VernSen showing the greatest improvements. Moreover, the root mean square of errors for grain number and yield was reduced by about 7% for CGM and 31% for CGM-WGP when using the sampled InitlSoilWCont values. Our results demonstrate the advantage of sampling missing environmental variables in CGM-WGP to improve prediction accuracy and increase the size of the reference population by enabling the utilisation of historical data that are missing environmental records.

Induction of innate immune memory via microRNA targeting of chromatin remodelling factors.

Prolonged exposure to microbial products such as lipopolysaccharide can induce a form of innate immune memory that blunts subsequent responses to unrelated pathogens, known as lipopolysaccharide tolerance. Sepsis is a dysregulated systemic immune response to disseminated infection that has a high mortality rate. In some patients, sepsis results in a period of immunosuppression (known as 'immunoparalysis')1 characterized by reduced inflammatory cytokine output2, increased secondary infection3 and an increased risk of organ failure and mortality4. Lipopolysaccharide tolerance recapitulates several key features of sepsis-associated immunosuppression5. Although various epigenetic changes have previously been observed in tolerized macrophages6-8, the molecular basis of tolerance, immunoparalysis and other forms of innate immune memory has remained unclear. Here we perform a screen for tolerance-associated microRNAs and identify miR-221 and miR-222 as regulators of the functional reprogramming of macrophages during lipopolysaccharide tolerization. Prolonged stimulation with lipopolysaccharide in mice leads to increased expression of miR-221 and mir-222, both of which regulate brahma-related gene 1 (Brg1, also known as Smarca4). This increased expression causes the transcriptional silencing of a subset of inflammatory genes that depend on chromatin remodelling mediated by SWI/SNF (switch/sucrose non-fermentable) and STAT (signal transducer and activator of transcription), which in turn promotes tolerance. In patients with sepsis, increased expression of miR-221 and miR-222 correlates with immunoparalysis and increased organ damage. Our results show that specific microRNAs can regulate macrophage tolerization and may serve as biomarkers of immunoparalysis and poor prognosis in patients with sepsis.

Crosstalk in NF-κB signaling pathways.

NF-κB transcription factors are critical regulators of immunity, stress responses, apoptosis and differentiation. A variety of stimuli coalesce on NF-κB activation, which can in turn mediate varied transcriptional programs. Consequently, NF-κB-dependent transcription is not only tightly controlled by positive and negative regulatory mechanisms but also closely coordinated with other signaling pathways. This intricate crosstalk is crucial to shaping the diverse biological functions of NF-κB into cell type- and context-specific responses.

Using genomic prediction with crop growth models enables the prediction of associated traits in wheat.

Crop growth models (CGM) can predict the performance of a cultivar in untested environments by sampling genotype-specific parameters. As they cannot predict the performance of new cultivars, it has been proposed to integrate CGMs with whole genome prediction (WGP) to combine the benefits of both models. Here, we used a CGM-WGP model to predict the performance of new wheat (Triticum aestivum) genotypes. The CGM was designed to predict phenology, nitrogen, and biomass traits. The CGM-WGP model simulated more heritable GSPs compared with the CGM and gave smaller errors for the observed phenotypes. The WGP model performed better when predicting yield, grain number, and grain protein content, but showed comparable performance to the CGM-WGP model for heading and physiological maturity dates. However, the CGM-WGP model was able to predict unobserved traits (for which there were no phenotypic records in the reference population). The CGM-WGP model also showed superior performance when predicting unrelated individuals that clustered separately from the reference population. Our results demonstrate new advantages for CGM-WGP modelling and suggest future efforts should focus on calibrating CGM-WGP models using high-throughput phenotypic measures that are cheaper and less laborious to collect.

Integrating biophysical crop growth models and whole genome prediction for their mutual benefit: a case study in wheat phenology.

Running crop growth models (CGM) coupled with whole genome prediction (WGP) as a CGM-WGP model introduces environmental information to WGP and genomic relatedness information to the genotype-specific parameters modelled through CGMs. Previous studies have primarily used CGM-WGP to infer prediction accuracy without exploring its potential to enhance CGM and WGP. Here, we implemented a heading and maturity date wheat phenology model within a CGM-WGP framework and compared it with CGM and WGP. The CGM-WGP resulted in more heritable genotype-specific parameters with more biologically realistic correlation structures between genotype-specific parameters and phenology traits compared with CGM-modelled genotype-specific parameters that reflected the correlation of measured phenotypes. Another advantage of CGM-WGP is the ability to infer accurate prediction with much smaller and less diverse reference data compared with that required for CGM. A genome-wide association analysis linked the genotype-specific parameters from the CGM-WGP model to nine significant phenology loci including Vrn-A1 and the three PPD1 genes, which were not detected for CGM-modelled genotype-specific parameters. Selection on genotype-specific parameters could be simpler than on observed phenotypes. For example, thermal time traits are theoretically more independent candidates, compared with the highly correlated heading and maturity dates, which could be used to achieve an environment-specific optimal flowering period. CGM-WGP combines the advantages of CGM and WGP to predict more accurate phenotypes for new genotypes under alternative or future environmental conditions.

PlmCas12e (CasX2) cleavage of CCR5: impact of guide RNA spacer length and PAM sequence on cleavage activity.

Gene editing using CRISPR/Cas (clustered regularly interspaced palindromic repeats/CRISPR-associated) is under development as a therapeutic tool for the modification of genes in eukaryotic cells. While much effort has focused on CRISPR/Cas9 systems from Streptococcus pyogenes and Staphylococcus aureus, alternative CRISPR systems have been identified from non-pathogenic microbes, including previously unknown class 2 systems, adding to a diverse toolbox of CRISPR/Cas enzymes. The Cas12e enzymes from non-pathogenic Deltaproteobacteria (CasX1, DpeCas12e) and Planctomycetes (CasX2, PlmCas12e) are smaller than Cas9, have a selective protospacer adjacent motif (PAM), and deliver a staggered cleavage cut with a 5-7 nucleotide overhang. We investigated the impact of guide RNA spacer length and alternative PAM sequences on cleavage activity to determine optimal conditions for PlmCas12e cleavage of the cellular gene CCR5 (CC-Chemokine receptor-5). CCR5 encodes the CCR5 coreceptor used by human immunodeficiency virus-type 1 (HIV-1) to infect target cells. A 32 base-pair deletion in CCR5 (CCR5-[Formula: see text]32) is responsible for HIV-1 resistance and reported cures following bone marrow transplantation. Consequently, CCR5 has been an important target for gene editing utilizing CRISPR/Cas. We determined that CCR5 cleavage activity varied with the target site, spacer length, and the fourth nucleotide in the previously described PAM sequence, TTCN. Our analyses demonstrated a PAM preference for purines (adenine, guanine) over pyrimidines (thymidine, cytosine) in the fourth position of the CasX2 PAM. This improved understanding of CasX2 cleavage requirements facilitates the development of therapeutic strategies to recreate the CCR5-[Formula: see text]32 mutation in haematopoietic stem cells.

PDK1 Is Required for Maintenance of CD4+ Foxp3+ Regulatory T Cell Function.

Regulatory T (Treg) cells have an essential role in maintaining immune homeostasis, in part by suppressing effector T cell functions. Phosphoinositide-dependent kinase 1 (PDK1) is a pleiotropic kinase that acts as a key effector downstream of PI3K in many cell types. In T cells, PDK1 has been shown to be critical for activation of NF-κB and AKT signaling upon TCR ligation and is therefore essential for effector T cell activation, proliferation, and cytokine production. Using Treg cell-specific conditional deletion, we now demonstrate that PDK1 is also essential for Treg cell suppressive activity in vivo. Ablation of Pdk1 specifically in Treg cells led to systemic, lethal, scurfy-like inflammation in mice. Genome-wide analysis confirmed that PDK1 is essential for the regulation of key Treg cell signature gene expression and, further, suggested that PDK1 acts primarily to control Treg cell gene expression through regulation of the canonical NF-κB pathway. Consistent with these results, the scurfy-like phenotype of mice lacking PDK1 in Treg cells was rescued by enforced activation of NF-κB downstream of PDK1. Therefore, PDK1-mediated activation of the NF-κB signaling pathway is essential for regulation of Treg cell signature gene expression and suppressor function.

Shared principles in NF-kappaB signaling.

The transcription factor NF-kappaB has served as a standard for inducible transcription factors for more than 20 years. The numerous stimuli that activate NF-kappaB, and the large number of genes regulated by NF-kappaB, ensure that this transcription factor is still the subject of intense research. Here, we attempt to synthesize some of the basic principles that have emerged from studies of NF-kappaB, and we aim to generate a more unified view of NF-kappaB regulation.

Physical Mapping of QTLs for Root Traits in a Population of Recombinant Inbred Lines of Hexaploid Wheat.

Root architecture is key in determining how effective plants are at intercepting and absorbing nutrients and water. Previously, the wheat (Triticum aestivum) cultivars Spica and Maringa were shown to have contrasting root morphologies. These cultivars were crossed to generate an F6:1 population of recombinant inbred lines (RILs) which was genotyped using a 90 K single nucleotide polymorphisms (SNP) chip. A total of 227 recombinant inbred lines (RILs) were grown in soil for 21 days in replicated trials under controlled conditions. At harvest, the plants were scored for seven root traits and two shoot traits. An average of 7.5 quantitative trait loci (QTL) were associated with each trait and, for each of these, physical locations of the flanking markers were identified using the Chinese Spring reference genome. We also compiled a list of genes from wheat and other monocotyledons that have previously been associated with root growth and morphology to determine their physical locations on the Chinese Spring reference genome. This allowed us to determine whether the QTL discovered in our study encompassed genes previously associated with root morphology in wheat or other monocotyledons. Furthermore, it allowed us to establish if the QTL were co-located with the QTL identified from previously published studies. The parental lines together with the genetic markers generated here will enable specific root traits to be introgressed into elite wheat lines. Moreover, the comprehensive list of genes associated with root development, and their physical locations, will be a useful resource for researchers investigating the genetics of root morphology in cereals.

Mining the Australian Grains Gene Bank for Rust Resistance in Barley.

Global barley production is threatened by plant pathogens, especially the rusts. In this study we used a targeted genotype-by-sequencing (GBS) assisted GWAS approach to identify rust resistance alleles in a collection of 287 genetically distinct diverse barley landraces and historical cultivars available in the Australian Grains Genebank (AGG) and originally sourced from Eastern Europe. The accessions were challenged with seven US-derived cereal rust pathogen races including Puccinia hordei (Ph-leaf rust) race 17VA12C, P. coronata var. hordei (Pch-crown rust) race 91NE9305 and five pathogenically diverse races of P. striiformis f. sp. hordei (Psh-stripe rust) (PSH-33, PSH-48, PSH-54, PSH-72 and PSH-100) and phenotyped quantitatively at the seedling stage. Novel resistance factors were identified on chromosomes 1H, 2H, 4H and 5H in response to Pch, whereas a race-specific QTL on 7HS was identified that was effective only to Psh isolates PSH-72 and PSH-100. A major effect QTL on chromosome 5HL conferred resistance to all Psh races including PSH-72, which is virulent on all 12 stripe rust differential tester lines. The same major effect QTL was also identified in response to leaf rust (17VA12C) suggesting this locus contains several pathogen specific rust resistance genes or the same gene is responsible for both leaf rust and stripe rust resistance. Twelve accessions were highly resistant to both leaf and stripe rust diseases and also carried the 5HL QTL. We subsequently surveyed the physical region at the 5HL locus for across the barley pan genome variation in the presence of known resistance gene candidates and identified a rich source of high confidence protein kinase and antifungal genes in the QTL region.

Ca2+-saturated calmodulin binds tightly to the N-terminal domain of A-type fibroblast growth factor homologous factors.

Voltage-gated sodium channels (Navs) are tightly regulated by multiple conserved auxiliary proteins, including the four fibroblast growth factor homologous factors (FGFs), which bind the Nav EF-hand like domain (EFL), and calmodulin (CaM), a multifunctional messenger protein that binds the NaV IQ motif. The EFL domain and IQ motif are contiguous regions of NaV cytosolic C-terminal domains (CTD), placing CaM and FGF in close proximity. However, whether the FGFs and CaM act independently, directly associate, or operate through allosteric interactions to regulate channel function is unknown. Titrations monitored by steady-state fluorescence spectroscopy, structural studies with solution NMR, and computational modeling demonstrated for the first time that both domains of (Ca2+)4-CaM (but not apo CaM) directly bind two sites in the N-terminal domain (NTD) of A-type FGF splice variants (FGF11A, FGF12A, FGF13A, and FGF14A) with high affinity. The weaker of the (Ca2+)4-CaM-binding sites was known via electrophysiology to have a role in long-term inactivation of the channel but not known to bind CaM. FGF12A binding to a complex of CaM associated with a fragment of the NaV1.2 CTD increased the Ca2+-binding affinity of both CaM domains, consistent with (Ca2+)4-CaM interacting preferentially with its higher-affinity site in the FGF12A NTD. Thus, A-type FGFs can compete with NaV IQ motifs for (Ca2+)4-CaM. During spikes in the cytosolic Ca2+ concentration that accompany an action potential, CaM may translocate from the NaV IQ motif to the FGF NTD, or the A-type FGF NTD may recruit a second molecule of CaM to the channel.

Adult plant stem rust resistance in durum wheat Glossy Huguenot: mapping, marker development and validation.

KEY MESSAGE: Adult plant stem rust resistance locus, QSrGH.cs-2AL, was identified in durum wheat Glossy Huguenot and mendelised as Sr63. Markers closely linked with Sr63 were developed. An F3 population from a Glossy Huguenot (GH)/Bansi cross used in a previous Australian study was advanced to F6 for molecular mapping of adult plant stem rust resistance. Maturity differences among F6 lines confounded assessments of stem rust response. GH was crossed with a stem rust susceptible F6 recombinant inbred line (RIL), GHB14 (M14), with similar maturity and an F6:7 population was developed through single seed descent method. F7 and F8 RILs were tested along with the parents at different locations. The F6 individual plants and both parents were genotyped using the 90 K single nucleotide polymorphism (SNP) wheat array. Stem rust resistance QTL on the long arms of chromosomes 1B (QSrGH.cs-1BL) and 2A (QSrGH.cs-2AL) were detected. QSrGH.cs-1BL and QSrGH.cs-2AL were both contributed by GH and explained 22% and 18% adult plant stem rust response variation, respectively, among GH/M14 RIL population. RILs carrying combinations of these QTL reduced more than 14% stem rust severity compared to those that possessed QSrGH.cs-1BL and QSrGH.cs-2AL individually. QSrGH.cs1BL was demonstrated to be the same as Sr58/Lr46/Yr29/Pm39 through marker genotyping. Lines lacking QSrGH.cs-1BL were used to Mendelise QSrGH.cs-2AL. Based on genomic locations of previously catalogued stem rust resistance genes and the QSrGH.cs-2AL map, it appeared to represent a new APR locus and was permanently named Sr63. SNP markers associated with Sr63 were converted to kompetetive allele-specific PCR (KASP) assays and were validated on a set of durum cultivars.

Recognition of profilin by Toll-like receptor 12 is critical for host resistance to Toxoplasma gondii.

Toll-like receptor 11 (TLR11) recognizes T. gondii profilin (TgPRF) and is required for interleukin-12 production and induction of immune responses that limit cyst burden in Toxoplasma gondii-infected mice. However, TLR11 only modestly affects survival of T. gondii-challenged mice. We report that TLR12, a previously uncharacterized TLR, also recognized TgPRF. TLR12 was sufficient for recognition of TgPRF by plasmacytoid dendritic cells (pDCs), whereas TLR11 and TLR12 were both required in macrophages and conventional DCs. In contrast to TLR11, TLR12-deficient mice succumb rapidly to T. gondii infection. TLR12-dependent induction of IL-12 and IFN-α in pDCs led to production of IFN-γ by NK cells. Consistent with this observation, the partial resistance of Tlr11(-/-) mice is lost upon pDC or NK cell depletion. Thus, TLR12 is critical for the innate immune response to T. gondii, and this TLR may promote host resistance by triggering pDC and NK cell function.

Genome-wide association studies dissect the G × E interaction for agronomic traits in a worldwide collection of safflowers (Carthamus tinctorius L.).

Genome-wide association studies were conducted using a globally diverse safflower (Carthamus tinctorius L.) Genebank collection for grain yield (YP), days to flowering (DF), plant height (PH), 500 seed weight (SW), seed oil content (OL), and crude protein content (PR) in four environments (sites) that differed in water availability. Phenotypic variation was observed for all traits. YP exhibited low overall genetic correlations (rGoverall) across sites, while SW and OL had high rGoverall and high pairwise genetic correlations (rGij) across all pairwise sites. In total, 92 marker-trait associations (MTAs) were identified using three methods, single locus genome-wide association studies (GWAS) using a mixed linear model (MLM), the Bayesian multi-locus method (BayesR), and meta-GWAS. MTAs with large effects across all sites were detected for OL, SW, and PR, and MTAs specific for the different water stress sites were identified for all traits. Five MTAs were associated with multiple traits; 4 of 5 MTAs were variously associated with the three traits of SW, OL, and PR. This study provided insights into the phenotypic variability and genetic architecture of important safflower agronomic traits under different environments. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-022-01295-8.

Mining the Wheat Grain Proteome.

Bread wheat is the most widely cultivated crop worldwide, used in the production of food products and a feed source for animals. Selection tools that can be applied early in the breeding cycle are needed to accelerate genetic gain for increased wheat production while maintaining or improving grain quality if demand from human population growth is to be fulfilled. Proteomics screening assays of wheat flour can assist breeders to select the best performing breeding lines and discard the worst lines. In this study, we optimised a robust LC-MS shotgun quantitative proteomics method to screen thousands of wheat genotypes. Using 6 cultivars and 4 replicates, we tested 3 resuspension ratios (50, 25, and 17 µL/mg), 2 extraction buffers (with urea or guanidine-hydrochloride), 3 sets of proteases (chymotrypsin, Glu-C, and trypsin/Lys-C), and multiple LC settings. Protein identifications by LC-MS/MS were used to select the best parameters. A total 8738 wheat proteins were identified. The best method was validated on an independent set of 96 cultivars and peptides quantities were normalised using sample weights, an internal standard, and quality controls. Data mining tools found particularly useful to explore the flour proteome are presented (UniProt Retrieve/ID mapping tool, KEGG, AgriGO, REVIGO, and Pathway Tools).