Intraspecific nucleotide variation in Anopheles gambiae: new insights into the biology of malaria vectors

The Anopheles gambiae genome sequence, together with the recent development of molecular tools for genome-wide analysis, promises new insights into the biology of the malaria vector. These insights should help define the best possible breakdown point for interrupting transmission in the mosquito vector. A survey of the intraspecific nucleotide diversity in coding regions of three different mosquito strains showed an average of one single nucleotide polymorphism (SNP) every 125 coding base pairs. High levels of nucleotide polymorphism were observed in mosquito immune-related genes and pathogen recognition receptors harbored higher replacement substitutions. Genotyping at SNP loci in natural populations of An. gambiae from three malaria foci showed contrasting patterns. The distribution of mutation Y443H in the thioester-containing protein 3 (TEP3) gene suggested this mutational event has occurred under selective constraints. Our results show that SNP-based studies will be valuable in identifying the sequence variation associated with phenotypic traits shaping vector competence.     

Candidate odorant receptors from the malaria vector mosquito Anopheles gambiae and evidence of down-regulation in response to blood feeding

Olfaction plays a major role in host preference and blood feeding, integral behaviors for disease transmission by the malaria vector mosquito Anopheles gambiae sensu stricto (henceforth A. gambiae). We have identified four genes encoding candidate odorant receptors from A. gambiae that are selectively expressed in olfactory organs, contain approximately seven transmembrane domains, and show significant similarity to several putative odorant receptors in Drosophila melanogaster. Furthermore, one of the putative A. gambiae odorant receptors exhibits female-specific antennal expression and is down-regulated 12 h after blood feeding, a period during which substantial reduction in olfactory responses to human odorants has been observed. Taken together, these data suggest these genes encode a family of odorant receptors in A. gambiae, whose further study may aid in the design of novel antimalarial programs.     

Genome expression analysis of Anopheles gambiae: responses to injury,bacterial challenge,and malaria infection

The complex gene expression responses of Anopheles gambiae to microbial and malaria challenges, injury, and oxidative stress (in the mosquito and/or a cultured cell line) were surveyed by using cDNA microarrays constructed from an EST-clone collection. The expression profiles were broadly subdivided into induced and down-regulated gene clusters. Gram+ and Gram- bacteria and microbial elicitors up-regulated a diverse set of genes, many belonging to the immunity class, and the response to malaria partially overlapped with this response. Oxidative stress activated a distinctive set of genes, mainly implicated in oxidoreductive processes. Injury up- and down-regulated gene clusters also were distinctive, prominently implicating glycolysis-related genes and citric acid cycle/oxidative phosphorylation/redox-mitochondrial functions, respectively. Cross-comparison of in vivo and in vitro responses indicated the existence of tightly coregulated gene groups that may correspond to gene pathways.     

Robust gut-specific gene expression in transgenic Aedes aegypti mosquitoes

Genetic modification of the vectorial capacity of mosquito vectors of human disease requires promoters capable of driving gene expression with appropriate tissue and stage specificity. We report on the characterization in transgenic Aedes aegypti of two mosquito gut-specific promoters. A 1.4-kb DNA fragment adjacent to the 5' end of the coding region of the Ae. aegypti carboxypeptidase (AeCP) gene and a corresponding 3.4-kb DNA fragment at the 5' end of the Anopheles gambiae carboxypeptidase (AgCP) gene were linked to a firefly luciferase reporter gene and introduced into the Ae. aegypti germ line by using Hermes and mariner (Mos1) transposons. Six independent transgenic lines were obtained with the AeCP construct and one with the AgCP construct. Luciferase mRNA and protein were abundantly expressed in the guts of transgenic mosquitoes in four of the six AeCP lines and in the AgCP line. Expression of the reporter gene was gut-specific and reached peak levels at about 24 h post-blood ingestion. The AeCP and AgCP promoters can be used to drive the expression of genes that hinder parasite development in the mosquito gut.     

Transitive functional annotation by shortest-path analysis of gene expression data

Current methods for the functional analysis of microarray gene expression data make the implicit assumption that genes with similar expression profiles have similar functions in cells. However, among genes involved in the same biological pathway, not all gene pairs show high expression similarity. Here, we propose that transitive expression similarity among genes can be used as an important attribute to link genes of the same biological pathway. Based on large-scale yeast microarray expression data, we use the shortest-path analysis to identify transitive genes between two given genes from the same biological process. We find that not only functionally related genes with correlated expression profiles are identified but also those without. In the latter case, we compare our method to hierarchical clustering, and show that our method can reveal functional relationships among genes in a more precise manner. Finally, we show that our method can be used to reliably predict the function of unknown genes from known genes lying on the same shortest path. We assigned functions for 146 yeast genes that are considered as unknown by the Saccharomyces Genome Database and by the Yeast Proteome Database. These genes constitute around 5% of the unknown yeast ORFome.     

Decoupled evolution of coding region and mRNA expression patterns after gene duplication: implications for the neutralist-selectionist debate

The neutralist perspective on molecular evolution maintains that the vast majority of mutations affecting gene function are neutral or deleterious. After a gene duplication where both genes are retained, it predicts that original and duplicate genes diverge at clock-like rates. This prediction is usually tested for coding sequences, but can also be applied to another important aspect of gene function, the genes' expression pattern. Moreover, if both sequence and expression pattern diverge at clock-like rates, a correlation between divergence in sequence and divergence in expression patterns is expected. Duplicate gene pairs with more highly diverged sequences should also show more highly diverged expression patterns. This prediction is tested for a large sample of duplicated genes in the yeast Saccharomyces cerevisiae, using both genome sequence and microarray expression data. Only a weak correlation is observed, suggesting that coding sequence and mRNA expression patterns of duplicate gene pairs evolve independently and at vastly different rates. Implications of this finding for the neutralist-selectionist debate are discussed.     

Anopheles gambiae pilot gene discovery project: identification of mosquito innate immunity genes from expressed sequence tags generated from immune-competent cell lines

Together with AIDS and tuberculosis, malaria is at the top of the list of devastating infectious diseases. However, molecular genetic studies of its major vector, Anopheles gambiae, are still quite limited. We have conducted a pilot gene discovery project to accelerate progress in the molecular analysis of vector biology, with emphasis on the mosquito's antimalarial immune defense. A total of 5,925 expressed sequence tags were determined from normalized cDNA libraries derived from immune-responsive hemocyte-like cell lines. The 3,242 expressed sequence tag-containing cDNA clones were grouped into 2,380 clone clusters, potentially representing unique genes. Of these, 1,118 showed similarities to known genes from other organisms, but only 27 were identical to previously known mosquito genes. We identified 38 candidate genes, based on sequence similarity, that may be implicated in immune reactions including antimalarial defense; 19 of these were shown experimentally to be inducible by bacterial challenge, lending support to their proposed involvement in mosquito immunity.     

Delineating developmental and metabolic pathways in vivo by expression profiling using the RIKEN set of 18,816 full-length enriched mouse cDNA arrays

We have systematically characterized gene expression patterns in 49 adult and embryonic mouse tissues by using cDNA microarrays with 18,816 mouse cDNAs. Cluster analysis defined sets of genes that were expressed ubiquitously or in similar groups of tissues such as digestive organs and muscle. Clustering of expression profiles was observed in embryonic brain, postnatal cerebellum, and adult olfactory bulb, reflecting similarities in neurogenesis and remodeling. Finally, clustering genes coding for known enzymes into 78 metabolic pathways revealed a surprising coordination of expression within each pathway among different tissues. On the other hand, a more detailed examination of glycolysis revealed tissue-specific differences in profiles of key regulatory enzymes. Thus, by surveying global gene expression by using microarrays with a large number of elements, we provide insights into the commonality and diversity of pathways responsible for the development and maintenance of the mammalian body plan.     

Multivariate regression analysis of distance matrices for testing associations between gene expression patterns and related variables

A fundamental step in the analysis of gene expression and other high-dimensional genomic data is the calculation of the similarity or distance between pairs of individual samples in a study. If one has collected N total samples and assayed the expression level of G genes on those samples, then an N x N similarity matrix can be formed that reflects the correlation or similarity of the samples with respect to the expression values over the G genes. This matrix can then be examined for patterns via standard data reduction and cluster analysis techniques. We consider an alternative to conventional data reduction and cluster analyses of similarity matrices that is rooted in traditional linear models. This analysis method allows predictor variables collected on the samples to be related to variation in the pairwise similarity/distance values reflected in the matrix. The proposed multivariate method avoids the need for reducing the dimensions of a similarity matrix, can be used to assess relationships between the genes used to construct the matrix and additional information collected on the samples under study, and can be used to analyze individual genes or groups of genes identified in different ways. The technique can be used with any high-dimensional assay or data type and is ideally suited for testing subsets of genes defined by their participation in a biochemical pathway or other a priori grouping. We showcase the methodology using three published gene expression data sets.     

Functionalization of a protosynaptic gene expression network

Assembly of a functioning neuronal synapse requires the precisely coordinated synthesis of many proteins. To understand the evolution of this complex cellular machine, we tracked the developmental expression patterns of a core set of conserved synaptic genes across a representative sampling of the animal kingdom. Coregulation, as measured by correlation of gene expression over development, showed a marked increase as functional nervous systems emerged. In the earliest branching animal phyla (Porifera), in which a nearly complete set of synaptic genes exists in the absence of morphological synapses, these "protosynaptic" genes displayed a lack of global coregulation although small modules of coexpressed genes are readily detectable by using network analysis techniques. These findings suggest that functional synapses evolved by exapting preexisting cellular machines, likely through some modification of regulatory circuitry. Evolutionarily ancient modules continue to operate seamlessly within the synapses of modern animals. This work shows that the application of network techniques to emerging genomic and expression data can provide insights into the evolution of complex cellular machines such as the synapse.     

人心力衰竭心肌细胞基因表达谱分析方法初探

目的通过对多中心人类心力衰竭心肌细胞基因表达谱数据整理、分析,探索识别不同表达谱中相同结构蛋白基因的方法。方法选取14个已发表的原发病为扩张型心肌病和（或）缺血性心肌病的心力衰竭基因表达谱,筛选出各表达谱中差异表达≥1.4倍的心肌结构蛋白基因及功能蛋白基因。编写程序自动下载相关基因片段的序列,并建立数据库。根据表达谱中提供的识别号及相关注释信息,用BatchGenAna网站在染色体基因组上进行序列定位,通过定位图显示同一基因不同片段的定位信息。结果通过GenBank序号及芯片中克隆序号得到不同表达谱中相同序号的基因仅有23组,检出率低,漏检了相同基因的不同片段。通过染色体基因组比对、定位,检出不同表达谱中定位在同一染色体区域的相同基因的不同片段共51组,且集中定位在1、2、9、11、12、16号染色体上。通过定位明确显示此基因的编码区及调控区序列。结论基因组定位方法可识别不同表达谱中相同基因的不同片段,明确定位基因编码区及调控区,为后续的功能基因组研究提供可靠信息。     

A ribosomal RNA gene probe differentiates member species of the Anopheles gambiae complex

A 0.59 kilobase DNA fragment cloned from an rDNA cistron of the mosquito Anopheles gambiae can be used as a probe to differentiate between A. gambiae, A. arabiensis, and A. melas, three morphologically identical sibling species in the A. gambiae complex which otherwise can be reliably distinguished only by polytene chromosome banding patterns. Although all are important (and often sympatric) African malaria vectors, their relative roles in malaria transmission have thus far been difficult to assess. The probe, an EcoRI-SalI fragment from the 3' end of the 28S beta coding region of the cistron, is present in all three species, but the species differ uniquely with respect to the location of an EcoRI site in the nontranscribed spacer (NTS) downstream of the fragment. We have routinely used the probe to identify A. gambiae complex mosquitoes to species on the basis of genomic DNA extracted from individual air dried specimens. A single mosquito abdomen provides more than sufficient DNA for the assay, and neither eggs nor a bloodmeal in the abdomen interfere with DNA yield. Moreover, the DNA extraction procedure does not degrade the bloodmeal IgG, so the residual protein pellet can be used to identify the mosquito bloodmeal source. Since the rDNA cistron organization as detected by the probe does not differ between male and female mosquitoes, the probe can be used for either sex. Preliminary experiments show that the probe is equally useful for mosquito larvae and pupae.     

Genome-wide expression profiling of mid-gestation placenta and embryo using a 15,000 mouse developmental cDNA microarray

cDNA microarray technology has been increasingly used to monitor global gene expression patterns in various tissues and cell types. However, applications to mammalian development have been hampered by the lack of appropriate cDNA collections, particularly for early developmental stages. To overcome this problem, a PCR-based cDNA library construction method was used to derive 52,374 expressed sequence tags from pre- and peri-implantation embryos, embryonic day (E) 12.5 female gonad/mesonephros, and newborn ovary. From these cDNA collections, a microarray representing 15,264 unique genes (78% novel and 22% known) was assembled. In initial applications, the divergence of placental and embryonic gene expression profiles was assessed. At stage E12.5 of development, based on triplicate experiments, 720 genes (6.5%) displayed statistically significant differences in expression between placenta and embryo. Among 289 more highly expressed in placenta, 61 placenta-specific genes encoded, for example, a novel prolactin-like protein. The number of genes highly expressed (and frequently specific) for placenta has thereby been increased 5-fold over the total previously reported, illustrating the potential of the microarrays for tissue-specific gene discovery and analysis of mammalian developmental programs.     

Cross-species analysis of biological networks by Bayesian alignment

Complex interactions between genes or proteins contribute a substantial part to phenotypic evolution. Here we develop an evolutionarily grounded method for the cross-species analysis of interaction networks by alignment, which maps bona fide functional relationships between genes in different organisms. Network alignment is based on a scoring function measuring mutual similarities between networks, taking into account their interaction patterns as well as sequence similarities between their nodes. High-scoring alignments and optimal alignment parameters are inferred by a systematic Bayesian analysis. We apply this method to analyze the evolution of coexpression networks between humans and mice. We find evidence for significant conservation of gene expression clusters and give network-based predictions of gene function. We discuss examples where cross-species functional relationships between genes do not concur with sequence similarity.