Intron Phase Research Articles

Multiple calmodulin genes of the Pacific abalone, Haliotis discus hannai (Mollusca: Vetigastropoda: Haliotidae)

ABSTRACTIn this study, we identified four canonical calmodulin genes in the Pacific abalone, Haliotis discus hannai. Their full-length cDNAs were variable in the 5′ and 3′ untranslated regions, but highly similar (91–97%) in the coding region. Each of the genes encoded 149 amino acids, with 93–97% similarity among themselves and 94–98% similarity with human CAM I. There were 54 substitutions distributed unevenly throughout the coding regions, found mostly in the third codon position. Gene structure analysis revealed that each of the calmodulin genes comprised five exons and four introns. The intron positions and phases were identical and there were no introns in the fourth exon. The corresponding introns differed in their sequences and sizes. Expression profiles of nine tissues from abalone revealed that the calmodulin genes were transcribed in common in gill and mantle tissue, but differentially in the other tissues. A phylogenetic analysis based on the amino acid sequences revealed that calmodulin C was the most common isoform in Gastropoda and calmodulin was the most diverged isoform. An in silico analysis of the calmodulin genes identified paralogous genes in other Haliotis species, indicating that gene duplication might have occurred in the last common ancestor of Haliotis.Abbreviations: ORF: open reading frame; RACE: random amplification of cDNA end; TSA: transcriptome shotgun assembly; UTR: untranslated region

Structural and functional evolution of an auxin efflux carrier PIN1 and its functional characterization in common wheat.

Particularly PIN1, PIN protein-mediated rate-limiting auxin distribution plays a critical role in plant differentiation. Although well-characterized in Arabidopsis, little is known about the structural and functional relationship of the PIN1 gene among other plants. Here, we report that the gene structure remained conserved among bryophytes and angiosperms while the gene size varied by ~ 17%. Although the positions were conserved, highly variable intron phase suggests preference for specific regions in the gene sequence for independent events of intron insertion. Significant variation was observed across gene length for insertions and deletions that were mainly localized to the exonic regions flanking intron 1, possibly demarcating the sequences prone to deletions/duplications. The N and C-terminals showed a higher protein sequence similarity (~ 80%) compared to the central hydrophilic loop (~ 26%). In addition to the signature domains and motifs, we identified four novel uncharacterized motifs in the central divergent loop of PIN1 protein. Three different homo-loci, one each on chromosome groups 4, 6, and 7, were identified in wheat each showing dramatically different expression patterns during different plant developmental stages. Virus-induced gene silencing of the TaPIN1 gene resulted up to 26% reduction in plant height. Because of its direct role in controlling plant height along with a higher expression during stem elongation, the TaPIN1 gene can be manipulated to regulate plant height.

Genome-wide identification and expression analysis of the fatty acid desaturase genes in Medicago truncatula

Fatty acid desaturases (FADs) are of great importance and play critical roles in regulating plant fatty acid (FA) compositions. But to date, no reports about characterization of the FAD genes have been reported in the model dicotyledonous grass species Medicago truncatula. In this study, using database searches, 20 full-length FAD genes were identified in M. truncatula. These FAD genes were unevenly distributed on six chromosomes except the chromosome 6 and 8. Phylogenetic analysis showed the FAD genes in M. truncatula were clustered into six subfamilies and had similar exon number and intron phase in the same subfamily. Moreover, expression analysis based on qRT-PCR indicated these FAD genes were extensively involved in cold and heat responses. This study would provide an important foundation for future cloning and functional studies of FAD genes in M. truncatula and other related legume species.

Genome-Wide Identification and Characterization of the Potato bHLH Transcription Factor Family.

Plant basic/helix–loop–helix (bHLH) transcription factors participate in a number of biological processes, such as growth, development and abiotic stress responses. The bHLH family has been identified in many plants, and several bHLH transcription factors have been functionally characterized in Arabidopsis. However, no systematic identification of bHLH family members has been reported in potato (Solanum tuberosum). Here, 124 StbHLH genes were identified and named according to their chromosomal locations. The intron numbers varied from zero to seven. Most StbHLH proteins had the highly conserved intron phase 0, which accounted for 86.2% of the introns. According to the Neighbor-joining phylogenetic tree, 259 bHLH proteins acquired from Arabidopsis and potato were divided into 15 groups. All of the StbHLH genes were randomly distributed on 12 chromosomes, and 20 tandem duplicated genes and four pairs of duplicated gene segments were detected in the StbHLH family. The gene ontology (GO) analysis revealed that StbHLH mainly function in protein and DNA binding. Through the RNA-seq and quantitative real time PCR (qRT-PCR) analyses, StbHLH were found to be expressed in various tissues and to respond to abiotic stresses, including salt, drought and heat. StbHLH1, 41 and 60 were highly expressed in flower tissues, and were predicted to be involved in flower development by GO annotation. StbHLH45 was highly expressed in salt, drought and heat stress, which suggested its important role in abiotic stress response. The results provide comprehensive information for further analyses of the molecular functions of the StbHLH gene family.

Conservation and divergence of Starch Synthase III genes of monocots and dicots.

Starch Synthase (SS) plays an important role in extending the α-1,4 glucan chains during starch biosynthesis by catalyzing the transfer of the glucosyl moiety from ADP-glucose to the non-reducing end of a pre-existing glucan chain. SS has five distinct isoforms of which SSIII is involved in the formation of longer glucan chain length. Here we report identification and detailed characterization of ‘true’ orthologs of the well-characterized maize SSIII (ZmSSIII), among six monocots and two dicot species. ZmSSIII orthologs have nucleotide sequence similarity ranging from 56–81%. Variation in gene size among various orthologs ranged from 5.49 kb in Arabidopsis to 11.62 kb in Brachypodium and the variation was mainly due to intron size and indels present in the exons 1 and 3. Number of exons and introns were highly conserved among all orthologs however. While the intron number was conserved, intron phase showed variation at group, genera and species level except for intron 1 and 5. Several species, genera, and class specific cis-acting regulatory elements were identified in the promoter region. The predicted protein size of the SSIII orthologs ranged from 1094 amino acid (aa) in Arabidopsis to 1688 aa in Brachypodium with sequence identity ranging from 60%-89%. The N-terminal region of the protein was highly variable whereas the C-terminal region containing the Glycosyltransferase domain was conserved with >80% sequence similarity among the orthologs. In addition to confirming the known motifs, eleven novel motifs possibly providing species, genera and group specific functions, were identified in the three carbohydrate binding domains. Despite of significant sequence variation among orthologs, most of the motifs and their relative distances are highly conserved among the orthologs. The 3-D structure of catalytic region of SSIII orthologs superimposed with higher confidence confirming the presence of similar binding sites with five unidentified conserved regions in the catalytic (glycosyltransferase) domain including the pockets involved in catalysis and binding of ligands. Homeologs of wheat SSIII gene showed tissue and developmental stage specific expression pattern with the highest expression recorded in developing grains.

Слайдинг и вариабельность длины интронов для генов, обогащенных длинными интронами фазы 1

Due to high mutagenesis of intron sequences, intron evolution is usually considered in terms of evolution of exon-intron structures (EIS). The shifting of intron over short distances (rare evolutionary event called intron sliding) could lead to the change of intron phase, i.e. the intron position relative to the open reading frame. Here we analyze the EIS from four datasets of eukaryotic orthologues in order to find out the preferable choice of intron phase during sliding and to study the correlation between orthologous intron lengths. To identify the orthologous introns we have constructed the alignments of EIS of orthologous genes. Several sliding events with intron phase change were revealed from the analysis; however, our initial hypothesis that in the process of sliding introns prefer to change its phase to 0 more frequently, was not been confirmed. Nevertheless, it is necessary to expand the analysis on a larger dataset for making a proper conclusions. Despite high variability of intron length, some taxonomic groups share the similar length values. Moreover, some length conservation could be observed if instead of intron length L we consider a normalized length N = (L-A)/A, where A is an average length within an orthologous intron group. E.g. for ptprd genes of birds (28 species) the normalized value is in the interval (-0.15, 0.15) for 85.2 % of introns what is significantly higher than the values for random lengths set in accordance with the intron lengths distribution. That length “conservation” leads us to the question what intron length was in the ancient introns.

Genome-Wide Survey and Characterization of Fatty Acid Desaturase Gene Family in Brassica napus and Its Parental Species.

Rapeseed (Brassica napus) is an important oilseed crop worldwide, and fatty acid (FA) compositions determine the nutritional and economic value of its seed oil. Fatty acid desaturases (FADs) play a pivotal role in regulating FA compositions, but to date, no comprehensive genome-wide analysis of FAD gene family in rapeseed and its parent species has been reported. In this study, using homology searches, 84, 45, and 44 FAD genes were identified in rapeseed, Brassica rapa, and Brassica oleracea genomes, respectively. These FAD genes were unevenly located in 17 chromosomes and 2 scaffolds of rapeseed, 9 chromosomes and 1 scaffold of B. rapa, and all the chromosomes of B. oleracea. Phylogenetic analysis showed that the soluble and membrane-bound FADs in the three Brassica species were divided into four and six subfamilies, respectively. Generally, the soluble FADs contained two conserved histidine boxes, while three highly conserved histidine boxes were harbored in membrane-bound FADs. Exon-intron structure, intron phase, and motif composition and position were highly conserved in each FAD subfamily. Putative subcellular locations of FAD proteins in three Brassica species were consistent with those of corresponding known FADs. In total, 25 of simple sequence repeat (SSR) loci were found in FAD genes of the three Brassica species. Transcripts of selected FAD genes in the three species were examined in various organs/tissues or stress treatments from NCBI expressed sequence tag (EST) database. This study provides a critical molecular basis for quality improvement of rapeseed oil and facilitates our understanding of key roles of FAD genes in plant growth and development and stress response.

SnRK2 Homologs in Gossypium and GhSnRK2.6 Improved Salt Tolerance in Transgenic Upland Cotton and Arabidopsis

SnRK2s are a large family of plant-specific protein kinases, which play important roles in multiple abiotic stress responses in various plant species. But the family in Gossypium has not been well studied. Here, we identified 13, 10, and 13 members of the SnRK2 family from Gossypium raimondii, Gossypium arboreum, and Gossypium hirsutum, respectively, and analyzed the locations of SnRK2 homologs in chromosomes based on genome data of cotton species. Phylogenetic tree analysis of SnRK2 proteins showed that these families were classified into three groups. All SnRK2 genes were comprised of nine exons and eight introns, and the exon distributions and the intron phase of homolog genes among different cotton species were analogous. Moreover, GhSnRK2.6 was overexpressed in Arabidopsis and upland cotton, respectively. Under salt treatment, overexpressed Arabidopsis could maintain higher biomass accumulation than wild-type plants, and GhSnRK2.6 overexpression in cotton exhibited higher germination rate than the control. So, the gene GhSnRK2.6 could be utilized in cotton breeding for salt tolerance.

Mechanism for DNA transposons to generate introns on genomic scales.

Discovered four decades ago, the existence of introns was one of the most unexpected findings in molecular biology1. Introns are sequences interrupting genes that must be removed as part of mRNA production. Genome sequencing projects have documented that most eukaryotic genes contain at least one and frequently many introns2,3. Comparison of these genomes reveals a history of long evolutionary periods with little intron gain punctuated by episodes of rapid, extensive gain2,3. However, no detailed mechanism for such episodic intron generation has been empirically supported on a sufficient scale, despite several proposals4–8. Here we show how short non-autonomous DNA transposons independently generated hundreds to thousands of introns in the prasinophyte Micromonas pusilla and the pelagophyte Aureococcus anophagefferens. Each transposon carries one splice site. The other splice site is co-opted from gene sequence duplicated upon transposon insertion, allowing perfect splicing out of RNA. The distributions of sequences that can be co-opted are biased with respect to codons, and phasing of transposon-generated introns is similarly biased. These transposons insert between preexisting nucleosomes, so that multiple nearby insertions generate nucleosome-sized intervening segments. Thus, transposon insertion and sequence co-option may explain the intron phase biases2 and prevalence of nucleosome-sized exons9 observed in eukaryotes. Overall, the two independent examples of proliferating elements illustrate a general DNA transposon mechanism plausibly accounting for episodes of rapid, extensive intron gain during eukaryotic evolution2,3.

Genome-Wide Identification, Evolution, and Co-expression Network Analysis of Mitogen-Activated Protein Kinase Kinase Kinases in Brachypodium distachyon

Mitogen-activated protein kinase (MAPK) cascades are the conserved and universal signal transduction modules in all eukaryotes, which play the vital roles in plant growth, development, and in response to multiple stresses. In this study, we used bioinformatics methods to identify 86 MAPKKK protein encoded by 73 MAPKKK genes in Brachypodium. Phylogenetic analysis of MAPKKK family from Arabidopsis, rice, and Brachypodium has classified them into three subfamilies, of which 28 belonged to MEKK, 52 to Raf, and 6 to ZIK subfamily, respectively. Conserved protein motif, exon-intron organization, and splicing intron phase in kinase domains supported the evolutionary relationships inferred from the phylogenetic analysis. And gene duplication analysis suggested the chromosomal segment duplication happened before the divergence of the rice and Brachypodium, while all of three tandem duplicated gene pairs happened after their divergence. We further demonstrated that the MAPKKKs have evolved under strong purifying selection, implying the conservation of them. The splicing transcripts expression analysis showed that the splicesome translating longest protein tended to be adopted. Furthermore, the expression analysis of BdMAPKKKs in different organs and development stages as well as heat, virus and drought stresses revealed that the MAPKKK genes were involved in various signaling pathways. And the circadian analysis suggested there were 41 MAPKKK genes in Brachypodium showing cycled expression in at least one condition, of which seven MAPKKK genes expressed in all conditions and the promoter analysis indicated these genes possessed many cis-acting regulatory elements involved in circadian and light response. Finally, the co-expression network of MAPK, MAPKK, and MAPKKK in Brachypodium was constructed using 144 microarray and RNA-seq datasets, and ten potential MAPK cascades pathway were predicted. To conclude, our study provided the important information for evolutionary and functional characterization of MAPKKK family in Brachypodium, which will facilitate the functional analysis of BdMAPKKK genes, and also will facilitate better understanding the MAPK signal pathway in Brachypodium and beyond.

Comparative and evolutionary analysis of α-amylase gene across monocots and dicots

α-amylase is an important enzyme involved in starch degradation to provide energy to the germinating seedling. The present study was conducted to reveal structural and functional evolution of this gene among higher plants. Discounting polyploidy, most plant species showed only a single copy of the gene making multiple isoforms in different tissues and developmental stages. Genomic length of the gene ranged from 1472bp in wheat to 2369bp in soybean, and the size variation was mainly due to differences in the number and size of introns. In spite of this variation, the intron phase distribution and insertion sites were mostly conserved. The predicted protein size ranged from 414 amino acid (aa) in soybean to 449aa in Brachypodium. Overall, the protein sequence similarity among orthologs ranged from 56.4 to 97.4%. Key motifs and domains along with their relative distances were conserved among plants although several species, genera, and class specific motifs were identified. The glycosyl hydrolase superfamily domain length varied from 342aa in soybean to 384aa in maize and sorghum while length of the C-terminal β-sheet domain was highly conserved with 61aa in all monocots and Arabidopsis but was 59aa in soybean and Medicago. Compared to rice, 3D structure of the proteins showed 89.8 to 91.3% similarity among the monocots and 72.7 to 75.8% among the dicots. Sequence and relative location of the five key aa required for the ligand binding were highly conserved in all species except rice.

Dramatic Number Variation of R Genes in Solanaceae Species Accounted for by a Few R Gene Subfamilies.

Most disease resistance genes encode nucleotide-binding-site (NBS) and leucine-rich-repeat (LRR) domains, and the NBS-LRR encoding genes are often referred to as R genes. Using newly developed approach, 478, 485, 1,194, 1,665, 2,042 and 374 R genes were identified from the genomes of tomato Heinz1706, wild tomato LA716, potato DM1-3, pepper Zunla-1 and wild pepper Chiltepin and tobacco TN90, respectively. The majority of R genes from Solanaceae were grouped into 87 subfamilies, including 16 TIR-NBS-LRR (TNL) and 71 non-TNL subfamilies. Each subfamily was annotated manually, including identification of intron/exon structure and intron phase. Interestingly, TNL subfamilies have similar intron phase patterns, while the non-TNL subfamilies have diverse intron phase due to frequent gain of introns. Prevalent presence/absence polymorphic R gene loci were found among Solanaceae species, and an integrated map with 427 R loci was constructed. The pepper genome (2,042 in Chiltepin) has at least four times of R genes as in tomato (478 in Heinz1706). The high number of R genes in pepper genome is due to the amplification of R genes in a few subfamilies, such as the Rpi-blb2 and BS2 subfamilies. The mechanism underlying the variation of R gene number among different plant genomes is discussed.

Identification, duplication, evolution and expression analyses of caleosins in Brassica plants and Arabidopsis subspecies.

Caleosins are a class of Ca(2+) binding proteins that appear to be ubiquitous in plants. Some of the main proteins embedded in the lipid monolayer of lipid droplets, caleosins, play critical roles in the degradation of storage lipids during germination and in lipid trafficking. Some of them have been shown to have histidine-dependent peroxygenase activity, which is believed to participate in stress responses in Arabidopsis. In the model plant Arabidopsis thaliana, caleosins have been examined extensively. However, little is known on a genome-wide scale about these proteins in other members of the Brassicaceae. In this study, 51 caleosins in Brassica plants and Arabidopsis lyrata were investigated and analyzed in silico. Among them, 31 caleosins, including 7 in A. lyrata, 11 in Brassica oleracea and 13 in Brassica napus, are herein identified for the first time. Segmental duplication was the main form of gene expansion. Alignment, motif and phylogenetic analyses showed that Brassica caleosins belong to either the H-family or the L-family with different motif structures and physicochemical properties. Our findings strongly suggest that L-caleosins are evolved from H-caleosins. Predicted phosphorylation sites were differentially conserved in H-caleosin and L-caleosins, respectively. 'RY-repeat' elements and phytohormone-related cis-elements were identified in different caleosins, which suggest diverse physiological functions. Gene structure analysis indicated that most caleosins (38 out of 44) contained six exons and five introns and their intron phases were highly conserved. Structurally integrated caleosins, such as BrCLO3-3 and BrCLO4-2, showed high expression levels and may have important roles. Some caleosins, such as BrCLO2 and BoCLO8-2, lost motifs of the calcium binding domain, proline knot, potential phosphorylation sites and haem-binding sites. Combined with their low expression, it is suggested that these caleosins may have lost function.

Sequence and Expression Analysis of Interferon Regulatory Factor 10 (IRF10) in Three Diverse Teleost Fish Reveals Its Role in Antiviral Defense.

BackgroundInterferon regulatory factor (IRF) 10 was first found in birds and is present in the genome of other tetrapods (but not humans and mice), as well as in teleost fish. The functional role of IRF10 in vertebrate immunity is relatively unknown compared to IRF1-9. The target of this research was to clone and characterize the IRF10 genes in three economically important fish species that will facilitate future evaluation of this molecule in fish innate and adaptive immunity.Molecular Characterization of IRF10 in Three Fish SpeciesIn the present study, a single IRF10 gene was cloned in grass carp Ctenopharyngodon idella and Asian swamp eel Monopterus albus, and two, named IRF10a and IRF10b, in rainbow trout Oncorhynchus mykiss. The fish IRF10 molecules share highest identities to other vertebrate IRF10s, and have a well conserved DNA binding domain, IRF-associated domain, and an 8 exon/7 intron structure with conserved intron phase. The presence of an upstream ATG or open reading frame (ORF) in the 5’-untranslated region of different fish IRF10 cDNA sequences suggests potential regulation at the translational level, and this has been verified by in vitro transcription/translation experiments of the trout IRF10a cDNA, but would still need to be validated in fish cells.Expression Analysis of IRF10 In Vivo and In VitroBoth trout IRF10 paralogues are highly expressed in thymus, blood and spleen but are relatively low in head kidney and caudal kidney. Trout IRF10b expression is significantly higher than IRF10a in integumentary tissues i.e. gills, scales, skin, intestine, adipose fin and tail fins, suggesting that IRF10b may be more important in mucosal immunity. The expression of both trout IRF10 paralogues is up-regulated by recombinant IFN-γ. The expression of the IRF10 genes is highly induced by Poly I:C in vitro and in vivo, and by viral infection, but is less responsive to peptidoglycan and bacterial infection, suggesting an important role of fish IRF10 in antiviral defense.

Alternative splicing acting as a bridge in evolution.

Alternative splicing (AS) regulates diverse cellular and developmental functions through alternative protein structures of different isoforms. Alternative exons dominate AS in vertebrates; however, very little is known about the extent and function of AS in lower eukaryotes. To understand the role of introns in gene evolution, we examined AS from a green algal and five fungal genomes using a novel EST-based gene-modeling algorithm (COMBEST). AS from each genome was classified with COMBEST that maps EST sequences to genomes to build gene models. Various aspects of AS were analyzed through statistical methods. The interplay of intron 3n length, phase, coding property, and intron retention (RI) were examined with Chi-square testing. With 3 to 834 times EST coverage, we identified up to 73% of AS in intron-containing genes and found preponderance of RI among 11 types of AS. The number of exons, expression level, and maximum intron length correlated with number of AS per gene (NAG), and intron-rich genes suppressed AS. Genes with AS were more ancient, and AS was conserved among fungal genomes. Among stopless introns, non-retained introns (NRI) avoided, but major RI preferred 3n length. In contrast, stop-containing introns showed uniform distribution among 3n, 3n+1, and 3n+2 lengths. We found a clue to the intron phase enigma: it was the coding function of introns involved in AS that dictates the intron phase bias. Majority of AS is non-functional, and the extent of AS is suppressed for intron-rich genes. RI through 3n length, stop codon, and phase bias bridges the transition from functionless to functional alternative isoforms.

GenePlots: Large-Scale Gene Structures Visualization and Comparison Based on GenBank CDS Annotations

Drawing and comparing structures of homologous genes is an important way to study the gene evolution. In recent years, several programs and web tools have been proposed to draw gene structures according to genome sequences or GFF (General Feature Format) gene data. It is unfitted for them to draw a large number of gene structures due to time-consuming preparation of input data. The visualization tools for comparing the structures of hundreds or even thousands of genes are particularly lacking. We developed GenePlots, a comparative and interactive viewer for large scale gene structures based on CDS annotations in GenBank. With CDS annotations formatted like FASTA, GenePlots can draw all exon-intron structures for numerous genes in a unified scale. The most important assets of GenePlots are: (i) using GenBank CDS annotations as input data; (ii) the capability for analyzing thousands of genes simultaneously; (iii) the exhibition of intron phase; (iv) the flexibility allowing users to change the resulting diagram. GenePlots provides a simple and high efficient way to produce a large number of gene structures for large-scale genes comparison. Keywords: Exon and intron, intron phase, GenBank CDS annotation, GenePlots, gene structure visualization, software.

Recombination of chl-fus gene (Plastid Origin) downstream of hop: a locus of chromosomal instability.

BackgroundThe co-chaperone Hop [heat shock protein (HSP) organizing protein] has been shown to act as an adaptor for protein folding and maturation, in concert with Hsp70 and Hsp90. The hop gene is of eukaryotic origin. Likewise, the chloroplast elongation factor G (cEF-G) catalyzes the translocation step in chloroplast protein synthesis. The chl-fus gene, which encodes the cEF-G protein, is of plastid origin. Both proteins, Hop and cEF-G, derived from domain duplications. It was demonstrated that the nuclear chl-fus gene locates in opposite orientation to a hop gene in Glycine max. We explored 53 available plant genomes from Chlorophyta to higher plants, to determine whether the chl-fus gene was transferred directly downstream of the primordial hop in the proto-eukaryote host cell. Since both genes came from exon/module duplication events, we wanted to explore the involvement of introns in the early origin and the ensuing evolutionary changes in gene structure.ResultsWe reconstructed the evolutionary history of the two convergent plant genes, on the basis of their gene structure, microsynteny and microcolinearity, from 53 plant nuclear genomes. Despite a high degree (72 %) of microcolinearity among vascular plants, our results demonstrate that their adjacency was a product of chromosomal rearrangements. Based on predicted exon − intron structures, we inferred the molecular events giving rise to the current form of genes. Therefore, we propose a simple model of exon/module shuffling by intronic recombinations in which phase-0 introns were essential for domain duplication, and a phase-1 intron for transit peptide recruiting. Finally, we demonstrate a natural susceptibility of the intergenic region to recombine or delete, seriously threatening the integrity of the chl-fus gene for the future.ConclusionsOur results are consistent with the interpretation that the chl-fus gene was transferred from the chloroplast to a chromosome different from that of hop, in the primitive photosynthetic eukaryote, and much later before the appearance of angiosperms, it was recombined downstream of hop. Exon/module shuffling mediated by symmetric intron phases (i.e., phase-0 introns) was essential for gene evolution. The intergenic region is prone to recombine, risking the integrity of both genes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1780-1) contains supplementary material, which is available to authorized users.

Characterization and expression profiling of glutathione S-transferases in the diamondback moth, Plutella xylostella (L.).

BackgroundGlutathione S-transferases (GSTs) are multifunctional detoxification enzymes that play important roles in insects. The completion of several insect genome projects has enabled the identification and characterization of GST genes over recent years. This study presents a genome-wide investigation of the diamondback moth (DBM), Plutella xylostella, a species in which the GSTs are of special importance because this pest is highly resistant to many insecticides.ResultsA total of 22 putative cytosolic GSTs were identified from a published P. xylostella genome and grouped into 6 subclasses (with two unclassified). Delta, Epsilon and Omega GSTs were numerically superior with 5 genes for each of the subclasses. The resulting phylogenetic tree showed that the P. xylostella GSTs were all clustered into Lepidoptera-specific branches. Intron sites and phases as well as GSH binding sites were strongly conserved within each of the subclasses in the GSTs of P. xylostella. Transcriptome-, RNA-seq- and qRT-PCR-based analyses showed that the GST genes were developmental stage- and strain-specifically expressed. Most of the highly expressed genes in insecticide resistant strains were also predominantly expressed in the Malpighian tubules, midgut or epidermis.ConclusionsTo date, this is the most comprehensive study on genome-wide identification, characterization and expression profiling of the GST family in P. xylostella. The diversified features and expression patterns of the GSTs are inferred to be associated with the capacity of this species to develop resistance to a wide range of pesticides and biological toxins. Our findings provide a base for functional research on specific GST genes, a better understanding of the evolution of insecticide resistance, and strategies for more sustainable management of the pest.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1343-5) contains supplementary material, which is available to authorized users.

Identification and analysis of the TIFY gene family in Gossypium raimondii.

The highly conserved TIFY domain is included in the TIFY protein family of transcription factors, which is important in plant development. Here, 28 TIFY family genes were identified in the Gossypium raimondii genome and classified into JAZ (15 genes), ZML (8), PPD (3), and TIFY (2). The normal (TIF[F/Y]XG) motif was dominant in the TIFY family, excluding the ZML subfamily, in which TLSFXG was prevalent. TIFY family genes were unevenly distributed in the G. raimondii genome, with TIFY clusters present on chromosome 9. Phylogenetic analysis indicated abundant variations in the G. raimondii TIFY family, which were most closely related to those in Theobroma cacao among 5 species. Exon-intron organization and intron phases were homologous within each subfamily, correlating with their phylogeny. Intra-species synteny analyses indicated that genomic duplication contributed to the expansion of the TIFY family. Inter-species synteny analyses indicated that synteny regions involved in G. raimondii TIFY family genes were also present in the comparison of G. raimondii vs Arabidopsis thaliana or T. cacao, signifying that these genes had common ancestors and play the same or similar roles in biological processes. Greater synteny was present in the comparison of G. raimondii vs T. cacao than of G. raimondii vs A. thaliana. The expression patterns of TIFY family genes were characterized and most TIFY family genes were indicated to be involved in fiber development. Our study provides new data related to the evolution of TIFYs and their role as important regulators of transcription; these data can be useful for fiber development.

Genomewide identification and expression analysis of the ARF gene family in apple

Auxin response factors (ARF) are transcription factors that regulate auxin responses in plants. Although the genomewide analysis of this family has been performed in some species, little is known regarding ARF genes in apple (Malus domestica). In this study, 31 putative apple ARF genes have been identified and located within the apple genome. The phylogenetic analysis revealed that MdARFs could be divided into three subfamilies (groups I, II and III). The predicted MdARFs were distributed across 15 of 17 chromosomes with different densities. In addition, the analysis of exon-intron junctions and of the intron phase inside the predicted coding region of each candidate gene has revealed high levels of conservation within and between phylogenetic groups. Expression profile analyses of MdARF genes were performed in different tissues (root, stem, leaf, flower and fruit), and all the selected genes were expressed in at least one of the tissues that were tested, which indicated that MdARFs are involved in various aspects of physiological and developmental processes of apple. To our knowledge, this report is the first to provide a genomewide analysis of the apple ARF gene family. This study provides valuable information for understanding the classification and putative functions of the ARF signal in apple.