RPS4Y gene family evolution in primates

Background: The RPS4 gene codifies for ribosomal protein S4, a very well-conserved protein present in all kingdoms. In primates, RPS4 is codified by two functional genes located on both sex chromosomes: the RPS4X and RPS4Y genes. In humans, RPS4Y is duplicated and the Y chromosome therefore carries a third functional paralog: RPS4Y2, which presents a testis-specific expression pattern. Results: DNA sequence analysis of the intronic and cDNA regions of RPS4Y genes from species covering the entire primate phylogeny showed that the duplication event leading to the second Y-linked copy occurred after the divergence of New World monkeys, about 35 million years ago. Maximum likelihood analyses of the synonymous and non-synonymous substitutions revealed that positive selection was acting on RPS4Y2 gene in the human lineage, which represents the first evidence of positive selection on a ribosomal protein gene. Putative positive amino acid replacements affected the three domains of the protein: one of these changes is located in the KOW protein domain and affects the unique invariable position of this motif, and might thus have a dramatic effect on the protein function.Conclusion: Here, we shed new light on the evolutionary history of RPS4Y gene family, especially on that of RPS4Y2. The results point that the RPS4Y1 gene might be maintained to compensate gene dosage between sexes, while RPS4Y2 might have acquired a new function, at least in the lineage leading to humans.

EGASP: the human ENCODE Genome Annotation Assessment Project

Background: Non-long terminal repeat (non-LTR) retrotransposons have contributed to shaping the structure and function of genomes. In silico and experimental approaches have been used to identify the non-LTR elements of the urochordate Ciona intestinalis. Knowledge of the types and abundance of non-LTR elements in urochordates is a key step in understanding their contribution to the structure and function of vertebrate genomes. Results: Consensus elements phylogenetically related to the I, LINE1, LINE2, LOA and R2 elements of the 14 eukaryotic non-LTR clades are described from C. intestinalis. The ascidian elements showed conservation of both the reverse transcriptase coding sequence and the overall structural organization seen in each clade. The apurinic/apyrimidinic endonuclease and nucleic-acid-binding domains encoded upstream of the reverse transcriptase, and the RNase H and the restriction enzyme-like endonuclease motifs encoded downstream of the reverse transcriptase were identified in the corresponding Ciona families. Conclusions: The genome of C. intestinalis harbors representatives of at least five clades of non-LTR retrotransposons. The copy number per haploid genome of each element is low, less than 100, far below the values reported for vertebrate counterparts but within the range for protostomes. Genomic and sequence analysis shows that the ascidian non-LTR elements are unmethylated and flanked by genomic segments with a gene density lower than average for the genome. The analysis provides valuable data for understanding the evolution of early chordate genomes and enlarges the view on the distribution of the non-LTR retrotransposons in eukaryotes.

Hundreds of variants clustered in genomic loci and biological pathways affect human height.

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

Hereditary hepatic and systemic amyloidosis caused by a new deletion/insertion mutation in the apolipoprotein AI gene.

We report a Spanish family with autosomal-dominant non-neuropathic hereditary amyloidosis with a unique hepatic presentation and death from liver failure, usually by the sixth decade. The disease is caused by a previously unreported deletion/insertion mutation in exon 4 of the apolipoprotein AI (apoAI) gene encoding loss of residues 60-71 of normal mature apoAI and insertion at that position of two new residues, ValThr. Affected individuals are heterozygous for this mutation and have both normal apoAI and variant molecules bearing one extra positive charge, as predicted from the DNA sequence. The amyloid fibrils are composed exclusively of NH2-terminal fragments of the variant, ending mainly at positions corresponding to residues 83 and 92 in the mature wild-type sequence. Amyloid fibrils derived from the other three known amyloidogenic apoAI variants are also composed of similar NH2-terminal fragments. All known amyloidogenic apoAI variants carry one extra positive charge in this region, suggesting that it may be responsible for their enhanced amyloidogenicity. In addition to causing a new phenotype, this is the first deletion mutation to be described in association with hereditary amyloidosis and it significantly extends the value of the apoAI model for investigation of molecular mechanisms of amyloid fibrillogenesis.

RTS,S/AS02A malaria vaccine does not induce parasite CSP T cell epitope selection and reduces multiplicity of infection

Objective: The candidate malaria vaccine RTS,S/AS02A is a recombinant protein containing part of the circumsporozoite protein (CSP) sequence of Plasmodium falciparum, linked to the hepatitis B surface antigen and formulated in the proprietary adjuvant system AS02A. In a recent trial conducted in children younger than age five in southern Mozambique, the vaccinedemonstrated significant and sustained efficacy against both infection and clinical disease. In a follow-up study to the main trial, breakthrough infections identified in the trial were examined to determine whether the distribution of csp sequences was affected by the vaccine and to measure the multiplicity of infecting parasite genotypes. Design: P. falciparum DNA from isolates collected during the trial was used for genotype studies. Setting: The main trial was carried out in the Manhiça district, Maputo province, Mozambique, between April 2003 and May 2004. Participants: Children from the two cohorts of the main trial provided parasite isolates as follows: children from Cohort 1 who were admitted to hospital with clinical malaria; children from Cohort 1 who were parasite-positive in a cross-sectional survey at study month 8.5; children from Cohort 2 identified as parasite-positive during follow-up by active detection of infection. Outcome: Divergence of DNA sequence encoding the CSP T cell-epitope region sequence from that of the vaccine sequence was measured in 521 isolates. The number of distinct P. falciparum genotypes was also determined. Results: We found no evidence that parasite genotypes from children in the RTS,S/AS02A arm were more divergent than those receiving control vaccines. For Cohort 1 (survey at studymonth 8.5) and Cohort 2, infections in the vaccine group contained significantly fewer genotypes than those in the control group, (p 1/4 0.035, p 1/4 0.006), respectively, for the two cohorts. This was not the case for children in Cohort 1 who were admitted to hospital (p 1/4 0.478). Conclusions: RTS,S/AS02A did not select for genotypes encoding divergent T cell epitopes in the C-terminal region of CSP in this trial. In both cohorts, there was a modest reduction in the mean number of parasite genotypes harboured by vaccinated children compared with controls, but only among those with asymptomatic infections.

Molecular evidence of interhuman transmission in an outbreak of Pneumocystis jirovecii pneumonia among renal transplant recipients.

S. Gianella, L. Haeberli, B. Joos, B. Ledergerber, R.P. Wüthrich, R. Weber, H. Kuster, P.M. Hauser, T. Fehr, N.J. Mueller. Molecular evidence of interhuman transmission in an outbreak of Pneumocystis jirovecii pneumonia among renal transplant recipients. Transpl Infect Dis 2009. All rights reserved Abstract: Pneumocystis jirovecii pneumonia (PCP) remains an important cause of morbidity and mortality in immunocompromised individuals. The epidemiology and pathogenesis of this infection are poorly understood, and the exact mode of transmission remains unclear. Recent studies reported clusters of PCP among immunocompromised patients, raising the suspicion of interhuman transmission. An unexpected increase of the incidence of PCP cases in our nephrology outpatient clinic prompted us to conduct a detailed analysis. Genotyping of 7 available specimens obtained from renal transplant recipients was performed using multi-locus DNA sequence typing (MLST). Fragments of 4 variable regions of the P. jirovecii genome (ITS1, 26S, mt26S, beta-tubulin) were sequenced and compared with those of 4 independent control patients. MLST analysis revealed identical sequences of the 4 regions among all 7 renal allograft recipients with available samples, indicating an infection with the same P. jirovecii genotype. We observed that all but 1 of the 19 PCP-infected transplant recipients had at least 1 concomitant visit with another PCP-infected patient within a common waiting area. This study provides evidence that nosocomial transmission among immunocompromised patients may have occurred in our nephrology outpatient clinic. Our findings have epidemiological implications and suggest that prolonged chemoprophylaxis for PCP may be warranted in an era of more intense immunosuppression.

Effects of season, sex and body size on the feeding ecology of turtle-headed seasnakes (Emydocephalus annulatus) on IndoPacific inshore coral reefs

In terrestrial snakes, many cases of intraspecific shifts in dietary habits as a function of predator sex and body size are driven by gape-limitation - and hence, are most common in species that feed on relatively large prey, and exhibit a wide body-size range. Our data on seasnakes reveal an alternative mechanism for intraspecific niche partitioning, based on sex-specific seasonal anorexia induced by reproductive activities. Turtle-headed seasnakes (Emydocephalus annulatus) on coral reefs in the New Caledonian Lagoon feed entirely on the eggs of demersal-spawning fishes. DNA sequence data (cytochrome b gene) on eggs that we palpated from stomachs of 37 snakes showed that despite this ontogenetic-stage specialization, the prey come from a taxonomically diverse array of species including damselfish (41% of samples, at least 5 species), blennies (41%, 4 species) and gobies (19%, 5 species). The composition of snake diets shifted seasonally (with damselfish dominating in winter but not summer), presumably reflecting seasonality of fish reproduction. That seasonal shift affects male and female snakes differently, because reproduction is incompatible with foraging. Adult female seasnakes ceased feeding when they became heavily distended with developing embryos in late summer, and males ceased feeding while they were mate-searching in winter. The sex divergence in foraging habits may be amplified by sexual size dimorphism; females grow larger than males, and larger snakes (of both sexes) feed more on damselfish (which often lay their eggs in exposed sites) than on blennies and gobies (whose eggs are hidden within narrow crevices). Specific features of reproductive biology of coral-reef fish (seasonality and nest type) have generated intraspecific niche partitioning in these seasnakes, by mechanisms different from those that apply to terrestrial snakes.

Differences in the evolutionary history of disease genes affected by dominant or recessive mutations

Background: Global analyses of human disease genes by computational methods have yielded important advances in the understanding of human diseases. Generally these studies have treated the group of disease genes uniformly, thus ignoring the type of disease-causing mutations (dominant or recessive). In this report we present a comprehensive study of the evolutionary history of autosomal disease genes separated by mode of inheritance.Results: We examine differences in protein and coding sequence conservation between dominant and recessive human disease genes. Our analysis shows that disease genes affected by dominant mutations are more conserved than those affected by recessive mutations. This could be a consequence of the fact that recessive mutations remain hidden from selection while heterozygous. Furthermore, we employ functional annotation analysis and investigations into disease severity to support this hypothesis. Conclusion: This study elucidates important differences between dominantly- and recessively-acting disease genes in terms of protein and DNA sequence conservation, paralogy and essentiality. We propose that the division of disease genes by mode of inheritance will enhance both understanding of the disease process and prediction of candidate disease genes in the future.

The evolution of the thyroid hormone distributor protein transthyretin in the order insectivora, class mammalia.

Thyroid hormones are involved in the regulation of growth and metabolism in all vertebrates. Transthyretin is one of the extracellular proteins with high affinity for thyroid hormones which determine the partitioning of these hormones between extracellular compartments and intracellular lipids. During vertebrate evolution, both the tissue pattern of expression and the structure of the gene for transthyretin underwent characteristic changes. The purpose of this study was to characterize the position of Insectivora in the evolution of transthyretin in eutherians, a subclass of Mammalia. Transthyretin was identified by thyroxine binding and Western analysis in the blood of adult shrews, hedgehogs, and moles. Transthyretin is synthesized in the liver and secreted into the bloodstream, similar to the situation for other adult eutherians, birds, and diprotodont marsupials, but different from that for adult fish, amphibians, reptiles, monotremes, and Australian polyprotodont marsupials. For the characterization of the structure of the gene and the processing of mRNA for transthyretin, cDNA libraries were prepared from RNA from hedgehog and shrew livers, and full-length cDNA clones were isolated and sequenced. Sections of genomic DNA in the regions coding for the splice sites between exons 1 and 2 were synthesized by polymerase chain reaction and sequenced. The location of splicing was deduced from comparison of genomic with cDNA nucleotide sequences. Changes in the nucleotide sequence of the transthyretin gene during evolution are most pronounced in the region coding for the N-terminal region of the protein. Both the derived overall amino sequences and the N-terminal regions of the transthyretins in Insectivora were found to be very similar to those in other eutherians but differed from those found in marsupials, birds, reptiles, amphibians, and fish. Also, the pattern of transthyretin precursor mRNA splicing in Insectivora was more similar to that in other eutherians than to that in marsupials, reptiles, and birds. Thus, in contrast to the marsupials, with a different pattern of transthyretin gene expression in the evolutionarily "older" polyprotodonts compared with the evolutionarily "younger" diprotodonts, no separate lineages of transthyretin evolution could be identified in eutherians. We conclude that transthyretin gene expression in the liver of adult eutherians probably appeared before the branching of the lineages leading to modern eutherian species.

International Union of Pharmacology. LXI. Peroxisome proliferator-activated receptors.

The three peroxisome proliferator-activated receptors (PPARs) are ligand-activated transcription factors of the nuclear hormone receptor superfamily. They share a high degree of structural homology with all members of the superfamily, particularly in the DNA-binding domain and ligand- and cofactor-binding domain. Many cellular and systemic roles have been attributed to these receptors, reaching far beyond the stimulation of peroxisome proliferation in rodents after which they were initially named. PPARs exhibit broad, isotype-specific tissue expression patterns. PPARalpha is expressed at high levels in organs with significant catabolism of fatty acids. PPARbeta/delta has the broadest expression pattern, and the levels of expression in certain tissues depend on the extent of cell proliferation and differentiation. PPARgamma is expressed as two isoforms, of which PPARgamma2 is found at high levels in the adipose tissues, whereas PPARgamma1 has a broader expression pattern. Transcriptional regulation by PPARs requires heterodimerization with the retinoid X receptor (RXR). When activated by a ligand, the dimer modulates transcription via binding to a specific DNA sequence element called a peroxisome proliferator response element (PPRE) in the promoter region of target genes. A wide variety of natural or synthetic compounds was identified as PPAR ligands. Among the synthetic ligands, the lipid-lowering drugs, fibrates, and the insulin sensitizers, thiazolidinediones, are PPARalpha and PPARgamma agonists, respectively, which underscores the important role of PPARs as therapeutic targets. Transcriptional control by PPAR/RXR heterodimers also requires interaction with coregulator complexes. Thus, selective action of PPARs in vivo results from the interplay at a given time point between expression levels of each of the three PPAR and RXR isotypes, affinity for a specific promoter PPRE, and ligand and cofactor availabilities.

Leishmania major: histone H1 gene expression from the sw3 locus.

Histone H1 in the parasitic protozoan Leishmania is a developmentally regulated protein encoded by the sw3 gene. Here we report that histone H1 variants exist in different Leishmania species and strains of L. major and that they are encoded by polymorphic genes. Amplification of the sw3 gene from the genome of three strains of L. major gave rise to different products in each strain, suggesting the presence of a multicopy gene family. In L. major, these genes were all restricted to a 50-kb Bg/II fragment found on a chromosomal band of 1.3 Mb (chromosome 27). The detection of RFLPs in this locus demonstrated its heterogeneity within several species and strains of Leishmania. Two different copies of sw3 (sw3.0 and sw3.1) were identified after screening a cosmid library containing L. major strain Friedlin genomic DNA. They were identical in their 5' UTRs and open reading frames, but differed in their 3' UTRs. With respect to the originally cloned copy of sw3 from L. major strain LV39, their open reading frames lacked a repeat unit of 9 amino acids. Immunoblots of L. guyanensis parasites transfected with these cosmids revealed that both copies could give rise to the histone H1 protein. The characterization of this locus will now make possible a detailed analysis of the function of histone H1 in Leishmania, as well as permit the dissection of the molecular mechanisms governing the developmental regulation of the sw3 gene.

Biosynthesis of pyochelin and dihydroaeruginoic acid requires the iron-regulated pchDCBA operon in Pseudomonas aeruginosa.

The high-affinity siderophore salicylate is an intermediate in the biosynthetic pathway of pyochelin, another siderophore and chelator of transition metal ions, in Pseudomonas aeruginosa. The 2.5-kb region upstream of the salicylate biosynthetic genes pchBA was sequenced and found to contain two additional, contiguous genes, pchD and pchC, having the same orientation. The deduced amino acid sequence of the 60-kDa PchD protein was similar to those of the EntE protein (2,3-dihydroxybenzoate-AMP ligase) of Escherichia coli and other adenylate-forming enzymes, suggesting that salicylate might be adenylated at the carboxyl group by PchD. The 28-kDa PchC protein showed similarities to thioesterases of prokaryotic and eukaryotic origin and might participate in the release of the product(s) formed from activated salicylate. One potential product, dihydroaeruginoate (Dha), was identified in culture supernatants of iron-limited P. aeruginosa cells. The antifungal antibiotic Dha is thought to arise from the reaction of salicylate with cysteine, followed by cyclization of cysteine. Inactivation of the chromosomal pchD gene by insertion of the transcription and translation stop element omega Sm/Sp abolished the production of Dha and pyochelin, implying that PchD-mediated activation of salicylate may be a common first step in the synthesis of both metabolites. Furthermore, the pchD::omega Sm/Sp mutation had a strong polar effect on the expression of the pchBA genes, i.e., on salicylate synthesis, indicating that the pchDCBA genes constitute a transcriptional unit. A full-length pchDCBA transcript of ca. 4.4 kb could be detected in iron-deprived, growing cells of P. aeruginosa. Transcription of pchD started at tandemly arranged promoters, which overlapped with two Fur boxes (binding sites for the ferric uptake regulator) and the promoter of the divergently transcribed pchR gene encoding an activator of pyochelin biosynthesis. This promoter arrangement allows tight iron-mediated repression of the pchDCBA operon.

Trichophyton rubrum secreted and membrane-associated carboxypeptidases

Dermatophytes are the most common agents of superficial mycoses, and exclusively infect stratum corneum, nails or hair. Therefore, secreted proteolytic activity is considered a virulence trait of these fungi. In a medium containing protein as a sole nitrogen and carbon source Trichophyton rubrum secretes a metallocarboxypeptidase (TruMcpA) of the M14 family according to the MEROPS proteolytic enzyme database. TruMcpA is homologous to human pancreatic carboxypeptidase A, and is synthesized as a precursor in a preproprotein form. The propeptide is removed to generate the mature active enzyme alternatively by either one of two subtilisins which are concomitantly secreted by the fungus. In addition, T. rubrum was shown to possess two genes (TruSCPA and TruSCPB) encoding serine carboxypeptidases of the S10 family which are homologues of the previously characterized Aspergillus and Penicillium secreted acid carboxypeptidases. However, in contrast to the Aspergillus and Penicillium homologues, TruScpA and TruScpB enzymes are not secreted into the environment, but are membrane-associated with a glycosylphosphatidylinositol (GPI) anchor. During infection, T. rubrum secreted and GPI-anchored carboxypeptidases may contribute to fungal virulence by cooperating with previously characterized endoproteases and aminopeptidases in the degradation of compact keratinized tissues into assimilable amino acids and short peptides.

GENCODE: producing a reference annotation for ENCODE.

BACKGROUND: The GENCODE consortium was formed to identify and map all protein-coding genes within the ENCODE regions. This was achieved by a combination of initial manual annotation by the HAVANA team, experimental validation by the GENCODE consortium and a refinement of the annotation based on these experimental results. RESULTS: The GENCODE gene features are divided into eight different categories of which only the first two (known and novel coding sequence) are confidently predicted to be protein-coding genes. 5' rapid amplification of cDNA ends (RACE) and RT-PCR were used to experimentally verify the initial annotation. Of the 420 coding loci tested, 229 RACE products have been sequenced. They supported 5' extensions of 30 loci and new splice variants in 50 loci. In addition, 46 loci without evidence for a coding sequence were validated, consisting of 31 novel and 15 putative transcripts. We assessed the comprehensiveness of the GENCODE annotation by attempting to validate all the predicted exon boundaries outside the GENCODE annotation. Out of 1,215 tested in a subset of the ENCODE regions, 14 novel exon pairs were validated, only two of them in intergenic regions. CONCLUSION: In total, 487 loci, of which 434 are coding, have been annotated as part of the GENCODE reference set available from the UCSC browser. Comparison of GENCODE annotation with RefSeq and ENSEMBL show only 40% of GENCODE exons are contained within the two sets, which is a reflection of the high number of alternative splice forms with unique exons annotated. Over 50% of coding loci have been experimentally verified by 5' RACE for EGASP and the GENCODE collaboration is continuing to refine its annotation of 1% human genome with the aid of experimental validation.

The cell cycle–regulated genes of schizosaccharomyces pombe

Many genes are regulated as an innate part of the eukaryotic cell cycle, and a complex transcriptional network helps enable the cyclic behavior of dividing cells. This transcriptional network has been studied in Saccharomyces cerevisiae (budding yeast) and elsewhere. To provide more perspective on these regulatory mechanisms, we have used microarrays to measure gene expression through the cell cycle of Schizosaccharomyces pombe (fission yeast). The 750 genes with the most significant oscillations were identified and analyzed. There were two broad waves of cell cycle transcription, one in early/mid G2 phase, and the other near the G2/M transition. The early/mid G2 wave included many genes involved in ribosome biogenesis, possibly explaining the cell cycle oscillation in protein synthesis in S.pombe. The G2/M wave included at least three distinctly regulated clusters of genes: one large cluster including mitosis, mitotic exit, and cell separation functions, one small cluster dedicated to DNA replication, and another small cluster dedicated to cytokinesis and division. S. pombe cell cycle genes have relatively long, complex promoters containing groups of multiple DNA sequence motifs, often of two, three, or more different kinds. Many of the genes, transcription factors, and regulatory mechanisms are conserved between S. pombe and S. cerevisiae. Finally, we found preliminary evidence for a nearly genome-wide oscillation in gene expression: 2,000 or more genes undergo slight oscillations in expression as a function of the cell cycle, although whether this is adaptive, or incidental to other events in the cell, such as chromatin condensation, we do not know.