Showing papers in "Molecular Genetics and Genomics in 2016"

repRNA: a web server for generating various feature vectors of RNA sequences.

Abstract: With the rapid growth of RNA sequences generated in the postgenomic age, it is highly desired to develop a flexible method that can generate various kinds of vectors to represent these sequences by focusing on their different features. This is because nearly all the existing machine-learning methods, such as SVM (support vector machine) and KNN (k-nearest neighbor), can only handle vectors but not sequences. To meet the increasing demands and speed up the genome analyses, we have developed a new web server, called "representations of RNA sequences" (repRNA). Compared with the existing methods, repRNA is much more comprehensive, flexible and powerful, as reflected by the following facts: (1) it can generate 11 different modes of feature vectors for users to choose according to their investigation purposes; (2) it allows users to select the features from 22 built-in physicochemical properties and even those defined by users' own; (3) the resultant feature vectors and the secondary structures of the corresponding RNA sequences can be visualized. The repRNA web server is freely accessible to the public at http://bioinformatics.hitsz.edu.cn/repRNA/ .

AtMYB12 regulates flavonoids accumulation and abiotic stress tolerance in transgenic Arabidopsis thaliana.

Abstract: In plants, transcriptional regulation is the most important tool for modulating flavonoid biosynthesis. The AtMYB12 gene from Arabidopsis thaliana has been shown to regulate the expression of key enzyme genes involved in flavonoid biosynthesis, leading to the increased accumulation of flavonoids. In this study, the codon-optimized AtMYB12 gene was chemically synthesized. Subcellular localization analysis in onion epidermal cells indicated that AtMYB12 was localized to the nucleus. Its overexpression significantly increased accumulation of flavonoids and enhanced salt and drought tolerance in transgenic Arabidopsis plants. Real-time quantitative PCR (qRT-PCR) analysis showed that overexpression of AtMYB12 resulted in the up-regulation of genes involved in flavonoid biosynthesis, abscisic acid (ABA) biosynthesis, proline biosynthesis, stress responses and ROS scavenging under salt and drought stresses. Further analyses under salt and drought stresses showed significant increases of ABA, proline content, superoxide dismutase (SOD) and peroxidase (POD) activities, as well as significant reduction of H2O2 and malonaldehyde (MDA) content. The results demonstrate the explicit role of AtMYB12 in conferring salt and drought tolerance by increasing the levels of flavonoids and ABA in transgenic Arabidopsis. The AtMYB12 gene has the potential to be used to enhance tolerance to abiotic stresses in plants.

iRSpot-GAEnsC: identifing recombination spots via ensemble classifier and extending the concept of Chou's PseAAC to formulate DNA samples.

Abstract: Meiotic recombination is vital for maintaining the sequence diversity in human genome. Meiosis and recombination are considered the essential phases of cell division. In meiosis, the genome is divided into equal parts for sexual reproduction whereas in recombination, the diverse genomes are combined to form new combination of genetic variations. Recombination process does not occur randomly across the genomes, it targets specific areas called recombination "hotspots" and "coldspots". Owing to huge exploration of polygenetic sequences in data banks, it is impossible to recognize the sequences through conventional methods. Looking at the significance of recombination spots, it is indispensable to develop an accurate, fast, robust, and high-throughput automated computational model. In this model, the numerical descriptors are extracted using two sequence representation schemes namely: dinucleotide composition and trinucleotide composition. The performances of seven classification algorithms were investigated. Finally, the predicted outcomes of individual classifiers are fused to form ensemble classification, which is formed through majority voting and genetic algorithm (GA). The performance of GA-based ensemble model is quite promising compared to individual classifiers and majority voting-based ensemble model. iRSpot-GAEnsC has achieved 84.46 % accuracy. The empirical results revealed that the performance of iRSpot-GAEnsC is not only higher than the examined algorithms but also better than existing methods in the literature developed so far. It is anticipated that the proposed model might be helpful for research community, academia and for drug discovery.

Unrevealed mosaicism in the next-generation sequencing era

Abstract: Mosaicism refers to the presence in an individual of normal and abnormal cells that are genotypically distinct and are derived from a single zygote. The incidence of mosaicism events in the human body is underestimated as the genotypes in the mosaic ratio, especially in the low-grade mosaicism, stay unrevealed. This review summarizes various research outcomes and diagnostic questions in relation to different types of mosaicism. The impact of both tested biological material and applied method on the mosaicism detection rate is especially highlighted. As next-generation sequencing technologies constitute a promising methodological solution in mosaicism detection in the coming years, revisions in current diagnostic protocols are necessary to increase the detection rate of the unrevealed mosaicism events. Since mosaicism identification is a complex process, numerous examples of multistep mosaicism investigations are presented and discussed.

Genome-wide identification and characterization of novel lncRNAs in Populus under nitrogen deficiency.

Abstract: Long non-coding RNAs (lncRNAs) have been identified as important regulatory factors of gene expression in eukaryotic species, such as Homo sapiens, Arabidopsis thaliana, and Oryza sativa. However, the systematic identification of potential lncRNAs in trees is comparatively rare. In particular, the characteristics, expression, and regulatory roles of lncRNAs in trees under nutrient stress remain largely unknown. A genome-wide strategy was used in this investigation to identify and characterize novel and low-nitrogen (N)-responsive lncRNAs in Populus tomentosa; 388 unique lncRNA candidates belonging to 380 gene loci were detected and only seven lncRNAs were found to belong to seven conserved non-coding RNA families indicating the majority of P. tomentosa lncRNAs are species-specific. In total, 126 lncRNAs were significantly altered under low-N stress; 8 were repressed, and 118 were induced. Furthermore, 9 and 5 lncRNAs were detected as precursors of 11 known and 14 novel Populus miRNAs, respectively, whereas 4 lncRNAs were targeted by 29 miRNAs belonging to 5 families, including 22 conserved and 7 non-conserved miRNAs. In addition, 15 antisense lncRNAs were identified to be generated from opposite strands of 14 corresponding protein-coding genes. In total, 111 protein-coding genes with regions complementary to 38 lncRNAs were also predicted with some lncRNAs corresponding to multiple genes and vice versa, and their functions were annotated, which further demonstrated the complex regulatory relationship between lncRNAs and protein-coding genes in plants. Moreover, an interaction network among lncRNAs, miRNAs, and mRNAs was investigated. These findings enrich our understanding of lncRNAs in Populus, expand the methods of miRNA identification. Our results present the first global characterization of lncRNAs and their potential target genes in response to nitrogen stress in trees, which provides more information on low-nutrition adaptation mechanisms in woody plants.

Quantitative trait loci affecting pathogen resistance and ripening of grapevines

Abstract: Grapevines (Vitis vinifera L.) form the basis of viticulture, and are susceptible to diseases such as downy mildew (Plasmopara viticola) and powdery mildew (Erysiphe necator). Therefore, successful viticulture programs require the use of pesticides. Breeding for resistance is the only eco-friendly solution. Marker-assisted selection is currently widely used for grapevine breeding. Consequently, traits of interest must be tagged with molecular markers linked to quantitative trait loci (QTL). We herein present our findings regarding genetic mapping and QTL analysis of resistance to downy and powdery mildew diseases in the progenies of the GF.GA-47-42 ('Bacchus' × 'Seyval') × 'Villard blanc' cross. Simple sequence repeats and single nucleotide polymorphisms of 151 individuals were analyzed. A map consisting of 543 loci was screened for QTL analyses based on phenotypic variations observed in plants grown in the field or under controlled conditions. A major QTL for downy mildew resistance was detected on chromosome 18. For powdery mildew resistance, a QTL was identified on chromosome 15. This QTL was replaced by a novel QTL on chromosome 18 in 2003 (abnormally high temperatures) and 2004. Subsequently, both QTLs functioned together. Additionally, variations in the timing of the onset of veraison, which is a crucial step during grape ripening, were studied to identify genomic regions affecting this trait. A major QTL was detected on linkage group 16, which was supplemented by a minor QTL on linkage group 18. This study provides useful information regarding novel QTL-linked markers relevant for the breeding of disease-resistant grapevines adapted to current climatic conditions.

Primordial dwarfism: overview of clinical and genetic aspects

Abstract: Primordial dwarfism is a group of genetic disorders which include Seckel Syndrome, Silver-Russell Syndrome, Microcephalic Osteodysplastic Primordial Dwarfism types I/III, II and Meier-Gorlin Syndrome. This genetic disorder group is characterized by intra-uterine growth retardation and post-natal growth abnormalities which occur as a result of disorganized molecular and genomic changes in embryonic stage and, thus, it represents a unique area to study growth and developmental abnormalities. Lot of research has been carried out on different aspects; however, a consolidated review that discusses an overall spectrum of this disorder is not accessible. Recent research in this area points toward important molecular and cellular mechanisms in human body that regulate the complexity of growth process. Studies have emerged that have clearly associated with a number of abnormal chromosomal, genetic and epigenetic alterations that can predispose an embryo to develop PD-associated developmental defects. Finding and associating such fundamental changes to its subtypes will help in re-examination of alleged functions at both cellular and developmental levels and thus reveal the intrinsic mechanism that leads to a balanced growth. Although such findings have unraveled a subtle understanding of growth process, we further require active research in terms of identification of reliable biomarkers for different subtypes as an immediate requirement for clinical utilization. It is hoped that further study will advance the understanding of basic mechanisms regulating growth relevant to human health. Therefore, this review has been written with an aim to present an overview of chromosomal, molecular and epigenetic modifications reported to be associated with different subtypes of this heterogenous disorder. Further, latest findings with respect to clinical and molecular genetics research have been summarized to aid the medical fraternity in their clinical utility, for diagnosing disorders where there are overlapping physical attributes and simultaneously inform about the latest developments in PD biology.

Genome-wide characterization and expression analysis of common bean bHLH transcription factors in response to excess salt concentration

Abstract: Members of basic helix-loop-helix (bHLH) gene family found in all eukaryotes play crucial roles in response to stress. Though, most eukaryotes carry the proteins of this family, biological functions of the most bHLH family members are not deeply evaluated in plants. In this study, we conducted a comprehensive genome-wide analysis of bHLH transcription factors in salt tolerant common bean. We identified 155 bHLH protein-encoding genes (PvbHLH) by using in silico comparative genomics tools. Based on the phylogenetic tree, PvbHLH genes were classified into 8 main groups with 21 subfamilies. Exon–intron analysis indicated that proteins belonging to same main groups exhibited a closely related gene structure. While, the PvbHLH gene family has been mainly expanded through segmental duplications, a total of 11 tandem duplication were detected. Genome-wide expression analysis of bHLH genes showed that 63 PvbHLH genes were differentially expressed in at least one tissue. Three of them displayed higher expression values in both leaf and root tissues. The in silico micro-RNA target transcript analyses revealed that totally 100 PvHLH genes targeted by 86 plant miRNAs. The most abundant transcripts, which were targeted by all 18 plant miRNA, were belonging to PvHLH-22 and PvHLH-44 genes. The expression of 16 PvbHLH genes in the root and leaf tissues of salt-stressed common bean was evaluated using qRT-PCR. Among them, two of PvbHLHs, PvbHLH-54, PvbHLH-148, were found to be up-regulated in both tissues in correlation with RNA-seq measurements. The results of this study could help improve understanding of biological functions of common bean bHLH family under salt stress. Additionally, it may provide basic resources for analyzing bHLH protein function for improving economic, agronomic and ecological benefit in common bean and other species.

NAC transcription factors play an important role in ethylene biosynthesis, reception and signaling of tomato fruit ripening.

Abstract: NAC proteins comprise a large family of transcription factors that play important roles in diverse physiological processes during development. To explore the role of NAC transcription factors in the ripening of fruits, we predicted the secondary and tertiary structure as well as regulative function of the SNAC4 (SlNAC48, Accession number: NM 001279348.2) and SNAC9 (SlNAC19, Accession number: XM 004236996.2) transcription factors in tomato. We found that the tertiary structure of SNAC9 was similar to that of ATNAP, which played an important role in the fruit senescence and was required for ethylene stimulation. Likewise, the bio-function prediction results indicated that SNAC4 and SNAC9 participated in various plant hormone signaling and senescence processes. More information about SNACs was obtained by the application of VIGS (virus-induced gene silencing). The silencing of SNAC4 and SNAC9 dramatically repressed the LeACS2, LeACS4 and LeACO1 expression, which consequently led to the inhibition of the ripening process. The silencing of SNACs down-regulated the mRNA levels of the ethylene perception genes and, at the same time, suppressed the expression of ethylene signaling-related genes except for LeERF2 which was induced by the silencing of SNAC4. The expressions of LeRIN were different in two silenced fruits. In addition, the silencing of SNAC4 reduced its mRNA level, while the silencing of SNAC9 induced its expression. Furthermore, the silencing of LeACS4, LeACO1 and LeERF2 reduced the expression of SNAC4 and SNAC9, while the silencing of NR induced the expression of all of them. In particular, these results indicate that SNAC transcription factors bind to the promoter of the ethylene synthesis genes in vitro. This experimental evidence demonstrates that SNAC4 and SNAC9 could positively regulate the tomato fruit ripening process by functioning upstream of ethylene synthesis genes. These outcomes will be helpful to provide a theoretical foundation for further exploring the tomato fruit ripening and senescence mechanism.

Transcriptome-wide identification of Camellia sinensis WRKY transcription factors in response to temperature stress.

Abstract: Tea plant [Camellia sinensis (L.) O. Kuntze] is a leaf-type healthy non-alcoholic beverage crop, which has been widely introduced worldwide. Tea is rich in various secondary metabolites, which are important for human health. However, varied climate and complex geography have posed challenges for tea plant survival. The WRKY gene family in plants is a large transcription factor family that is involved in biological processes related to stress defenses, development, and metabolite synthesis. Therefore, identification and analysis of WRKY family transcription factors in tea plant have a profound significance. In the present study, 50 putative C. sinensis WRKY proteins (CsWRKYs) with complete WRKY domain were identified and divided into three Groups (Group I-III) on the basis of phylogenetic analysis results. The distribution of WRKY family transcription factors among plantae, fungi, and protozoa showed that the number of WRKY genes increased in higher plant, whereas the number of these genes did not correspond to the evolutionary relationships of different species. Structural feature and annotation analysis results showed that CsWRKY proteins contained WRKYGQK/WRKYGKK domains and C2H2/C2HC-type zinc-finger structure: D-X18-R-X1-Y-X2-C-X4-7-C-X23-H motif; CsWRKY proteins may be associated with the biological processes of abiotic and biotic stresses, tissue development, and hormone and secondary metabolite biosynthesis. Temperature stresses suggested that the candidate CsWRKY genes were involved in responses to extreme temperatures. The current study established an extensive overview of the WRKY family transcription factors in tea plant. This study also provided a global survey of CsWRKY transcription factors and a foundation of future functional identification and molecular breeding.

Calcium-dependent protein kinase (CDPK) and CDPK-related kinase (CRK) gene families in tomato: genome-wide identification and functional analyses in disease resistance

Abstract: Calcium-dependent protein kinases (CDPKs) and CDPK-related kinases (CRKs) play multiple roles in plant. Nevertheless, genome-wide identification of these two families is limited to several plant species, and role of CRKs in disease resistance remains unclear. In this study, we identified the CDPK and CRK gene families in genome of the economically important crop tomato (Solanum lycopersicum L.) and analyzed their function in resistance to various pathogens. Twenty-nine CDPK and six CRK genes were identified in tomato genome. Both SlCDPK and SlCRK proteins harbored an STKc_CAMK type protein kinase domain, while only SlCDPKs contained EF-hand type Ca(2+) binding domain(s). Phylogenetic analysis revealed that plant CRK family diverged early from CDPKs, and shared a common ancestor gene with subgroup IV CDPKs. Subgroup IV SlCDPK proteins were basic and their genes contained 11 introns, which were distinguished from other subgroups but similar to CRKs. Subgroup I SlCDPKs generally did not carry an N-terminal myristoylation motif while those of the remaining subgroups and SlCRKs universally did. SlCDPK and SlCRK genes were differently responsive to pathogenic stimuli. Furthermore, silencing analyses demonstrated that SlCDPK18 and SlCDPK10 positively regulated nonhost resistance to Xanthomonas oryzae pv. oryzae and host resistance to Pseudomonas syringae pv. tomato (Pst) DC3000, respectively, while SlCRK6 positively regulated resistance to both Pst DC3000 and Sclerotinia sclerotiorum in tomato. In conclusion, CRKs apparently evolved from CDPK lineage, SlCDPK and SlCRK genes regulate a wide range of resistance and SlCRK6 is the first CRK gene proved to function in plant disease resistance.

Characters, functions and clinical perspectives of long non-coding RNAs.

Abstract: It is well established that most of the human genome and those of other mammals and plants are transcribed into RNA without protein-coding capacity, which we define as non-coding RNA. From siRNA to microRNA, whose functions and features have been well characterized, non-coding RNAs have been a popular topic in life science research over the last decade. Long non-coding RNAs (lncRNAs), however, as a novel class of transcripts, are distinguished from these other small RNAs. Recent studies have revealed a diverse population of lncRNAs with different sizes and functions across different species. These populations are expressed dynamically and act as important regulators in a variety of biological processes, especially in gene expression. Nevertheless, the functions and mechanisms of most lncRNAs remain unclear. In this review, we present recent progress in the identification of lncRNAs, their functions and molecular mechanisms, their roles in human diseases, their potential diagnostic and therapeutic applications as well as newer technologies for identifying deregulated lncRNAs in disease tissues.

Genome-wide identification and comparative analysis of grafting-responsive mRNA in watermelon grafted onto bottle gourd and squash rootstocks by high-throughput sequencing

Abstract: Grafting is an important agricultural technique widely used to improve plant growth, yield, and adaptation to either biotic or abiotic stresses. However, the molecular mechanisms underlying grafting-induced physiological processes remain unclear. Watermelon (Citrullus lanatus L.) is an important horticultural crop worldwide. Grafting technique is commonly used in watermelon production for improving its tolerance to stresses, especially to the soil-borne fusarium wilt disease. In the present study, we used high-throughput sequencing to perform a genome-wide transcript analysis of scions from watermelon grafted onto bottle gourd and squash rootstocks. Our transcriptome and digital gene expression (DGE) profiling data provided insights into the molecular aspects of gene regulation in grafted watermelon. Compared with self-grafted watermelon, there were 787 and 3485 genes differentially expressed in watermelon grafted onto bottle gourd and squash rootstocks, respectively. These genes were associated with primary and secondary metabolism, hormone signaling, transcription factors, transporters, and response to stimuli. Grafting led to changes in expression of these genes, suggesting that they may play important roles in mediating the physiological processes of grafted seedlings. The potential roles of the grafting-responsive mRNAs in diverse biological and metabolic processes were discussed. Obviously, the data obtained in this study provide an excellent resource for unraveling the mechanisms of candidate genes function in diverse biological processes and in environmental adaptation in a graft system.

Identifying N 6 -methyladenosine sites in the Arabidopsis thaliana transcriptome

Abstract: N 6-Methyladenosine (m6A) plays important roles in many biological processes. The knowledge of the distribution of m6A is helpful for understanding its regulatory roles. Although the experimental methods have been proposed to detect m6A, the resolutions of these methods are still unsatisfying especially for Arabidopsis thaliana. Benefitting from the experimental data, in the current work, a support vector machine-based method was proposed to identify m6A sites in A. thaliana transcriptome. The proposed method was validated on a benchmark dataset using jackknife test and was also validated by identifying strain-specific m6A sites in A. thaliana. The obtained predictive results indicate that the proposed method is quite promising. For the convenience of experimental biologists, an online webserver for the proposed method was built, which is freely available at http://lin.uestc.edu.cn/server/M6ATH . These results indicate that the proposed method holds a potential to become an elegant tool in identifying m6A site in A. thaliana.

SmARF8, a transcription factor involved in parthenocarpy in eggplant

Abstract: Parthenocarpic fruit is a very attractive trait for consumers and especially in eggplants where seeds can lead to browning of the flesh and bitterness. However, the molecular mechanisms underlying parthenocarpy in eggplant still remain unknown. Some auxin response factors have been previously shown in model species, such as Arabidopsis and tomato, to play an important role in such a process. Here, we have identified a natural parthenocarpic mutant and showed that ARF8 from eggplant (SmARF8), is down-regulated in buds compared to wild-type plants. Further characterization of SmARF8 showed that it is a nuclear protein and an active transcriptional regulator. We determined that amino acids 629–773 of SmARF8 act as the transcriptional activation domain, the C terminus of SmARF8 is the protein-binding domain, and that SmARF8 might form homodimers. Expression analysis in eggplant showed that SmARF8 is expressed ubiquitously in all tissues and organs and is responsive to auxin. Eggplant transgenic lines harboring RNA interference of SmARF8 exhibited parthenocarpy in unfertilized flowers, suggesting that SmARF8 negatively regulates fruit initiation. Interestingly, SmARF8-overexpressing Arabidopsis lines also induced parthenocarpy. These results indicate that SmARF8 could affect the dimerization of auxin/indole acetic acid repressors with SmARF8 via domains III and IV and thus induce fruit development. Furthermore, the introduction of SmARF8 full-length cDNA could partially complement the parthenocarpic phenotypes in Arabidopsisarf8-1 and arf8-4 mutants. Collectively, our results demonstrate that SmARF8 may act as a key negative regulator involved in parthenocarpic fruit development of eggplant. These findings give more insights into the conserved mechanisms leading to parthenocarpy in which auxin signaling plays a pivotal role, and provide potential target for eggplant breeding.

Identification of compound-protein interactions through the analysis of gene ontology, KEGG enrichment for proteins and molecular fragments of compounds.

Abstract: Compound–protein interactions play important roles in every cell via the recognition and regulation of specific functional proteins. The correct identification of compound–protein interactions can lead to a good comprehension of this complicated system and provide useful input for the investigation of various attributes of compounds and proteins. In this study, we attempted to understand this system by extracting properties from both proteins and compounds, in which proteins were represented by gene ontology and KEGG pathway enrichment scores and compounds were represented by molecular fragments. Advanced feature selection methods, including minimum redundancy maximum relevance, incremental feature selection, and the basic machine learning algorithm random forest, were used to analyze these properties and extract core factors for the determination of actual compound–protein interactions. Compound–protein interactions reported in The Binding Databases were used as positive samples. To improve the reliability of the results, the analytic procedure was executed five times using different negative samples. Simultaneously, five optimal prediction methods based on a random forest and yielding maximum MCCs of approximately 77.55 % were constructed and may be useful tools for the prediction of compound–protein interactions. This work provides new clues to understanding the system of compound–protein interactions by analyzing extracted core features. Our results indicate that compound–protein interactions are related to biological processes involving immune, developmental and hormone-associated pathways.

Dynamic transcriptome analysis reveals AP2/ERF transcription factors responsible for cold stress in rapeseed ( Brassica napus L.)

Abstract: The APETALA2/ethylene response factor (AP2/ERF) transcription factor (TF) superfamily plays an important regulatory role in signal transduction of the plant responses to various stresses including low temperature. Significant progress has been made in understanding the mechanism of cold resistance in Brassica napus, an important oilseed crop. However, comprehensive studies on the induction and activity of these TFs under low temperature have been lacking. In this study, 132 AP2/ERF genes were identified by transcriptome sequencing of rapeseed leaves exposed to 0, 2, 6, 12, and 24 h of low (4 °C) temperature stress. The genes were classified into 4 subfamilies (AP2, DREB, ERF, and RAV) and 13 subgroups, among which the DREB subfamily and ERF subfamily contained 114 genes, no genes were assigned to soloist or DREB A3 subgroups. One hundred and eighteen genes were located on chromosomes A1 to C9. GO functional analysis and promoter sequence analysis revealed that these genes are involved in many molecular pathways that may enhance cold resistance in plants, such as the low-temperature responsiveness, methyl jasmonate, abscisic acid, and ethylene-responsiveness pathways. Their expression patterns revealed dynamic control at different times following initiation of cold stress; the RAV and DREB subfamilies were expressed at the early stage of cold stress, whereas the AP2 subfamily was expressed later. Quantitative PCR analyses of 13 cold-induced AP2/ERF TFs confirmed the accuracy of above results. This study is the first dynamic analysis of the AP2/ERF TFs responsible for cold stress in rapeseed. These findings will serve as a reference for future functional research on transcription in rapeseed.

GDSL esterase/lipase genes in Brassica rapa L.: genome-wide identification and expression analysis.

Abstract: GDSL esterase/lipase proteins (GELPs), a very large subfamily of lipolytic enzymes, have been identified in microbes and many plants, but only a few have been characterized with respect to their roles in growth, development, and stress responses. In Brassica crops, as in many other species, genome-wide systematic analysis and functional studies of these genes are still lacking. As a first step to study their function in B. rapa ssp. pekinensis (Chinese cabbage), we comprehensively identified all GELP genes in the genome. We found a total of 121 Brassica rapa GDSL esterase/lipase protein genes (BrGELPs), forming three clades in the phylogenetic analysis (two major and one minor), with an asymmetrical chromosomal distribution. Most BrGELPs possess four strictly conserved residues (Ser-Gly-Asn-His) in four separate conserved regions, along with short conserved and clade-specific blocks, suggesting functional diversification of these proteins. Detailed expression profiling revealed that BrGELPs were expressed in various tissues, including floral organs, implying that BrGELPs play diverse roles in various tissues and during development. Ten percent of BrGELPs were specifically expressed in fertile buds, rather than male-sterile buds, implying their involvement in pollen development. Analyses of EXL6 (extracellular lipase 6) expression and its co-expressed genes in both B. rapa and Arabidopsis, as well as knockdown of this gene in Arabidopsis, revealed that this gene plays an important role in pollen development in both species. The data described in this study will facilitate future investigations of other BrGELP functions.

Evolution of the KCS gene family in plants: the history of gene duplication, sub/neofunctionalization and redundancy

Abstract: Very long-chain fatty acids (VLCFAs) play an important role in the survival and development of plants, and VLCFA synthesis is regulated by β-ketoacyl-CoA synthases (KCSs), which catalyze the condensation of an acyl-CoA with malonyl-CoA. Here, we present a genome-wide survey of the genes encoding these enzymes, KCS genes, in 28 species (26 genomes and two transcriptomes), which represents a large phylogenetic scale, and also reconstruct the evolutionary history of this gene family. KCS genes were initially single-copy genes in the green plant lineage; duplication resulted in five ancestral copies in land plants, forming five fundamental monophyletic groups in the phylogenetic tree. Subsequently, KCS genes duplicated to generate 11 genes of angiosperm origin, expanding up to 20-30 members in further-diverged angiosperm species. During this process, tandem duplications had only a small contribution, whereas polyploidy events and large-scale segmental duplications appear to be the main driving force. Accompanying this expansion were variations that led to the sub- and neofunctionalization of different members, resulting in specificity that is likely determined by the 3-D protein structure. Novel functions involved in other physiological processes emerged as well, though redundancy is also observed, largely among recent duplications. Conserved sites and variable sites of KCS proteins are also identified by statistical analysis. The variable sites are likely to be involved in the emergence of product specificity and catalytic power, and conserved sites are possibly responsible for the preservation of fundamental function.

Selection of reference genes for analysis of stress-responsive genes after challenge with viruses and temperature changes in the silkworm Bombyx mori.

Abstract: Viruses and high temperature (HT) are the primary threats to silkworms. Changes in the expression of stress-response genes can be measured using quantitative polymerase chain reaction (qPCR) after exposure to viruses or HT. However, appropriate reference genes (RGs) for qPCR data normalization have not been established in this organism. In this study, we summarized the RGs used in the previous silkworm studies after infection with Bombyx mori nucleopolyhedrovirus (BmNPV), B. mori cytoplasmic polyhedrosis virus (BmCPV), or B. mori densovirus (BmDNV) or after HT treatment. The expression levels of these RGs were extracted from silkworm transcriptome data to screen for candidate RGs that were unaffected by the experimental conditions. Actin-1 (A1), actin-3 (A3), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and translation initiation factor 4a (TIF-4A) were selected for further qPCR verification. The results of RNA-seq and qPCR showed that GAPDH and TIF-4A were suitable RGs after BmNPV challenge or HT stress, whereas TIF-4A was an appropriate RG for BmCPV or BmDNV-Z challenge in silkworms. These results suggested that TIF-4A may be the most appropriate RG for gene expression analysis after challenge with viruses or HT in silkworms.

Genome-wide identification of the expansin gene family in tobacco (Nicotiana tabacum).

Abstract: Expansins are pH-dependent cell wall loosening proteins which form a large family in plants. They have been shown to be involved in various developmental processes and been implicated in enabling plants' ability to absorb nutrients from the soil as well as conferring biotic and abiotic stress resistances. It is therefore clear that they can be potential targets in genetic engineering for crop improvement. Tobacco (Nicotiana tabacum) is a major crop species as well as a model organism. Considering that only a few tobacco expansins have been studied, a genome-wide analysis of the tobacco expansin gene family is necessary. In this study, we identified 52 expansins in tobacco, which were classified into four subfamilies: 36 NtEXPAs, 6 NtEXPBs, 3 NtEXLAs and 7 NtEXLBs. Compared to other species, the NtEXLB subfamily size was relatively larger. Phylogenetic analysis showed that the 52 tobacco expansins were divided into 13 subgroups. Gene structure analysis revealed that genes within subfamilies/subgroups exhibited similar characteristics such as gene structure and protein motif arrangement. Whole-genome duplication and tandem duplication events may have played important roles in the expanding of tobacco expansins. Cis-Acting element analysis revealed that each expansin gene was regulated or several expansin genes were co-regulated by both internal and environmental factors. 35 of these genes were identified as being expressed according to a microarray analysis. In contrast to most NtEXPAs which had higher expression levels in young organs, NtEXLAs and NtEXLBs were preferentially expressed in mature or senescent tissues, suggesting that they might play different roles in different organs or at different developmental stages. As the first step towards genome-wide analysis of the tobacco expansin gene family, our work provides solid background information related to structure, evolution and expression as well as regulatory cis-acting elements of the tobacco expansins. This information will provide a strong foundation for cloning and functional exploration of expansin genes in tobacco.

Genome-wide identification and expression analysis of the expansin gene family in tomato.

Abstract: Plant expansins are capable of inducing pH-dependent cell wall extension and stress relaxation. They may be useful as targets for crop improvement to enhance fruit development and stress resistance. Tomato is a major agricultural crop and a model plant for studying fruit development. Because only some tomato expansins have been studied, a genome-wide analysis of the tomato expansin family is necessary. In this study, we identified 25 SlEXPAs, eight SlEXPBs, one SlEXLA, four SlEXLBs, and five short homologs in the tomato genome. 25 of these genes were identified as being expressed. Bioinformatic analysis showed that although tomato expansins share similarities with those from other plants, they also exhibit specific features regarding genetic structure and amino acid sequences, which indicates a unique evolutionary process. Segmental and tandem duplication events have played important roles in expanding the tomato expansin family. Additionally, the 3-exon/2-intron structure may form the basic organization of expansin genes. We identified new expansin genes preferentially expressed in fruits (SlEXPA8, SlEXPB8, and SlEXLB1), roots (SlEXPA9, SlEXLB2, and SlEXLB4), and floral organs. Among the analyzed genes those that were inducible by hormone or stress treatments, including SlEXPA3, SlEXPA7, SlEXPB1-B2, SlEXPB8, SlEXLB1-LB2, and SlEXLB4. Our findings may further clarify the biological activities of tomato expansins, especially those related to fruit development and stress resistance, and contribute to the genetic modification of tomato plants to improve crop quality and yield.

Genome-wide identification of WRKY family genes in peach and analysis of WRKY expression during bud dormancy

Abstract: Bud dormancy in deciduous fruit trees is an important adaptive mechanism for their survival in cold climates. The WRKY genes participate in several developmental and physiological processes, including dormancy. However, the dormancy mechanisms of WRKY genes have not been studied in detail. We conducted a genome-wide analysis and identified 58 WRKY genes in peach. These putative genes were located on all eight chromosomes. In bioinformatics analyses, we compared the sequences of WRKY genes from peach, rice, and Arabidopsis. In a cluster analysis, the gene sequences formed three groups, of which group II was further divided into five subgroups. Gene structure was highly conserved within each group, especially in groups IId and III. Gene expression analyses by qRT-PCR showed that WRKY genes showed different expression patterns in peach buds during dormancy. The mean expression levels of six WRKY genes (Prupe.6G286000, Prupe.1G393000, Prupe.1G114800, Prupe.1G071400, Prupe.2G185100, and Prupe.2G307400) increased during endodormancy and decreased during ecodormancy, indicating that these six WRKY genes may play a role in dormancy in a perennial fruit tree. This information will be useful for selecting fruit trees with desirable dormancy characteristics or for manipulating dormancy in genetic engineering programs.

Genome-wide analysis and expression patterns of ZF-HD transcription factors under different developmental tissues and abiotic stresses in Chinese cabbage

Abstract: The ZF-HD gene family plays an important role in plant developmental processes and stress responses. However, the function of the ZF-HD genes in Chinese cabbage remains largely unknown. Chinese cabbage (Brassica rapa ssp. pekinensis) is a member of one of the most important leaf vegetables grown worldwide. The entire Chinese cabbage genome sequence has been determined, and more than forty thousand proteins have been identified to date. In this study, 31 ZF-HD genes were identified in Chinese cabbage. We show here that the BraZF-HD genes could be categorized into ZHD and MIF subfamilies. Among them, ZHD genes are plant-specific, nearly all intronless, and related to MINI ZINC FINGER genes that possess only the zinc finger. Phylogenetic analysis suggested that ZHDs have expanded considerably during angiosperm evolution. In addition, the ZHD group has 24 members, which is twice as much as the Arabidopsis ZHD group, indicating that the Chinese cabbage ZHD genes have been retained more frequently than other group genes. Real-time PCR analysis showed that most of BraZF-HD genes are preferentially expressed in flower. Furthermore, most of these genes are significantly induced under photoperiod or vernalization conditions, as well as abiotic stresses. Thereby implying that they may play important roles in these processes. This study provides insight into the evolution of ZF-HD genes in Chinese cabbage genome and may aid efforts to further characterize the function of these predicted ZF-HD genes in flowering and resistance.

De novo assembly and characterization of the leaf, bud, and fruit transcriptome from the vulnerable tree Juglans mandshurica for the development of 20 new microsatellite markers using Illumina sequencing

Abstract: Manchurian walnut (Juglans mandshurica Maxim.) is a vulnerable, temperate deciduous tree valued for its wood and nut, but transcriptomic and genomic data for the species are very limited. Next generation sequencing (NGS) has made it possible to develop molecular markers for this species rapidly and efficiently. Our goal is to use transcriptome information from RNA-Seq to understand development in J. mandshurica and develop polymorphic simple sequence repeats (SSRs, microsatellites) to understand the species' population genetics. In this study, more than 47.7 million clean reads were generated using Illumina sequencing technology. De novo assembly yielded 99,869 unigenes with an average length of 747 bp. Based on sequence similarity search with known proteins, a total of 39,708 (42.32 %) genes were identified. Searching against the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG) identified 15,903 (16.9 %) unigenes. Further, we identified and characterized 63 new transcriptome-derived microsatellite markers. By testing the markers on 4 to 14 individuals from four populations, we found that 20 were polymorphic and easily amplified. The number of alleles per locus ranged from 2 to 8. The observed and expected heterozygosity per locus ranged from 0.209 to 0.813 and 0.335 to 0.842, respectively. These twenty microsatellite markers will be useful for studies of population genetics, diversity, and genetic structure, and they will undoubtedly benefit future breeding studies of this walnut species. Moreover, the information uncovered in this research will also serve as a useful genetic resource for understanding the transcriptome and development of J. mandshurica and other Juglans species.

Comparison of skeletal muscle miRNA and mRNA profiles among three pig breeds

Abstract: The pig is an important source of animal protein, and is also an ideal model for human disease. There are significant differences in growth rate, muscle mass, and meat quality between different breeds. To understand the molecular mechanisms underlying porcine skeletal muscle phenotypes, we performed mRNA and miRNA profiling of muscle from three different breeds of pig, Landrace (lean-type), Tongcheng (obese-type), and Wuzhishan (mini-type) by Solexa sequencing. Forty-three genes and 106 miRNAs were differentially expressed between Landrace and Tongcheng pigs, 92 genes and 151 miRNAs were differentially expressed between Tongcheng and Wuzhishan pigs, and 145 genes and 156 miRNAs were differential expressed between Landrace and Wuzhishan pigs. Gene ontology analysis suggested that genes differentially expressed between Landrace and Tongcheng pigs were mainly involved in the biological processes of oxidative stress and muscle organ development. Meanwhile, for Tongcheng vs Wuzhishan and Landrace vs Wuzhishan pigs, the differentially expressed genes were involved in fatty acid metabolism, oxidative stress, muscle contraction, and muscle organ development, processes that are closely related to meat quality. To investigate the molecular mechanisms underlying meat quality diversity based on differentially expressed genes and miRNAs, interaction networks were constructed, according to target prediction results and integration analysis of up-regulated genes with down-regulated miRNAs or down-regulated genes with up-regulated miRNAs. Our findings identify candidate genes and miRNAs associated with muscle development and indicate their potential roles in muscle phenotype variance between different pig breeds. These results serve as a foundation for further studies on muscle development and molecular breeding.

SNP-based discovery of salinity-tolerant QTLs in a bi-parental population of rice (Oryza sativa).

Abstract: Breeding for salt tolerance is the most promising approach to enhance the productivity of saline prone areas. However, polygenic inheritance of salt tolerance in rice acts as a bottleneck in conventional breeding for salt tolerance. Hence, we set our goals to construct a single nucleotide polymorphism (SNP)-based molecular map employing high-throughput SNP marker technology and to investigate salinity tolerant QTLs with closest flanking markers using an elite rice background. Seedling stage salinity responses were assessed in a population of 281 recombinant inbred lines (RILs) derived from the cross between At354 (salt tolerant) and Bg352 (salt susceptible), by 11 morpho-physiological indices under a hydroponic system. Selected extreme 94 RILs were genotyped using Illumina Infinium rice 6K SNP array and densely saturated molecular map spanning 1460.81 cM of the rice genome with an average interval of 1.29 cM between marker loci was constructed using 1135 polymorphic SNP markers. The results revealed 83 significant QTLs for 11 salt responsive traits explaining 12.5–46.7 % of phenotypic variation in respective traits. Of them, 72 QTLs responsible for 10 traits were co-localized together forming 14 QTL hotspots at 14 different genomic regions. The all QTL hotspots were flanked less than 1 Mb intervals and therefore the SNP loci associated with these QTL hotspots would be important in candidate gene discovery for salt tolerance.

Genetic analysis for rice grain quality traits in the YVB stable variant line using RAD-seq

Abstract: The future of rice breeding will likely be built on the basis of the further utilization of heterosis between elite cultivars and genetic resources from distant subspecies of rice. Previous studies have proved that exogenous genomic DNA transformation methods can be used to transfer genetic information from distant relatives (donor) into cultivated rice (recipient). However, the mechanism underlying this form of genetic transfer is poorly characterized, and the genes that cause the phenotypic changes in these variants are typically difficult to identify. This study examined YVB, a stable variant line with greatly improved grain quality traits that was derived from an indica variety (V20B) by transferring genomic DNA of O.minuta through the "spike-stalk injection method (SIM)". We used restriction-site associated DNA sequencing technology (RAD-seq) to evaluate a population of BC1F5 backcross lines (YVB × V20B); the RAD-seq data were used to construct a genetic linkage map with high-density SNPs for use in association analysis exploring genotype-phenotype relationships at the whole-genome level. A total of 17 quantitative trait loci (QTLs) for rice quality traits were mapped to chromosomes 3, 5, 6, 8, and 9. 8 major QTLs controlling different phenotypic variations were mapped to the same region of chromosome 5. This region contained the GS5 gene for grain weight and the qSW5/GW5 gene for grain width. This study provides new resources and insights into the molecular mechanisms of grain trait phenotypic variation and the transmission of genetic information via the introduction of genomic DNA to a distantly related crop relative species.

Comparative transcriptome analyses of deltamethrin-susceptible and -resistant Culex pipiens pallens by RNA-seq

Abstract: The widespread and improper use of pyrethroid insecticides, such as deltamethrin, has resulted in the evolution of resistance in many mosquito species, including Culex pipiens pallens. With the development of high-throughput sequencing, it is possible to massively screen pyrethroid resistance-associated gene. In this study, we used Illumina-Solexa transcriptome sequencing to identify genes that are expressed differently in deltamethrin-susceptible and -resistant strains of Culex pipiens pallens as a critical knowledge base for further studies. A total of 4,961,197,620 base pairs and 55,124,418 reads were sequenced, mapped to the Culex quinquefasciatus genome and assembled into 17,679 known genes. We recorded 1826 significantly differentially expressed genes (DEGs). Among them, 1078 genes were up-regulated and 748 genes were down-regulated in the deltamethrin-resistant strain compared to -susceptible strain. These DEGs contained cytochrome P450 s, cuticle proteins, UDP-glucuronosyltransferases, lipases, serine proteases, heat shock proteins, esterases and others. Among the 1826 DEGs, we found that the transcriptional levels of CYP6AA9 in the laboratory populations was elevated as the levels of deltamethrin resistance increased. Moreover, the expression levels of the CYP6AA9 were significantly higher in the resistant strains than the susceptible strains in three different field populations. We further confirmed the association between the CYP6AA9 gene and deltamethrin resistance in mosquitoes by RNA interfering (RNAi). Altogether, we explored massive potential pyrethroid resistance-associated genes and demonstrated that CYP6AA9 participated in the pyrethroid resistance in mosquitoes.

Identification and characterization of CBL and CIPK gene families in eggplant (Solanum melongena L.)

Abstract: Eggplant (Solanum melongena L.) is an edible vegetable cultivated and consumed worldwide. But the production of eggplant is significantly limited by the soil salinization in greenhouse cultivation system. The main ions are Na(+), Ca(2+), Mg(2+), K(+), Cl(-), and SO4 (2-) in the salty soils. Calcineurin B-like proteins (CBLs) are calcium sensors and control the affinities and activities of numerous ion transporters with CBL-interacting protein kinases (CIPKs). In this study, a total of 5 CBL and 15 CIPK genes from eggplant were identified first. The yeast two-hybrid (Y2H) assay and bimolecular fluorescence complementation (BiFC) assay demonstrated the interaction network between SmCBLs and SmCIPKs. Strikingly, some new CBL-CIPK complexes were found which have never been discovered in any other plant species. The expression level of each SmCBL or SmCIPK under 200 mM NaCl, low potassium (LK; 100 μM), high Mg with 20 mM MgCl2 and MgSO4 stresses were examined by quantitative real-time PCR (qRT-PCR) assay and these CBL and CIPK genes were found to respond to the four ion stresses differently. Interestingly, the differential expression level of SmCIPK3, -24 or -25 to Mg(2+) is higher than Na(+), and Cl(-) is higher than SO4 (2-). In addition, different magnesium salt can induce different response mechanisms in eggplant. In summary, this study provides insight into the characterization of CBLs and CIPKs in eggplant. It may be used in a novel biotechnological breeding program strategy to create new eggplant cultivars, leading to enhance different ion tolerance.