Deep phenotyping of post-infectious myalgic encephalomyelitis/chronic fatigue syndrome.

Post-infectious myalgic encephalomyelitis/chronic fatigue syndrome (PI-ME/CFS) is a disabling disorder, yet the clinical phenotype is poorly defined, the pathophysiology is unknown, and no disease-modifying treatments are available. We used rigorous criteria to recruit PI-ME/CFS participants with matched controls to conduct deep phenotyping. Among the many physical and cognitive complaints, one defining feature of PI-ME/CFS was an alteration of effort preference, rather than physical or central fatigue, due to dysfunction of integrative brain regions potentially associated with central catechol pathway dysregulation, with consequences on autonomic functioning and physical conditioning. Immune profiling suggested chronic antigenic stimulation with increase in naïve and decrease in switched memory B-cells. Alterations in gene expression profiles of peripheral blood mononuclear cells and metabolic pathways were consistent with cellular phenotypic studies and demonstrated differences according to sex. Together these clinical abnormalities and biomarker differences provide unique insight into the underlying pathophysiology of PI-ME/CFS, which may guide future intervention.

Adaptive immune responses to SARS-CoV-2 persist in the pharyngeal lymphoid tissue of children.

Most studies of adaptive immunity to SARS-CoV-2 infection focus on peripheral blood, which may not fully reflect immune responses at the site of infection. Using samples from 110 children undergoing tonsillectomy and adenoidectomy during the COVID-19 pandemic, we identified 24 samples with evidence of previous SARS-CoV-2 infection, including neutralizing antibodies in serum and SARS-CoV-2-specific germinal center and memory B cells in the tonsils and adenoids. Single-cell B cell receptor (BCR) sequencing indicated virus-specific BCRs were class-switched and somatically hypermutated, with overlapping clones in the two tissues. Expanded T cell clonotypes were found in tonsils, adenoids and blood post-COVID-19, some with CDR3 sequences identical to previously reported SARS-CoV-2-reactive T cell receptors (TCRs). Pharyngeal tissues from COVID-19-convalescent children showed persistent expansion of germinal center and antiviral lymphocyte populations associated with interferon (IFN)-γ-type responses, particularly in the adenoids, and viral RNA in both tissues. Our results provide evidence for persistent tissue-specific immunity to SARS-CoV-2 in the upper respiratory tract of children after infection.

Robust, persistent adaptive immune responses to SARS-CoV-2 in the oropharyngeal lymphoid tissue of children.

SARS-CoV-2 infection triggers adaptive immune responses from both T and B cells. However, most studies focus on peripheral blood, which may not fully reflect immune responses in lymphoid tissues at the site of infection. To evaluate both local and systemic adaptive immune responses to SARS-CoV-2, we collected peripheral blood, tonsils, and adenoids from 110 children undergoing tonsillectomy/adenoidectomy during the COVID-19 pandemic and found 24 with evidence of prior SARS-CoV-2 infection, including detectable neutralizing antibodies against multiple viral variants. We identified SARS-CoV-2-specific germinal center (GC) and memory B cells; single cell BCR sequencing showed that these virus-specific B cells were class-switched and somatically hypermutated, with overlapping clones in the adenoids and tonsils. Oropharyngeal tissues from COVID-19-convalescent children showed persistent expansion of GC and anti-viral lymphocyte populations associated with an IFN-γ-type response, with particularly prominent changes in the adenoids, as well as evidence of persistent viral RNA in both tonsil and adenoid tissues of many participants. Our results show robust, tissue-specific adaptive immune responses to SARS-CoV-2 in the upper respiratory tract of children weeks to months after acute infection, providing evidence of persistent localized immunity to this respiratory virus.

Ethnic Pharmacogenomic Differences in the Management of Asian Patients with Metastatic Prostate Cancer.

Progression to metastatic disease occurs in about half of all men who develop prostate cancer (PC), one of the most common cancers in men worldwide. Androgen deprivation therapy has been the mainstay therapy for patients with metastatic PC (mPC) since the 1940s. In the last decade, there has been unprecedented advancement in systemic therapies, e.g., taxane, androgen-signalling pathway inhibitors, and biomarker-driven targeted therapies for various stages of disease, resulting in overall survival improvement. Adding to ongoing controversies over how best to treat these patients is the recognition that ethnicity may influence prognosis and outcomes. This review discusses recent evidence for the impacts of Asian ethnicity specifically, which includes environmental, sociocultural, and genetic factors, on the approach to pharmacological management of mPC. Clear inter-ethnic differences in drug tolerability, serious adverse events (AEs), and genetic heterogeneity must all be considered when dosing and scheduling for treatment, as well as designing future precision studies in PC.

Colchicine's effects on metabolic and inflammatory molecules in adults with obesity and metabolic syndrome: results from a pilot randomized controlled trial.

OBJECTIVE: Recent clinical trials have demonstrated that colchicine may have metabolic and cardiovascular and benefits in at-risk patients; however, the mechanisms through which colchicine may improve outcomes are still unclear. We sought to examine colchicine's effects on circulating inflammatory and metabolic molecules in adults with obesity and metabolic syndrome (MetS). METHODS: Blood samples were collected pre- and post-intervention during a double-blind randomized controlled trial in which 40 adults with obesity and MetS were randomized to colchicine 0.6 mg or placebo twice-daily for 3 months. Serum samples were analyzed for 1305 circulating factors using the SomaScan Platform. The Benjamini-Hochberg procedure was used to adjust the false discovery rate (FDR) for multiple testing. RESULTS: At baseline, age (48.0 ± 13.8 vs. 44.7 ± 10.3 years) and BMI (39.8 ± 6.4 vs. 41.8 ± 8.2 kg/m2) were not different between groups. After controlling for the FDR, 34 molecules were significantly changed by colchicine. Colchicine decreased concentrations of multiple inflammatory molecules, including C-reactive protein, interleukin 6, and resistin, in addition to vascular-related proteins (e.g., oxidized low-density lipoprotein receptor, phosphodiesterase 5A). Conversely, relative to placebo, colchicine significantly increased concentrations of eight molecules including secreted factors associated with metabolism and anti-thrombosis. CONCLUSIONS: In adults with obesity, colchicine significantly affected concentrations of proteins involved in the innate immune system, endothelial function and atherosclerosis, uncovering new mechanisms behind its cardiometabolic effects. Further research is warranted to investigate whether colchicine's IL-6 suppressive effects may be beneficial in COVID-19.

Broad immune activation underlies shared set point signatures for vaccine responsiveness in healthy individuals and disease activity in patients with lupus.

Responses to vaccination and to diseases vary widely across individuals, which may be partly due to baseline immune variations. Identifying such baseline predictors of immune responses and their biological basis is of broad interest, given their potential importance for cancer immunotherapy, disease outcomes, vaccination and infection responses. Here we uncover baseline blood transcriptional signatures predictive of antibody responses to both influenza and yellow fever vaccinations in healthy subjects. These same signatures evaluated at clinical quiescence are correlated with disease activity in patients with systemic lupus erythematosus with plasmablast-associated flares. CITE-seq profiling of 82 surface proteins and transcriptomes of 53,201 single cells from healthy high and low influenza vaccination responders revealed that our signatures reflect the extent of activation in a plasmacytoid dendritic cell-type I IFN-T/B lymphocyte network. Our findings raise the prospect that modulating such immune baseline states may improve vaccine responsiveness and mitigate undesirable autoimmune disease activity.

Highly multiplexed proteomic assessment of human bone marrow in acute myeloid leukemia.

Acute myeloid leukemia (AML) is a genetically heterogeneous disease that is characterized by abnormal clonal proliferation of myeloid progenitor cells found predominantly within the bone marrow (BM) and blood. Recent studies suggest that genetic and phenotypic alterations in the BM microenvironment support leukemogenesis and allow leukemic cells to survive and evade chemotherapy-induced death. However, despite substantial evidence indicating the role of tumor-host interactions in AML pathogenesis, little is known about the complex microenvironment of the BM. To address this, we performed novel proteomic profiling of the noncellular compartment of the BM microenvironment in patients with AML (n = 10) and age- and sex-matched healthy control subjects (n = 10) using an aptamer-based, highly multiplexed, affinity proteomics platform (SOMAscan). We show that proteomic assessment of blood or RNA-sequencing of BM are suboptimal alternate screening strategies to determine the true proteomic composition of the extracellular soluble compartment of AML patient BM. Proteomic analysis revealed that 168 proteins significantly differed in abundance, with 91 upregulated and 77 downregulated in leukemic BM. A highly connected signaling network of cytokines and chemokines, including IL-8, was found to be the most prominent proteomic signature associated with AML in the BM microenvironment. We report the first description of significantly elevated levels of the myelosuppressive chemokine CCL23 (myeloid progenitor inhibitory factor-1) in both AML and myelodysplastic syndrome patients and perform functional experiments supportive of a role in the suppression of normal hematopoiesis. This unique paired RNA-sequencing and proteomics data set provides innovative mechanistic insights into AML and healthy aging and should serve as a useful public resource.

Circulating S100A8 and S100A9 protein levels in plasma of patients with acquired aplastic anemia and myelodysplastic syndromes.

The alarmin family members S100A8 and S100A9 are acute phase inflammation proteins, but they also have been proposed as biomarkers in many malignant and non-malignant diseases. In this study, circulating S100A8 and S100A9 homodimers and S100A8/A9 heterodimers in plasma were systematically investigated by ELISA in aplastic anemia (AA) and myelodysplastic syndromes (MDS). Plasma was obtained from 58 severe AA (SAA) and 30 MDS patients, and from 47 age- and sex-matched healthy donors. In 40 out of the 58 AA subjects, S100A protein levels were measured before and 6 months after immunosuppressive therapy (IST). No differences were observed in AA patients at diagnosis compared to healthy controls for circulating S100A homodimers and heterodimers. After therapy, SAA-responders showed significantly increased circulating S100A8. Non-responding patients had significantly higher levels of circulating S100A8/A9 compared to responders and healthy controls, but without variations of S100A8 and S100A9 homodimers. In MDS patients, circulating S100A8 was significantly elevated compared to those of AA and/or healthy controls. By Pearson correlation analysis of protein levels and blood counts, multiple correlations were found. However, as S100A8 and S100A9 are abundantly present in white blood cells and platelets, correlations with blood counts likely mirror the higher number of cells in the blood of some patients. In conclusion, our findings indicate that circulating S100A8 is increased in MDS but not in AA, and that may be useful to distinguish these diseases in the differential diagnosis of bone marrow failure syndromes. Clinicaltrials.gov identifiers: NCT00260689, NCT00604201, NCT01328587, NCT01623167, NCT00001620, NCT00001397.

Impaired B cell immunity in acute myeloid leukemia patients after chemotherapy.

BACKGROUND: Changes in adaptive immune cells after chemotherapy in adult acute myeloid leukemia (AML) may have implications for the success of immunotherapy. This study was designed to determine the functional capacity of the immune system in adult patients with AML who have completed chemotherapy and are potential candidates for immunotherapy. METHODS: We used the response to seasonal influenza vaccination as a surrogate for the robustness of the immune system in 10 AML patients in a complete remission post-chemotherapy and performed genetic, phenotypic, and functional characterization of adaptive immune cell subsets. RESULTS: Only 2 patients generated protective titers in response to vaccination, and a majority of patients had abnormal frequencies of transitional and memory B-cells. B-cell receptor sequencing showed a B-cell repertoire with little evidence of somatic hypermutation in most patients. Conversely, frequencies of T-cell populations were similar to those seen in healthy controls, and cytotoxic T-cells demonstrated antigen-specific activity after vaccination. Effector T-cells had increased PD-1 expression in AML patients least removed from chemotherapy. CONCLUSION: Our results suggest that while some aspects of cellular immunity recover quickly, humoral immunity is incompletely reconstituted in the year following intensive cytotoxic chemotherapy for AML. The observed B-cell abnormalities may explain the poor response to vaccination often seen in AML patients after chemotherapy. Furthermore, the uncoupled recovery of B-cell and T-cell immunity and increased PD-1 expression shortly after chemotherapy might have implications for the success of several modalities of immunotherapy.

Effects of Systemically Administered Hydrocortisone on the Human Immunome.

Corticosteroids have been used for decades to modulate inflammation therapeutically, yet there is a paucity of data on their effects in humans. We examined the changes in cellular and molecular immune system parameters, or "immunome", in healthy humans after systemic corticosteroid administration. We used multiplexed techniques to query the immunome in 20 volunteers at baseline, and after intravenous hydrocortisone (HC) administered at moderate (250 mg) and low (50 mg) doses, to provide insight into how corticosteroids exert their effects. We performed comprehensive phenotyping of 120 lymphocyte subsets by high dimensional flow cytometry, and observed a decline in circulating specific B and T cell subsets, which reached their nadir 4-8 hours after administration of HC. However, B and T cells rebounded above baseline 24 hours after HC infusion, while NK cell numbers remained stable. Whole transcriptome profiling revealed down regulation of NF-κB signaling, apoptosis, and cell death signaling transcripts that preceded lymphocyte population changes, with activation of NK cell and glucocorticoid receptor signaling transcripts. Our study is the first to systematically characterize the effects of corticosteroids on the human immunome, and we demonstrate that HC exerts differential effects on B and T lymphocytes and natural killer cells in humans.

Ultra-low dose interleukin-2 promotes immune-modulating function of regulatory T cells and natural killer cells in healthy volunteers.

Low-dose interleukin-2 (IL-2) expands regulatory T cells (Tregs) and natural killer (NK) cells after stem cell transplantation (SCT) and may reduce graft-versus-host disease (GVHD). We hypothesized that ultra-low dose (ULD) IL-2 could serve as an immune-modulating agent for stem cell donors to prevent GVHD following SCT. However, the safety, dose level, and immune signatures of ULD IL-2 in immune-competent healthy subjects remain unknown. Here, we have characterized the phenotype and function of Tregs and NK cells as well as the gene expression and cytokine profiles of 21 healthy volunteers receiving 50,000 to 200,000 units/m(2)/day IL-2 for 5 days. ULD IL-2 was well tolerated and induced a significant increase in the frequency of Tregs with increased suppressive function. There was a marked expansion of CD56(bright) NK cells with enhanced interferon-γ (IFN-γ) production. Serum cytokine profiling demonstrated increase of IFN-γ induced protein 10 (IP-10). Gene expression analysis revealed significant changes in a highly restricted set of genes, including FOXP3, IL-2RA, and CISH. This is the first study to evaluate global immune-modulating function of ULD IL-2 in healthy subjects and to support the future studies administrating ULD IL-2 to stem cell donors.

In silico and in vitro characterization of phospholipase A2 isoforms from soybean (Glycine max).

At the present, no secreted phospholipase A2 (sPLA2) from soybean (Glycine max) was investigated in detail. In this work we identified five sequences of putative secreted sPLA2 from soybean after a BLAST search in G. max database. Sequence analysis showed a conserved PA2c domain bearing the Ca²âº binding loop and the active site motif. All the five mature proteins contain 12 cysteine residues, which are commonly conserved in plant sPLA2s. We propose a phylogenetic tree based on sequence alignment of reported plant sPLA2s including the novel enzymes from G. max. According to PLA2 superfamily, two of G. max sPLA2s are grouped as XIA and the rest of sequences as XIB, on the basis of differences found in their molecular weights and deviating sequences especially in the N- and C-terminal regions of the isoenzymes. Furthermore, we report the cloning, expression and purification of one of the putative isoenzyme denoted as GmsPLA2-XIA-1. We demonstrate that this mature sPLA2 of 114 residues had PLA2 activity on Triton:phospholipid mixed micelles and determine the kinetic parameters for this system. We generate a model based on the known crystal structure of sPLA2 from rice (isoform II), giving first insights into the three-dimensional structure of folded GmsPLA2-XIA-1. Besides describing the spatial arrangement of highly conserved pair HIS-49/ASP-50 and the Ca⁺² loop domains, we propose the putative amino acids involved in the interfacial recognition surface. Additionally, molecular dynamics simulations indicate that calcium ion, besides its key function in the catalytic cycle, plays an important role in the overall stability of GmsPLA2-XIA-1 structure.

A newly-developed community microarray resource for transcriptome profiling in Brassica species enables the confirmation of Brassica-specific expressed sequences.

BACKGROUND: The Brassica species include an important group of crops and provide opportunities for studying the evolutionary consequences of polyploidy. They are related to Arabidopsis thaliana, for which the first complete plant genome sequence was obtained and their genomes show extensive, although imperfect, conserved synteny with that of A. thaliana. A large number of EST sequences, derived from a range of different Brassica species, are available in the public database, but no public microarray resource has so far been developed for these species. RESULTS: We assembled unigenes using approximately 800,000 EST sequences, mainly from three species: B. napus, B. rapa and B. oleracea. The assembly was conducted with the aim of co-assembling ESTs of orthologous genes (including homoeologous pairs of genes in B. napus from each of the A and C genomes), but resolving assemblies of paralogous, or paleo-homoeologous, genes (i.e. the genes related by the ancestral genome triplication observed in diploid Brassica species). 90,864 unique sequence assemblies were developed. These were incorporated into the BAC sequence annotation for the Brassica rapa Genome Sequencing Project, enabling the identification of cognate genomic sequences for a proportion of them. A 60-mer oligo microarray comprising 94,558 probes was developed using the unigene sequences. Gene expression was analysed in reciprocal resynthesised B. napus lines and the B. oleracea and B. rapa lines used to produce them. The analysis showed that significant expression could consistently be detected in leaf tissue for 35,386 unigenes. Expression was detected across all four genotypes for 27,355 unigenes, genome-specific expression patterns were observed for 7,851 unigenes and 180 unigenes displayed other classes of expression pattern. Principal component analysis (PCA) clearly resolved the individual microarray datasets for B. rapa, B. oleracea and resynthesised B. napus. Quantitative differences in expression were observed between the resynthesised B. napus lines for 98 unigenes, most of which could be classified into non-additive expression patterns, including 17 that showed cytoplasm-specific patterns. We further characterized the unigenes for which A genome-specific expression was observed and cognate genomic sequences could be identified. Ten of these unigenes were found to be Brassica-specific sequences, including two that originate from complex loci comprising gene clusters. CONCLUSION: We succeeded in developing a Brassica community microarray resource. Although expression can be measured for the majority of unigenes across species, there were numerous probes that reported in a genome-specific manner. We anticipate that some proportion of these will represent species-specific transcripts and the remainder will be the consequence of variation of sequences within the regions represented by the array probes. Our studies demonstrated that the datasets obtained from the arrays can be used for typical analyses, including PCA and the analysis of differential expression. We have also demonstrated that Brassica-specific transcripts identified in silico in the sequence assembly of public EST database accessions are indeed reported by the array. These would not be detectable using arrays designed using A. thaliana sequences.

Analysis of the Pythium ultimum transcriptome using Sanger and Pyrosequencing approaches.

BACKGROUND: Pythium species are an agriculturally important genus of plant pathogens, yet are not understood well at the molecular, genetic, or genomic level. They are closely related to other oomycete plant pathogens such as Phytophthora species and are ubiquitous in their geographic distribution and host rage. To gain a better understanding of its gene complement, we generated Expressed Sequence Tags (ESTs) from the transcriptome of Pythium ultimum DAOM BR144 (= ATCC 200006 = CBS 805.95) using two high throughput sequencing methods, Sanger-based chain termination sequencing and pyrosequencing-based sequencing-by-synthesis. RESULTS: A single half-plate pyrosequencing (454 FLX) run on adapter-ligated cDNA from a normalized cDNA population generated 90,664 reads with an average read length of 190 nucleotides following cleaning and removal of sequences shorter than 100 base pairs. After clustering and assembly, a total of 35,507 unique sequences were generated. In parallel, 9,578 reads were generated from a library constructed from the same normalized cDNA population using dideoxy chain termination Sanger sequencing, which upon clustering and assembly generated 4,689 unique sequences. A hybrid assembly of both Sanger- and pyrosequencing-derived ESTs resulted in 34,495 unique sequences with 1,110 sequences (3.2%) that were solely derived from Sanger sequencing alone. A high degree of similarity was seen between P. ultimum sequences and other sequenced plant pathogenic oomycetes with 91% of the hybrid assembly derived sequences > 500 bp having similarity to sequences from plant pathogenic Phytophthora species. An analysis of Gene Ontology assignments revealed a similar representation of molecular function ontologies in the hybrid assembly in comparison to the predicted proteomes of three Phytophthora species, suggesting a broad representation of the P. ultimum transcriptome was present in the normalized cDNA population. P. ultimum sequences with similarity to oomycete RXLR and Crinkler effectors, Kazal-like and cystatin-like protease inhibitors, and elicitins were identified. Sequences with similarity to thiamine biosynthesis enzymes that are lacking in the genome sequences of three Phytophthora species and one downy mildew were identified and could serve as useful phylogenetic markers. Furthermore, we identified 179 candidate simple sequence repeats that can be used for genotyping strains of P. ultimum. CONCLUSION: Through these two technologies, we were able to generate a robust set (approximately 10 Mb) of transcribed sequences for P. ultimum. We were able to identify known sequences present in oomycetes as well as identify novel sequences. An ample number of candidate polymorphic markers were identified in the dataset providing resources for phylogenetic and diagnostic marker development for this species. On a technical level, in spite of the depth possible with 454 FLX platform, the Sanger and pyro-based sequencing methodologies were complementary as each method generated sequences unique to each platform.

Insights into the Musa genome: syntenic relationships to rice and between Musa species.

BACKGROUND: Musa species (Zingiberaceae, Zingiberales) including bananas and plantains are collectively the fourth most important crop in developing countries. Knowledge concerning Musa genome structure and the origin of distinct cultivars has greatly increased over the last few years. Until now, however, no large-scale analyses of Musa genomic sequence have been conducted. This study compares genomic sequence in two Musa species with orthologous regions in the rice genome. RESULTS: We produced 1.4 Mb of Musa sequence from 13 BAC clones, annotated and analyzed them along with 4 previously sequenced BACs. The 443 predicted genes revealed that Zingiberales genes share GC content and distribution characteristics with eudicot and Poaceae genomes. Comparison with rice revealed microsynteny regions that have persisted since the divergence of the Commelinid orders Poales and Zingiberales at least 117 Mya. The previously hypothesized large-scale duplication event in the common ancestor of major cereal lineages within the Poaceae was verified. The divergence time distributions for Musa-Zingiber (Zingiberaceae, Zingiberales) orthologs and paralogs provide strong evidence for a large-scale duplication event in the Musa lineage after its divergence from the Zingiberaceae approximately 61 Mya. Comparisons of genomic regions from M. acuminata and M. balbisiana revealed highly conserved genome structure, and indicated that these genomes diverged circa 4.6 Mya. CONCLUSION: These results point to the utility of comparative analyses between distantly-related monocot species such as rice and Musa for improving our understanding of monocot genome evolution. Sequencing the genome of M. acuminata would provide a strong foundation for comparative genomics in the monocots. In addition a genome sequence would aid genomic and genetic analyses of cultivated Musa polyploid genotypes in research aimed at localizing and cloning genes controlling important agronomic traits for breeding purposes.

Phosphorus stress in common bean: root transcript and metabolic responses.

Phosphorus (P) is an essential element for plant growth. Crop production of common bean (Phaseolus vulgaris), the most important legume for human consumption, is often limited by low P in the soil. Functional genomics were used to investigate global gene expression and metabolic responses of bean plants grown under P-deficient and P-sufficient conditions. P-deficient plants showed enhanced root to shoot ratio accompanied by reduced leaf area and net photosynthesis rates. Transcript profiling was performed through hybridization of nylon filter arrays spotted with cDNAs of 2,212 unigenes from a P deficiency root cDNA library. A total of 126 genes, representing different functional categories, showed significant differential expression in response to P: 62% of these were induced in P-deficient roots. A set of 372 bean transcription factor (TF) genes, coding for proteins with Inter-Pro domains characteristic or diagnostic for TF, were identified from The Institute of Genomic Research/Dana Farber Cancer Institute Common Bean Gene Index. Using real-time reverse transcription-polymerase chain reaction analysis, 17 TF genes were differentially expressed in P-deficient roots; four TF genes, including MYB TFs, were induced. Nonbiased metabolite profiling was used to assess the degree to which changes in gene expression in P-deficient roots affect overall metabolism. Stress-related metabolites such as polyols accumulated in P-deficient roots as well as sugars, which are known to be essential for P stress gene induction. Candidate genes have been identified that may contribute to root adaptation to P deficiency and be useful for improvement of common bean.

Accumulation of genome-specific transcripts, transcription factors and phytohormonal regulators during early stages of fiber cell development in allotetraploid cotton.

Gene expression during the early stages of fiber cell development and in allopolyploid crops is poorly understood. Here we report computational and expression analyses of 32 789 high-quality ESTs derived from Gossypium hirsutum L. Texas Marker-1 (TM-1) immature ovules (GH_TMO). The ESTs were assembled into 8540 unique sequences including 4036 tentative consensus sequences (TCs) and 4504 singletons, representing approximately 15% of the unique sequences in the cotton EST collection. Compared with approximately 178 000 existing ESTs derived from elongating fibers and non-fiber tissues, GH_TMO ESTs showed a significant increase in the percentage of genes encoding putative transcription factors such as MYB and WRKY and genes encoding predicted proteins involved in auxin, brassinosteroid (BR), gibberellic acid (GA), abscisic acid (ABA) and ethylene signaling pathways. Cotton homologs related to MIXTA, MYB5, GL2 and eight genes in the auxin, BR, GA and ethylene pathways were induced during fiber cell initiation but repressed in the naked seed mutant (N1N1) that is impaired in fiber formation. The data agree with the known roles of MYB and WRKY transcription factors in Arabidopsis leaf trichome development and the well-documented phytohormonal effects on fiber cell development in immature cotton ovules cultured in vitro. Moreover, the phytohormonal pathway-related genes were induced prior to the activation of MYB-like genes, suggesting an important role of phytohormones in cell fate determination. Significantly, AA sub-genome ESTs of all functional classifications including cell-cycle control and transcription factor activity were selectively enriched in G. hirsutum L., an allotetraploid derived from polyploidization between AA and DD genome species, a result consistent with the production of long lint fibers in AA genome species. These results suggest general roles for genome-specific, phytohormonal and transcriptional gene regulation during the early stages of fiber cell development in cotton allopolyploids.

Comparative genomics of Brassica oleracea and Arabidopsis thaliana reveal gene loss, fragmentation, and dispersal after polyploidy.

We sequenced 2.2 Mb representing triplicated genome segments of Brassica oleracea, which are each paralogous with one another and homologous with a segmentally duplicated region of the Arabidopsis thaliana genome. Sequence annotation identified 177 conserved collinear genes in the B. oleracea genome segments. Analysis of synonymous base substitution rates indicated that the triplicated Brassica genome segments diverged from a common ancestor soon after divergence of the Arabidopsis and Brassica lineages. This conclusion was corroborated by phylogenetic analysis of protein families. Using A. thaliana as an outgroup, 35% of the genes inferred to be present when genome triplication occurred in the Brassica lineage have been lost, most likely via a deletion mechanism, in an interspersed pattern. Genes encoding proteins involved in signal transduction or transcription were not found to be significantly more extensively retained than those encoding proteins classified with other functions, but putative proteins predicted in the A. thaliana genome were underrepresented in B. oleracea. We identified one example of gene loss from the Arabidopsis lineage. We found evidence for the frequent insertion of gene fragments of nuclear genomic origin and identified four apparently intact genes in noncollinear positions in the B. oleracea and A. thaliana genomes.

Comparative analysis of 87,000 expressed sequence tags from the fumonisin-producing fungus Fusarium verticillioides.

Fusarium verticillioides (teleomorph Gibberella moniliformis) is a pathogen of maize worldwide and produces fumonisins, a family of mycotoxins that have been associated with several animal diseases as well as cancer in humans. In this study, we sought to identify fungal genes that affect fumonisin production and/or the plant-fungal interaction. We generated over 87,000 expressed sequence tags from nine different cDNA libraries that correspond to 11,119 unique sequences and are estimated to represent 80% of the genomic complement of genes. A comparative analysis of the libraries showed that all 15 genes in the fumonisin gene cluster were differentially expressed. In addition, nine candidate fumonisin regulatory genes and a number of genes that may play a role in plant-fungal interaction were identified. Analysis of over 700 FUM gene transcripts from five different libraries provided evidence for transcripts with unspliced introns and spliced introns with alternative 3' splice sites. The abundance of the alternative splice forms and the frequency with which they were found for genes involved in the biosynthesis of a single family of metabolites as well as their differential expression suggest they may have a biological function. Finally, analysis of an EST that aligns to genomic sequence between FUM12 and FUM13 provided evidence for a previously unidentified gene (FUM20) in the FUM gene cluster.

PrP (prion) gene expression in sheep may be modulated by alternative polyadenylation of its messenger RNA.

Scrapie-associated fibrils and their major protein component, PrP or prion protein, accumulate in the brains and some other tissues of all species affected by transmissible spongiform encephalopathies or prion diseases. To investigate the role of PrP gene expression in the hosts of these diseases, we have analysed some characteristics of PrP gene RNA transcripts in sheep and cattle tissues and made comparisons with PrP RNA transcripts in human and mouse tissues. Two PrP messenger RNAs of 4.6 kb and 2.1 kb, the result of alternative polyadenylation, were found first in sheep peripheral tissues and also occurred at low levels in sheep brain and bovine tissues, but not in human and mouse tissues. Our results from transfection assays of murine neuroblastoma cells with constructs expressing different regions of ovine PrP messenger RNA revealed the presence of sequences in the 3' untranslated region of the gene that modulate protein synthesis.