A combination of the F-box motif and kelch repeats defines a large Arabidopsis family of F-box proteins

In the sequences released by the Arabidopsis Genome Initiative (AGI), we have discovered a new large gene family (48 genes as of July 2000). A detailed computational and biochemical analysis of the predicted gene products reveals a novel family of plant F-box proteins, where the amino (N)-terminal F-box motif is followed by four kelch repeats and a characteristic carboxy-terminal domain. F-box proteins are an expanding family of eukaryotic proteins, which have been shown in some cases to be critical for the controlled degradation of cellular regulatory proteins via the ubiquitin pathway. The F-box motif of the At5g48990 gene product, a member of the family, was shown to be functionally active by its ability to mediate the in vitro interaction between At5g48990 and ASK1 proteins. F-box proteins specifically recruit the targets to be ubiquitinated, mainly through protein-protein interaction modules such as WD-40 domains or leucine-rich repeats (LRRs). The kelch repeats of the family described here form a potential protein-protein interaction domain, as molecular modelling of the kelch repeats according to the galactose oxidase crystal structure (the only solved structure containing kelch repeats) predicts a -propeller. The identification of this family of F-box proteins greatly expands the field of plant F-box proteins and suggests that controlled degradation of cellular proteins via the ubiquitin pathway could play a critical role in multiple plant cellular processes.     

EST-based in silico identification and in vitro test of antimicrobial peptides in Brassica napus

Background

Brassica napus is the third leading source of vegetable oil in the world after soybean and oil palm. The accumulation of gene sequences, especially expressed sequence tags (ESTs) from plant cDNA libraries, has provided a rich resource for genes discovery including potential antimicrobial peptides (AMPs). In this study, we used ESTs including those generated from B. napus cDNA libraries of seeds, pathogen-challenged leaves and deposited in the public databases, as a model, to perform in silico identification and consequently in vitro confirmation of putative AMP activities through a highly efficient system of recombinant AMP prokaryotic expression.

Results

In total, 35,788 were generated from cDNA libraries of pathogen-challenged leaves and 187,272 ESTs from seeds of B. napus, and the 644,998 ESTs of B. napus were downloaded from the EST database of PlantGDB. They formed 201,200 unigenes. First, all the known AMPs from the AMP databank (APD2 database) were individually queried against all the unigenes using the BLASTX program. A total of 972 unigenes that matched the 27 known AMP sequences in APD2 database were extracted and annotated using Blast2GO program. Among these unigenes, 237 unigenes from B. napus pathogen-challenged leaves had the highest ratio (1.15 %) in this unigene dataset, which is 13 times that of the unigene datasets of B. napus seeds (0.09 %) and 2.3 times that of the public EST dataset. About 87 % of each EST library was lipid-transfer protein (LTP) (32 % of total unigenes), defensin, histone, endochitinase, and gibberellin-regulated proteins. The most abundant unigenes in the leaf library were endochitinase and defensin, and LTP and histone in the pub EST library. After masking of the repeat sequence, 606 peptides that were orthologous matched to different AMP families were found. The phylogeny and conserved structural motifs of seven AMPs families were also analysed. To investigate the antimicrobial activities of the predicted peptides, 31 potential AMP genes belonging to different AMP families were selected to test their antimicrobial activities after bioinformatics identification. The AMP genes were all optimized according to Escherichia coli codon usage and synthetized through one-step polymerase chain reaction method. The results showed that 28 recombinant AMPs displayed expected antimicrobial activities against E. coli and Micrococcus luteus and Sclerotinia sclerotiorum strains.

Conclusion

The study not only significantly expanded the number of known/predicted peptides, but also contributed to long-term plant genetic improvement for increased resistance to diverse pathogens of B.napus. These results proved that the high-throughput method developed that combined an in silico procedure with a recombinant AMP prokaryotic expression system is considerably efficient for identification of new AMPs from genome or EST sequence databases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1849-x) contains supplementary material, which is available to authorized users.     

Genome-wide association study reveals the genetic architecture of flowering time in rapeseed (Brassica napus L.)

Flowering time adaptation is a major breeding goal in the allopolyploid species Brassica napus. To investigate the genetic architecture of flowering time, a genome-wide association study (GWAS) of flowering time was conducted with a diversity panel comprising 523 B. napus cultivars and inbred lines grown in eight different environments. Genotyping was performed with a Brassica 60K Illumina Infinium SNP array. A total of 41 single-nucleotide polymorphisms (SNPs) distributed on 14 chromosomes were found to be associated with flowering time, and 12 SNPs located in the confidence intervals of quantitative trait loci (QTL) identified in previous researches based on linkage analyses. Twenty-five candidate genes were orthologous to Arabidopsis thaliana flowering genes. To further our understanding of the genetic factors influencing flowering time in different environments, GWAS was performed on two derived traits, environment sensitivity and temperature sensitivity. The most significant SNPs were found near Bn-scaff_16362_1-p380982, just 13 kb away from BnaC09g41990D, which is orthologous to A. thaliana CONSTANS (CO), an important gene in the photoperiod flowering pathway. These results provide new insights into the genetic control of flowering time in B. napus and indicate that GWAS is an effective method by which to reveal natural variations of complex traits in B. napus.     

Diversity array technology markers: genetic diversity analyses and linkage map construction in rapeseed (Brassica napus L.)

We developed Diversity Array Technology (DArT) markers for application in genetic studies of Brassica napus and other Brassica species with A or C genomes. Genomic representation from 107 diverse genotypes of B. napus L. var. oleifera (rapeseed, AACC genomes) and B. rapa (AA genome) was used to develop a DArT array comprising 11 520 clones generated using PstI/BanII and PstI/BstN1 complexity reduction methods. In total, 1547 polymorphic DArT markers of high technical quality were identified and used to assess molecular diversity among 89 accessions of B. napus, B. rapa, B. juncea, and B. carinata collected from different parts of the world. Hierarchical cluster and principal component analyses based on genetic distance matrices identified distinct populations clustering mainly according to their origin/pedigrees. DArT markers were also mapped in a new doubled haploid population comprising 131 lines from a cross between spring rapeseed lines 'Lynx-037DH' and 'Monty-028DH'. Linkage groups were assigned on the basis of previously mapped simple sequence repeat (SSRs), intron polymorphism (IP), and gene-based markers. The map consisted of 437 DArT, 135 SSR, 6 IP, and 6 gene-based markers and spanned 2288 cM. Our results demonstrate that DArT markers are suitable for genetic diversity analysis and linkage map construction in rapeseed.     

Computational identification of novel microRNAs and targets in Brassica napus

MicroRNAs (miRNAs) are a newly discovered class of non-protein-coding small RNAs with roughly 22 nucleotide-long. Increasing evidence has shown that miRNAs play multiple roles in biological processes, including development, cell proliferation and apoptosis and stress responses. In this research, several approaches were combined to make computational prediction of potential miRNAs and their targets in Brassica napus. We used previously known miRNAs from Arabidopsis, rice and other plant species against both expressed sequence tags (EST) and genomic survey sequence (GSS) databases to search for potential miRNAs in B. napus. A total of 21 potential miRNAs were detected following a range of strict filtering criteria. Using these potential miRNA sequences, we could further blast the mRNA database and found 67 potential targets in this species. According to the mRNA target information provided by NCBI (http://www.ncbi.nlm.nih.gov/), most of the target mRNAs appeared to be involved in plant growth, development and stress responses. To validate the prediction of miRNAs in B. napus, we performed a RT-PCR based assay of mature miRNA expression. Five miRNAs were identified in response to auxin, cadmium stress and phosphate starvation. So far, little is known about experimental or computational identification of miRNA in B. napus species. To improve efficiency for blast search, we developed an implementation (miRNAassist) that can identify homologs of miRNAs and their targets, with high sensitivity and specificity. The program is allowed to be run on Windows Operation System platform. miRNAassist is freely available if required.     

Development and storage-protein synthesis in Brassica napus L. embryos in vivo and in vitro

Immature embryos of Brassica napus were cultured in vitro with and without various concentrations of germination inhibitors, and the progress of embryogeny was monitored by comparing accumulation of storage proteins in culture with the normal accumulation in seeds. The two major B. napus storage proteins (12S and 1.7S) were purified from seed extracts and analyzed by rocket immunoelectrophoresis (12S protein) or by sodium lauryl sulfate polyacrylamide gel electrophoresis (1.7S protein). During embryo development within seeds both the 12S and 1.7S proteins were first detected when the cotyledons were well developed (embryo dry weight, 0.4 mg), and each storage protein accumulated at an average rate of 26 g d^-1 during maximum deposition. Accumulation of the 1.7S protein stopped when the water content of the embryo began to decline (embryo DW, 2.7 mg), but accumulation of the 12S protein continued until seed maturity (embryo DW, 3.6 mg). At the end of embryo development the 12S and the 1.7S proteins comprised approx. 60 and 20% of the total salt-soluble protein, respectively. When embryos were removed from seeds at day 27, just as storage protein was starting to accumulate, and placed in culture on a basal medium, they precociously germinated within 3d, and incorporation of amino acids into the 12S storage protein dropped from 3% of total incorporation to less than 1%. If 10^-6 M abscisic acid (ABA) was included in the medium, amino-acid incorporation into the 12S protein increased from 3% of total incorporation when embryos were placed into culture to 18%, 5d later, and the accumulation rate (27.1±2.6 g embryo^-1 d^-1) matched the maximum rate observed in the seed. High osmotica, such as 0.29 M sucrose or mannitol, added to the basal medium, also inhibited precocious germination, but there was a lag period before 12S-protein synthesis rates equaled the rates on ABA media. These results indicate that some factor in the seed environment is necessary for storage-protein synthesis to proceed, and that ABA is a possible candidate.Abbreviations ABA abscisic acid - PAGE polyacrylamide gel electrophoresis - PMSF phenylmethylsulfonylfluoride - SDS sodium lauryl sulfate     

Phylogenic analysis of the genus Leishmania by cytochrome b gene sequencing

In a previous report (Luyo-Acero et al., 2004), we demonstrated that cytochrome b (Cyt b) gene analysis is an effective method for classifying several isolates of the genus Leishmania; hence, we have further applied this method to other Leishmania species in an effort to enhance the accuracy of the procedure and to construct a new phylogenic tree. In this study, a total of 30 Leishmania and Endotrypanum WHO reference strains, clinical isolates from our patients assigned to 28 strains (human and non-human pathogenic species) and two species of the genus Endotrypanum were analyzed. The Cyt b gene in each sample was amplified by PCR, and was then sequenced by several primers, as reported previously. The phylogenic tree was constructed based on the results obtained by the computer software MEGA v3.1 and PAUP* v4.0 Beta. The present phylogenic tree was almost identical to the traditional method of classification proposed by Lainson and Shaw (1987). However, it produces the following suggestions: (1) exclusion of L. (Leishmania) major from the L. (L.) tropica complex; (2) placement of L.tarentolae in the genus Sauroleishmania; (3) L. (L.) hertigi complex and L. (V.) equatorensis close to the genus Endotrypanum; (4) L. (L.) enrietti, defined as L. (L.) mexicana complex, placed in another position; and (5) L. (L.) turanica and L. (L.) arabica are located in an area far from human pathogenic Leishmania strains. Cyt b gene analysis is thus applicable to the analyzing phylogeny of the genus Leishmania and may be useful for separating non-human pathogenic species from human pathogenic species.     

Genome-wide and molecular evolution analysis of the subtilase gene family in Vitis vinifera

Background

Vitis vinifera (grape) is one of the most economically significant fruit crops in the world. The availability of the recently released grape genome sequence offers an opportunity to identify and analyze some important gene families in this species. Subtilases are a group of subtilisin-like serine proteases that are involved in many biological processes in plants. However, no comprehensive study incorporating phylogeny, chromosomal location and gene duplication, gene organization, functional divergence, selective pressure and expression profiling has been reported so far for the grape.

Results

In the present study, a comprehensive analysis of the subtilase gene family in V. vinifera was performed. Eighty subtilase genes were identified. Phylogenetic analyses indicated that these subtilase genes comprised eight groups. The gene organization is considerably conserved among the groups. Distribution of the subtilase genes is non-random across the chromosomes. A high proportion of these genes are preferentially clustered, indicating that tandem duplications may have contributed significantly to the expansion of the subtilase gene family. Analyses of divergence and adaptive evolution show that while purifying selection may have been the main force driving the evolution of grape subtilases, some of the critical sites responsible for the divergence may have been under positive selection. Further analyses of real-time PCR data suggested that many subtilase genes might be important in the stress response and functional development of plants.

Conclusions

Tandem duplications as well as purifying and positive selections have contributed to the functional divergence of subtilase genes in V. vinifera. The data may contribute to a better understanding of the grape subtilase gene family.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1116) contains supplementary material, which is available to authorized users.     

Genome-wide analysis of alternative splicing of pre-mRNA under salt stress in Arabidopsis

Background

Alternative splicing (AS) of precursor mRNA (pre-mRNA) is an important gene regulation process that potentially regulates many physiological processes in plants, including the response to abiotic stresses such as salt stress.

Results

To analyze global changes in AS under salt stress, we obtained high-coverage (~200 times) RNA sequencing data from Arabidopsis thaliana seedlings that were treated with different concentrations of NaCl. We detected that ~49% of all intron-containing genes were alternatively spliced under salt stress, 10% of which experienced significant differential alternative splicing (DAS). Furthermore, AS increased significantly under salt stress compared with under unstressed conditions. We demonstrated that most DAS genes were not differentially regulated by salt stress, suggesting that AS may represent an independent layer of gene regulation in response to stress. Our analysis of functional categories suggested that DAS genes were associated with specific functional pathways, such as the pathways for the responses to stresses and RNA splicing. We revealed that serine/arginine-rich (SR) splicing factors were frequently and specifically regulated in AS under salt stresses, suggesting a complex loop in AS regulation for stress adaptation. We also showed that alternative splicing site selection (SS) occurred most frequently at 4 nucleotides upstream or downstream of the dominant sites and that exon skipping tended to link with alternative SS.

Conclusions

Our study provided a comprehensive view of AS under salt stress and revealed novel insights into the potential roles of AS in plant response to salt stress.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-431) contains supplementary material, which is available to authorized users.     

Identification of genome-wide single nucleotide polymorphisms in allopolyploid crop Brassica napus

Background

Single nucleotide polymorphisms (SNPs) are the most common type of genetic variation. Identification of large numbers of SNPs is helpful for genetic diversity analysis, map-based cloning, genome-wide association analyses and marker-assisted breeding. Recently, identifying genome-wide SNPs in allopolyploid Brassica napus (rapeseed, canola) by resequencing many accessions has become feasible, due to the availability of reference genomes of Brassica rapa (2n = AA) and Brassica oleracea (2n = CC), which are the progenitor species of B. napus (2n = AACC). Although many SNPs in B. napus have been released, the objective in the present study was to produce a larger, more informative set of SNPs for large-scale and efficient genotypic screening. Hence, short-read genome sequencing was conducted on ten elite B. napus accessions for SNP discovery. A subset of these SNPs was randomly selected for sequence validation and for genotyping efficiency testing using the Illumina GoldenGate assay.

Results

A total of 892,536 bi-allelic SNPs were discovered throughout the B. napus genome. A total of 36,458 putative amino acid variants were located in 13,552 protein-coding genes, which were predicted to have enriched binding and catalytic activity as a result. Using the GoldenGate genotyping platform, 94 of 96 SNPs sampled could effectively distinguish genotypes of 130 lines from two mapping populations, with an average call rate of 92%.

Conclusions

Despite the polyploid nature of B. napus, nearly 900,000 simple SNPs were identified by whole genome resequencing. These SNPs were predicted to be effective in high-throughput genotyping assays (51% polymorphic SNPs, 92% average call rate using the GoldenGate assay, leading to an estimated >450 000 useful SNPs). Hence, the development of a much larger genotyping array of informative SNPs is feasible. SNPs identified in this study to cause non-synonymous amino acid substitutions can also be utilized to directly identify causal genes in association studies.     

The Fate of Chromosomes and Alleles in an Allohexaploid Brassica Population

Production of allohexaploid Brassica (2n = AABBCC) is a promising goal for plant breeders due to the potential for hybrid heterosis and useful allelic contributions from all three of the Brassica genomes present in the cultivated diploid (2n = AA, 2n = BB, 2n = CC) and allotetraploid (2n = AABB, 2n = AACC, and 2n = BBCC) crop species (canola, cabbages, mustards). We used high-throughput SNP molecular marker assays, flow cytometry, and fluorescent in situ hybridization (FISH) to characterize a population of putative allohexaploids derived from self-pollination of a hybrid from the novel cross (B. napus × B. carinata) × B. juncea to investigate whether fertile, stable allohexaploid Brassica can be produced. Allelic segregation in the A and C genomes generally followed Mendelian expectations for an F₂ population, with minimal nonhomologous chromosome pairing. However, we detected no strong selection for complete 2n = AABBCC chromosome complements, with weak correlations between DNA content and fertility (r² = 0.11) and no correlation between missing chromosomes or chromosome segments and fertility. Investigation of next-generation progeny resulting from one highly fertile F₂ plant using FISH revealed general maintenance of high chromosome numbers but severe distortions in karyotype, as evidenced by recombinant chromosomes and putative loss/duplication of A- and C-genome chromosome pairs. Our results show promise for the development of meiotically stable allohexaploid lines, but highlight the necessity of selection for 2n = AABBCC karyotypes.