ERG Functionally Overlaps with Other Ets Proteins in Promoting TH9 Cell Expression of Il9 during Allergic Lung Inflammation.

CD4+ TH cells develop into subsets that are specialized in the secretion of particular cytokines to mediate restricted types of inflammation and immune responses. Among the subsets that promote development of allergic inflammatory responses, IL-9-producing TH9 cells are regulated by a number of transcription factors. We have previously shown that the E26 transformation-specific (Ets) family members PU.1 and Ets translocation variant 5 (ETV5) function in parallel to regulate IL-9. In this study we identified a third member of the Ets family of transcription factors, Ets-related gene (ERG), that mediates IL-9 production in TH9 cells in the absence of PU.1 and ETV5. Chromatin immunoprecipitation assays revealed that ERG interaction at the Il9 promoter region is restricted to the TH9 lineage and is sustained during murine TH9 polarization. Knockdown or knockout of ERG during murine or human TH9 polarization in vitro led to a decrease in IL-9 production in TH9 cells. Deletion of ERG in vivo had modest effects on IL-9 production in vitro or in vivo. However, in the absence of PU.1 and ETV5, ERG was required for residual IL-9 production in vitro and for IL-9 production by lung-derived CD4 T cells in a mouse model of chronic allergic airway disease. Thus, ERG contributes to IL-9 regulation in TH9 cells.

An adaptive teosinte mexicana introgression modulates phosphatidylcholine levels and is associated with maize flowering time.

Native Americans domesticated maize (Zea mays ssp. mays) from lowland teosinte parviglumis (Zea mays ssp. parviglumis) in the warm Mexican southwest and brought it to the highlands of Mexico and South America where it was exposed to lower temperatures that imposed strong selection on flowering time. Phospholipids are important metabolites in plant responses to low-temperature and phosphorus availability and have been suggested to influence flowering time. Here, we combined linkage mapping with genome scans to identify High PhosphatidylCholine 1 (HPC1), a gene that encodes a phospholipase A1 enzyme, as a major driver of phospholipid variation in highland maize. Common garden experiments demonstrated strong genotype-by-environment interactions associated with variation at HPC1, with the highland HPC1 allele leading to higher fitness in highlands, possibly by hastening flowering. The highland maize HPC1 variant resulted in impaired function of the encoded protein due to a polymorphism in a highly conserved sequence. A meta-analysis across HPC1 orthologs indicated a strong association between the identity of the amino acid at this position and optimal growth in prokaryotes. Mutagenesis of HPC1 via genome editing validated its role in regulating phospholipid metabolism. Finally, we showed that the highland HPC1 allele entered cultivated maize by introgression from the wild highland teosinte Zea mays ssp. mexicana and has been maintained in maize breeding lines from the Northern United States, Canada, and Europe. Thus, HPC1 introgressed from teosinte mexicana underlies a large metabolic QTL that modulates phosphatidylcholine levels and has an adaptive effect at least in part via induction of early flowering time.

Unraveling Prevalence and Effects of Deleterious Mutations in Maize Elite Lines across Decades of Modern Breeding.

Future breeding is likely to involve the detection and removal of deleterious alleles, which are mutations that negatively affect crop fitness. However, little is known about the prevalence of such mutations and their effects on phenotypic traits in the context of modern crop breeding. To address this, we examined the number and frequency of deleterious mutations in 350 elite maize inbred lines developed over the past few decades in China and the United States. Our findings reveal an accumulation of weakly deleterious mutations and a decrease in strongly deleterious mutations, indicating the dominant effects of genetic drift and purifying selection for the two types of mutations, respectively. We also discovered that slightly deleterious mutations, when at lower frequencies, were more likely to be heterozygous in the developed hybrids. This is consistent with complementation as a potential explanation for heterosis. Subsequently, we found that deleterious mutations accounted for more of the variation in phenotypic traits than nondeleterious mutations with matched minor allele frequencies, especially for traits related to leaf angle and flowering time. Moreover, we detected fewer deleterious mutations in the promoter and gene body regions of differentially expressed genes across breeding eras than in nondifferentially expressed genes. Overall, our results provide a comprehensive assessment of the prevalence and impact of deleterious mutations in modern maize breeding and establish a useful baseline for future maize improvement efforts.

Conserved noncoding sequences provide insights into regulatory sequence and loss of gene expression in maize.

Thousands of species will be sequenced in the next few years; however, understanding how their genomes work, without an unlimited budget, requires both molecular and novel evolutionary approaches. We developed a sensitive sequence alignment pipeline to identify conserved noncoding sequences (CNSs) in the Andropogoneae tribe (multiple crop species descended from a common ancestor ∼18 million years ago). The Andropogoneae share similar physiology while being tremendously genomically diverse, harboring a broad range of ploidy levels, structural variation, and transposons. These contribute to the potential of Andropogoneae as a powerful system for studying CNSs and are factors we leverage to understand the function of maize CNSs. We found that 86% of CNSs were comprised of annotated features, including introns, UTRs, putative cis-regulatory elements, chromatin loop anchors, noncoding RNA (ncRNA) genes, and several transposable element superfamilies. CNSs were enriched in active regions of DNA replication in the early S phase of the mitotic cell cycle and showed different DNA methylation ratios compared to the genome-wide background. More than half of putative cis-regulatory sequences (identified via other methods) overlapped with CNSs detected in this study. Variants in CNSs were associated with gene expression levels, and CNS absence contributed to loss of gene expression. Furthermore, the evolution of CNSs was associated with the functional diversification of duplicated genes in the context of maize subgenomes. Our results provide a quantitative understanding of the molecular processes governing the evolution of CNSs in maize.

Fine mapping of major QTL qshgd1 for spontaneous haploid genome doubling in maize (Zea mays L.).

KEY MESSAGE: A large-effect QTL was fine mapped, which revealed 79 gene models, with 10 promising candidate genes, along with a novel inversion. In commercial maize breeding, doubled haploid (DH) technology is arguably the most efficient resource for rapidly developing novel, completely homozygous lines. However, the DH strategy, using in vivo haploid induction, currently requires the use of mutagenic agents which can be not only hazardous, but laborious. This study focuses on an alternative approach to develop DH lines-spontaneous haploid genome duplication (SHGD) via naturally restored haploid male fertility (HMF). Inbred lines A427 and Wf9, the former with high HMF and the latter with low HMF, were selected to fine-map a large-effect QTL associated with SHGD-qshgd1. SHGD alleles were derived from A427, with novel haploid recombinant groups having varying levels of the A427 chromosomal region recovered. The chromosomal region of interest is composed of 45 megabases (Mb) of genetic information on chromosome 5. Significant differences between haploid recombinant groups for HMF were identified, signaling the possibility of mapping the QTL more closely. Due to suppression of recombination from the proximity of the centromere, and a newly discovered inversion region, the associated QTL was only confined to a 25 Mb region, within which only a single recombinant was observed among ca. 9,000 BC1 individuals. Nevertheless, 79 gene models were identified within this 25 Mb region. Additionally, 10 promising candidate genes, based on RNA-seq data, are described for future evaluation, while the narrowed down genome region is accessible for straightforward introgression into elite germplasm by BC methods.

Allele-specific Expression Reveals Multiple Paths to Highland Adaptation in Maize.

Maize is a staple food of smallholder farmers living in highland regions up to 4,000 m above sea level worldwide. Mexican and South American highlands are two major highland maize growing regions, and population genetic data suggest the maize's adaptation to these regions occurred largely independently, providing a case study for convergent evolution. To better understand the mechanistic basis of highland adaptation, we crossed maize landraces from 108 highland and lowland sites of Mexico and South America with the inbred line B73 to produce F1 hybrids and grew them in both highland and lowland sites in Mexico. We identified thousands of genes with divergent expression between highland and lowland populations. Hundreds of these genes show patterns of convergent evolution between Mexico and South America. To dissect the genetic architecture of the divergent gene expression, we developed a novel allele-specific expression analysis pipeline to detect genes with divergent functional cis-regulatory variation between highland and lowland populations. We identified hundreds of genes with divergent cis-regulation between highland and lowland landrace alleles, with 20 in common between regions, further suggesting convergence in the genes underlying highland adaptation. Further analyses suggest multiple mechanisms contribute to this convergence in gene regulation. Although the vast majority of evolutionary changes associated with highland adaptation were region specific, our findings highlight an important role for convergence at the gene expression and gene regulation levels as well.

Molecular Parallelism Underlies Convergent Highland Adaptation of Maize Landraces.

Convergent phenotypic evolution provides some of the strongest evidence for adaptation. However, the extent to which recurrent phenotypic adaptation has arisen via parallelism at the molecular level remains unresolved, as does the evolutionary origin of alleles underlying such adaptation. Here, we investigate genetic mechanisms of convergent highland adaptation in maize landrace populations and evaluate the genetic sources of recurrently selected alleles. Population branch excess statistics reveal substantial evidence of parallel adaptation at the level of individual single-nucleotide polymorphism (SNPs), genes, and pathways in four independent highland maize populations. The majority of convergently selected SNPs originated via migration from a single population, most likely in the Mesoamerican highlands, while standing variation introduced by ancient gene flow was also a contributor. Polygenic adaptation analyses of quantitative traits reveal that alleles affecting flowering time are significantly associated with elevation, indicating the flowering time pathway was targeted by highland adaptation. In addition, repeatedly selected genes were significantly enriched in the flowering time pathway, indicating their significance in adapting to highland conditions. Overall, our study system represents a promising model to study convergent evolution in plants with potential applications to crop adaptation across environmental gradients.

Major locus for spontaneous haploid genome doubling detected by a case-control GWAS in exotic maize germplasm.

KEY MESSAGE: A major locus for spontaneous haploid genome doubling was detected by a case-control GWAS in an exotic maize germplasm. The combination of double haploid breeding method with this locus leads to segregation distortion on genomic regions of chromosome five. Temperate maize (Zea mays L.) breeding programs often rely on limited genetic diversity, which can be expanded by incorporating exotic germplasm. The aims of this study were to perform characterization of inbred lines derived from the tropical BS39 population using different breeding methods, to identify genomic regions showing segregation distortion in lines derived by the DH process using spontaneous haploid genome doubling (SHGD), and use case-control association mapping to identify loci controlling SHGD. Four different sets were used: BS39_DH and BS39_SSD were derived from the BS39 population by DH and single-seed descendent (SSD) methods, and BS39 × A427_DH and BS39 × A427_SSD from the cross between BS39 and A427. A total of 663 inbred lines were genotyped. The analyses of gene diversity and genetic differentiation for the DH sets provided evidence of the presence of a SHGD locus near the centromere of chromosome 5. The case-control GWAS for the DH set also pinpointed this locus. Haplotype sharing analysis showed almost 100% exclusive contribution of the A427 genome in the same region on chromosome 5 of BS39 × A427_DH, presumably due to an allele in this region affecting SHGD. This locus enables DH line production in exotic populations without colchicine or other artificial haploid genome doubling.

SequelTools: a suite of tools for working with PacBio Sequel raw sequence data.

BACKGROUND: PacBio sequencing is an incredibly valuable third-generation DNA sequencing method due to very long read lengths, ability to detect methylated bases, and its real-time sequencing methodology. Yet, hitherto no tool was available for analyzing the quality of, subsampling, and filtering PacBio data. RESULTS: Here we present SequelTools, a command-line program containing three tools: Quality Control, Read Subsampling, and Read Filtering. The Quality Control tool quickly processes PacBio Sequel raw sequence data from multiple SMRTcells producing multiple statistics and publication-quality plots describing the quality of the data including N50, read length and count statistics, PSR, and ZOR. The Read Subsampling tool allows the user to subsample reads by one or more of the following criteria: longest subreads per CLR or random CLR selection. The Read Filtering tool provides options for normalizing data by filtering out certain low-quality scraps reads and/or by minimum CLR length. SequelTools is implemented in bash, R, and Python using only standard libraries and packages and is platform independent. CONCLUSIONS: SequelTools is a program that provides the only free, fast, and easy-to-use quality control tool, and the only program providing this kind of read subsampling and read filtering for PacBio Sequel raw sequence data, and is available at https://github.com/ISUgenomics/SequelTools .

GenomeQC: a quality assessment tool for genome assemblies and gene structure annotations.

BACKGROUND: Genome assemblies are foundational for understanding the biology of a species. They provide a physical framework for mapping additional sequences, thereby enabling characterization of, for example, genomic diversity and differences in gene expression across individuals and tissue types. Quality metrics for genome assemblies gauge both the completeness and contiguity of an assembly and help provide confidence in downstream biological insights. To compare quality across multiple assemblies, a set of common metrics are typically calculated and then compared to one or more gold standard reference genomes. While several tools exist for calculating individual metrics, applications providing comprehensive evaluations of multiple assembly features are, perhaps surprisingly, lacking. Here, we describe a new toolkit that integrates multiple metrics to characterize both assembly and gene annotation quality in a way that enables comparison across multiple assemblies and assembly types. RESULTS: Our application, named GenomeQC, is an easy-to-use and interactive web framework that integrates various quantitative measures to characterize genome assemblies and annotations. GenomeQC provides researchers with a comprehensive summary of these statistics and allows for benchmarking against gold standard reference assemblies. CONCLUSIONS: The GenomeQC web application is implemented in R/Shiny version 1.5.9 and Python 3.6 and is freely available at https://genomeqc.maizegdb.org/ under the GPL license. All source code and a containerized version of the GenomeQC pipeline is available in the GitHub repository https://github.com/HuffordLab/GenomeQC.

Haplotype structure in commercial maize breeding programs in relation to key founder lines.

KEY MESSAGE: High-density haplotype analysis revealed significant haplotype sharing between ex-PVPs registered from 1976 to 1992 and key maize founders, and uncovered similarities and differences in haplotype sharing patterns by company and heterotic group. Proprietary inbreds developed by the private seed industry have been the major source for driving genetic gain in successful North American maize hybrids for decades. Much of the history of industry germplasm can be traced back to key founder lines, some of which were pivotal in the development of prominent heterotic groups. Previous studies have summarized pedigree-based relationships, genetic diversity and population structure among commercial inbreds with expired Plant Variety Protection (ex-PVP). However, less is known about the extent of haplotype sharing between historical founders and ex-PVPs. A better understanding of the relationships between founders and ex-PVPs provides insight into the haplotype and heterotic group structure among industry germplasm. We performed high-density haplotype analysis with 11.3 million SNPs on 212 maize inbreds, which included 157 ex-PVPs registered 1976-1992 and 55 public lines relevant to PVPs. Among these lines were 12 key founders identified in literature review: 207, A632, B14, B37, B73, LH123HT, LH82, Mo17, Oh43, OH7, PHG39 and Wf9. Our results revealed that, on average, 81.6% of an ex-PVP's genome is shared with at least 1 of these 12 founder lines and more than half when limited to B73, Mo17 and 207. Quantifiable similarities and contrasts among heterotic groups and major US seed industry companies were also observed. The results from this study provide high-resolution haplotype data on ex-PVP germplasm, confirm founder relationship trends observed in previous studies, uncover region-specific haplotype structure differences and demonstrate how haplotype sharing analysis can be used as a tool to explore germplasm diversity.

Evidence for a Unique DNA-Dependent RNA Polymerase in Cereal Crops.

Gene duplication is an important driver for the evolution of new genes and protein functions. Duplication of DNA-dependent RNA polymerase (Pol) II subunits within plants led to the emergence of RNA Pol IV and V complexes, each of which possess unique functions necessary for RNA-directed DNA Methylation. Comprehensive identification of Pol V subunit orthologs across the monocot radiation revealed a duplication of the largest two subunits within the grasses (Poaceae), including critical cereal crops. These paralogous Pol subunits display sequence conservation within catalytic domains, but their carboxy terminal domains differ in length and character of the Ago-binding platform, suggesting unique functional interactions. Phylogenetic analysis of the catalytic region indicates positive selection on one paralog following duplication, consistent with retention via neofunctionalization. Positive selection on residue pairs that are predicted to interact between subunits suggests that paralogous subunits have evolved specific assembly partners. Additional Pol subunits as well as Pol-interacting proteins also possess grass-specific paralogs, supporting the hypothesis that a novel Pol complex with distinct function has evolved in the grass family, Poaceae.

The contribution of joint and skin improvements to the health-related quality of life of patients with psoriatic arthritis: a post hoc analysis of two randomised controlled studies.

OBJECTIVE: Determine the contribution of joint and skin improvements to health-related quality of life (HRQoL) in patients with psoriatic arthritis (PsA). METHODS: SPIRIT-P1 and SPIRIT-P2 are phase 3 trials investigating ixekizumab, an interleukin-17A antagonist, in the treatment of patients with active PsA. Patients were randomised to ixekizumab or placebo. Outcomes included the Disease Activity Index for Psoriatic Arthritis (DAPSA), the Psoriasis Area and Severity Index (PASI), the European Quality of Life-Five Dimensions (EQ-5D) Visual Analogue Score (VAS), the 36-Item Short-Form Health Survey (SF-36) and the Work Productivity and Activity Impairment (WPAI) Questionnaire. The contribution of joint and skin improvements to HRQoL was modelled using a smoothing spline method and depicted with response surface graphics. RESULTS: In this integrated analysis, 402 patients with PsA had baseline psoriasis of ≥3% of body surface area. We applied response surface modelling to this patient data set to investigate the relationship between DAPSA, PASI and HRQoL improvements at week 24. The greatest improvement in EQ-5D VAS was associated with the largest per cent improvements in both DAPSA and PASI together, rather than DAPSA or PASI alone. Similar observations were made in domains of SF-36 and WPAI. CONCLUSION: Optimal improvements in patients' HRQoL were dependent on successful treatment of both joint and skin symptoms.

The extent of adaptive wild introgression in crops.

The study of crop evolution has focused primarily on the process of initial domestication. Post-domestication adaptation during the expansion of crops from their centers of origin has received considerably less attention. Recent research has revealed that, in at least some instances, crops have received introgression from their wild relatives that has facilitated adaptation to novel conditions encountered during expansion. Such adaptive introgression could have an important impact on the basic study of domestication, affecting estimates of several evolutionary processes of interest (e.g. the strength of the domestication bottleneck, the timing of domestication, the targets of selection during domestication). Identification of haplotypes introgressed from the wild may also help in the identification of alleles that are beneficial under particular environmental conditions. Here we review mounting evidence for substantial adaptive wild introgression in several crops and consider the implications of such gene flow to our understanding of crop histories.

Technological advances in maize breeding: past, present and future.

Maize has for many decades been both one of the most important crops worldwide and one of the primary genetic model organisms. More recently, maize breeding has been impacted by rapid technological advances in sequencing and genotyping technology, transformation including genome editing, doubled haploid technology, parallelled by progress in data sciences and the development of novel breeding approaches utilizing genomic information. Herein, we report on past, current and future developments relevant for maize breeding with regard to (1) genome analysis, (2) germplasm diversity characterization and utilization, (3) manipulation of genetic diversity by transformation and genome editing, (4) inbred line development and hybrid seed production, (5) understanding and prediction of hybrid performance, (6) breeding methodology and (7) synthesis of opportunities and challenges for future maize breeding.

Etv5 Regulates IL-10 Production in Th Cells.

IL-10 is an immunoregulatory cytokine that has broad effects across the immune system. In Th cell subsets, Th2 cells produce considerable amounts of IL-10. The transcription factors that regulate IL-10 production in Th2 cells are still incompletely described. In this study, we demonstrate that the ETS family transcription factor ETS variant (Etv)5 regulates IL-10 production in Th2 cells. T cell-specific Etv5-deficient and littermate control mice demonstrated that IL-10 production and gene expression were significantly decreased in the absence of Etv5. In an Aspergillus fumigatus extract-induced inflammation model, IL-10-producing CD4+ T cells in bronchoalveolar lavage and lung were significantly decreased in mice that lacked Etv5 in T cells, compared with control mice. We showed that Etv5 directly binds to the Il10 locus conserved noncoding sequence 3 site and that it activates gene expression in a luciferase reporter assay and following retroviral transduction. Etv5 deficiency did not affect the expression of other transcription factors known to be important for expression of IL-10, including Jun family members, GATA3, E4BP4, and IFN regulatory factor 4. However, in the absence of Etv5, binding of these transcription factors to the Il10 locus was dramatically reduced. Ectopic Etv5 expression in Th2 cells that lack Etv5 restored IL-10 production and the binding of IL-10-inducing transcription factors including E4BP4, IFN regulatory factor 4, and GATA3. Taken together, we conclude that Etv5 plays a crucial role in regulating IL-10 production in Th2 cells by facilitating the binding of IL-10-inducing transcription factors at the Il10 locus.

STAT3 Impairs STAT5 Activation in the Development of IL-9-Secreting T Cells.

Th cell subsets develop in response to multiple activating signals, including the cytokine environment. IL-9-secreting T cells develop in response to the combination of IL-4 and TGF-ß, although they clearly require other cytokine signals, leading to the activation of transcription factors including STAT5. In Th17 cells, there is a molecular antagonism of STAT5 with STAT3 signaling, although whether this paradigm exists in other Th subsets is not clear. In this paper, we demonstrate that STAT3 attenuates the ability of STAT5 to promote the development of IL-9-secreting T cells. We demonstrate that production of IL-9 is increased in the absence of STAT3 and cytokines that result in a sustained activation of STAT3, including IL-6, have the greatest potency in repressing IL-9 production in a STAT3-dependent manner. Increased IL-9 production in the absence of STAT3 correlates with increased endogenous IL-2 production and STAT5 activation, and blocking IL-2 responses eliminates the difference in IL-9 production between wild-type and STAT3-deficient T cells. Moreover, transduction of developing Th9 cells with a constitutively active STAT5 eliminates the ability of IL-6 to reduce IL-9 production. Thus, STAT3 functions as a negative regulator of IL-9 production through attenuation of STAT5 activation and function.

The ETS Family Transcription Factors Etv5 and PU.1 Function in Parallel To Promote Th9 Cell Development.

The IL-9-secreting Th9 subset of CD4 Th cells develop in response to an environment containing IL-4 and TGF-ß, promoting allergic disease, autoimmunity, and resistance to pathogens. We previously identified a requirement for the ETS family transcription factor PU.1 in Th9 development. In this report, we demonstrate that the ETS transcription factor ETS variant 5 (ETV5) promotes IL-9 production in Th9 cells by binding and recruiting histone acetyltransferases to the Il9 locus at sites distinct from PU.1. In cells that are deficient in both PU.1 and ETV5 there is lower IL-9 production than in cells lacking either factor alone. In vivo loss of PU.1 and ETV5 in T cells results in distinct effects on allergic inflammation in the lung, suggesting that these factors function in parallel. Together, these data define a role for ETV5 in Th9 development and extend the paradigm of related transcription factors having complementary functions during differentiation.

PU.1 Expression in T Follicular Helper Cells Limits CD40L-Dependent Germinal Center B Cell Development.

PU.1 is an ETS family transcription factor that is important for the development of multiple hematopoietic cell lineages. Previous work demonstrated a critical role for PU.1 in promoting Th9 development and in limiting Th2 cytokine production. Whether PU.1 has functions in other Th lineages is not clear. In this study, we examined the effects of ectopic expression of PU.1 in CD4(+) T cells and observed decreased expression of genes involved with the function of T follicular helper (Tfh) cells, including Il21 and Tnfsf5 (encoding CD40L). T cells from conditional mutant mice that lack expression of PU.1 in T cells (Sfpi1(lck-/-)) demonstrated increased production of CD40L and IL-21 in vitro. Following adjuvant-dependent or adjuvant-independent immunization, we observed that Sfpi1(lck-/-) mice had increased numbers of Tfh cells, increased germinal center B cells (GCB cells), and increased Ab production in vivo. This correlated with increased expression of IL-21 and CD40L in Tfh cells from Sfpi1(lck-/-) mice compared with control mice. Finally, although blockade of IL-21 did not affect GCB cells in Sfpi1(lck-/-) mice, anti-CD40L treatment of immunized Sfpi1(lck-/-) mice decreased GCB cell numbers and Ag-specific Ig concentrations. Together, these data indicate an inhibitory role for PU.1 in the function of Tfh cells, germinal centers, and Tfh-dependent humoral immunity.