Extensive natural variation for cellular hydrogen peroxide release is genetically controlled.

Natural variation in DNA sequence contributes to individual differences in quantitative traits. While multiple studies have shown genetic control over gene expression variation, few additional cellular traits have been investigated. Here, we investigated the natural variation of NADPH oxidase-dependent hydrogen peroxide (H(2)O(2) release), which is the joint effect of reactive oxygen species (ROS) production, superoxide metabolism and degradation, and is related to a number of human disorders. We assessed the normal variation of H(2)O(2) release in lymphoblastoid cell lines (LCL) in a family-based 3-generation cohort (CEPH-HapMap), and in 3 population-based cohorts (KORA, GenCord, HapMap). Substantial individual variation was observed, 45% of which were associated with heritability in the CEPH-HapMap cohort. We identified 2 genome-wide significant loci of Hsa12 and Hsa15 in genome-wide linkage analysis. Next, we performed genome-wide association study (GWAS) for the combined KORA-GenCord cohorts (n = 279) using enhanced marker resolution by imputation (>1.4 million SNPs). We found 5 significant associations (p<5.00×10-8) and 54 suggestive associations (p<1.00×10-5), one of which confirmed the linked region on Hsa15. To replicate our findings, we performed GWAS using 58 HapMap individuals and ∼2.1 million SNPs. We identified 40 genome-wide significant and 302 suggestive SNPs, and confirmed genome signals on Hsa1, Hsa12, and Hsa15. Genetic loci within 900 kb from the known candidate gene p67phox on Hsa1 were identified in GWAS in both cohorts. We did not find replication of SNPs across all cohorts, but replication within the same genomic region. Finally, a highly significant decrease in H(2)O(2) release was observed in Down Syndrome (DS) individuals (p<2.88×10-12). Taken together, our results show strong evidence of genetic control of H(2)O(2) in LCL of healthy and DS cohorts and suggest that cellular phenotypes, which themselves are also complex, may be used as proxies for dissection of complex disorders.

Comparative analyses of classical phenotypic method and ribosomal RNA gene sequencing for identification of medically relevant Candida species

As the distribution of Candida species and their susceptibility to antifungal agents have changed, a new means of accurately and rapidly identifying these species is necessary for the successful early resolution of infection and the subsequent reduction of morbidity and mortality. The current work aimed to evaluate ribosomal RNA gene sequencing for the identification of medically relevant Candida species in comparison with a standard phenotypic method. Eighteen reference strains (RSs), 69 phenotypically identified isolates and 20 inconclusively identified isolates were examined. Internal transcribed spaces (ITSs) and D1/D2 of the 26S ribosomal RNA gene regions were used as targets for sequencing. Additionally, the sequences of the ITS regions were used to establish evolutionary relationships. The sequencing of the ITS regions was successful for 88% (94/107) of the RS and isolates, whereas 100% of the remaining 12% (13/107) of the samples were successfully analysed by sequencing the D1/D2 region. Similarly, genotypic analysis identified all of the RS and isolates, including the 20 isolates that were not phenotypically identified. Phenotypic analysis, however, misidentified 10% (7/69) of the isolates. Phylogenetic analysis allowed the confirmation of the relationships between evolutionarily close species. Currently, the use of genotypic methods is necessary for the correct identification of Candida species.

The gamma subunit is a specific component of the Na,K-ATPase and modulates its transport function.

The role of small, hydrophobic peptides that are associated with ion pumps or channels is still poorly understood. By using the Xenopus oocyte as an expression system, we have characterized the structural and functional properties of the gamma peptide which co-purifies with Na,K-ATPase. Immuno-radiolabeling of epitope-tagged gamma subunits in intact oocytes and protease protection assays show that the gamma peptide is a type I membrane protein lacking a signal sequence and exposing the N-terminus to the extracytoplasmic side. Co-expression of the rat or Xenopus gamma subunit with various proteins in the oocyte reveals that it specifically associates only with isozymes of Na,K-ATPase. The gamma peptide does not influence the formation and cell surface expression of functional Na,K-ATPase alpha-beta complexes. On the other hand, the gamma peptide itself needs association with Na,K-ATPase in order to be stably expressed in the oocyte and to be transported efficiently to the plasma membrane. Gamma subunits do not associate with individual alpha or beta subunits but only interact with assembled, transport-competent alpha-beta complexes. Finally, electrophysiological measurements indicate that the gamma peptide modulates the K+ activation of Na,K pumps. These data document for the first time the membrane topology, the specificity of association and a potential functional role for the gamma subunit of Na,K-ATPase.

Novel Chlamydiales strains isolated from a water treatment plant.

Chlamydiae are obligate intracellular bacteria infecting free-living amoebae, vertebrates and some invertebrates. Novel members are regularly discovered, and there is accumulating evidence supporting a very important diversity of chlamydiae in the environment. In this study, we investigated the presence of chlamydiae in a drinking water treatment plant. Samples were used to inoculate Acanthamoeba monolayers (Acanthamoeba co-culture), and to recover autochthonous amoebae onto non-nutritive agar. Chlamydiae were searched for by a pan-chlamydia 16S rRNA gene PCR from both Acanthamoeba co-cultures and autochthonous amoebae, and phylotypes determined by 16S rRNA gene sequencing. Autochthonous amoebae also were identified by 18S rRNA gene amplification and sequencing. From a total of 79 samples, we recovered eight chlamydial strains by Acanthamoeba co-culture, but only one of 28 amoebae harboured a chlamydia. Sequencing results and phylogenetic analysis showed our strains belonging to four distinct chlamydial lineages. Four strains, including the strain recovered within its natural host, belonged to the Parachlamydiaceae; two closely related strains belonged to the Criblamydiaceae; two distinct strains clustered with Rhabdochlamydia spp.; one strain clustered only with uncultured environmental clones. Our results confirmed the usefulness of amoeba co-culture to recover novel chlamydial strains from complex samples and demonstrated the huge diversity of chlamydiae in the environment, by identifying several new species including one representing the first strain of a new family.

Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.

Knockout confirmation for Hurries: rapid genotype identification of Trypanosoma cruzi transfectants by polymerase chain reaction directly from liquid culture

Gene knockout is a widely used approach to evaluate loss-of-function phenotypes and it can be facilitated by the incorporation of a DNA cassette having a drug-selectable marker. Confirmation of the correct knockout cassette insertion is an important step in gene removal validation and has generally been performed by polymerase chain reaction (PCR) assays following a time-consuming DNA extraction step. Here, we show a rapid procedure for the identification of Trypanosoma cruzi transfectants by PCR directly from liquid culture - without prior DNA extraction. This simple approach enabled us to generate PCR amplifications from different cultures varying from 106-108 cells/mL. We also show that it is possible to combine different primer pairs in a multiplex detection reaction and even to achieve knockout confirmation with an extremely simple interpretation of a real-time PCR result. Using the “culture PCR” approach, we show for the first time that we can assess different DNA sequence combinations by PCR directly from liquid culture, saving time in several tasks for T. cruzi genotype interrogation.

Identification of structural variation in mouse genomes.

Structural variation is variation in structure of DNA regions affecting DNA sequence length and/or orientation. It generally includes deletions, insertions, copy-number gains, inversions, and transposable elements. Traditionally, the identification of structural variation in genomes has been challenging. However, with the recent advances in high-throughput DNA sequencing and paired-end mapping (PEM) methods, the ability to identify structural variation and their respective association to human diseases has improved considerably. In this review, we describe our current knowledge of structural variation in the mouse, one of the prime model systems for studying human diseases and mammalian biology. We further present the evolutionary implications of structural variation on transposable elements. We conclude with future directions on the study of structural variation in mouse genomes that will increase our understanding of molecular architecture and functional consequences of structural variation.

Polyphasic characterisation of Burkholderia cepaciacomplex species isolated from children with cystic fibrosis

Cystic fibrosis (CF) patients with Burkholderia cepacia complex (Bcc) pulmonary infections have high morbidity and mortality. The aim of this study was to compare different methods for identification of Bcc species isolated from paediatric CF patients. Oropharyngeal swabs from children with CF were used to obtain isolates of Bcc samples to evaluate six different tests for strain identification. Conventional (CPT) and automatised (APT) phenotypic tests, polymerase chain reaction (PCR)-recA, restriction fragment length polymorphism-recA, recAsequencing, and matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) were applied. Bacterial isolates were also tested for antimicrobial susceptibility. PCR-recA analysis showed that 36 out of the 54 isolates were Bcc. Kappa index data indicated almost perfect agreement between CPT and APT, CPT and PCR-recA, and APT and PCR-recA to identify Bcc, and MALDI-TOF and recAsequencing to identify Bcc species. The recAsequencing data and the MALDI-TOF data agreed in 97.2% of the isolates. Based on recA sequencing, the most common species identified were Burkholderia cenocepacia IIIA (33.4%),Burkholderia vietnamiensis (30.6%), B. cenocepaciaIIIB (27.8%), Burkholderia multivorans (5.5%), and B. cepacia (2.7%). MALDI-TOF proved to be a useful tool for identification of Bcc species obtained from CF patients, although it was not able to identify B. cenocepacia subtypes.

Epigenetic mechanisms in non-alcoholic fatty liver disease: An emerging field.

Non-alcoholic fatty liver disease (NAFLD) is an emerging health concern in both developed and non-developed world, encompassing from simple steatosis to non-alcoholic steatohepatitis (NASH), cirrhosis and liver cancer. Incidence and prevalence of this disease are increasing due to the socioeconomic transition and change to harmful diet. Currently, gold standard method in NAFLD diagnosis is liver biopsy, despite complications and lack of accuracy due to sampling error. Further, pathogenesis of NAFLD is not fully understood, but is well-known that obesity, diabetes and metabolic derangements played a major role in disease development and progression. Besides, gut microbioma and host genetic and epigenetic background could explain considerable interindividual variability. Knowledge that epigenetics, heritable events not caused by changes in DNA sequence, contribute to development of diseases has been a revolution in the last few years. Recently, evidences are accumulating revealing the important role of epigenetics in NAFLD pathogenesis and in NASH genesis. Histone modifications, changes in DNA methylation and aberrant profiles or microRNAs could boost development of NAFLD and transition into clinical relevant status. PNPLA3 genotype GG has been associated with a more progressive disease and epigenetics could modulate this effect. The impact of epigenetic on NAFLD progression could deserve further applications on therapeutic targets together with future non-invasive methods useful for the diagnosis and staging of NAFLD.

Organization of lin genes and IS6100 among different strains of hexachlorocyclohexane-degrading Sphingomonas paucimobilis: evidence for horizontal gene transfer.

The organization of lin genes and IS6100 was studied in three strains of Sphingomonas paucimobilis (B90A, Sp+, and UT26) which degraded hexachlorocyclohexane (HCH) isomers but which had been isolated at different geographical locations. DNA-DNA hybridization data revealed that most of the lin genes in these strains were associated with IS6100, an insertion sequence classified in the IS6 family and initially found in Mycobacterium fortuitum. Eleven, six, and five copies of IS6100 were detected in B90A, Sp+, and UT26, respectively. IS6100 elements in B90A were sequenced from five, one, and one regions of the genomes of B90A, Sp+, and UT26, respectively, and were found to be identical. DNA-DNA hybridization and DNA sequencing of cosmid clones also revealed that S. paucimobilis B90A contains three and two copies of linX and linA, respectively, compared to only one copy of these genes in strains Sp+ and UT26. Although the copy number and the sequence of the remaining genes of the HCH degradative pathway (linB, linC, linD, and linE) were nearly the same in all strains, there were striking differences in the organization of the linA genes as a result of replacement of portions of DNA sequences by IS6100, which gave them a strange mosaic configuration. Spontaneous deletion of linD and linE from B90A and of linA from Sp+ occurred and was associated either with deletion of a copy of IS6100 or changes in IS6100 profiles. The evidence gathered in this study, coupled with the observation that the G+C contents of the linA genes are lower than that of the remaining DNA sequence of S. paucimobilis, strongly suggests that all these strains acquired the linA gene through horizontal gene transfer mediated by IS6100. The association of IS6100 with the rest of the lin genes further suggests that IS6100 played a role in shaping the current lin gene organization.

Structure-function analyses point to a polynucleotide-accommodating groove essential for APOBEC3A restriction activities.

Members of the human APOBEC3 family of editing enzymes can inhibit various mobile genetic elements. APOBEC3A (A3A) can block the retrotransposon LINE-1 and the parvovirus adeno-associated virus type 2 (AAV-2) but does not inhibit retroviruses. In contrast, APOBEC3G (A3G) can block retroviruses but has only limited effects on AAV-2 or LINE-1. What dictates this differential target specificity remains largely undefined. Here, we modeled the structure of A3A based on its homology with the C-terminal domain of A3G and further compared the sequence of human A3A to those of 11 nonhuman primate orthologues. We then used these data to perform a mutational analysis of A3A, examining its ability to restrict LINE-1, AAV-2, and foreign plasmid DNA and to edit a single-stranded DNA substrate. The results revealed an essential functional role for the predicted single-stranded DNA-docking groove located around the A3A catalytic site. Within this region, amino acid differences between A3A and A3G are predicted to affect the shape of the polynucleotide-binding groove. Correspondingly, transferring some of these A3A residues to A3G endows the latter protein with the ability to block LINE-1 and AAV-2. These results suggest that the target specificity of APOBEC3 family members is partly defined by structural features influencing their interaction with polynucleotide substrates.

The radiation of the clownfishes has two geographical replicates

AimThe study of adaptive radiations provides an evolutionary perspective on the interactions between organisms and their environment, and is necessary to understand global biodiversity. Adaptive radiations can sometimes be replicated over several disjunct geographical entities, but most examples are found on island or in lakes. Here, we investigated the biogeographical history of the clownfishes, a clade of coral reef fish with ranges that now span most of the Indo-Pacific Ocean, in order to explore the geographical structure of an unusual adaptive radiation. LocationIndian Ocean, Indo-Australian Archipelago (IAA) and Central Pacific Ocean. MethodsWe generated DNA sequence data comprising seven nuclear markers for 27 of the 30 clownfish species. We then inferred a Bayesian phylogeny and reconstructed the biogeographical history of the group using three different methods. Finally, we applied a biogeographical model of diversification to assess whether diversification patterns differ between the Indian and Pacific Oceans. ResultsThe phylogenetic tree is highly supported and allows reconstruction of the biogeographical history of the clade. While most species arose in the IAA, one clade colonized the eastern shores of Africa and diversified there. We found that the diversification rate of clownfishes does not differ between the main radiation and the African clade. Main conclusionsThe clownfishes first appeared and diversified in the IAA. Following a colonization event, a geographically independent radiation occurred in the Indian Ocean off East Africa. This rare example of replicated adaptive radiation in the marine realm provides intriguing possibilities for further research on ecological speciation in the sea.

Point Mutations in the Monocarboxylate Transporter SLC16A12 Lead to Juvenile and Age-Related Cataract

Purpose: Previously we reported on a premature termination mutation in SLC16A12 that leads to dominant juvenile cataract and renal glucosuria. To assess the mutation rate and genotype-phenotype correlations of SLC16A12 in juvenile or age-related forms of cataract, we performed a mutation screen in cataract patients. Methods: Clinical data of approximately 660 patients were collected, genomic DNA was isolated and analyzed. Exons 3 to 8 including flanking intron sequences of SLC16A12 were PCR amplified and DNA sequence was determined. Selected mutations were tested by cell culture assays, in silico analysis and RT-PCR. Results: We found sequence alterations at a rate of approximately 1/75 patients. None of them was found in 360 control alleles. Alterations affect splice site and regulatory region but most mutations caused an amino acid substitution. The majority of the coding region mutations maps to trans-membrane domains. One mutation located to the 5'UTR. It affects translational efficiency of SLC16A12. In addition, we identified a cataract-predisposing SNP in the non-coding region that causes allele-specific splicing of the 5'UTR region. Conclusions: Altered translational efficiency of the solute carrier SLC16A12 and its allele-specific splicing strongly support a model of challenged homeostasis to cause various forms of cataract. In addition, the pathogenic property of the here reported sequence alterations is supported by the lack of known sequence variations within the coding region of SLC16A12. Due to the relatively high mutation rate, we suggest to include SLC16A12 in diagnostic cataract screening. Generally, our data recommend the assessment of regulatory sequences for diagnostic purposes.

Host-parasite relationships during a biologic invasion: 75 years postinvasion, cane toads and sympatric Australian frogs retain separate lungworm faunas.

Invasive species may carry with them parasites from their native range, differing from parasite taxa found in the invaded range. Host switching by parasites (either from the invader to native fauna or from native fauna to the invader) may have important consequences for the viability of either type of host (e.g., their survivorship, fecundity, dispersal ability, or geographic distribution). Rhabdias pseudosphaerocephala (Nematoda) is a common parasite of cane toads (Rhinella marina) in the toad's native range (South and Central America) and also in its introduced Australian range. This lungworm can depress host viability and is capable of infecting Australian frogs in laboratory trials. Despite syntopy between toads and frogs for up to 75 yr, our analyses, based on DNA sequence data of lungworms from 80 frogs and 56 toads, collected from 2008 to 2011, did not reveal any cases of host switching in nature: toads and native frogs retain entirely different lungworm faunas. All lungworms in cane toads were the South and Central American species Rhabdias pseudosphaerocephala, whereas Australian frogs contained at least four taxa (mostly undescribed and currently lumped under the name Rhabdias cf. hylae). General patterns of prevalence and intensity, based on the dissection of 1,315 frogs collected between 1989 and 2011 across the toads' Australian range, show that these Australian endemic Rhabdias spp. are widely distributed geographically and across host taxa but are more common in some frog species (especially, large-bodied species) than they are in others.

Genomic organization and gene expression in a chromosomal region of Leishmania major.

Little is known about the relation between the genome organization and gene expression in Leishmania. Bioinformatic analysis can be used to predict genes and find homologies with known proteins. A model was proposed, in which genes are organized into large clusters and transcribed from only one strand, in the form of large polycistronic primary transcripts. To verify the validity of this model, we studied gene expression at the transcriptional, post-transcriptional and translational levels in a unique locus of 34kb located on chr27 and represented by cosmid L979. Sequence analysis revealed 115 ORFs on either DNA strand. Using computer programs developed for Leishmania genes, only nine of these ORFs, localized on the same strand, were predicted to code for proteins, some of which show homologies with known proteins. Additionally, one pseudogene, was identified. We verified the biological relevance of these predictions. mRNAs from nine predicted genes and proteins from seven were detected. Nuclear run-on analyses confirmed that the top strand is transcribed by RNA polymerase II and suggested that there is no polymerase entry site. Low levels of transcription were detected in regions of the bottom strand and stable transcripts were identified for four ORFs on this strand not predicted to be protein-coding. In conclusion, the transcriptional organization of the Leishmania genome is complex, raising the possibility that computer predictions may not be comprehensive.