Whole-genome resequencing reveals Brassica napus origin and genetic loci involved in its improvement.

Brassica napus (2n = 4x = 38, AACC) is an important allopolyploid crop derived from interspecific crosses between Brassica rapa (2n = 2x = 20, AA) and Brassica oleracea (2n = 2x = 18, CC). However, no truly wild B. napus populations are known; its origin and improvement processes remain unclear. Here, we resequence 588 B. napus accessions. We uncover that the A subgenome may evolve from the ancestor of European turnip and the C subgenome may evolve from the common ancestor of kohlrabi, cauliflower, broccoli, and Chinese kale. Additionally, winter oilseed may be the original form of B. napus. Subgenome-specific selection of defense-response genes has contributed to environmental adaptation after formation of the species, whereas asymmetrical subgenomic selection has led to ecotype change. By integrating genome-wide association studies, selection signals, and transcriptome analyses, we identify genes associated with improved stress tolerance, oil content, seed quality, and ecotype improvement. They are candidates for further functional characterization and genetic improvement of B. napus.

Genome-Wide Identification and Characterization of NODULE-INCEPTION-Like Protein (NLP) Family Genes in Brassica napus.

NODULE-INCEPTION-like proteins (NLPs) are conserved, plant-specific transcription factors that play crucial roles in responses to nitrogen deficiency. However, the evolutionary relationships and characteristics of NLP family genes in Brassica napus are unclear. In this study, we identified 31 NLP genes in B. napus, including 16 genes located in the A subgenome and 15 in the C subgenome. Subcellular localization predictions indicated that most BnaNLP proteins are localized to the nucleus. Phylogenetic analysis suggested that the NLP gene family could be divided into three groups and that at least three ancient copies of NLP genes existed in the ancestor of both monocots and dicots prior to their divergence. The ancestor of group III NLP genes may have experienced duplication more than once in the Brassicaceae species. Three-dimensional structural analysis suggested that 14 amino acids in BnaNLP7-1 protein are involved in DNA binding, whereas no binding sites were identified in the two RWP-RK and PB1 domains conserved in BnaNLP proteins. Expression profile analysis indicated that BnaNLP genes are expressed in most organs but tend to be highly expressed in a single organ. For example, BnaNLP6 subfamily members are primarily expressed in roots, while the four BnaNLP7 subfamily members are highly expressed in leaves. BnaNLP genes also showed different expression patterns in response to nitrogen-deficient conditions. Under nitrogen deficiency, all members of the BnaNLP1/4/5/9 subfamilies were upregulated, all BnaNLP2/6 subfamily members were downregulated, and BnaNLP7/8 subfamily members showed various expression patterns in different organs. These results provide a comprehensive evolutionary history of NLP genes in B. napus, and insight into the biological functions of BnaNLP genes in response to nitrogen deficiency.

Genome-Wide Association and Transcriptome Analyses Reveal Candidate Genes Underlying Yield-determining Traits in Brassica napus.

Yield is one of the most important yet complex crop traits. To improve our understanding of the genetic basis of yield establishment, and to identify candidate genes responsible for yield improvement in Brassica napus, we performed genome-wide association studies (GWAS) for seven yield-determining traits [main inflorescence pod number (MIPN), branch pod number (BPN), pod number per plant (PNP), seed number per pod (SPP), thousand seed weight, main inflorescence yield (MIY), and branch yield], using data from 520 diverse B. napus accessions from two different yield environments. In total, we detected 128 significant single nucleotide polymorphisms (SNPs), 93 of which were revealed as novel by integrative analysis. A combination of GWAS and transcriptome sequencing on 21 haplotype blocks from samples pooled by four extremely high-yielding or low-yielding accessions revealed the differential expression of 14 crucial candiate genes (such as Bna.MYB83, Bna.SPL5, and Bna.ROP3) associated with multiple traits or containing multiple SNPs associated with the same trait. Functional annotation and expression pattern analyses further demonstrated that these 14 candiate genes might be important in developmental processes and biomass accumulation, thus affecting the yield establishment of B. napus. These results provide valuable information for understanding the genetic mechanisms underlying the establishment of high yield in B. napus, and lay the foundation for developing high-yielding B. napus varieties.

A combination of genome-wide association and transcriptome analysis reveals candidate genes controlling harvest index-related traits in Brassica napus.

Harvest index (HI), the ratio of seed mass to total biomass of the aboveground plant parts, is an important trait for harvestable yield of crops. Unfortunately, HI of Brassica napus is lower than that of other economically important crops. To identify candidate genes associated with high HI, a genome-wide association study of HI and four HI-related traits was conducted with 520 B. napus accessions cultivated in both Yunnan and Chongqing. We detected 294 single nucleotide polymorphisms significantly associated with the abovementioned traits, including 79 SNPs that affected two or more traits. Differentially expressed genes between extremely high- and low-HI accessions were identified in 8 tissues at two cultivated regions. Combination of linkage disequilibrium and transcriptome analyses revealed 33 functional candidate genes located within the confidence intervals of significant SNPs associated with more than one trait, such as SHOOT GRAVITROPISM 5 (Bna.SGR5), ATP-CITRATE LYASE A-3 (Bna.ACLA-3) and CAROTENOID CLEAVAGE DIOXYGENASE 1 (Bna.CCD1), their orthologs in the Arabidopsis thaliana have been shown to play key roles in photosynthesis, inflorescence, and silique development. Our results provide insight into the molecular mechanisms underlying establishment of high-HI B. napus and lay a foundation for characterization of candidate genes aimed at developing high-HI B. napus varieties.

Solution NMR and calorimetric analysis of Rem2 binding to the Ca2+ channel ß4 subunit: a low affinity interaction is required for inhibition of Cav2.1 Ca2+ currents.

Rem, Rad, Kir/Gem (RGK) proteins, including Rem2, mediate profound inhibition of high-voltage activated Ca(2+) channels containing intracellular regulatory ß subunits. All RGK proteins bind to voltage-gated Ca(2+) channel ß subunit (Cavß) subunits in vitro, but the necessity of the interaction for current inhibition remains controversial. This study applies NMR and calorimetric techniques to map the binding site for Rem2 on human Cavß4a and measure its binding affinity. Our experiments revealed 2 binding surfaces on the ß4 guanylate kinase domain contributing to a 156 ± 18 µM Kd interaction: a hydrophobic pocket lined by 4 critical residues (L173, N261, H262, and V303), mutation of any of which completely disrupted binding, and a nearby surface containing 3 residues (D206, L209, and D258) that when individually mutated decreased affinity. Voltage-gated Ca(2+) channel α1A subunit (Cav2.1) Ca(2+) currents were completely inhibited by Rem2 when co-expressed with wild-type Cavß4a, but were unaffected by Rem2 when coexpressed with a Cavß4a site 1 (L173A/V303A) or site 2 (D258A) mutant. These results provide direct evidence for a low-affinity Rem2/Cavß4 interaction and show definitively that the interaction is required for Cav2.1 inhibition.

¹H, ¹³C, and ¹5N backbone resonance assignments of the 37 kDa voltage-gated Ca²âº channel ß4 subunit core SH3-GK domains.

The ß subunit of the voltage-gated Ca(2+) channel (α1, α2δ, and ß subunits) is a member of the MAGUK family of proteins and plays an essential role in regulating Ca(2+) channel trafficking and gating. It also serves as a central interaction partner for various Ca(2+) channel regulatory proteins. We report here the nearly complete (1)H, (13)C, and (15)N backbone resonance assignments of the 37 kDa core SH3-GK domains of the ß4 subunit. This is the first report of solution assignments for ß subunits, and as such will lay the foundation for future investigations of interaction site mapping, functional dynamics, and protein complex structure determination.

Selective modulation of nuclear factor of activated T-cell function in restenosis by a potent bipartite peptide inhibitor.

RATIONALE: Nuclear factor of activated T-cells (NFAT) is importantly implicated in pathological cardiac remodeling and vascular lesion formation. NFAT functionality is mainly regulated by calcineurin, a Ca(2+)-dependent multi-effector phosphatase. Calcineurin inhibitors such as cyclosporine A (CsA) were shown to be effective in the treatment of restenosis and vascular inflammation but with adverse side effects. OBJECTIVE: This prompted the design of more selective inhibitors such as VIVIT and inhibitors of NFAT-calcineurin association, which unfortunately have a poor potency precluding clinical use. METHODS AND RESULTS: Here, we describe the rational design of a potent bipartite inhibitor of NFAT-calcineurin interaction, MCV1, which targets two separate calcineurin docking motifs. Modeling, site-directed mutagenesis, and functional studies demonstrated that MCV1 acts by allosteric modulation of calcineurin. Comparable to CsA, MCV1 prevents NFAT activation at nanomolar potency without impairing calcineurin phosphatase activity, nuclear factor-κB nuclear import, and general cell signaling. In contrast, CsA but not MCV1-activated basal level extracellular signal-regulated kinases activity and prevented nuclear import of calcineurin, independent of NFAT activation. In vivo MCV1 abrogated NFAT-mediated T-cell activation in a model of PMA-elicited peritonitis, whereas topical application of MCV1 markedly reduced neointima formation in a mouse model of restenosis. CONCLUSIONS: We designed a bipartite NFAT inhibitor that is more potent than VIVIT and more selective than CsA. MCV1 constitutes not only a powerful tool to unravel NFAT function but also a potential drug candidate for the treatment of diseases implicating NFAT activation.

The Ca2+ channel beta4c subunit interacts with heterochromatin protein 1 via a PXVXL binding motif.

The ß subunits of voltage-gated Ca(2+) channels are best known for their roles in regulating surface expression and gating of voltage-gated Ca(2+) channel α(1) subunits. Recent evidence, however, indicates that these proteins have a variety of Ca(2+) channel-independent functions. For example, on the molecular level, they regulate gene expression, and on the whole animal level, they regulate early cell movements in zebrafish development. In the present study, an alternatively spliced, truncated ß4 subunit (ß4c) is identified in the human brain and shown to be highly expressed in nuclei of vestibular neurons. Pull-down assays, nuclear magnetic resonance, and isothermal titration calorimetry demonstrate that the protein interacts with the chromo shadow domain (CSD) of heterochromatin protein 1γ. Site-directed mutagenesis reveals that the primary CSD interaction occurs through a ß4c C-terminal PXVXL consensus motif, adding the ß4c subunit to a growing PXVXL protein family with epigenetic responsibilities. These proteins have multiple nuclear functions, including transcription regulation (TIF1α) and nucleosome assembly (CAF1). An NMR-based two-site docking model of ß4c in complex with dimerized CSD is presented. Possible roles for the interaction are discussed.

Structural and functional characterization of the ga-substituted ferredoxin from Synechocystis sp. PCC6803, a mimic of the native protein.

In photosynthetic organisms, ferredoxin (Fd) interacts with many proteins, acting as a shuttle for electrons from Photosystem I to a group of enzymes involved in NADP(+) reduction, sulfur and nitrogen assimilation, and the regulation of carbon assimilation. The study of the dynamic interactions between ferredoxin and these enzymes by nuclear magnetic resonance is severely hindered by the paramagnetic [2Fe-2S] cluster of a ferredoxin. To establish whether ferredoxin in which the cluster has been replaced by Ga is a suitable diamagnetic mimic, the solution structure of Synechocystis Ga-substituted ferredoxin has been determined and compared with the structure of the native protein. The ensemble of 10 structures with the lowest energies has an average root-mean-square deviation of 0.30 +/- 0.05 A for backbone atoms and 0.65 +/- 0.04 A for all heavy atoms. Comparison of the NMR structure of GaFd with the crystal structure of the native Fd indicates that the general structural fold found for the native, iron-containing ferredoxin is conserved in GaFd. The ferredoxin contains a single gallium and no inorganic sulfide. The distortion of the metal binding loop caused by the single gallium substitution is small. The binding site on Fd for binding ferredoxin:NADP(+) reductase in solution, determined using GaFd, includes the metal binding loop and its surroundings, consistent with the crystal structures of related complexes. The results provide a structural justification for the use of the gallium-substituted analogue in interaction studies.

Ternary protein complex of ferredoxin, ferredoxin:thioredoxin reductase, and thioredoxin studied by paramagnetic NMR spectroscopy.

In oxygenic photosynthetic cells, carbon metabolism is regulated by a light-dependent redox signaling pathway through which the light signal is transmitted in the form of electrons via a redox chain comprising ferredoxin (Fd), ferredoxin:thioredoxin reductase (FTR), and thioredoxin (Trx). Trx affects the activity of a variety of enzymes via dithiol oxidation and reduction reactions. FTR reduces an intramolecular disulfide bridge of Trx, and Trx reduction involves a transient cross-link with FTR. NMR spectroscopy was used to investigate the interaction of Fd, FTR, and an m-type Trx. NMR titration experiments indicate that FTR uses distinct sites to bind Fd and Trx simultaneously to form a noncovalent ternary complex. The orientation of Trx-m relative to FTR was determined from the intermolecular paramagnetic broadening caused by the [4Fe-4S] cluster of FTR. Two models of the noncovalent binary complex of FTR/Trx-m based on the paramagnetic distance restraints were obtained. The models suggest that either a modest or major rotational movement of Trx must take place when the noncovalent binary complex proceeds to the covalent complex. This study demonstrates the complementarity of paramagnetic NMR and X-ray diffraction of crystals in the elucidation of dynamics in a transient protein complex.

Intermolecular dynamics studied by paramagnetic tagging.

Yeast cytochrome c and bovine adrenodoxin form a dynamic electron transfer complex, which is a pure encounter complex. It is demonstrated that the dynamic nature of the interaction can readily be probed by using a rigid lanthanide tag attached to cytochrome c. The tag, Caged Lanthanide NMR Probe 5, induces pseudocontact shifts and residual dipolar couplings and does not perturb the binding interface. Due to the dynamics in the complex, residual dipolar couplings in adrenodoxin are very small. Simulation shows that cytochrome c needs to sample a large part of the surface of adrenodoxin to explain the small degree of alignment observed for adrenodoxin. The applied method provides a simple and straightforward way to observe dynamics in protein complexes or domain-domain mobility without the need for external alignment media.

Dynamics in a pure encounter complex of two proteins studied by solution scattering and paramagnetic NMR spectroscopy.

In the general view of protein-complex formation, a transient and dynamic encounter complex proceeds to form a more stable, well-defined, and active form. In weak protein complexes, however, the encounter state can represent a significant population of the complex. The redox proteins adrenodoxin (Adx) and cytochrome c (C c) associate to form such a weak and short-lived complex, which is nevertheless active in electron transfer. To study the conformational freedom within the protein complex, the native complex has been compared to a cross-linked counterpart by using solution scattering and NMR spectroscopy. Oligomerization behavior of the native complex in solution revealed by small-angle X-ray scattering indicates a stochastic nature of complex formation. For the cross-linked complex, interprotein paramagnetic effects are observed, whereas for the native complex, extensive averaging occurs, consistent with multiple orientations of the proteins within the complex. Simulations show that C c samples about half of the surface area of adrenodoxin. It is concluded that the complex of Adx/C c is entirely dynamic and can be considered as a pure encounter complex.

Ferredoxin/ferredoxin-thioredoxin reductase complex: Complete NMR mapping of the interaction site on ferredoxin by gallium substitution.

The reduction of ferredoxin-thioredoxin reductase (FTR) by plant-type ferredoxin plays an important role in redox regulation in plants and cyanobacteria. Nuclear magnetic resonance (NMR) was used to map the binding sites on Synechocystis ferredoxin for FTR. A gallium-substituted structural analog of this [2Fe-2S] ferredoxin was obtained by reconstituting the apoprotein in a refolding buffer containing gallium. For the first time, the complete interaction interface of a [2Fe-2S] ferredoxin with a target enzyme has been mapped by NMR chemical shift perturbation with this diamagnetic structural analog.

NMR structure of a type IVb pilin from Salmonella typhi and its assembly into pilus.

The structure of the N-terminal-truncated Type IVb structural pilin (t-PilS) from Salmonella typhi was determined by NMR. Although topologically similar to the recently determined x-ray structure of pilin from Vibrio cholerae toxin-coregulated pilus, the only Type IVb pilin with known structure, t-PilS contains many distinct structural features. The protein contains an extra pair of beta-strands in the N-terminal alphabeta loop that align with the major beta-strands to form a continuous 7-stranded antiparallel beta-sheet. The C-terminal disulfide-bonded region of t-PilS is only half the length of that of toxin-coregulated pilus pilin. A model of S. typhi pilus has been proposed and mutagenesis studies suggested that residues on both the alphabeta loop and the C-terminal disulfide-bonded region of PilS might be involved in binding specificity of the pilus. This model structure reveals an exposed surface between adjacent subunits of PilS that could be a potential binding site for the cystic fibrosis transmembrane conductance regulator.