Complex history of a chromosomal paralogy region: insights from amphioxus aromatic amino acid hydroxylase genes and insulin-related genes

Aromatic amino acid hydroxylase (AAAH) genes and insulin-like genes form part of an extensive paralogy region shared by human chromosomes 11 and 12, thought to have arisen by tetraploidy in early vertebrate evolution. Cloning of a complementary DNA (cDNA) for an amphioxus (Branchiostoma floridae) hydroxylase gene (AmphiPAH) allowed us to investigate the ancestry of the human chromosome 11/12 paralogy region. Molecular phylogenetic evidence reveals that AmphiPAH is orthologous to vertebrate phenylalanine (PAH) genes; the implication is that all three vertebrate AAAH genes arose early in metazoan evolution, predating vertebrates. In contrast, our phylogenetic analysis of amphioxus and vertebrate insulin-related gene sequences is consistent with duplication of these genes during early chordate ancestry. The conclusion is that two tightly linked gene families on human chromosomes 11 and 12 were not duplicated coincidentally. We rationalize this paradox by invoking gene loss in the AAAH gene family and conclude that paralogous genes shared by paralogous chromosomes need not have identical evolutionary histories.

Evolutionary genomics of the Fox genes: origin of gene families and the ancestry of gene clusters

Over the past decade genomic approaches have begun to revolutionise the study of animal diversity. In particular, genome sequencing programmes have spread beyond the traditional model species to encompass an increasing diversity of animals from many different phyla, as well as unicellular eukaryotes that are closely related to the animals. Whole genome sequences allow researchers to establish, with reasonable confidence, the full complement of any particular family of genes in a genome. Comparison of gene complements from appropriate genomes can reveal the evolutionary history of gene families, indicating when both gene diversification and gene loss have occurred. More than that, however, assembled genomes allow the genomic environment in which individual genes are found to be analysed and compared between species. This can reveal how gene diversification occurred. Here, we focus on the Fox genes, drawing from multiple animal genomes to develop an evolutionary framework explaining the timing and mechanism of origin of the diversity of animal Fox genes. Ancient linkages between genes are a prominent feature of the Fox genes, depicting a history of gene clusters, some of which may be relevant to understanding Fox gene function.

Miniature transposable sequences are frequently mobilized in the bacterial plant pathogen Pseudomonas syringae pv. phaseolicola

Mobile genetic elements are widespread in Pseudomonas syringae, and often associate with virulence genes. Genome reannotation of the model bean pathogen P. syringae pv. phaseolicola 1448A identified seventeen types of insertion sequences and two miniature inverted-repeat transposable elements (MITEs) with a biased distribution, representing 2.8% of the chromosome, 25.8% of the 132-kb virulence plasmid and 2.7% of the 52-kb plasmid. Employing an entrapment vector containing sacB, we estimated that transposition frequency oscillated between 2.661025 and 1.161026, depending on the clone, although it was stable for each clone after consecutive transfers in culture media. Transposition frequency was similar for bacteria grown in rich or minimal media, and from cells recovered from compatible and incompatible plant hosts, indicating that growth conditions do not influence transposition in strain 1448A. Most of the entrapped insertions contained a full-length IS801 element, with the remaining insertions corresponding to sequences smaller than any transposable element identified in strain 1448A, and collectively identified as miniature sequences. From these, fragments of 229, 360 and 679-nt of the right end of IS801 ended in a consensus tetranucleotide and likely resulted from one-ended transposition of IS801. An average 0.7% of the insertions analyzed consisted of IS801 carrying a fragment of variable size from gene PSPPH_0008/PSPPH_0017, showing that IS801 can mobilize DNA in vivo. Retrospective analysis of complete plasmids and genomes of P. syringae suggests, however, that most fragments of IS801 are likely the result of reorganizations rather than one-ended transpositions, and that this element might preferentially contribute to genome flexibility by generating homologous regions of recombination. A further miniature sequence previously found to affect host range specificity and virulence, designated MITEPsy1 (100-nt), represented an average 2.4% of the total number of insertions entrapped in sacB, demonstrating for the first time the mobilization of a MITE in bacteria.

Characterization of the complex locus of bean encoding polygalacturonase-inhibiting proteins reveals subfunctionalization for defense against fungi and insects.

Polygalacturonase-inhibiting proteins (PGIPs) are extracellular plant inhibitors of fungal endopolygalacturonases (PGs) that belong to the superfamily of Leu-rich repeat proteins. We have characterized the full complement of pgip genes in the bean (Phaseolus vulgaris) genotype BAT93. This comprises four clustered members that span a 50-kb region and, based on their similarity, form two pairs (Pvpgip1/Pvpgip2 and Pvpgip3/Pvpgip4). Characterization of the encoded products revealed both partial redundancy and subfunctionalization against fungal-derived PGs. Notably, the pair PvPGIP3/PvPGIP4 also inhibited PGs of two mirid bugs (Lygus rugulipennis and Adelphocoris lineolatus). Characterization of Pvpgip genes of Pinto bean showed variations limited to single synonymous substitutions or small deletions. A three-amino acid deletion encompassing a residue previously identified as crucial for recognition of PG of Fusarium moniliforme was responsible for the inability of BAT93 PvPGIP2 to inhibit this enzyme. Consistent with the large variations observed in the promoter sequences, reverse transcription-PCR expression analysis revealed that the different family members differentially respond to elicitors, wounding, and salicylic acid. We conclude that both biochemical and regulatory redundancy and subfunctionalization of pgip genes are important for the adaptation of plants to pathogenic fungi and phytophagous insects.

Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Gene expression is a quantitative trait that can be mapped genetically in structured populations to identify expression quantitative trait loci (eQTL). Genes and regulatory networks underlying complex traits can subsequently be inferred. Using a recently released genome sequence, we have defined cis- and trans-eQTL and their environmental response to low phosphorus (P) availability within a complex plant genome and found hotspots of trans-eQTL within the genome. Interval mapping, using P supply as a covariate, revealed 18,876 eQTL. trans-eQTL hotspots occurred on chromosomes A06 and A01 within Brassica rapa; these were enriched with P metabolism-related Gene Ontology terms (A06) as well as chloroplast-and photosynthesis-related terms (A01). We have also attributed heritability components to measures of gene expression across environments, allowing the identification of novel gene expression markers and gene expression changes associated with low P availability. Informative gene expression markers were used to map eQTL and P use efficiency-related QTL. Genes responsive to P supply had large environmental and heritable variance components. Regulatory loci and genes associated with P use efficiency identified through eQTL analysis are potential targets for further characterization and may have potential for crop improvement.

Genes driving potato tuber initiation and growth: identification based on transcriptional changes using the POCI array

The increasing amount of available expressed gene sequence data makes whole-transcriptome analysis of certain crop species possible. Potato currently has the second largest number of publicly available expressed sequence tag (EST) sequences among the Solanaceae. Most of these ESTs, plus other proprietary sequences, were combined and used to generate a unigene assembly. The set of 246,182 sequences produced 46,345 unigenes, which were used to design a 44K 60-mer oligo array (Potato Oligo Chip Initiative: POCI). In this study, we attempt to identify genes controlling and driving the process of tuber initiation and growth by implementing large-scale transcriptional changes using the newly developed POCI array. Major gene expression profiles could be identified exhibiting differential expression at key developmental stages. These profiles were associated with functional roles in cell division and growth. A subset of genes involved in the regulation of the cell cycle, based on their Gene Ontology classification, exhibit a clear transient upregulation at tuber onset indicating increased cell division during these stages. The POCI array allows the study of potato gene expression on a much broader level than previously possible and will greatly enhance analysis of transcriptional control mechanisms in a wide range of potato research areas. POCI sequence and annotation data are publicly available through the POCI database (http://pgrc.ipk-gatersleben.de/poci).

Distribution of mosquitoes in the South East of Argentina and first report on the analysis based on 18S rDNA and COI sequences

Although Mar del Plata is the most important city on the Atlantic coast of Argentina, mosquitoes inhabiting such area are almost uncharacterized. To increase our knowledge in their distribution, we sampled specimens of natural populations. After the morphological identification based on taxonomic keys, sequences of DNA from small ribosomal subunit (18S rDNA) and cytochrome c oxidase I (COI) genes were obtained from native species and the phylogenetic analysis of these sequences were done. Fourteen species from the genera Uranotaenia, Culex, Ochlerotatus and Psorophora were found and identified. Our 18S rDNA and COI-based analysis indicates the relationships among groups at the supra-species level in concordance with mosquito taxonomy. The introduction and spread of vectors and diseases carried by them are not known in Mar del Plata, but some of the species found in this study were reported as pathogen vectors.

Transcriptional dysregulation of coding and non-coding genes in cellular models of Huntington's disease

HD (Huntington's disease) is a late onset heritable neurodegenerative disorder that is characterized by neuronal dysfunction and death, particularly in the cerebral cortex and medium spiny neurons of the striatum. This is followed by progressive chorea, dementia and emotional dysfunction, eventually resulting in death. HD is caused by an expanded CAG repeat in the first exon of the HD gene that results in an abnormally elongated polyQ (polyglutamine) tract in its protein product, Htt (Huntingtin). Wild-type Htt is largely cytoplasmic; however, in HD, proteolytic N-terminal fragments of Htt form insoluble deposits in both the cytoplasm and nucleus, provoking the idea that mutHtt (mutant Htt) causes transcriptional dysfunction. While a number of specific transcription factors and co-factors have been proposed as mediators of mutHtt toxicity, the causal relationship between these Htt/transcription factor interactions and HD pathology remains unknown. Previous work has highlighted REST [RE1 (repressor element 1)-silencing transcription factor] as one such transcription factor. REST is a master regulator of neuronal genes, repressing their expression. Many of its direct target genes are known or suspected to have a role in HD pathogenesis, including BDNF (brain-derived neurotrophic factor). Recent evidence has also shown that REST regulates transcription of regulatory miRNAs (microRNAs), many of which are known to regulate neuronal gene expression and are dysregulated in HD. Thus repression of miRNAs constitutes a second, indirect mechanism by which REST can alter the neuronal transcriptome in HD. We will describe the evidence that disruption to the REST regulon brought about by a loss of interaction between REST and mutHtt may be a key contributory factor in the widespread dysregulation of gene expression in HD.

Variation in genes related to hepatic lipid metabolism and changes in waist circumference and body weight.

We analysed single nucleotide polymorphisms (SNPs) tagging the genetic variability of six candidate genes (ATF6, FABP1, LPIN2, LPIN3, MLXIPL and MTTP) involved in the regulation of hepatic lipid metabolism, an important regulatory site of energy balance for associations with body mass index (BMI) and changes in weight and waist circumference. We also investigated effect modification by sex and dietary intake. Data of 6,287 individuals participating in the European prospective investigation into cancer and nutrition were included in the analyses. Data on weight and waist circumference were followed up for 6.9 ± 2.5 years. Association of 69 tagSNPs with baseline BMI and annual changes in weight as well as waist circumference were investigated using linear regression analysis. Interactions with sex, GI and intake of carbohydrates, fat as well as saturated, monounsaturated and polyunsaturated fatty acids were examined by including multiplicative SNP-covariate terms into the regression model. Neither baseline BMI nor annual weight or waist circumference changes were significantly associated with variation in the selected genes in the entire study population after correction for multiple testing. One SNP (rs1164) in LPIN2 appeared to be significantly interacting with sex (p = 0.0003) and was associated with greater annual weight gain in men (56.8 ± 23.7 g/year per allele, p = 0.02) than in women (-25.5 ± 19.8 g/year per allele, p = 0.2). With respect to gene-nutrient interaction, we could not detect any significant interactions when accounting for multiple testing. Therefore, out of our six candidate genes, LPIN2 may be considered as a candidate for further studies.

Phylogenetic relationships among species of Anopheles (Nyssorhynchus) (Diptera, Culicidae) based on nuclear and mitochondrial gene sequences

Phylogenetic relationships among 21 species of mosquitoes in subgenus Nyssorhynchus were inferred from the nuclear white and mitochondrial NADH dehydrogenase subunit 6 (ND6) genes. Bayestan phylogenetic methods found that none of the three Sections within Nyssorhynchus (Albimanus, Argyritarsis, Myzorhynchella) were supported in all analyses, although Myzorhynchella was found to be monophyletic at the combined genes Within the Albimanus Section the monophyly of the Stroder Subgroup was strongly supported and within the Myzorhynchella Section Anopheles anrunesi and An lutzu formed a strongly supported monophyletic group The epidemiologically significant Albitarsis Complex showed evidence of paraphyly (relative to An lanet-Myzorhynchella) and discordance across gene trees, and the previously synonomized species of An. dunhami and An goeldii were recovered as sister species Finally, there was evidence of complexes in several species, including An antunesi, An deaneorum, and An. strodei (c) 2010 Elsevier B.V. All rights reserved

Expression pattern of four storage xyloglucan mobilization-related genes during seedling development of the rain forest tree Hymenaea courbaril L.

During seedling establishment, cotyledons of the rain forest tree Hymenaea courbaril mobilize storage cell wall xyloglucan to sustain growth. The polysaccharide is degraded and its products are transported to growing sink tissues. Auxin from the shoot controls the level of xyloglucan hydrolytic enzymes. It is not yet known how important the expression of these genes is for the control of storage xyloglucan degradation. In this work, partial cDNAs of the genes xyloglucan transglycosylase hydrolase (HcXTH1) and beta-galactosidase (HcBGAL1), both related to xyloglucan degradation, and two other genes related to sucrose metabolism [alkaline invertase (HcAlkIN1) and sucrose synthase (HcSUS1)], were isolated. The partial sequences were characterized by comparison with sequences available in the literature, and phylogenetic trees were assembled. Gene expression was evaluated at intervals of 6 h during 24 h in cotyledons, hypocotyl, roots, and leaves, using 45-d-old plantlets. HcXTH1 and HcBGAL1 were correlated to xyloglucan degradation and responded to auxin and light, being down-regulated when transport of auxin was prevented by N-1-naphthylphthalamic acid (NPA) and stimulated by constant light. Genes related to sucrose metabolism, HcAlkIN1 and HcSUS1, responded to inhibition of auxin transport in consonance with storage mobilization in the cotyledons. A model is proposed suggesting that auxin and light are involved in the control of the expression of genes related to storage xyloglucan mobilization in seedlings of H. courbaril. It is concluded that gene expression plays a role in the control of the intercommunication system of the source-sink relationship during seeding growth, favouring its establishment in the shaded environment of the rain forest understorey.

Genetic structure, phylogeny, and biogeography of Brazilian eyelid-less lizards of genera Calyptommatus and Nothobachia (Squamata, Gymnophthalmidae) as inferred from mitochondrial DNA sequences

Calyptommatus and Nothobachia genera of gymnophthalmid lizards are restricted to sandy open habitats on Sao Francisco River margins, northeastern Brazil. Phylogenetic relationships and geographic distribution of the four recognized species of Calyptommatus were analyzed from partial mitochondrial cyt b, 12S, and 16S rRNA genes sequencing, taking allopatric populations of the monotypic Nothobachia ablephara as the outgroup. In Calyptommatus a basal split separated C. sinebrachiatus, a species restricted to the eastern bank of the river, from the three other species. In this clade, C. confusionibus, found on western margin, was recovered as the sister group of the two other species, C. leiolepis and C. nicterus, from opposite margins. According to approximate date estimations, C. sinebrachiatus would have separated from the other congeneric species by 4.4-6.5 my, and C. nicterus, also from eastern bank, would be diverging by 1.8-2.6 my from C. leiolepis, the sister species on the opposite margin. C. confusionibus and C. leiolepis, both from western sandy areas, would be differentiating by 2.8-5.0 my. Divergence times of about 3.0-4.0 my were estimated for allopatric populations of Nothobachia restricted to western margin. Significant differences in 16S rRNA secondary structure relatively to other vertebrates are reported. Distinct evolutionary patterns are proposed for different taxa in those sandy areas, probably related to historical changes in the course of Sao Francisco River. (C) 2010 Elsevier Inc. All rights reserved.

Phylogeny of the Neotropical genus Acestrorhynchus (Ostariophysi: Characiformes) based on nuclear and mitochondrial gene sequences and morphology: A total evidence approach

Acestrorhynchus is the sole genus of the family Acestrorhynchidae which includes 14 species currently recognized as valid. Species of Acestrorhynchus comprise small-to-medium sized piscivorous fishes and have been traditionally grouped on the basis of well-defined color patterns. A recent phylogeny, based on morphological characters, could not resolve the phylogenetic affinities of A. heterolepis and the relationships among the species of the clade formed by A. abbreviatus, A. altus, A. falcatus, A. lacustris, and A. pantaneiro. The simultaneous analysis of two mitochondrial genes (16S and ATP synthase subunits 6 and 8) and one nuclear intron (S7) was able to resolve the latter clade, but the position of A. heterolepis remained unresolved. The combination of the molecular and morphological data sets in a total evidence analysis resulted in a well-resolved hypothesis regarding the phylogenetic relationships of Acestrorhynchus species. (C) 2009 Elsevier Inc. All rights reserved.

Divergence Involving Global Regulatory Gene Mutations in an Escherichia coli Population Evolving under Phosphate Limitation

Many of the important changes in evolution are regulatory in nature. Sequenced bacterial genomes point to flexibility in regulatory circuits but we do not know how regulation is remodeled in evolving bacteria. Here, we study the regulatory changes that emerge in populations evolving under controlled conditions during experimental evolution of Escherichia coli in a phosphate-limited chemostat culture. Genomes were sequenced from five clones with different combinations of phenotypic properties that coexisted in a population after 37 days. Each of the distinct isolates contained a different mutation in 1 of 3 highly pleiotropic regulatory genes (hfq, spoT, or rpoS). The mutations resulted in dissimilar proteomic changes, consistent with the documented effects of hfq, spoT, and rpoS mutations. The different mutations do share a common benefit, however, in that the mutations each redirect cellular resources away from stress responses that are redundant in a constant selection environment. The hfq mutation lowers several individual stress responses as well the small RNA-dependent activation of rpoS translation and hence general stress resistance. The spoT mutation reduces ppGpp levels, decreasing the stringent response as well as rpoS expression. The mutations in and upstream of rpoS resulted in partial or complete loss of general stress resistance. Our observations suggest that the degeneracy at the core of bacterial stress regulation provides alternative solutions to a common evolutionary challenge. These results can explain phenotypic divergence in a constant environment and also how evolutionary jumps and adaptive radiations involve altered gene regulation.

Characterization of the alpha-haemolysin determinant from the human enteropathogenic Escherichia coli O26 plasmid pEO5

The 157-kb conjugative plasmid pEO5 encoding alpha-haemolysin in strains of human enteropathogenic Escherichia coli (EPEC) O26 was investigated for its relationship with EHEC-haemolysin-encoding plasmids of enterohaemorrhagic E. coli (EHEC) O26 and O157 strains. Plasmid pEO5 was found to be compatible with EHEC-virulence plasmids and did not hybridize in Southern blots with plasmid pO157 from the EHEC O157:H7 strain EDL933, indicating that both plasmids were unrelated. A 9227-bp stretch of pEO5 DNA encompassing the entire alpha-hlyCABD operon was sequenced and compared for similarity to plasmid and chromosomally inherited alpha-hly determinants. The alpha-hly determinant of pEO5 (7252 bp) and its upstream region was most similar to corresponding sequences of the murine E. coli alpha-hly plasmid pHly152, in particular, the structural alpha-hlyCABD genes (99.2% identity) and the regulatory hlyR regions (98.8% identity). pEO5 and alpha-hly plasmids of EPEC O26 strains from humans and cattle were very similar for the regions encompassing the structural alpha-hlyCABD genes. The major difference found between the hly regions of pHly152 and pEO5 is caused by the insertion of an IS2 element upstream of the hlyC gene in pHly152. The presence of transposon-like structures at both ends of the alpha-hly sequence indicates that this pEO5 virulence factor was probably acquired by horizontal gene transfer.