Heme acts through the Bach1b/Nrf2a-MafK pathway to regulate exocrine peptidase precursor genes in porphyric zebrafish

Using a zebrafish model of hepatoerythropoietic porphyria (HEP), we identify a previously unknown mechanism underlying heme-mediated regulation of exocrine zymogens. Zebrafish bach1b, nrf2a and mafK are all expressed in the zebrafish exocrine pancreas. Overexpression of bach1b or knockdown of nrf2a result in the downregulation of the expression of the exocrine zymogens, whereas overexpression of nrf2a or knockdown of bach1b cause their upregulation. In vitro luciferase assays demonstrate that heme activates the zymogens in a dosage-dependent manner and that the zymogen promoter activities require the integral Maf recognition element (MARE) motif. The Bach1b-MafK heterodimer represses the zymogen promoters, whereas the Nrf2a-MafK heterodimer activates them. Furthermore, chromatin immunoprecipitation (ChIP) assays show that MafK binds to the MARE sites in the 5′ regulatory regions of the zymogens. Taken together, these data indicate that heme stimulates the exchange of Bach1b for Nrf2a at MafK-occupied MARE sites and that, particularly in heme-deficient porphyria, the repressive Bach1b-MafK heterodimer dominates, which can be exchanged for the activating Nrf2a-MafK heterodimer upon treatment with hemin. These results provide novel insights into the regulation of exocrine function, as well as the pathogenesis of porphyria, and should be useful for designing new therapies for both types of disease.KEY WORDS: Bach1b, Nrf2a, MafK, Heme, Porphyria, Zymogens, Zebrafish     

Comparative genomic analysis of teleost fish <Emphasis Type="Italic">bmal</Emphasis> genes

Bmal1 (Brain and muscle ARNT like 1) gene is a key circadian clock gene. Tetrapods also have the second Bmal gene, Bmal2. Fruit fly has only one bmal1/cycle gene. Interrogation of the five teleost fish genome sequences coupled with phylogenetic and splice site analyses found that zebrafish have two bmal1 genes, bmal1a and bmal1b, and bmal2a; Japanese pufferfish (fugu), green spotted pufferfish (tetraodon) and Japanese medaka fish each have two bmal2 genes, bmal2a and bmal2b, and bmal1a; and three-spine stickleback have bmal1a and bmal2b. Syntenic analysis further indicated that zebrafish bmal1a/bmal1b, and fugu, tetraodon and medaka bmal2a/bmal2b are ancient duplicates. Although the dN/dS ratios of these four fish bmal duplicates are all <1, implicating they have been under purifying selection, the Tajima relative rate test showed that fugu, tetraodon and medaka bmal2a/bmal2b have asymmetric evolutionary rates, suggesting that one of these duplicates have been subject to positive selection or relaxed functional constraint. These results support the notion that teleost fish bmal genes were derived from the fish-specific genome duplication (FSGD), divergent resolution following the duplication led to retaining different ancient bmal duplicates in different fishes, which could have shaped the evolution of the complex teleost fish timekeeping mechanisms. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

SUMO E3 Ligases GmSIZ1a and GmSIZ1b regulate vegetative growth in soybean

SIZ1 is a small ubiquitin‐related modifier (SUMO) E3 ligase that mediates post‐translational SUMO modification of target proteins and thereby regulates developmental processes and hormonal and environmental stress responses in Arabidopsis. However, the role of SUMO E3 ligases in crop plants is largely unknown. Here, we identified and characterized two Glycine max (soybean) SUMO E3 ligases, GmSIZ1a and GmSIZ1b. Expression of GmSIZ1a and GmSIZ1b was induced in response to salicylic acid (SA), heat, and dehydration treatment, but not in response to cold, abscisic acid (ABA), and NaCl treatment. Although GmSIZ1a was expressed at higher levels than GmSIZ1b, both genes encoded proteins with SUMO E3 ligase activity in vivo. Heterologous expression of GmSIZ1a or GmSIZ1b rescued the mutant phenotype of Arabidopsis siz1‐2, including dwarfism, constitutively activated expression of pathogen‐related genes, and ABA‐sensitive seed germination. Simultaneous downregulation of GmSIZ1a and GmSIZ1b (GmSIZ1a/b) using RNA interference (RNAi)‐mediated gene silencing decreased heat shock‐induced SUMO conjugation in soybean. Moreover, GmSIZ1RNAi plants exhibited reduced plant height and leaf size. However, unlike Arabidopsis siz1‐2 mutant plants, flowering time and SA levels were not significantly altered in GmSIZ1RNAi plants. Taken together, our results indicate that GmSIZ1a and GmSIZ1b mediate SUMO modification and positively regulate vegetative growth in soybean.     

Plasticity and innovation of regulatory mechanisms underlying seed oil content mediated by duplicated genes in the palaeopolyploid soybean

Many plants have undergone whole genome duplication (WGD). However, how regulatory networks underlying a particular trait are reshaped in polyploids has not been experimentally investigated. Here we show that the regulatory pathways modulating seed oil content, which involve WRINKLED1 (WRI1), LEAFY COTYLEDON1 (LEC1), and LEC2 in Arabidopsis, have been modified in the palaeopolyploid soybean. Such modifications include functional reduction of GmWRI1b of the GmWRI1a/GmWRI1b homoeologous pair relevant to WRI1, complementary non‐allelic dosage effects of the GmLEC1a/GmLEC1b homoeologous pair relevant to LEC1, pseudogenization of the singleton GmLEC2 relevant to LEC2, and the rise of the LEC2‐like function of GmABI3b, contrasting to its homoeolog GmABI3a, which maintains the ABSCISIC ACID INSENSITIVE 3 (ABI3)‐like function in modulating seed maturation and dormancy. The function of GmABI3b in modulating seed oil biosynthesis was fulfilled by direct binding to a RY (CATGCA) cis‐regulatory element in the GmWRI1a promoter, which was absent in the GmWRI1b promoter, resulting in reduction of the GmWRI1b expression. Nevertheless, the three regulators each exhibited similar intensities of purifying selection to their respective duplicates since these pairs were formed by a WGD event that is proposed to have occurred approximately 13 million years ago (mya), suggesting that the differentiation in spatiotemporal expression between the duplicated genes is more likely to be the outcome of neutral variation in regulatory sequences. This study thus exemplifies the plasticity, dynamics, and novelty of regulatory networks mediated by WGD.     

Mutation of Oryza sativa CORONATINE INSENSITIVE 1b (OsCOI1b) delays leaf senescence

Jasmonic acid (JA) functions in plant development, including senescence and immunity. Arabidopsis thaliana CORONATINE INSENSITIVE 1 encodes a JA receptor and functions in the JA‐responsive signaling pathway. The Arabidopsis genome harbors a single COI gene, but the rice (Oryza sativa) genome harbors three COI homologs, OsCOI1a, OsCOI1b, and OsCOI2. Thus, it remains unclear whether each OsCOI has distinct, additive, synergistic, or redundant functions in development. Here, we use the oscoi1b‐1 knockout mutants to show that OsCOI1b mainly affects leaf senescence under senescence‐promoting conditions. oscoi1b‐1 mutants stayed green during dark‐induced and natural senescence, with substantial retention of chlorophylls and photosynthetic capacity. Furthermore, several senescence‐associated genes were downregulated in oscoi1b‐1 mutants, including homologs of Arabidopsis thaliana ETHYLENE INSENSITIVE 3 and ORESARA 1, important regulators of leaf senescence. These results suggest that crosstalk between JA signaling and ethylene signaling affects leaf senescence. The Arabidopsis coi1‐1 plants containing 35S:OsCOI1a or 35S:OsCOI1b rescued the delayed leaf senescence during dark incubation, suggesting that both OsCOI1a and OsCOI1b are required for promoting leaf senescence in rice. oscoi1b‐1 mutants showed significant decreases in spikelet fertility and grain weight, leading to severe reduction of grain yield, indicating that OsCOI1‐mediated JA signaling affects spikelet fertility and grain filling.     

Structure, organization and expression of the metallothionein gene family inArabidopsis

Metallothioneins (MTs) are cysteine-rich proteins required for heavy metal tolerance in animals and fungi. Recent results indicate that plants also possess functional metallothionein genes. Here we report the cloning and characterization of five metallothionein genes fromArabidopsis thaliana. The position of the single intron in each gene is conserved. The proteins encoded by these genes can be divided into two groups (MT1 and MT2) based on the presence or absence of a central domain separating two cysteine-rich domains. Four of the MT genes (MT1a,MT1c,MT2a andMT2b) are transcribed inArabidopsis. Several lines of evidence suggest that the fifth gene,MT1b, is inactive. There is differential regulation of the MT gene family. MT1 mRNA is expressed highly in roots, moderately in leaves and is barely detected in inflorescences and siliques. MT2a and MT2b mRNAs are more abundant in leaves, inflorescences and in roots from mature plants, but are also detected in roots of young plants, and in siliques. MT2a mRNA is strongly induced in seedlings by CUSO₄, whereas MT2b mRNA is relatively abundant in this tissue and levels increase only slightly upon exposure to copper.MT1a andMT1c are located within 2 kb of each other and have been mapped to chromosome 1.MT1b andMT2b map to separate loci on chromosome V, andMT2a is located on chromosome III. The locations of these MT genes are different from that ofCAD1, a gene involved in cadmium tolerance inArabidopsis.     

Phylogenetic and evolutionary relationships and developmental expression patterns of the zebrafish <Emphasis Type="Italic">twist</Emphasis> gene family

Four members of the twist gene family (twist1a, 1b, 2, and 3) are found in the zebrafish, and they are thought to have arisen through three rounds of gene duplication, two of which occurred prior to the tetrapod-fish split. Phylogenetic analysis groups most of the vertebrate Twist1 peptides into clade I, except for the Twist1b proteins of the acanthopterygian fish (medaka, pufferfish, stickleback), which clustered within clade III. Paralogies and orthologies among the zebrafish, medaka, and human twist genes were determined using comparative synteny analysis of the chromosomal regions flanking these genes. Comparative nucleotide substitution analyses also revealed a faster rate of nucleotide mutation/substitution in the acanthopterygian twist1b compared to the zebrafish twist1b, thus accounting for their anomalous phylogenetic clustering. We also observed minimal expression overlap among the four twist genes, suggesting that despite their significant peptide similarity, their regulatory controls have diverged considerably, with minimal functional redundancy between them. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Length–weight relationships of 40 freshwater fish species from two decades of monitoring in Flanders (Belgium)

Nearly 263 000 individual length–weight (L/W) data, collected during 2839 fish stock assessments between 1992 and 2009, were used to calculate L/W relationships (W = aL^b) of 40 freshwater fish species from Flanders (Belgium). Roach Rutilus rutilus, perch Perca fluviatilis, three‐spine stickleback Gasterosteus aculeatus, gudgeon Gobio gobio and European eel Anguilla anguilla were the most abundant species in our surveys. Mean a and b values were in accordance with those available in FishBase. Mean b across all examined species is 3.094 (95% CL = 3.041 to 3.148) indicating a tendency towards slightly positive allometric growth in most fishes (thirty‐three species with b > 3). On the basis of a form factor equation, we characterized 2 species as ‘eel‐like’, 7 as ‘elongated’, 22 as ‘fusiform’ and 9 species as ‘short and deep’. This paper adds a substantial dataset to the yet scarcely available a and b values of European populations of freshwater fish species.     

Accumulation of the NON‐YELLOW COLORING 1 protein of the chlorophyll cycle requires chlorophyll b in Arabidopsis thaliana

Chlorophyll a and chlorophyll b are interconverted in the chlorophyll cycle. The initial step in the conversion of chlorophyll b to chlorophyll a is catalyzed by the chlorophyll b reductases NON‐YELLOW COLORING 1 (NYC1) and NYC1‐like (NOL), which convert chlorophyll b to 7‐hydroxymethyl chlorophyll a. This step is also the first stage in the degradation of the light‐harvesting chlorophyll a/b protein complex (LHC). In this study, we examined the effect of chlorophyll b on the level of NYC1. NYC1 mRNA and NYC1 protein were in low abundance in green leaves, but their levels increased in response to dark‐induced senescence. When the level of chlorophyll b was enhanced by the introduction of a truncated chlorophyllide a oxygenase gene and the leaves were incubated in the dark, the amount of NYC1 was greatly increased compared with that of the wild type; however, the amount of NYC1 mRNA was the same as in the wild type. In contrast, NYC1 did not accumulate in the mutant without chlorophyll b, even though the NYC1 mRNA level was high after incubation in the dark. Quantification of the LHC protein showed no strong correlation between the levels of NYC1 and LHC proteins. However, the level of chlorophyll fluorescence of the dark adapted plant (F_o) was closely related to the accumulation of NYC1, suggesting that the NYC1 level is related to the energetically uncoupled LHC. These results and previous reports on the degradation of chlorophyllide a oxygenase suggest that the a feedforward and feedback network is included in chlorophyll cycle.