Statistical analysis of an RNA titration series evaluates microarray precision and sensitivity on a whole-array basis

Background

It is one of the ultimate goals for modern biological research to fully elucidate the intricate interplays and the regulations of the molecular determinants that propel and characterize the progression of versatile life phenomena, to name a few, cell cycling, developmental biology, aging, and the progressive and recurrent pathogenesis of complex diseases. The vast amount of large-scale and genome-wide time-resolved data is becoming increasing available, which provides the golden opportunity to unravel the challenging reverse-engineering problem of time-delayed gene regulatory networks.

Results

In particular, this methodological paper aims to reconstruct regulatory networks from temporal gene expression data by using delayed correlations between genes, i.e., pairwise overlaps of expression levels shifted in time relative each other. We have thus developed a novel model-free computational toolbox termed TdGRN (Time-delayed Gene Regulatory Network) to address the underlying regulations of genes that can span any unit(s) of time intervals. This bioinformatics toolbox has provided a unified approach to uncovering time trends of gene regulations through decision analysis of the newly designed time-delayed gene expression matrix. We have applied the proposed method to yeast cell cycling and human HeLa cell cycling and have discovered most of the underlying time-delayed regulations that are supported by multiple lines of experimental evidence and that are remarkably consistent with the current knowledge on phase characteristics for the cell cyclings.

Conclusion

We established a usable and powerful model-free approach to dissecting high-order dynamic trends of gene-gene interactions. We have carefully validated the proposed algorithm by applying it to two publicly available cell cycling datasets. In addition to uncovering the time trends of gene regulations for cell cycling, this unified approach can also be used to study the complex gene regulations related to the development, aging and progressive pathogenesis of a complex disease where potential dependences between different experiment units might occurs.     

GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists

Background

Although expression microarrays have become a standard tool used by biologists, analysis of data produced by microarray experiments may still present challenges. Comparison of data from different platforms, organisms, and labs may involve complicated data processing, and inferring relationships between genes remains difficult.

Results

S TAR N ET 2 is a new web-based tool that allows post hoc visual analysis of correlations that are derived from expression microarray data. S TAR N ET 2 facilitates user discovery of putative gene regulatory networks in a variety of species (human, rat, mouse, chicken, zebrafish, Drosophila, C. elegans, S. cerevisiae, Arabidopsis and rice) by graphing networks of genes that are closely co-expressed across a large heterogeneous set of preselected microarray experiments. For each of the represented organisms, raw microarray data were retrieved from NCBI's Gene Expression Omnibus for a selected Affymetrix platform. All pairwise Pearson correlation coefficients were computed for expression profiles measured on each platform, respectively. These precompiled results were stored in a MySQL database, and supplemented by additional data retrieved from NCBI. A web-based tool allows user-specified queries of the database, centered at a gene of interest. The result of a query includes graphs of correlation networks, graphs of known interactions involving genes and gene products that are present in the correlation networks, and initial statistical analyses. Two analyses may be performed in parallel to compare networks, which is facilitated by the new H EAT S EEKER module.

Conclusion

S TAR N ET 2 is a useful tool for developing new hypotheses about regulatory relationships between genes and gene products, and has coverage for 10 species. Interpretation of the correlation networks is supported with a database of previously documented interactions, a test for enrichment of Gene Ontology terms, and heat maps of correlation distances that may be used to compare two networks. The list of genes in a S TAR N ET network may be useful in developing a list of candidate genes to use for the inference of causal networks. The tool is freely available at http://vanburenlab.medicine.tamhsc.edu/starnet2.html, and does not require user registration.     

Mapping and identification of a Cicer arietinum NSP2 gene involved in nodulation pathway

Key message

For the first time the putative NSP2 gene in chickpea has been identified using pairs of NILs differing for the Rn1 / rn1 nodulation gene that was located in LG5 of chickpea genetic map.

Abstract

An intraspecific cross between the mutant non-nodulating genotype PM233, carrying the recessive gene rn1, and the wild-type CA2139 was used to develop two pairs of near-isogenic lines (NILs) for nodulation in chickpea. These pairs of NILs were characterized using sequence tagged microsatellite site (STMS) markers distributed across different linkage groups (LGs) of the chickpea genetic map leading to the detection of polymorphic markers located in LG5. Using this information, together with the genome annotation in Medicago truncatula, a candidate gene (NSP2) known to be involved in nodulation pathway was selected for mapping in chickpea. The full length sequence obtained in chickpea wild-type (CaNSP2) was 1,503 bp. Linkage analysis in an F₃ population of 118 plants derived from the cross between the pair of NILS NIL7-2A (nod) × NIL7-2B (non-nod) revealed a co-localization between CaNSP2 and Rn1 gene. These data implicate the CaNSP2 gene as a candidate for identity to Rn1, and suggest that it could act in the nodulation signaling transduction pathway similarly to that in other legumes species.     

Hierarchical Cluster-based Partial Least Squares Regression (HC-PLSR) is an efficient tool for metamodelling of nonlinear dynamic models

Background

Network inference methods reconstruct mathematical models of molecular or genetic networks directly from experimental data sets. We have previously reported a mathematical method which is exclusively data-driven, does not involve any heuristic decisions within the reconstruction process, and deliveres all possible alternative minimal networks in terms of simple place/transition Petri nets that are consistent with a given discrete time series data set.

Results

We fundamentally extended the previously published algorithm to consider catalysis and inhibition of the reactions that occur in the underlying network. The results of the reconstruction algorithm are encoded in the form of an extended Petri net involving control arcs. This allows the consideration of processes involving mass flow and/or regulatory interactions. As a non-trivial test case, the phosphate regulatory network of enterobacteria was reconstructed using in silico-generated time-series data sets on wild-type and in silico mutants.

Conclusions

The new exact algorithm reconstructs extended Petri nets from time series data sets by finding all alternative minimal networks that are consistent with the data. It suggested alternative molecular mechanisms for certain reactions in the network. The algorithm is useful to combine data from wild-type and mutant cells and may potentially integrate physiological, biochemical, pharmacological, and genetic data in the form of a single model.     

Fitting hidden Markov models of protein domains to a target species: application to Plasmodium falciparum

Background

The identification of gene sets that are significantly impacted in a given condition based on microarray data is a crucial step in current life science research. Most gene set analysis methods treat genes equally, regardless how specific they are to a given gene set.

Results

In this work we propose a new gene set analysis method that computes a gene set score as the mean of absolute values of weighted moderated gene t-scores. The gene weights are designed to emphasize the genes appearing in few gene sets, versus genes that appear in many gene sets. We demonstrate the usefulness of the method when analyzing gene sets that correspond to the KEGG pathways, and hence we called our method P athway A nalysis with D own-weighting of O verlapping G enes (PADOG). Unlike most gene set analysis methods which are validated through the analysis of 2-3 data sets followed by a human interpretation of the results, the validation employed here uses 24 different data sets and a completely objective assessment scheme that makes minimal assumptions and eliminates the need for possibly biased human assessments of the analysis results.

Conclusions

PADOG significantly improves gene set ranking and boosts sensitivity of analysis using information already available in the gene expression profiles and the collection of gene sets to be analyzed. The advantages of PADOG over other existing approaches are shown to be stable to changes in the database of gene sets to be analyzed. PADOG was implemented as an R package available at: http://bioinformaticsprb.med.wayne.edu/PADOG/or http://www.bioconductor.org.     

Evaluation of methods for detecting conversion events in gene clusters

Background

Several platforms for the analysis of genome-wide association data are available. However, these platforms focus on the evaluation of the genotype inherited by affected (i.e. case) individuals, whereas for some conditions (e.g. birth defects) the genotype of the mothers of affected individuals may also contribute to risk. For such conditions, it is critical to evaluate associations with both the maternal and the inherited (i.e. case) genotype. When genotype data are available for case-parent triads, a likelihood-based approach using log-linear modeling can be used to assess both the maternal and inherited genotypes. However, available software packages for log-linear analyses are not well suited to the analysis of typical genome-wide association data (e.g. including missing data).

Results

An integrated platform, Maternal and Inherited Analyses for Genome-wide Association Studies (MI-GWAS) for log-linear analyses of maternal and inherited genetic effects in large, genome-wide datasets, is described. MI-GWAS uses SAS and LEM software in combination to appropriately format data, perform the log-linear analyses and summarize the results. This platform was evaluated using existing genome-wide data and was shown to perform accurately and relatively efficiently.

Conclusions

The MI-GWAS platform provides a valuable tool for the analysis of association of a phenotype or condition with maternal and inherited genotypes using genome-wide data from case-parent triads. The source code for this platform is freely available at http://www.sph.uth.tmc.edu/sbrr/mi-gwas.htm.     

QTL mapping for salt tolerance and domestication-related traits in Vigna marina subsp. oblonga,a halophytic species

Key message

QTL mapping in F ₂ population [ V. luteola × V. marina subsp. oblonga ] revealed that the salt tolerance in V. marina subsp. oblonga is controlled by a single major QTL.

Abstract

The habitats of beach cowpea (Vigna marina) are sandy beaches in tropical and subtropical regions. As a species that grows closest to the sea, it has potential to be a gene source for breeding salt-tolerant crops. We reported here for the first time, quantitative trait loci (QTLs) mapping for salt tolerance in V. marina. A genetic linkage map was constructed from an F₂ population of 120 plants derived from an interspecific cross between V. luteola and V. marina subsp. oblonga. The map comprised 150 SSR markers. The markers were clustered into 11 linkage groups spanning 777.6 cM in length with a mean distance between the adjacent markers of 5.59 cM. The F_2:3 population was evaluated for salt tolerance under hydroponic conditions at the seedling and developmental stages. Segregation analysis indicated that salt tolerance in V. marina is controlled by a few genes. Multiple interval mapping consistently identified one major QTL which can explain about 50 % of phenotypic variance. The flanking markers may facilitate transfer of the salt tolerance allele from V. marina subsp. oblonga into related Vigna crops. The QTL for domestication-related traits from V. marina are also discussed.     

5′ Regulatory region of ubiquitin 2 gene from Porteresia coarctata makes efficient promoters for transgene expression in monocots and dicots

Key message

Porteresia ubiquitin 5′ regulatory region drives transgene expression in monocots and dicots.

Abstract

Ubiquitin promoters are promising candidates for constitutive transgene expression in plants. In this study, we isolated and characterized a novel 5′ regulatory sequence of a ubiquitin gene from Porteresia coarctata, a stress-tolerant wild grass species. Through functional analysis in heterologous plant systems, we have demonstrated that full length (Port Ubi2.3) or truncated sequence (PD2) of the isolated regulatory fragment can drive constitutive expression of GUS in monocots and/or dicots. In silico analysis of Port Ubi2.3 has revealed the presence of a 640 bp core promoter region followed by two exons and two introns with numerous putative cis-acting sites scattered throughout the regulatory region. Transformation and expression studies of six different deletion constructs in rice, tobacco and sugarcane revealed that the proximal intron has an enhancing effect on the activity of the core promoter in both monocots and dicots, whereas, Port Ubi2.3 was able to render strong expression only in monocots. This regulatory sequence is quite distinct from the other reported ubiquitin promoters in structure and performs better in monocots compared to other commonly used promoters—maize Ubi1 and Cauliflower Mosaic Virus 35S.     

Generation of selectable marker-free transgenic eggplant resistant to Alternaria solani using the R/RS site-specific recombination system

Key message

Marker-free transgenic eggplants, exhibiting enhanced resistance to Alternaria solani , can be generated on plant growth regulators (PGRs)- and antibiotic-free MS medium employing the multi-auto-transformation (MAT) vector, pMAT21 - wasabi defensin , wherein isopentenyl transferase ( ipt ) gene is used as a positive selection marker.

Abstract

Use of the selection marker genes conferring antibiotic or herbicide resistance in transgenic plants has been considered a serious problem for environment and the public. Multi-auto-transformation (MAT) vector system has been one of the tools to excise the selection marker gene and produce marker-free transgenic plants. Ipt gene was used as a selection marker gene. Wasabi defensin gene, isolated from Wasabia japonica (a Japanese horseradish which has been a potential source of antimicrobial proteins), was used as a gene of interest. Wasabi defensin gene was cloned from the binary vector, pEKH-WD, to an ipt-type MAT vector, pMAT21, by gateway cloning technology and transferred to Agrobacterium tumefaciens strain EHA105. Infected cotyledon explants of eggplant were cultured on PGRs- and antibiotic-free MS medium. Extreme shooty phenotype/ipt shoots were produced by the explants infected with the pMAT21-wasabi defensin (WD). The same PGRs- and antibiotic-free MS medium was used in subcultures of the ipt shoots. Subsequently, morphologically normal shoots emerged from the Ipt shoots. Molecular analyses of genomic DNA from transgenic plants confirmed the integration of the WD gene and excision of the selection marker (ipt gene). Expression of the WD gene was confirmed by RT-PCR and Northern blot analyses. In vitro whole plant and detached leaf assay of the marker-free transgenic plants exhibited enhanced resistance against Alternaria solani.     

Purple-leaved Ficus lyrata plants produced by overexpressing a grapevine VvMybA1 gene

Key message

This study established an efficient method of regenerating plants of Ficus lyrata and producing purple-leaved F. lyrata plants through genetic transformation using a VvMybA1 gene of grapevine.

Abstract

Ficus lyrata, a species with unique violin- or guitar-shaped leaves, was regenerated from leaf-derived calli cultured on Murashige and Skoog (MS) basal medium supplemented with 4.5 μM N-phenyl-N’-1, 2, 3-thiadiazol-5-yl urea (TDZ) and 0.5 μM α-naphthalene acetic acid (NAA). Leaf discs were inoculated with Agrobacterium tumefaciens strain EHA 105 harboring a binary vector DEAT that contains the VvMybA1 gene and neomycin phosphotransferase (npt II) gene and subsequently cultured on the established regeneration medium supplemented with 100 mg l^?1 kanamycin. Results showed that 87.5 % of the leaf discs produced kanamycin-resistant callus, and 68.8 % of them produced adventitious shoots. Transgenic plants with three leaf colors including green, green-purple, and purple were produced. Regular and quantitative real-time PCR analyses confirmed the integration of transgenes into the host genome. Semi-quantitative RT-PCR analysis indicated that the VvMybA1 gene was responsible for the purple-colored phenotype. Purple-leaved plants with strong color stability grew vigorously in a greenhouse. This study illustrated the feasibility of using a genetically engineered VvMybA1 gene for drastic modification of leaf color of an important woody ornamental plant.     

The Vaccinium corymbosum FLOWERING LOCUS T-like gene (VcFT): a flowering activator reverses photoperiodic and chilling requirements in blueberry

Key message

The blueberry FLOWERING LOCUS T ( FT )-like gene ( VcFT ) cloned from the cDNA of a tetraploid, northern highbush blueberry ( Vaccinium corymbosum L.) is able to reverse the photoperiodic and chilling requirements and drive early and continuous flowering.

Abstract

Blueberry is a woody perennial bush with a longer juvenile period than annual crops, requiring vernalization to flower normally. Few studies have been reported on the molecular mechanism of flowering in blueberry or other woody plants. Because FLOWERING LOCUS T (FT) from Arabidopsis thaliana plays a multifaceted role in generating mobile molecular signals to regulate plant flowering time, isolation and functional analysis of the blueberry (Vaccinium corymbosum L.) FT-like gene (VcFT) will facilitate the elucidation of molecular mechanisms of flowering in woody plants. Based on EST sequences, a 525-bpVcFT was identified and cloned from the cDNA of a tetraploid, northern highbush blueberry cultivar, Bluecrop. Ectopic expression of 35S:VcFT in tobacco induced flowering an average of 28 days earlier than wild-type plants. Expression of the 35S:VcFT in the blueberry cultivar Aurora resulted in an extremely early flowering phenotype, which flowered not only during in vitro culture, a growth stage when nontransgenic shoots had not yet flowered, but also in 6–10-week old, soil-grown transgenic plants, in contrast to the fact that at least 1 year and 800 chilling hours are required for the appearance of the first flower of both nontransgenic ‘Aurora’ and transgenic controls with the gusA. These results demonstrate that the VcFT is a functional floral activator and overexpression of the VcFT is able to reverse the photoperiodic and chilling requirements and drive early and continuous flowering.     

Estimating the absorptive root area in Norway spruce by using the common direct and indirect earth impedance methods

Aims

Estimates of root absorption magnitude are needed for the balanced management of forest ecosystems, but no methods able to work on the whole tree and stand level were available. Modified earth impedance method was developed recently and here it was tested, by comparing the results with those obtained by combination of several classical methods.

Methods

We used direct (soil cores, scanning and microscopy) and indirect (sap flow patterns and modified earth impedance) methods in an attempt to estimate the absorptive root area indexes (RAI) at two sites of about 25 and 40-years-old Norway spruce. We considered the geometric surfaces of all scanned fine roots to be equal to the fine root absorptive area (RAI _scan ). To estimate the potentially physically permeable area of fine roots, we microscopically evaluated the point of secondary xylem appearance and calculated the geometric area of root portions with primary structure (RAI _micro ). We termed the area of electrically conductive root surface as the active (ion) absorptive area (RAI _mei ) and measured its extent by the modified earth impedance (MEI) method.

Results

The highest values for absorptive root areas at the two experimental sites we obtained with the scanning method (RAI _scan was considered to be 100%), followed by the RAI _micro (51%) and RAI _mei (32%). RAImei reached about 2/3 of RAImicro. The surface area of the ectomycorrhizal hyphae was an order of magnitude larger than that of all fine roots, but the MEI did not measure such increase.

Conclusions

We showed that the absorptive root area, indirectly estimated by the MEI, provides consistent results that approach the values obtained for fine roots with a primary structure estimated by traditional direct methods. The similar range of the values for the absorptive root surface area obtained by microscopy and by the MEI method indicates that this method is feasible and that it could be used to determine the extent of active absorptive root surface areas in forests.     

Fine mapping and candidate gene analysis of the novel thermo-sensitive genic male sterility tms9-1 gene in rice

Key message

Fine mapping of the novel thermo-sensitive genic male sterility locus tms9 - 1 in the traditional TGMS line HengnongS-1 revealed that the MALE STERILITY1 homolog OsMS1 is the candidate gene.

Abstract

Photoperiod-thermo-sensitive genic male sterility (P/TGMS) has been widely used in the two-line hybrid rice breeding system. HengnongS-1 is one of the oldest TGMS lines and is often used in indica two-line breeding programs in China. In this study, our genetic analysis showed that the TGMS gene in HengnongS-1 was controlled by a single recessive gene that was non-allelic with the other TGMS loci identified, including C815S, Zhu1S and Y58S. Using SSR markers and bulked segregant analysis, we located the TGMS locus on chromosome 9 and named the gene tms9-1. Fine mapping further narrowed the tms9-1 loci to a 162 kb interval between two dCAPS markers. Sequence analysis revealed that a T to C substitution results in an amino acid change in the tms9-1 candidate gene (Os09g27620) in HengnongS-1 as compared to Minghui63. Sequencing of other rice accessions, including six P/TGMS lines, seven indica varieties and nine japonica varieties, showed that this SNP was exclusive to HengnongS-1. With multiple sequence alignment and expression pattern analyses, the rice MALE STERILITY1 homolog OsMS1 gene was identified as the candidate gene for tms9-1. Therefore, our study identified a novel TGMS locus and will facilitate the functional identification of the tms9-1 gene. Moreover, the markers linked to the tms9-1 gene will provide useful tools for the development of new TGMS lines by marker-assisted selection in two-line hybrid rice breeding programs.