Publications

Jennifer Doucet, Christina Truong, Elizabeth Frank-Webb, Hyun Kyung Lee, Anna Daneva, Zhen Gao, Moritz K. Nowack & Daphne R. Goring  

Plant Reproduction 2019 32(3):307-322. 10.1007/s00497-019-00372-x. PMID:31069543

Abstract

We describe a function for a novel Arabidopsis gene, E6-like 1 (E6L1), that was identified as a highly expressed gene in the stigma and plays a role in early post-pollination stages. In Arabidopsis, successful pollen-stigma interactions are dependent on rapid recognition of compatible pollen by the stigmatic papillae located on the surface of the pistil and the subsequent regulation of pollen hydration and germination, and followed by the growth of pollen tubes through the stigma surface. Here we have described the function of a novel gene, E6-like 1 (E6L1), that was identified through the analysis of transcriptome datasets, as one of highest expressed genes in the stigma, and furthermore, its expression was largely restricted to the stigma and trichomes. The first E6 gene was initially identified as a highly expressed gene during cotton fiber development, and related E6-like predicted proteins are found throughout the Angiosperms. To date, no orthologous genes have been assigned a biological function. Both the Arabidopsis E6L1 and cotton E6 proteins are predicted to be secreted, and this was confirmed using an E6L1:RFP fusion construct. To further investigate E6L1’s function, one T-DNA and two independent CRISPR-generated mutants were analyzed for compatible pollen-stigma interactions, and pollen hydration, pollen adhesion, and seed set were mildly impaired for the e6l1 mutants. This work identifies E6L1 as a novel stigmatic factor that plays a role during the early post-pollination stages in Arabidopsis.

Jasmina Uzunović, Emily B Josephs, John R Stinchcombe, Stephen I Wright

Molecular Biology and Evolution 2019 36(8):1734-1745. 10.1093/molbev/msz098. PMID: 31028401

Abstract

Transposable elements (TEs) make up a significant portion of eukaryotic genomes and are important drivers of genome evolution. However, the extent to which TEs affect gene expression variation on a genome-wide scale in comparison with other types of variants is still unclear. We characterized TE insertion polymorphisms and their association with gene expression in 124 whole-genome sequences from a single population of Capsella grandiflora, and contrasted this with the effects of single nucleotide polymorphisms (SNPs). Population frequency of insertions was negatively correlated with distance to genes, as well as density of conserved noncoding elements, suggesting that the negative effects of TEs on gene regulation are important in limiting their abundance. Rare TE variants strongly influence gene expression variation, predominantly through downregulation. In contrast, rare SNPs contribute equally to up- and down-regulation, but have a weaker individual effect than TEs. An expression quantitative trait loci (eQTL) analysis shows that a greater proportion of common TEs are eQTLs as opposed to common SNPs, and a third of the genes with TE eQTLs do not have SNP eQTLs. In contrast with rare TE insertions, common insertions are more likely to increase expression, consistent with recent models of cis-regulatory evolution favoring enhancer alleles. Taken together, these results imply that TEs are a significant contributor to gene expression variation and are individually more likely than rare SNPs to cause extreme changes in gene expression.

Yi F, Gu W, Chen J, Song N, Gao X, Zhang X, Zhou Y, Ma X, Song W, Zhao H, Esteban E, Pasha A, Provart NJ, Lai J.

Plant Cell 2019 31(5):974-992. 10.1105/tpc.18.00961 PMID: 30914497

Abstract

The early maize (Zea mays) seed undergoes several developmental stages after double fertilization to become fully differentiated within a short period of time, but the genetic control of this highly dynamic and complex developmental process remains largely unknown. Here, we report a high temporal-resolution investigation of transcriptomes using 31 samples collected at an interval of 4 or 6 h within the first six days of seed development. These time-course transcriptomes were clearly separated into four distinct groups corresponding to the stages of double fertilization, coenocyte formation, cellularization, and differentiation. A total of 22,790 expressed genes including 1415 transcription factors (TFs) were detected in early stages of maize seed development. In particular, 1093 genes including 110 TFs were specifically expressed in the seed and displayed high temporal specificity by expressing only in particular period of early seed development. There were 160, 22, 112, and 569 seed-specific genes predominantly expressed in the first 16 h after pollination, coenocyte formation, cellularization, and differentiation stage, respectively. In addition, network analysis predicted 31,256 interactions among 1317 TFs and 14,540 genes. The high temporal-resolution transcriptome atlas reported here provides an important resource for future functional study to unravel the genetic control of seed development.

Dong S, Lau V, Song R, Ierullo M, Esteban E, Wu Y, Sivieng T, Nahal H, Gaudinier A, Pasha A, Oughtred R, Dolinski K, Tyers M, Brady SM, Grene R, Usadel B, Provart NJ.

Plant Physiol 2019 179(4):1893-1907. 10.1104/pp.18.01216 PMID: 30679268

Abstract

Determining the complete Arabidopsis (Arabidopsis thaliana) protein-protein interaction network is essential for understanding the functional organization of the proteome. Numerous small-scale studies and a couple of large-scale ones have elucidated a fraction of the estimated 300,000 binary protein-protein interactions in Arabidopsis. In this study, we provide evidence that a docking algorithm has the ability to identify real interactions using both experimentally determined and predicted protein structures. We ranked 0.91 million interactions generated by all possible pairwise combinations of 1,346 predicted structure models from an Arabidopsis predicted “structure-ome” and found a significant enrichment of real interactions for the top-ranking predicted interactions, as shown by cosubcellular enrichment analysis and yeast two-hybrid validation. Our success rate for computationally predicted, structure-based interactions was 63% of the success rate for published interactions naively tested using the yeast two-hybrid system and 2.7 times better than for randomly picked pairs of proteins. This study provides another perspective in interactome exploration and biological network reconstruction using protein structural information. We have made these interactions freely accessible through an improved Arabidopsis Interactions Viewer and have created community tools for accessing these and ∼2.8 million other protein-protein and protein-DNA interactions for hypothesis generation by researchers worldwide. The Arabidopsis Interactions Viewer is freely available at http://bar.utoronto.ca/interactions2/.

Zou A, Nadeau K, Wang PW, Lee JY, Guttman DS, Sharif S, Korver DR, Brumell JH, Parkinson J.

BMC Genet. 2020 Jan 17;21(1):5. doi: 10.1186/s12863-020-0807-z.

PMID: 31952471

Abstract

To satisfy an increasing demand for dietary protein, the poultry industry has employed genetic selection to increase the growth rate of broilers by over 400% in the past 50 years. Although modern broilers reach a marketable weight of ~ 2 kg in a short span of 35 days, a speed twice as fast as a broiler 50 years ago, the expedited growth has been associated with several negative detrimental consequences. Aside from heart and musculoskeletal problems, which are direct consequences of additional weight, the immune response is also thought to be altered in modern broilers.

Given that identifying the underlying genetic basis responsible for a less sensitive innate immune response would be economically beneficial for poultry breeding, we decided to compare the genomes of two unselected meat control strains that are representative of broilers from 1957 and 1978, and a current commercial broiler line. Through analysis of genetic variants, we developed a custom prioritization strategy to identify genes and pathways that have accumulated genetic changes and are biologically relevant to immune response and growth performance. Our results highlight two genes, TLR3 and PLIN3, with genetic variants that are predicted to enhance growth performance at the expense of immune function.

Placing these new genomes in the context of other chicken lines, reveal genetic changes that have specifically arisen in selective breeding programs that were implemented in the last 50 years.

Drall KM, Tun HM, Morales-Lizcano NP, Konya TB, Guttman DS, Field CJ, Mandal R, Wishart DS, Becker AB, Azad MB, Lefebvre DL, Mandhane PJ, Moraes TJ, Sears MR, Turvey SE, Subbarao P, Scott JA, Kozyrskyj AL.

Front Immunol. 2019 Dec 11;10:2866. doi: 10.3389/fimmu.2019.02866

PMID: 31921134

Abstract

Colonization with Clostridioides difficile occurs in up to half of infants under the age of 3 months, is strongly influenced by feeding modality and is largely asymptomatic. In spite of this, C. difficile‘s presence has been associated with susceptibility to chronic disease later in childhood, perhaps by promoting or benefiting from changes in infant gut microbiome development, including colonization with pathogenic bacteria and disrupted production of microbial bioactive metabolites and proteins. In this study, the microbiomes of 1554 infants from the CHILD Cohort Study were described according to C. difficile colonization status and feeding mode at 3-4 months of age. C. difficile colonization was associated with a different gut microbiome profile in exclusively breastfed (EBF) vs. exclusively formula fed (EFF) infants. EBF infants colonized with C. difficile had an increased relative abundance of Firmicutes and Proteobacteria, decreased relative abundance of Bifidobacteriaceae, greater microbiota alpha-diversity, greater detectable fecal short chain fatty acids (SCFA), and lower detectable fecal secretory Immunoglobulin A (sIgA) than those not colonized. Similar but less pronounced differences were seen among partially breastfed infants (PBF) but EFF infants did not possess these differences in the gut microbiome according to colonization status. Thus, breastfed infants colonized with C. difficile appear to possess a gut microbiome that differs from non-colonized infants and resembles that of EFF infants, but the driving force and direction of this association remains unknown. Understanding these compositional differences as drivers of C. difficile colonization may be important to ensure future childhood health.

Bastedo DP, Khan M, Martel A, Seto D, Kireeva I, Zhang J, Masud W, Millar D, Lee JY, Lee AH, Gong Y, Santos-Severino A, Guttman DS, Desveaux D.

PLoS Pathog. 2019 Jul 3;15(7):e1007900. doi: 10.1371/journal.ppat.1007900

PMID: 31269090

Abstract

The Pseudomonas syringae acetyltransferase HopZ1a is delivered into host cells by the type III secretion system to promote bacterial growth. However, in the model plant host Arabidopsis thaliana, HopZ1a activity results in an effector-triggered immune response (ETI) that limits bacterial proliferation. HopZ1a-triggered immunity requires the nucleotide-binding, leucine-rich repeat domain (NLR) protein, ZAR1, and the pseudokinase, ZED1. Here we demonstrate that HopZ1a can acetylate members of a family of ‘receptor-like cytoplasmic kinases’ (RLCK family VII; also known as PBS1-like kinases, or PBLs) and promote their interaction with ZED1 and ZAR1 to form a ZAR1-ZED1-PBL ternary complex. Interactions between ZED1 and PBL kinases are determined by the pseudokinase features of ZED1, and mutants designed to restore ZED1 kinase motifs can (1) bind to PBLs, (2) recruit ZAR1, and (3) trigger ZAR1-dependent immunity in planta, all independently of HopZ1a. A ZED1 mutant that mimics acetylation by HopZ1a also triggers immunity in planta, providing evidence that effector-induced perturbations of ZED1 also activate ZAR1. Overall, our results suggest that interactions between these two RLCK families are promoted by perturbations of structural features that distinguish active from inactive kinase domain conformations. We propose that effector-induced interactions between ZED1/ZRK pseudokinases (RLCK family XII) and PBL kinases (RLCK family VII) provide a sensitive mechanism for detecting perturbations of either kinase family to activate ZAR1-mediated ETI.

Doucet J, Lee HK, Udugama N, Xu J, Qi B, Goring DR

BMC Plant Biol. 2019 Dec 11;19(1):549. doi: 10.1186/s12870-019-2160-9.

PMID:31829135

Abstract

BACKGROUND:

In the Brassicaceae, the early stages of compatible pollen-stigma interactions are tightly controlled with early checkpoints regulating pollen adhesion, hydration and germination, and pollen tube entry into the stigmatic surface. However, the early signalling events in the stigma which trigger these compatible interactions remain unknown.

RESULTS:

A set of stigma-expressed pseudokinase genes, termed BRASSIKINs (BKNs), were identified and found to be present in only core Brassicaceae genomes. In Arabidopsis thaliana Col-0, BKN1 displayed stigma-specific expression while the BKN2 gene was expressed in other tissues as well. CRISPR deletion mutations were generated for the two tandemly linked BKNs, and very mild hydration defects were observed for wild-type Col-0 pollen when placed on the bkn1/2 mutant stigmas. In further analyses, the predominant transcript for the stigma-specific BKN1 was found to have a premature stop codon in the Col-0 ecotype, but a survey of the 1001 Arabidopsis genomes uncovered three ecotypes that encoded a full-length BKN1 protein. Furthermore, phylogenetic analyses identified intact BKN1 orthologues in the closely related outcrossing Arabidopsis species, A. lyrata and A. halleri. Finally, the BKN pseudokinases were found to be plasma-membrane localized through the dual lipid modification of myristoylation and palmitoylation, and this localization would be consistent with a role in signaling complexes.

CONCLUSION:

In this study, we have characterized the novel Brassicaceae-specific family of BKN pseudokinase genes, and examined the function of BKN1 and BKN2 in the context of pollen-stigma interactions in A. thaliana Col-0. Additionally, premature stop codons were identified in the predicted stigma specific BKN1 gene in a number of the 1001 A. thaliana ecotype genomes, and this was in contrast to the out-crossing Arabidopsis species which carried intact copies of BKN1. Thus, understanding the function of BKN1 in other Brassicaceae species will be a key direction for future studies.

Yu X, Xu G, Li B, de Souza Vespoli L, Liu H, Moeder W, Chen S, de Oliveira MVV, Ariádina de Souza S, Shao W, Rodrigues B, Ma Y, Chhajed S, Xue S, Berkowitz GA, Yoshioka K, He P, Shan L.

Curr Biol. 2019 Oct 23. pii: S0960-9822(19)31178-9. doi: 10.1016/j.cub.2019.09.018.

PMID: 31679931

Abstract:

Cell death is a vital and ubiquitous process that is tightly controlled in all organisms. However, the mechanisms underlying precise cell death control remain fragmented. As an important shared module in plant growth, development, and immunity, Arabidopsis thaliana BRASSINOSTEROID INSENSITIVE 1-associated receptor kinase 1 (BAK1) and somatic embryogenesis receptor kinase 4 (SERK4) redundantly and negatively regulate plant cell death. By deploying an RNAi-based genetic screen for bak1/serk4 cell death suppressors, we revealed that cyclic nucleotide-gated channel 20 (CNGC20) functions as a hyperpolarization-activated Ca²⁺-permeable channel specifically regulating bak1/serk4 cell death. BAK1 directly interacts with and phosphorylates CNGC20 at specific sites in the C-terminal cytosolic domain, which in turn regulates CNGC20 stability. CNGC19, the closest homolog of CNGC20 with a low abundance compared with CNGC20, makes a quantitative genetic contribution to bak1/serk4 cell death only in the absence of CNGC20, supporting the biochemical data showing homo- and heteromeric assembly of the CNGC20 and CNGC19 channel complexes. Transcripts of CNGC20 and CNGC19 are elevated in bak1/serk4 compared with wild-type plants, further substantiating a critical role of homeostasis of CNGC20 and CNGC19 in cell death control. Our studies not only uncover a unique regulation of ion channel stability by cell-surface-resident receptor kinase-mediated phosphorylation but also provide evidence for fine-tuning Ca²⁺ channel functions in maintaining cellular homeostasis by the formation of homo- and heterotetrameric complexes.

Robertson SJ, Lemire P, Maughan H, Goethel A, Turpin W, Bedrani L, Guttman DS, Croitoru K, Girardin SE, Philpott DJ

Cell Rep. 2019 May 7;27(6):1910-1919.e2. doi: 10.1016/j.celrep.2019.04.023

PMID: 31067473

Abstract:

The intestinal microbiota is a fundamental factor that broadly influences physiology. Thus, studies using transgenic animals should be designed to limit the confounding effects of microbiota variation between strains. Here, we report the impact on intestinal microbiota of co-housed versus F2-generation littermates, two commonly used techniques to standardize microbiota in animal models. Our results establish that while fecal microbiota is partially normalized by extended co-housing, mucosal communities associated with the proximal colon and terminal ileum remain stable and distinct. In contrast, strain inter-crossing to generate F2 littermates allows robust microbiota standardization in fecal, colon, and ileum sampling locations. Using reciprocal inter-crosses of P1 parents, we identify dissymmetry in F2 community structures caused by maternal transmission, in particular of the Verrucomicrobiaceae. Thus, F2 littermate animals from a unidirectional P1 cross should be used as a standard method to minimize the influence of the microbiota in genotype-phenotype studies.