Analysis of the conservation of synteny between Fugu and human chromosome 12

Background

The pufferfish Fugu rubripes (Fugu) with its compact genome is increasingly recognized as an important vertebrate model for comparative genomic studies. In particular, large regions of conserved synteny between human and Fugu genomes indicate its utility to identify disease-causing genes. The human chromosome 12p12 is frequently deleted in various hematological malignancies and solid tumors, but the actual tumor suppressor gene remains unidentified.

Results

We investigated approximately 200 kb of the genomic region surrounding the ETV6 locus in Fugu (fETV6) in order to find conserved functional features, such as genes or regulatory regions, that could give insight into the nature of the genes targeted by deletions in human cancer cells. Seven genes were identified near the fETV6 locus. We found that the synteny with human chromosome 12 was conserved, but extensive genomic rearrangements occurred between the Fugu and human ETV6 loci.

Conclusion

This comparative analysis led to the identification of previously uncharacterized genes in the human genome and some potentially important regulatory sequences as well. This is a good indication that the analysis of the compact Fugu genome will be valuable to identify functional features that have been conserved throughout the evolution of vertebrates.

     

Functional characterization of transposon-tagged abiotic stress-responsive rice genes and their molecular polymorphisms among various stress-tolerant genotypes

Dissociation (Ds) insertional mutagenesis has been regarded as an efficient tool to generate insertion mutants for functional genomics and molecular breeding. However, little is known about the application of the tool on exploring biological functions of abiotic stress-related genes and their molecular breeding experience. In this study, a total of 833 Ds insertion lines have been obtained, which showed significantly higher tolerance or sensitivity to high salinity, drought or cold stress, by screening around 20,000 Ds lines. Analysis of Ds flanking sequence tags revealed that 165 rice genes were tagged by Ds insertion. Gene Ontology (GO) and gene set enrichment analysis showed that over-represented Ds-tagged genes might function in the response to exogenous stimuli. These Ds-tagged genes showed expression divergence among five high salinity and five drought tolerant rice lines under either high salinity or drought stress. Higher percentages of Ds-tagged genes were down- or up-regulated by these abiotic stresses. These Ds-tagged genes were also frequently reduced or suppressed by various phytohormones including abscisic acid and jasmonate. On the other hand, we have also detected single nucleotide polymorphisms (SNPs) and 1–10 base pairs of insertion and deletions (indels) of these Ds-tagged genes among ten high salinity/drought tolerant rice lines by comparing with the reference genome Nipponbare. Our data showed that SNPs were detected among 102 out of 165 genes and indels were identified in 39 out of 165 genes. All the data provided additional information to further explore the biological functions of these genes or to carry out molecular breeding.     

Predicting gene function using few positive examples and unlabeled ones

Background

A large amount of functional genomic data have provided enough knowledge in predicting gene function computationally, which uses known functional annotations and relationship between unknown genes and known ones to map unknown genes to GO functional terms. The prediction procedure is usually formulated as binary classification problem. Training binary classifier needs both positive examples and negative ones that have almost the same size. However, from various annotation database, we can only obtain few positive genes annotation for most offunctional terms, that is, there are only few positive examples for training classifier, which makes predicting directly gene function infeasible.

Results

We propose a novel approach SPE_RNE to train classifier for each functional term. Firstly, positive examples set is enlarged by creating synthetic positive examples. Secondly, representative negative examples are selected by training SVM(support vector machine) iteratively to move classification hyperplane to a appropriate place. Lastly, an optimal SVM classifier are trained by using grid search technique. On combined kernel ofYeast protein sequence, microarray expression, protein-protein interaction and GO functional annotation data, we compare SPE_RNE with other three typical methods in three classical performance measures recall R, precise P and their combination F: twoclass considers all unlabeled genes as negative examples, twoclassbal selects randomly same number negative examples from unlabeled gene, PSoL selects a negative examples set that are far from positive examples and far from each other.

Conclusions

In test data and unknown genes data, we compute average and variant of measure F. The experiments showthat our approach has better generalized performance and practical prediction capacity. In addition, our method can also be used for other organisms such as human.

     

Comparative Genomics and Synteny Analysis of <Emphasis Type="BoldItalic">KCS17</Emphasis>-<Emphasis Type="BoldItalic">KCS18</Emphasis> Cluster Across Different Genomes and Sub-genomes of Brassicaceae for Analysis of Its Evolutionary History

Comparative genomics-based synteny analysis has proved to be an effective strategy to understand evolution of genomic regions spanning a single gene (micro-unit) to large segments encompassing hundreds of kilobases to megabases. Brassicaceae is in a unique position to contribute to understanding genome and trait evolution through comparative genomics because whole genome sequences from as many as nine species have been completed and are available for analysis. In the present work, we compared genomic loci surrounding the KCS17-KCS18 cluster across these nine genomes. KCS18 or FAE 1 gene encodes beta-ketoacyl synthase, (β-KCS) a membrane-bound enzyme that catalyses the key rate-limiting step during synthesis of VLCFAs such as erucic acid (C22) present in seed oil in Brassicaceae by elongating carbon chain from C18 to C22; knowledge on role of KCS17 in plant development is however lacking. Synteny across the genomic segments harbouring FAE1 showed variable levels of gene retention ranging between 26% (Arabidopsis thaliana and Brassica napus C03) and 89% (between A. thaliana and Brassica rapa A01), and gene density ranged between 1 gene/2.86 kb and 1 gene/4.88 kb. Interestingly, in diploid Brassica species, FAE1 was retained in only one of the sub-genomes in spite of the presence of three sub-genomes created as a result of genome triplication; in contrast, FAE1 was present at three loci, with four copies in Camellina sativa which is also known to have experienced a recent genome triplication revealing contrasting fates upon duplication. The organization of KCS17 and KCS18 as head-to-tail cluster was conserved across most of the species, except the C genome containing Brassicas, namely B. oleracea and B. napus, where disruptions because of other genes were observed. Even in the conserved blocks, the distance between KCS17 and KCS18 varied; the functional implication of the organization of KCS17-KCS18 as a cluster vis-à-vis fatty acid biosynthesis needs to be dissected, as the cis-regulatory region is expected to be present in the intergenic space. Phylogenetic analysis of KCS gene family along with KCS17-KCS18 from members of Brassicaceae reveals their ancestral relationship with KCS8-KCS9 block. Further comparative functional analysis between KCS8, KCS9, KCS16, KCS17 and KCS18 across evolutionary time-scale will be required to understand the conservation or diversification of roles of these members of KCS family in fatty acid biosynthesis during course of evolution.     

Resolving neutral and deterministic contributions to genomic structure in <Emphasis Type="Italic">Syntrichia ruralis</Emphasis> (Bryophyta,Pottiaceae) informs propagule sourcing for dryland restoration

Syntrichia ruralis is a cosmopolitan moss that occupies steep environmental gradients. In arid to semi-arid regions of the world it is a key component of biological soil crusts, which are fundamental to healthy dryland ecosystem processes. As such, S. ruralis has attracted the attention of conservationists seeking to restore degraded biological soil crust communities and their associated vascular flora. Here, we generate genomic data for S. ruralis populations that span climatic gradients across the Colorado Plateau of the southwestern USA to investigate the contributions of neutral and deterministic processes to the partitioning of genomic structure. Although S. ruralis appears to be highly dispersible, geographic proximity significantly predicts genomic similarity. In addition, even when taking into account apparently high migration rates among populations and spatial autocorrelation of allele frequencies, some genomic variation is explained by environmental gradients correlated with elevation and latitude. Consequently, efforts to restore dryland ecosystems by establishing S. ruralis as a foundation should include strategies to ensure that propagule sources of this moss are environmentally stratified and targeted to the current/future climates of restoration sites.     

In silico identification and computational characterization of endogenous small interfering RNAs from diverse grapevine tissues and stages

Small interfering RNAs (siRNAs) are effectors of regulatory pathways underlying plant development, metabolism, and stress- and nutrient-signaling regulatory networks. The endogenous siRNAs are generally not conserved between plants; consequently, it is necessary and important to identify and characterize siRNAs from various plants. To address the nature and functions of siRNAs, and understand the biological roles of the huge siRNA population in grapevine (Vitis vinifera L.). The high-throughput sequencing technology was used to identify a large set of putative endogenous siRNAs from six grapevine tissues/organs. Subsequently, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed to classify the target genes of siRNA. In total, 520,519 candidate siRNAs were identified and their expression profiles exhibited typical temporal characters during grapevine development. In addition, we identified two grapevine trans-acting siRNA (TAS) gene homologs (VvTAS3 and VvTAS4) and the derived trans-acting siRNAs (tasiRNAs) that could target grapevine auxin response factor (ARF) and myeloblastosis (MYB) genes. Furthermore, the GO and KEGG analysis of target genes showed that most of them covered a broad range of functional categories, especially involving in disease-resistance process. The large-scale and completely genome-wide level identification and characterization of grapevine endogenous siRNAs from the diverse tissues by high throughput technology revealed the nature and functions of siRNAs in grapevine.     

Genes encoding members of 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) gene family from <Emphasis Type="Italic">Azadirachta indica</Emphasis> and correlation with azadirachtin biosynthesis

Azadirachta indica (A. Juss) commonly known as Neem is an important source of valuable natural products and occupies an important place in traditional healthcare system. Naturally, this plant synthesizes a number of tetranortriterpenoids utilizing isoprenoid as substrate flux. Although various phytochemical and pharmacological studies in A. indica have been carried out, but very limited information is available about the biosynthetic pathway as well as structural and regulatory genes involved in synthesis of bioactive molecules. In this study, we have cloned and characterized two genes, AiHMGR1 and AiHMGR2, encoding 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR; EC 1.1.1.34) catalyzing the rate limiting step of the isoprenoid biosynthesis. Two isoforms, AiHMGR1 and AiHMGR2, contain an open reading frame of 1707 and 1695 bp encoding polypeptides of 568 and 545 amino acid residues, respectively. The nucleotide and encoded amino acid sequence analyses suggest that both genes encode polypeptides with necessary structural domains present in other plant HMGRs, however, have different genomic organization. The relative expression analysis suggests that two genes express differentially in various tissues. Out of the two genes, expression of AiHMGR2 showed a direct correlation with azadirachtin accumulation in fruit tissue. The common as well as unique cis-regulatory elements present in both genes might be responsible for differential expression of both the genes in various tissues. The color complementation assay in Escherichia coli suggests that though both AiHMGR1 and AiHMGR2 encode functional proteins, AiHMGR2 is more active as compared to AiHMGR1.     

New resources for functional analysis of omics data for the genus <Emphasis Type="Italic">Aspergillus</Emphasis>

Background

Detailed and comprehensive genome annotation can be considered a prerequisite for effective analysis and interpretation of omics data. As such, Gene Ontology (GO) annotation has become a well accepted framework for functional annotation. The genus Aspergillus comprises fungal species that are important model organisms, plant and human pathogens as well as industrial workhorses. However, GO annotation based on both computational predictions and extended manual curation has so far only been available for one of its species, namely A. nidulans.

Results

Based on protein homology, we mapped 97% of the 3,498 GO annotated A. nidulans genes to at least one of seven other Aspergillus species: A. niger, A. fumigatus, A. flavus, A. clavatus, A. terreus, A. oryzae and Neosartorya fischeri. GO annotation files compatible with diverse publicly available tools have been generated and deposited online. To further improve their accessibility, we developed a web application for GO enrichment analysis named FetGOat and integrated GO annotations for all Aspergillus species with public genome sequences. Both the annotation files and the web application FetGOat are accessible via the Broad Institute's website (http://www.broadinstitute.org/fetgoat/index.html). To demonstrate the value of those new resources for functional analysis of omics data for the genus Aspergillus, we performed two case studies analyzing microarray data recently published for A. nidulans, A. niger and A. oryzae.

Conclusions

We mapped A. nidulans GO annotation to seven other Aspergilli. By depositing the newly mapped GO annotation online as well as integrating it into the web tool FetGOat, we provide new, valuable and easily accessible resources for omics data analysis and interpretation for the genus Aspergillus. Furthermore, we have given a general example of how a well annotated genome can help improving GO annotation of related species to subsequently facilitate the interpretation of omics data.

     

Genome-wide survey and expression analysis of the stress-associated protein gene family in desert poplar, <Emphasis Type="Italic">Populus euphratica</Emphasis>

Stress-associated proteins (SAPs) are a novel class of zinc finger proteins that extensively participate in abiotic stress responses. To date, no overall analysis and expression profiling of SAP genes in woody plants have been reported. Populus euphratica is distributed in desert regions and is extraordinarily adaptable to abiotic stresses. Thus, it is regarded as a promising candidate for studying abiotic stress resistance mechanisms of woody plants. In this study, 18 non-redundant SAP genes were identified from the genome of P. euphratica using basic local alignment search tool algorithms and functional domain verification. Among these 18 PeuSAP genes, 15 were intronless. To investigate the evolutionary relationships of SAP genes in P. euphratica and other Salicaceae plants, phylogenetic analyses were performed. Subsequently, the expression profiles of the 18 PeuSAP genes were analyzed in different tissues and under various stresses (drought, salt, heat, cold, and abscisic acid (ABA) treatment) using quantitative real-time PCR. Tissue expression analysis indicated that PeuSAPs showed no tissue specificity. PeuSAPs were induced by multiple abiotic stresses, especially drought, salt, and heat stresses, perhaps because of abundant cis-acting heat shock elements and drought-inducible elements in the promoter regions of the PeuSAPs. Moreover, single nucleotide polymorphisms (SNPs) variant analysis revealed many synonymous and non-synonymous SNPs in PeuSAP genes, but the zinc finger structure was conserved during evolution. These results provide an overview of the SAP gene family in P. euphratica and a reference for further functional research on PeuSAP genes.     

Genome sequencing and analysis of <Emphasis Type="Italic">Kloeckera apiculata</Emphasis> strain 34-9, a biocontrol agent against postharvest pathogens in citrus

Kloeckera apiculata, as the anamorphic state of Hanseniaspora uvarum from the Ascomycota phylum, plays an important role in the inhibition of fungal diseases in plants and spontaneous wine fermentation. This study was performed to sequence and analyze the whole genome of K. apiculata strain 34-9; This analysis provides further genomic features and assists functional research. The complete genome was determined using an Illumina HiSeq 2000 system applying paired-end and mate-pair methods to construct four reads libraries. The data assembly of all the reads resulted in a total genome size of 8.1 Mb, including 106 contigs, which were assembled into 41 scaffolds with a 31.95 % G+C content and a 234X sequence coverage. The performance of the gene prediction and functional annotation revealed that 2724 of 3786 protein-coding genes matched the KOG database, and 1127 genes were classified into GO categories. Further genome features analyses found 1066 microsatellite sites, 71 tRNAs, 3 rRNAs and 3 microRNAs in the genomic DNA. A prediction of the metabolic pathways identified potentially crucial genes for explaining the phenylalanine pathway involved in biocontrol. Comparisons with the typical yeasts Lachancea thermotolerans, Kluyveromyces lactis and Saccharomyces cerevisiae revealed the particularity and difference of K. apiculata strain 34-9. The genome alignments among Hanseniaspora vineae T02/19AF, K. apiculata DSM 2768 and 34-9 indicated numerous homologous regions distributed over the genomes between strain DSM2768 and 34-9. A SSR analysis identified that mono- and tri- nucleotide repeat types were more abundant in all six types, likely affecting the evolution of K. apiculata.     

Genome-wide bioinformatics analysis of Dof transcription factor gene family of chickpea and its comparative phylogenetic assessment with Arabidopsis and rice

The genome mining of chickpea (Cicer arietinum L.) revealed a total of 37 putative Dof genes using NCBI BLAST search against the genome with a highly conserved Dof domain. The translated Dof proteins possessed 150–493 amino acid residues with molecular weight ranging from 16.9 to 54.4 kD and pI varied from 4.98 to 9.64 as revealed by ExPASy server ProtParam. The exon–intron organization showed predominance of intronless Dof genes in chickpea. The predicted Dof genes were distributed among the eight chromosomes with a maximum of 9 Dof genes present on chromosome 7 and a single Dof gene was found on chromosome 8.The predominance of segmental gene duplication as compared to tandem duplication was observed which might be the prime cause of Dof gene family expansion in chickpea. The cis-regulatory element analysis revealed the presence of light-responsive, hormone-responsive, endosperm-specific, meristem-specific and stress-responsive elements. Comprehensive phylogenetic analyses of Dof genes of chickpea with Arabidopsis, rice, soybean and pigeonpea revealed several orthologs and paralogs assisting in understanding the putative functions of CaDof genes. The functional divergence and site-specific selective pressures of chickpea Dof genes have been investigated. The bioinformatics-based genome-wide assessment of Dof gene family of chickpea attempted in the present study could be a significant step for deciphering novel Dof genes based on genome-wide expression profiling.     

Genetic differentiation among natural populations of the lizard complex <Emphasis Type="Italic">Darevskia raddei</Emphasis> as inferred from genome microsatellite marking

The article presents the genetic parameters of the populations of lizards of the Darevskia raddei complex (D. raddei nairensis and D. raddei raddei) and the populations of D. valentini calculated on the basis of the analysis of variability of 50 allelic variants of the three nuclear genome microsatellite-containing loci of 83 individuals. It was demonstrated that the F_st genetic distances between the populations of D. raddei nairensis and D. raddei raddei were not statistically significantly different from the F_st genetic distances between the populations of different species, D. raddei and D. valentini. At the same time, these distances were statistically significantly higher than the F_st distances between the populations belonging to one species within the genus Darevskia. These data suggest deep divergence between the populations of D. raddei raddei and D. raddei nairensis of the D. raddei complex and there arises the question on considering them as separate species.     

Functional analysis of the two cyclophilin isoforms of <Emphasis Type="Italic">Sinorhizobium meliloti</Emphasis>

The nitrogen fixing Sinorhizobium meliloti possesses two genes, ppiA and ppiB, encoding two cyclophilin isoforms which belong to the superfamily of peptidyl prolyl cis/trans isomerases (PPIase, EC: 5.2.1.8). Here, we functionally characterize the two proteins and we demonstrate that both recombinant cyclophilins are able to isomerise the Suc-AAPF-pNA synthetic peptide but neither of them displays chaperone function in the citrate synthase thermal aggregation assay. Furthermore, we observe that the expression of both enzymes increases the viability of E. coli BL21 in the presence of abiotic stress conditions such as increased heat and salt concentration. Our results support and strengthen previous high-throughput studies implicating S. meliloti cyclophilins in various stress conditions.     

Gene loss/retention and evolutionary pattern of ascorbic acid biosynthesis and recycling genes in <Emphasis Type="Italic">Brassica rapa</Emphasis> following whole genome triplication

Ascorbic acid (AsA) is an inevitable antioxidant found abundantly in higher plants. Despite the importance of AsA in plants, how AsA biosynthesis (ABGs; d-mannose/l-galactose pathway) and AsA recycling genes (ARGs) evolved through polyploidization has not been addressed to date. Here, we evaluated the impacts of whole genome triplication (WGT) on ABGs and ARGs in Chinese cabbage (Brassica rapa ssp. pekinensis), which diverged from Arabidopsis thaliana before the WGT event. Twenty-three ABGs coded in 13 loci representing nine different enzyme classes and 29 ARGs coded in 19 loci representing five different enzyme classes were identified in the B. rapa genome by whole-genome screening through comparative genomic analyses. Five of these loci maintained three gene copies, 10 loci maintained two gene copies and the majority of the loci (n = 17) maintained single gene copies. Segmental (62 %) and tandem duplication (6 %), and fragment (21 %) and large-scale recombination (10 %) events accelerated the diversification of ABGs and ARGs. Thirteen of the 52 (25 %) identified genes experienced intron losses and two (4 %) experienced intron gains implying that intron losses outnumbered intron gains. The expansion and the retention of ABGs and ARGs were similar to the whole genome gene expansion and retention (P > 0.05). These findings provide new insights into the structural characteristics and evolutionary trends of ABGs and ARGs. In addition, our data could become a useful resource to further the functional characterization of these genes.     

Molecular mapping of powdery mildew resistance gene <Emphasis Type="Italic">PmSGD</Emphasis> in Chinese wheat landrace Shangeda using RNA-seq with bulk segregant analysis

Wheat powdery mildew, caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt), is one of the most devastating diseases of wheat in China and causes serious yield losses. Resistance genes are urgently needed by wheat breeding programs to combat this disease. In the present study, genetic analysis of powdery mildew resistance was conducted on segregated F₂ and F_2:3 populations derived from the cross of Shangeda (providing good resistance to powdery mildew) and Chancellor (susceptible to powdery mildew). The results showed that the resistance of Shangeda to E09 was controlled by a single recessive gene, tentatively designated as PmSGD. In addition, RNA sequencing of the parental lines Shangeda and Chancellor and the corresponding bulked pools derived from homozygous resistant or susceptible F_2:3 lines was implemented to identify single-nucleotide polymorphisms (SNPs). The PmSGD gene was estimated to be located in the 240–250-Mb region of chromosome 7B based on the characteristics of putative SNP loci distributed on 21 wheat chromosomes. Among the developed SNP markers, 17 (57%) markers were linked to PmSGD flanked by SNP2-57 and SNP2-46, with genetic distances of 0.4 and 0.8 cM, respectively. The reaction patterns of Shangeda and cultivars (lines) carrying the Pm5e, Pmhym, mlxbd, and PmTm4 genes to 22 Bgt isolates indicated that PmSGD may be allelic or very closely linked to those genes. All of the SNP loci linked to PmSGD were used to test 38 cultivars with known Pm gene(s), and the results suggested that these SNP loci are useful for pyramiding PmSGD by marker-assisted selection.