Strategies for genetic model specification in the screening of genome-wide meta-analysis signals for further replication

Background Meta-analysis is increasingly being employed as a screening procedure in large-scale association studies to select promising variants for follow-up studies. However, standard methods for meta-analysis require the assumption of an underlying genetic model, which is typically unknown a priori. This drawback can introduce model misspecifications, causing power to be suboptimal, or the evaluation of multiple genetic models, which augments the number of false-positive associations, ultimately leading to waste of resources with fruitless replication studies. We used simulated meta-analyses of large genetic association studies to investigate naive strategies of genetic model specification to optimize screenings of genome-wide meta-analysis signals for further replication. Methods Different methods, meta-analytical models and strategies were compared in terms of power and type-I error. Simulations were carried out for a binary trait in a wide range of true genetic models, genome-wide thresholds, minor allele frequencies (MAFs), odds ratios and between-study heterogeneity (tau(2)). Results Among the investigated strategies, a simple Bonferroni-corrected approach that fits both multiplicative and recessive models was found to be optimal in most examined scenarios, reducing the likelihood of false discoveries and enhancing power in scenarios with small MAFs either in the presence or in absence of heterogeneity. Nonetheless, this strategy is sensitive to tau(2) whenever the susceptibility allele is common (MAF epsilon 30%), resulting in an increased number of false-positive associations compared with an analysis that considers only the multiplicative model. Conclusion Invoking a simple Bonferroni adjustment and testing for both multiplicative and recessive models is fast and an optimal strategy in large meta-analysis-based screenings. However, care must be taken when examined variants are common, where specification of a multiplicative model alone may be preferable.

Phylogenetic Analysis of Hepatitis B Virus Genotype F Complete Genome Sequences From Chilean Patients With Chronic Infection

Molecular epidemiological data concerning the hepatitis B virus (HBV) in Chile are not known completely. Since the HBV genotype F is the most prevalent in the country, the goal of this study was to obtain full HBV genome sequences from patients infected chronically in order to determine their subgenotypes and the occurrence of resistance-associated mutations. Twenty-one serum samples from antiviral drug-naive patients with chronic hepatitis B were subjected to full-length PCR amplification, and both strands of the whole genomes were fully sequenced. Phylogenetic analyses were performed along with reference sequences available from GenBank (n = 290). The sequences were aligned using Clustal X and edited in the SE-AL software. Bayesian phylogenetic analyses were conducted by Markov Chain Monte Carlo simulations (MCMC) for 10 million generations in order to obtain the substitution tree using BEAST. The sequences were also analyzed for the presence of primary drug resistance mutations using CodonCode Aligner Software. The phylogenetic analyses indicated that all sequences were found to be the HBV subgenotype F1b, clustered into four different groups, suggesting that diverse lineages of this subgenotype may be circulating within this population of Chilean patients. J. Med. Virol. 83: 1530-1536, 2011. (C) 2011 Wiley-Liss, Inc.

In familial hyperaldosteronism type I, hybrid gene-induced aldosterone production dominates that induced by wild-type genes

We compared the aldosterone-producing potency of the angiotensin II-sensitive wild-type aldosterone synthase genes and the ACTH-sensitive hybrid 11 beta-hydroxylase/aldosterone synthase gene by examining aldosterone, PRA, and cortisol day-curves (2-hourly levels over 24 h) in patients with familial hyperaldosteronism type I, before and during long-term (0.8-13.5 yr) glucocorticoid treatment. In 8 untreated patients, PRA levels were usually suppressed, and aldosterone correlated strongly with cortisol (r = 0.69-0.99). Fourteen studies were performed on 10 patients receiving glucocorticoid treatment that corrected hypertension, hypokalemia, and PRA suppression in all. ACTH was markedly and continuously suppressed in 6 studies, 3 of which demonstrated strong correlations between aldosterone and PRA (r = 0.77-0.92), ACTH was only partially suppressed in the remaining 8 studies; aldosterone correlated strongly: 1) with cortisol alone in 5 (r = 0.71-0.98); 2) with cortisol (r = 0.90) and PRA (r = 0.74) in one; 3) with PRA only in one (r = 0.80); and 4) with neither PRA nor cortisol in one. Unless ACTH is markedly and continuously suppressed, aldosterone is more responsive to ACTH than to renin/angiotensin II, despite the latter being unsuppressed. This is consistent with the hybrid gene being more powerfully expressed than the wild-type aldosterone synthase genes in familial hyperaldosteronism type I.

Whole-Genome Characterization of Human and Simian Immunodeficiency Virus Intrahost Diversity by Ultradeep Pyrosequencing

Rapid evolution and high intrahost sequence diversity are hallmarks of human and simian immunodeficiency virus (HIV/SIV) infection. Minor viral variants have important implications for drug resistance, receptor tropism, and immune evasion. Here, we used ultradeep pyrosequencing to sequence complete HIV/SIV genomes, detecting variants present at a frequency as low as 1%. This approach provides a more complete characterization of the viral population than is possible with conventional methods, revealing low-level drug resistance and detecting previously hidden changes in the viral population. While this work applies pyrosequencing to immunodeficiency viruses, this approach could be applied to virtually any viral pathogen.

High frequency of submicroscopic chromosomal imbalances in patients with syndromic craniosynostosis detected by a combined approach of microsatellite segregation analysis, multiplex ligation-dependent probe amplification and array-based comparative genome hybridisation

We present the first comprehensive study, to our knowledge, on genomic chromosomal analysis in syndromic craniosynostosis. In total, 45 patients with craniosynostotic disorders were screened with a variety of methods including conventional karyotype, microsatellite segregation analysis, subtelomeric multiplex ligation-dependent probe amplification) and whole-genome array-based comparative genome hybridisation. Causative abnormalities were present in 42.2% (19/45) of the samples, and 27.8% (10/36) of the patients with normal conventional karyotype carried submicroscopic imbalances. Our results include a wide variety of imbalances and point to novel chromosomal regions associated with craniosynostosis. The high incidence of pure duplications or trisomies suggests that these are important mechanisms in craniosynostosis, particularly in cases involving the metopic suture.

Distinct patterns of somatic alterations in a lymphoblastoid and a tumor genome derived from the same individual

Although patterns of somatic alterations have been reported for tumor genomes, little is known on how they compare with alterations present in non-tumor genomes. A comparison of the two would be crucial to better characterize the genetic alterations driving tumorigenesis. We sequenced the genomes of a lymphoblastoid (HCC1954BL) and a breast tumor (HCC1954) cell line derived from the same patient and compared the somatic alterations present in both. The lymphoblastoid genome presents a comparable number and similar spectrum of nucleotide substitutions to that found in the tumor genome. However, a significant difference in the ratio of non-synonymous to synonymous substitutions was observed between both genomes (P = 0.031). Protein-protein interaction analysis revealed that mutations in the tumor genome preferentially affect hub-genes (P = 0.0017) and are co-selected to present synergistic functions (P < 0.0001). KEGG analysis showed that in the tumor genome most mutated genes were organized into signaling pathways related to tumorigenesis. No such organization or synergy was observed in the lymphoblastoid genome. Our results indicate that endogenous mutagens and replication errors can generate the overall number of mutations required to drive tumorigenesis and that it is the combination rather than the frequency of mutations that is crucial to complete tumorigenic transformation.

The genome of the blood fluke Schistosoma mansoni

Schistosoma mansoni is responsible for the neglected tropical disease schistosomiasis that affects 210 million people in 76 countries. Here we present analysis of the 363 megabase nuclear genome of the blood fluke. It encodes at least 11,809 genes, with an unusual intron size distribution, and new families of micro-exon genes that undergo frequent alternative splicing. As the first sequenced flatworm, and a representative of the Lophotrochozoa, it offers insights into early events in the evolution of the animals, including the development of a body pattern with bilateral symmetry, and the development of tissues into organs. Our analysis has been informed by the need to find new drug targets. The deficits in lipid metabolism that make schistosomes dependent on the host are revealed, and the identification of membrane receptors, ion channels and more than 300 proteases provide new insights into the biology of the life cycle and new targets. Bioinformatics approaches have identified metabolic chokepoints, and a chemogenomic screen has pinpointed schistosome proteins for which existing drugs may be active. The information generated provides an invaluable resource for the research community to develop much needed new control tools for the treatment and eradication of this important and neglected disease.

Occurrence of TRGV-BJ hybrid gene in SV40-transformed fibroblast cell lines

Illegitimate V(D)J-recombination in lymphoid malignancies involves rearrangements in immunoglobulin or T-cell receptor genes, and these rearrangements may play a role in oncogenic events. High frequencies of TRGV-BJ hybrid gene (rearrangement between the TRB and TRG loci at 7q35 and 7p14-15, respectively) have been detected in lymphocytes from patients with ataxia telangiectasia (AT), and also in patients with lymphoid malignancies. Although the TRGV-BJ gene has been described only in T-lymphocytes, we previously detected the presence of TRGV-BJ hybrid gene in the genomic DNA extracted from SV40-transformed AT5BIVA fibroblasts from an AT patient. Aiming to determine whether the AT phenotype or the SV40 transformation could be responsible for the production of the hybrid gene by illegitimate V(D)J-recombination, DNA samples were extracted from primary and SV40-transformed (normal and AT) cell lines, following Nested-PCR with TRGV- and TRBJ-specific primers. The hybrid gene was only detected in SV40-transformed fibroblasts (AT-5BIVA and MRC-5). Sequence alignment of the cloned PCR products using the BLAST program confirmed that the fragments corresponded to the TRGV-BJ hybrid gene. The present results indicate that the rearrangement can be produced in nonlymphoid cells, probably as a consequence of the genomic instability caused by the SV40-transformation, and independently of ATM gene mutation.

The Genome Sequence of Taurine Cattle: A Window to Ruminant Biology and Evolution

To understand the biology and evolution of ruminants, the cattle genome was sequenced to about sevenfold coverage. The cattle genome contains a minimum of 22,000 genes, with a core set of 14,345 orthologs shared among seven mammalian species of which 1217 are absent or undetected in noneutherian (marsupial or monotreme) genomes. Cattle-specific evolutionary breakpoint regions in chromosomes have a higher density of segmental duplications, enrichment of repetitive elements, and species-specific variations in genes associated with lactation and immune responsiveness. Genes involved in metabolism are generally highly conserved, although five metabolic genes are deleted or extensively diverged from their human orthologs. The cattle genome sequence thus provides a resource for understanding mammalian evolution and accelerating livestock genetic improvement for milk and meat production.

Complete genome analysis of a rabies virus isolate from Brazilian wild fox

The complete genome sequence of wild-type rabies virus (RABV) isolated from a wild Brazilian hoary fox (Dusicyon sp.), the BR-Pfx1 isolate, was determined and compared with fixed RABV strains. The genome structure and organization of the BR-Pfx1 isolate were composed of 11,924 nt and included the five standard genes of rhabdoviruses. Sequences of mRNA start and stop signals for transcription were highly conserved among all structural protein genes of the BR-Pfx1 isolate. All amino acid residues in the glycoprotein (G) gene associated with pathogenicity were retained in the BR-Pfx1 isolate, while unique amino acid substitutions were found in antigenic region I of the nucleoprotein gene and III of G. These results suggest that although the standard genome structure and organization of the RABV isolate are common between the BR-Pfx1 isolate and fixed RABV strains, the unique amino acid substitutions in functional sites of the BR-Pfx1 isolate may result in different biological characteristics from fixed RABV strains.

Micro-hardness evaluation of a micro-hybrid composite resin light cured with halogen light, light-emitting diode and argon ion laser

This in vitro study aimed to determine whether the micro-hardness of a composite resin is modified by the light units or by the thickness of the increment. Composite resin disks were divided into 15 groups (n = 5), according to the factors under study: composite resin thickness (0 mm, 1 mm, 2 mm , 3 mm and 4 mm) and light units. The light activation was performed with halogen light (HL) (40 s, 500 mW/cm(2)), argon ion laser (AL) (30 s, 600 mW/cm(2)) or light-emitting diode (LED) (30 s, 400 mW/cm(2)). Vickers micro-hardness tests were performed after 1 week and were carried out on the top surface (0 mm-control) and at different depths of the samples. Analysis of variance (ANOVA) and Tukey tests (P a parts per thousand currency signaEuro parts per thousand 0.05) revealed no statistically significant difference among the light units for the groups of 0 mm and 1 mm thickness. At 2 mm depth, the AL was not statistically different from the HL, but the latter showed higher micro-hardness values than the LED. In groups with 3 mm and 4 mm thickness, the HL also showed higher micro-hardness values than the groups activated by the AL and the LED. Only the HL presented satisfactory polymerization with 3 mm of thickness. With a 4 mm increment no light unit was able to promote satisfactory polymerization.

Effects of Nd:YAG Laser Irradiation on the Hybrid Layer of Different Adhesive Systems

Purpose: The aim of this in vitro study was to evaluate the microtensile bond strength (mu TBS) and hybrid layer morphology of different adhesive systems, either followed by treatment with Nd:YAG laser irradiation or not. Previous studies have shown the effects of Nd:YAG laser irradiation on the dentin surface at restoration margins, but there are few reports about the significance of the irradiation on the hybrid layer. Materials and Methods: The flattened coronal and root dentin samples of 24 bovine teeth were randomly divided into 8 groups, according to the adhesive system used - Scotchbond Multi Purpose (SBMP) or Clearfil SE Bond (CSEB) - and were either irradiated with Nd:YAG or not, with different parameters: 0.8 W/10 Hz, 0.8 W/20 Hz, 1.2 W/10 Hz, 1.2 W/20 Hz. The left sides of specimens were the control groups, and right sides were irradiated. A composite crown was built over bonded surfaces and stored in water (24 h at 37 degrees C). Specimens were sectioned vertically into slabs that were subjected to mu TBS testing and observed by SEM. Results: Control groups (27.81 +/- 1.38) showed statistically higher values than lased groups (21.37 +/- 0.99), and CSEB control group values (31.26 +/- 15.71) were statistically higher than those of SBMP (24.3 +/- 10.66). There were no significant differences between CSEB (20.34 +/- 10.01) and SBMP (22.43 +/- 9.82) lased groups. Among parameters tested, 0.8 W/10 Hz showed the highest value (25.54 +/- 11.74). Nd:YAG laser irradiation caused dentin to melt under the adhesive layer of both adhesive systems tested. Conclusion: With the parameters used in this study, Nd:YAG laser irradiation of the hybrid layer promoted morphological changes in dentin and negatively influenced the bond strength of both adhesive systems.

Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera)

The complete arrangement of genes in the mitochondrial (mt) genome is known for 12 species of insects, and part of the gene arrangement in the mt genome is known for over 300 other species of insects. The arrangement of genes in the mt genome is very conserved in insects studied, since all of the protein-coding and rRNA genes and most of the tRNA genes are arranged in the same way. We sequenced the entire mt genome of the wallaby louse, Heterodoxus macropus, which is 14,670 bp long and has the 37 genes typical of animals and some noncoding regions. The largest noncoding region is 73 bp long (93% A+T), and the second largest is 47 bp long (92% AST). Both of these noncoding regions seem to be able to form stem-loop structures. The arrangement of genes in the mt genome of this louse is unlike that of any other animal studied. All tRNA genes have moved and/or inverted relative to the ancestral gene arrangement of insects, which is present in the fruit fly Drosophila yakuba. At least nine protein-coding genes (atp6, atp8, cox2, cob, nad1-nad3, nad5, and nad6) have moved; moreover, four of these genes (atp6, atp8, nad1, and nad3) have inverted. The large number of gene rearrangements in the mt genome of H. macropus is unprecedented for an arthropod.

Familial varieties of primary aldosteronism

1. Improved approaches to screening and diagnosis have revealed primary aldosteronism (PAL) to be much more common than previously thought, with most patients normokalaemic. The spectrum of this disorder has been further broadened by the study of familial varieties. 2. Familial hyperaldosteronism type I (FH-I) is a glucocorticoid-remediable form of PAL caused by the inheritance of an adrenocorticotrophic hormone (ACTH)- regulated, hybrid CYP11B1/CYP11B2 gene. Diagnosis has been greatly facilitated by the advent of genetic testing. The severity of hypertension varies widely in FH-I, even among members of the same family, and has demonstrated relationships with gender, degree of biochemical disturbance and hybrid gene crossover point position. Hormone day curve studies show that the hybrid gene dominates over wild-type CYP11B2 in terms of aldosterone regulation. This may be due, in part, to a defect in wild-type CYP11B2-induced aldosterone production. Control of hypertension in FH-I requires only partial suppression of ACTH and much smaller glucocorticoid doses than previously recommended. 3. Familial hyperaldosteronism type II (FH-II) is not glucocorticoid remediable and is not associated with the hybrid gene mutation. Familial hyperaldosteronism type II is clinically, biochemically and morphologically indistinguishable from apparently non-familial PAL. Linkage studies in one informative family did not show segregation of FH-II with the CYP11B2, AT1 or MEN1 genes, but a genome-wide search has revealed linkage with a locus in chromosome 7. As has already occurred in FH-I, elucidation of causative mutations is likely to facilitate earlier detection of PAL.