Li Liu

ESPRIT: estimating species richness using large collections of 16S rRNA pyrosequences

Recent metagenomics studies of environmental samples suggested that microbial communities are much more diverse than previously reported, and deep sequencing will significantly increase the estimate of total species diversity. Massively parallel pyrosequencing technology enables ultra-deep sequencing of complex microbial populations rapidly and inexpensively. However, computational methods for analyzing large collections of 16S ribosomal sequences are limited. We proposed a new algorithm, referred to as ESPRIT, which addresses several computational issues with prior methods. We developed two versions of ESPRIT, one for personal computers (PCs) and one for computer clusters (CCs). The PC version is used for small- and medium-scale data sets and can process several tens of thousands of sequences within a few minutes, while the CC version is for large-scale problems and is able to analyze several hundreds of thousands of reads within one day. Large-scale experiments are presented that clearly demonstrate the effectiveness of the newly proposed algorithm. The source code and user guide are freely available at http://www.biotech.ufl.edu/people/sun/esprit.html.

Neuronal Transcriptome of Aplysia: Neuronal Compartments and Circuitry

Molecular analyses of Aplysia, a well-established model organism for cellular and systems neural science, have been seriously handicapped by a lack of adequate genomic information. By sequencing cDNA libraries from the central nervous system (CNS), we have identified over 175,000 expressed sequence tags (ESTs), of which 19,814 are unique neuronal gene products and represent 50%-70% of the total Aplysia neuronal transcriptome. We have characterized the transcriptome at three levels: (1) the central nervous system, (2) the elementary components of a simple behavior: the gill-withdrawal reflex-by analyzing sensory, motor, and serotonergic modulatory neurons, and (3) processes of individual neurons. In addition to increasing the amount of available gene sequences of Aplysia by two orders of magnitude, this collection represents the largest database available for any member of the Lophotrochozoa and therefore provides additional insights into evolutionary strategies used by this highly successful diversified lineage, one of the three proposed superclades of bilateral animals.

A global view of the selectivity of zinc deprivation and excess on genes expressed in human THP-1 mononuclear cells

Among the micronutrients required by humans, zinc has particularly divergent modes of action. cDNA microarray and quantitative PCR technologies were used to investigate the zinc responsiveness of known genes that influence zinc homeostasis and to identify, through global screening, genes that may relate to phenotypic outcomes of altered dietary zinc intake. Human monocytic/macrophage THP-1 cells were either acutely zinc depleted, using a cell-permeable zinc-specific chelator, or were supplemented with zinc to alter intracellular zinc concentrations. Initially, genes associated with zinc homeostasis were evaluated by quantitative PCR to establish ranges for fold changes in transcript abundance that might be expected with global screening. Zinc transporter-1 and zinc transporter-7 expression increased when cellular zinc increased, whereas Zip-2 expression, the most zinc-responsive gene examined, was markedly increased by zinc depletion. Microarrays composed of ≈22,000 elements were used to identify those genes responsive to either zinc depletion, zinc supplementation, or both conditions. Hierarchal clustering and ANOVA revealed that ≈5% or 1,045 genes were zinc responsive. Further sorting based on this pattern of the zinc responsiveness of these genes into seven groups revealed that 104 genes were linearly zinc responsive in a positive mode (i.e., increased expression as cellular zinc increases) and 86 genes that were linearly zinc responsive in a negative mode (i.e., decreased expression as cellular zinc increases). Expression of some genes was responsive to only zinc depletion or supplementation. Categorization by function revealed numerous genes needed for host defense were among those identified as zinc responsive, including cytokine receptors and genes associated with amplification of the Th1 immune response.

Improved breast cancer prognosis through the combination of clinical and genetic markers

MOTIVATION Accurate prognosis of breast cancer can spare a significant number of breast cancer patients from receiving unnecessary adjuvant systemic treatment and its related expensive medical costs. Recent studies have demonstrated the potential value of gene expression signatures in assessing the risk of post-surgical disease recurrence. However, these studies all attempt to develop genetic marker-based prognostic systems to replace the existing clinical criteria, while ignoring the rich information contained in established clinical markers. Given the complexity of breast cancer prognosis, a more practical strategy would be to utilize both clinical and genetic marker information that may be complementary. METHODS A computational study is performed on publicly available microarray data, which has spawned a 70-gene prognostic signature. The recently proposed I-RELIEF algorithm is used to identify a hybrid signature through the combination of both genetic and clinical markers. A rigorous experimental protocol is used to estimate the prognostic performance of the hybrid signature and other prognostic approaches. Survival data analyses is performed to compare different prognostic approaches. RESULTS The hybrid signature performs significantly better than other methods, including the 70-gene signature, clinical makers alone and the St. Gallen consensus criterion. At the 90% sensitivity level, the hybrid signature achieves 67% specificity, as compared to 47% for the 70-gene signature and 48% for the clinical makers. The odds ratio of the hybrid signature for developing distant metastases within five years between the patients with a good prognosis signature and the patients with a bad prognosis is 21.0 (95% CI:6.5-68.3), far higher than either genetic or clinical markers alone. AVAILABILITY The breast cancer dataset is available at www.nature.com and Matlab codes are available upon request.

Integrating omic technologies into aquatic ecological risk assessment and environmental monitoring: hurdles, achievements, and future outlook.

Background In this commentary we present the findings from an international consortium on fish toxicogenomics sponsored by the U.K. Natural Environment Research Council (Fish Toxicogenomics—Moving into Regulation and Monitoring, held 21–23 April 2008 at the Pacific Environmental Science Centre, Vancouver, BC, Canada). Objectives The consortium from government agencies, academia, and industry addressed three topics: progress in ecotoxicogenomics, regulatory perspectives on roadblocks for practical implementation of toxicogenomics into risk assessment, and dealing with variability in data sets. Discussion Participants noted that examples of successful application of omic technologies have been identified, but critical studies are needed to relate molecular changes to ecological adverse outcome. Participants made recommendations for the management of technical and biological variation. They also stressed the need for enhanced interdisciplinary training and communication as well as considerable investment into the generation and curation of appropriate reference omic data. Conclusions The participants concluded that, although there are hurdles to pass on the road to regulatory acceptance, omics technologies are already useful for elucidating modes of action of toxicants and can contribute to the risk assessment process as part of a weight-of-evidence approach.

Long-term retinal function and structure rescue using capsid mutant AAV8 vector in the rd10 mouse, a model of recessive retinitis pigmentosa.

The retinal degeneration 10 (rd10) mouse is a well-characterized model of autosomal recessive retinitis pigmentosa (RP), which carries a spontaneous mutation in the β subunit of rod cGMP-phosphodiesterase (PDEβ). Rd10 mouse exhibits photoreceptor dysfunction and rapid rod photoreceptor degeneration followed by cone degeneration and remodeling of the inner retina. Here, we evaluate whether gene replacement using the fast-acting tyrosine-capsid mutant AAV8 (Y733F) can provide long-term therapy in this model. AAV8 (Y733F)-smCBA-PDEβ was subretinally delivered to postnatal day 14 (P14) rd10 mice in one eye only. Six months after injection, spectral domain optical coherence tomography (SD-OCT), electroretinogram (ERG), optomotor behavior tests, and immunohistochemistry showed that AAV8 (Y733F)-mediated PDEβ expression restored retinal function and visual behavior and preserved retinal structure in treated rd10 eyes for at least 6 months. This is the first demonstration of long-term phenotypic rescue by gene therapy in an animal model of PDEβ-RP. It is also the first example of tyrosine-capsid mutant AAV8 (Y733F)-mediated correction of a retinal phenotype. These results lay the groundwork for the development of PDEβ-RP gene therapy trial and suggest that tyrosine-capsid mutant AAV vectors may be effective for treating other rapidly degenerating models of retinal degeneration.

Identification of Open Reading Frames Unique to a Select Agent: Ralstonia solanacearum Race 3 Biovar 2

An 8x draft genome was obtained and annotated for Ralstonia solanacearum race 3 biovar 2 (R3B2) strain UW551, a United States Department of Agriculture Select Agent isolated from geranium. The draft UW551 genome consisted of 80,169 reads resulting in 582 contigs containing 5,925,491 base pairs, with an average 64.5% GC content. Annotation revealed a predicted 4,454 protein coding open reading frames (ORFs), 43 tRNAs, and 5 rRNAs; 2,793 (or 62%) of the ORFs had a functional assignment. The UW551 genome was compared with the published genome of R. solanacearum race 1 biovar 3 tropical tomato strain GMI1000. The two phylogenetically distinct strains were at least 71% syntenic in gene organization. Most genes encoding known pathogenicity determinants, including predicted type III secreted effectors, appeared to be common to both strains. A total of 402 unique UW551 ORFs were identified, none of which had a best hit or >45% amino acid sequence identity with any R. solanacearum predicted protein; 16 had strong (E < 10(-13)) best hits to ORFs found in other bacterial plant pathogens. Many of the 402 unique genes were clustered, including 5 found in the hrp region and 38 contiguous, potential prophage genes. Conservation of some UW551 unique genes among R3B2 strains was examined by polymerase chain reaction among a group of 58 strains from different races and biovars, resulting in the identification of genes that may be potentially useful for diagnostic detection and identification of R3B2 strains. One 22-kb region that appears to be present in GMI1000 as a result of horizontal gene transfer is absent from UW551 and encodes enzymes that likely are essential for utilization of the three sugar alcohols that distinguish biovars 3 and 4 from biovars 1 and 2.

Reduced-Order Nonlinear Unsteady Aerodynamic Modeling Using a Surrogate-Based Recurrence Framework

A reduced-order nonlinear unsteady aerodynamic modeling approach suitable for analyzing pitching/plunging airfoils subject to fixed or time-varying freestream Mach numbers is described. The reduced-order model uses kriging surrogates to account for flow nonlinearities and recurrence solutions to account for time-history effects associated with unsteadiness. The resulting surrogate-based recurrence framework generates time-domain predictionsofunsteadylift,moment,anddragthataccuratelyapproximate computational fluiddynamicssolutions, but at a fraction of the computational cost. Results corresponding to transonic conditions demonstrate that the surrogate-based recurrence framework can mimic computational fluid dynamics predictions of unsteady aerodynamic responses when flow nonlinearities are present. For an unsteady aerodynamic modeling problem considered in this study, an accurate reduced-order model was generated by the surrogate-based recurrence framework approach with significantly fewer computational fluid dynamics evaluations compared to results reported in the literature for a similar problem in which a proper-orthogonal-decomposition-based approach was applied. Furthermore, the results show that the surrogate-based approach can accurately model time-varying freestream Mach number effects and is therefore applicable to rotary-wing applications in addition to fixed-wing applications.

AAV-mediated cone rescue in a naturally occurring mouse model of CNGA3-achromatopsia.

Achromatopsia is a rare autosomal recessive disorder which shows color blindness, severely impaired visual acuity, and extreme sensitivity to bright light. Mutations in the alpha subunits of the cone cyclic nucleotide-gated channels (CNGA3) are responsible for about 1/4 of achromatopsia in the U.S. and Europe. Here, we test whether gene replacement therapy using an AAV5 vector could restore cone-mediated function and arrest cone degeneration in the cpfl5 mouse, a naturally occurring mouse model of achromatopsia with a CNGA3 mutation. We show that gene therapy leads to significant rescue of cone-mediated ERGs, normal visual acuities and contrast sensitivities. Normal expression and outer segment localization of both M- and S-opsins were maintained in treated retinas. The therapeutic effect of treatment lasted for at least 5 months post-injection. This study is the first demonstration of substantial, relatively long-term restoration of cone-mediated light responsiveness and visual behavior in a naturally occurring mouse model of CNGA3 achromatopsia. The results provide the foundation for development of an AAV5-based gene therapy trial for human CNGA3 achromatopsia.

Higher-Harmonic-Control Algorithm for Helicopter Vibration Reduction Revisited

The higher-harmonic-control (HHC) algorithm is examined from a control theory perspective. A brief review of the history and variants of HHC is given, followed by a careful development of the algorithm. An analytic convergence and robustness analysis is then performed. Online identification with the adaptive variant of the algorithm is also addressed. A new version of the algorithm, relaxed HHC, is introduced and shown to have beneficial robustness properties. Some numerical results comparing these variants of the HHC algorithm applied to helicopter vibration reduction are also presented. The results presented unify and extend previous work on the higher-harmonic-control algorithm.