Detecting epistasis with restricted response patterns in pairs of biallelic loci.

Well-established examples of genetic epistasis between a pair of loci typically show characteristic patterns of phenotypic distributions in joint genotype tables. However, inferring epistasis given such data is difficult due to the lack of power in commonly used approaches, which decompose the epistatic patterns into main plus interaction effects followed by testing the interaction term. Testing additive-only or all terms may have more power, but they are sensitive to nonepistatic patterns. Alternatively, the epistatic patterns of interest can be enumerated and the best matching one is found by searching through the possibilities. Although this approach requires multiple testing correction over possible patterns, each pattern can be fitted with a regression model with just one degree of freedom and thus the overall power can still be high, if the number of possible patterns is limited. Here we compare the power of the linear decomposition and pattern search methods, by applying them to simulated data generated under several patterns of joint genotype effects with simple biological interpretations. Interaction-only tests are the least powerful; while pattern search approach is the most powerful if the range of possibilities is restricted, but still includes the true pattern.

Novel testis germ cell-specific transcript of the CREB gene contains an alternatively spliced exon with multiple in-frame stop codons.

The gene encoding the cAMP-responsive transcription factor CREB consists of multiple small exons some of which undergo alternative RNA splicing. We describe the finding of a novel transcript of the CREB gene expressed at high levels in the germ cells of the rat testis. The transcript contains an alternatively spliced exon inserted within the sequence encoding the transcriptional transactivation domain of CREB and this exon contains multiple in-frame stop codons. Furthermore, the exon is conserved in both rat and human genes (75% nucleotide identity). Although the function(s) of this RNA or the truncated CREB protein predicted to result from the translation of this unusual transcript is unknown, the high level of expression in the testicular germ cells and remarkable conservation of sequences in rat and human suggests that it may have a unique biological function in these cells.

Mutation in the gene for connexin 30.3 in a family with erythrokeratodermia variabilis.

Erythrokeratodermia variabilis (EKV) is an autosomal dominant keratinization disorder characterized by migratory erythematous lesions and fixed keratotic plaques. All families with EKV show mapping to chromosome 1p34-p35, and mutations in the gene for connexin 31 (Cx31) have been reported in some but not all families. We studied eight affected and three healthy subjects in an Israeli family, of Kurdish origin, with EKV. After having mapped the disorder to chromosome 1p34-p35, we found no mutations in the genes for Cx31, Cx31.1, and Cx37. Further investigation revealed a heterozygous T-->C transition leading to the missense mutation (F137L) in the human gene for Cx30.3 that colocalizes on chromosome 1p34-p35. This nucleotide change cosegregated with the disease and was not found in 200 alleles from normal individuals. This mutation concerns a highly conserved phenylalanine, in the third transmembrane region of the Cx30.3 molecule, known to be implicated in the wall formation of the gap-junction pore. Our results show that mutations in the gene for Cx30.3 can be causally involved in EKV and point to genetic heterogeneity of this disorder. Furthermore, we suggest that our family presents a new type of EKV because of the hitherto unreported association with erythema gyratum repens.

The NAD biosynthesis inhibitor APO866 has potent antitumor activity against hematologic malignancies.

APO866 inhibits nicotinamide phosphoribosyltransferase (NMPRTase), a key enzyme involved in nicotinamide adenine dinucleotide (NAD) biosynthesis from the natural precursor nicotinamide. Intracellular NAD is essential for cell survival, and NAD depletion resulting from APO866 treatment elicits tumor cell death. Here, we determine the in vitro and in vivo sensitivities of hematologic cancer cells to APO866 using a panel of cell lines (n = 45) and primary cells (n = 32). Most cancer cells (acute myeloid leukemia [AML], acute lymphoblastic leukemia [ALL], mantle cell lymphoma [MCL], chronic lymphocytic leukemia [CLL], and T-cell lymphoma), but not normal hematopoietic progenitor cells, were sensitive to low concentrations of APO866 as measured in cytotoxicity and clonogenic assays. Treatment with APO866 decreased intracellular NAD and adenosine triphosphate (ATP) at 24 hours and 48 to72 hours, respectively. The NAD depletion led to cell death. At 96 hours, APO866-mediated cell death occurred in a caspase-independent mode, and was associated with mitochondrial dysfunction and autophagy. Further, in vivo administration of APO866 as a single agent prevented and abrogated tumor growth in animal models of human AML, lymphoblastic lymphoma, and leukemia without significant toxicity to the animals. The results support the potential of APO866 for treating hematologic malignancies.

The eukaryotic promoter database (EPD).

The Eukaryotic Promoter Database (EPD) is an annotated non-redundant collection of eukaryotic POL II promoters for which the transcription start site has been determined experimentally. Access to promoter sequences is provided by pointers to positions in nucleotide sequence entries. The annotation part of an entry includes a description of the initiation site mapping data, exhaustive cross-references to the EMBL nucleotide sequence database, SWISS-PROT, TRANSFAC and other databases, as well as bibliographic references. EPD is structured in a way that facilitates dynamic extraction of biologically meaningful promoter subsets for comparative sequence analysis. WWW-based interfaces have been developed that enable the user to view EPD entries in different formats, to select and extract promoter sequences according to a variety of criteria, and to navigate to related databases exploiting different cross-references. The EPD web site also features yearly updated base frequency matrices for major eukaryotic promoter elements. EPD can be accessed at http://www.epd.isb-sib.ch

Transmission disequilibrium of small CNVs in simplex autism.

We searched for disruptive, genic rare copy-number variants (CNVs) among 411 families affected by sporadic autism spectrum disorder (ASD) from the Simons Simplex Collection by using available exome sequence data and CoNIFER (Copy Number Inference from Exome Reads). Compared to high-density SNP microarrays, our approach yielded ∼2× more smaller genic rare CNVs. We found that affected probands inherited more CNVs than did their siblings (453 versus 394, p = 0.004; odds ratio [OR] = 1.19) and that the probands' CNVs affected more genes (921 versus 726, p = 0.02; OR = 1.30). These smaller CNVs (median size 18 kb) were transmitted preferentially from the mother (136 maternal versus 100 paternal, p = 0.02), although this bias occurred irrespective of affected status. The excess burden of inherited CNVs among probands was driven primarily by sibling pairs with discordant social-behavior phenotypes (p < 0.0002, measured by Social Responsiveness Scale [SRS] score), which contrasts with families where the phenotypes were more closely matched or less extreme (p > 0.5). Finally, we found enrichment of brain-expressed genes unique to probands, especially in the SRS-discordant group (p = 0.0035). In a combined model, our inherited CNVs, de novo CNVs, and de novo single-nucleotide variants all independently contributed to the risk of autism (p < 0.05). Taken together, these results suggest that small transmitted rare CNVs play a role in the etiology of simplex autism. Importantly, the small size of these variants aids in the identification of specific genes as additional risk factors associated with ASD.

Thirty new loci for age at menarche identified by a meta-analysis of genome-wide association studies.

To identify loci for age at menarche, we performed a meta-analysis of 32 genome-wide association studies in 87,802 women of European descent, with replication in up to 14,731 women. In addition to the known loci at LIN28B (P = 5.4 × 10⁻⁶⁰) and 9q31.2 (P = 2.2 × 10⁻³³), we identified 30 new menarche loci (all P < 5 × 10⁻⁸) and found suggestive evidence for a further 10 loci (P < 1.9 × 10⁻⁶). The new loci included four previously associated with body mass index (in or near FTO, SEC16B, TRA2B and TMEM18), three in or near other genes implicated in energy homeostasis (BSX, CRTC1 and MCHR2) and three in or near genes implicated in hormonal regulation (INHBA, PCSK2 and RXRG). Ingenuity and gene-set enrichment pathway analyses identified coenzyme A and fatty acid biosynthesis as biological processes related to menarche timing.

Genetic architecture of ambulatory blood pressure in the general population: insights from cardiovascular gene-centric array.

Genetic determinants of blood pressure are poorly defined. We undertook a large-scale, gene-centric analysis to identify loci and pathways associated with ambulatory systolic and diastolic blood pressure. We measured 24-hour ambulatory blood pressure in 2020 individuals from 520 white European nuclear families (the Genetic Regulation of Arterial Pressure of Humans in the Community Study) and genotyped their DNA using the Illumina HumanCVD BeadChip array, which contains ≈50 000 single nucleotide polymorphisms in >2000 cardiovascular candidate loci. We found a strong association between rs13306560 polymorphism in the promoter region of MTHFR and CLCN6 and mean 24-hour diastolic blood pressure; each minor allele copy of rs13306560 was associated with 2.6 mm Hg lower mean 24-hour diastolic blood pressure (P=1.2×10(-8)). rs13306560 was also associated with clinic diastolic blood pressure in a combined analysis of 8129 subjects from the Genetic Regulation of Arterial Pressure of Humans in the Community Study, the CoLaus Study, and the Silesian Cardiovascular Study (P=5.4×10(-6)). Additional analysis of associations between variants in gene ontology-defined pathways and mean 24-hour blood pressure in the Genetic Regulation of Arterial Pressure of Humans in the Community Study showed that cell survival control signaling cascades could play a role in blood pressure regulation. There was also a significant overrepresentation of rare variants (minor allele frequency: <0.05) among polymorphisms showing at least nominal association with mean 24-hour blood pressure indicating that a considerable proportion of its heritability may be explained by uncommon alleles. Through a large-scale gene-centric analysis of ambulatory blood pressure, we identified an association of a novel variant at the MTHFR/CLNC6 locus with diastolic blood pressure and provided new insights into the genetic architecture of blood pressure.

Innate immunogenetics: a tool for exploring new frontiers of host defence.

The discovery of innate immune genes, such as those encoding Toll-like receptors (TLRs), nucleotide-binding oligomerisation domain-like receptors (NLRs), and related signal-transducing molecules, has led to a substantial improvement of our understanding of innate immunity. Recent immunogenetic studies have associated polymorphisms of the genes encoding TLRs, NLRs, and key signal-transducing molecules, such as interleukin-1 receptor-associated kinase 4 (IRAK4), with increased susceptibility to, or outcome of, infectious diseases. With the availability of high-throughput genotyping techniques, it is becoming increasingly evident that analyses of genetic polymorphisms of innate immune genes will further improve our knowledge of the host antimicrobial defence response and help in identifying individuals who are at increased risk of life-threatening infections. This is likely to open new perspectives for the development of diagnostic, predictive, and preventive management strategies to combat infectious diseases.

trEST, trGEN and Hits: access to databases of predicted protein sequences.

High throughput genome (HTG) and expressed sequence tag (EST) sequences are currently the most abundant nucleotide sequence classes in the public database. The large volume, high degree of fragmentation and lack of gene structure annotations prevent efficient and effective searches of HTG and EST data for protein sequence homologies by standard search methods. Here, we briefly describe three newly developed resources that should make discovery of interesting genes in these sequence classes easier in the future, especially to biologists not having access to a powerful local bioinformatics environment. trEST and trGEN are regularly regenerated databases of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Hits is a web-based data retrieval and analysis system providing access to precomputed matches between protein sequences (including sequences from trEST and trGEN) and patterns and profiles from Prosite and Pfam. The three resources can be accessed via the Hits home page (http://hits. isb-sib.ch).

Pro-opiomelanocortin gene and melanin-based colour polymorphism in a reptile

Colour polymorphism is widespread among vertebrates and plays important roles in prey-predator interactions, thermoregulation, social competition, and sexual selection. However, the genetic mechanisms involved in colour variation have been studied mainly in domestic mammals and birds, whereas information on wild animals remains scarce. Interestingly, the pro-opiomelanocortin gene (POMC) gives rise to melanocortin hormones that trigger melanogenesis (by binding the melanocortin-1-receptor; Mc1r) and other physiological and behavioural functions (by binding the melanocortin receptors Mc1-5rs). Owing to its pleiotropic effect, the POMC gene could therefore account for the numerous covariations between pigmentation and other phenotypic traits. We screened the POMC and Mc1r genes in 107 wild asp vipers (Vipera aspis) that can exhibit four discrete colour morphs (two unpatterned morphs: concolor or melanistic; two patterned morphs: blotched or lined) in a single population. Our study revealed a correlation between a single nucleotide polymorphism situated within the 3-untranslated region of the POMC gene and colour variation, whereas Mc1r was not found to be polymorphic. To the best of our knowledge, we disclose for the first time a relationship between a mutation at the POMC gene and coloration in a wild animal, as well as a correlation between a genetic marker and coloration in a snake species. Interestingly, similar mutations within the POMC 3-untranslated region are linked to human obesity and alcohol and drug dependence. Combined with our results, this suggests that the 3-untranslated region of the POMC gene may play a role in its regulation in distant vertebrates.

Genome-wide association and genetic functional studies identify autism susceptibility candidate 2 gene (AUTS2) in the regulation of alcohol consumption.

Alcohol consumption is a moderately heritable trait, but the genetic basis in humans is largely unknown, despite its clinical and societal importance. We report a genome-wide association study meta-analysis of ∼2.5 million directly genotyped or imputed SNPs with alcohol consumption (gram per day per kilogram body weight) among 12 population-based samples of European ancestry, comprising 26,316 individuals, with replication genotyping in an additional 21,185 individuals. SNP rs6943555 in autism susceptibility candidate 2 gene (AUTS2) was associated with alcohol consumption at genome-wide significance (P = 4 × 10(-8) to P = 4 × 10(-9)). We found a genotype-specific expression of AUTS2 in 96 human prefrontal cortex samples (P = 0.026) and significant (P < 0.017) differences in expression of AUTS2 in whole-brain extracts of mice selected for differences in voluntary alcohol consumption. Down-regulation of an AUTS2 homolog caused reduced alcohol sensitivity in Drosophila (P < 0.001). Our finding of a regulator of alcohol consumption adds knowledge to our understanding of genetic mechanisms influencing alcohol drinking behavior.

A genome-wide association study of early menopause and the combined impact of identified variants.

Early menopause (EM) affects up to 10% of the female population, reducing reproductive lifespan considerably. Currently, it constitutes the leading cause of infertility in the western world, affecting mainly those women who postpone their first pregnancy beyond the age of 30 years. The genetic aetiology of EM is largely unknown in the majority of cases. We have undertaken a meta-analysis of genome-wide association studies (GWASs) in 3493 EM cases and 13 598 controls from 10 independent studies. No novel genetic variants were discovered, but the 17 variants previously associated with normal age at natural menopause as a quantitative trait (QT) were also associated with EM and primary ovarian insufficiency (POI). Thus, EM has a genetic aetiology which overlaps variation in normal age at menopause and is at least partly explained by the additive effects of the same polygenic variants. The combined effect of the common variants captured by the single nucleotide polymorphism arrays was estimated to account for ∼30% of the variance in EM. The association between the combined 17 variants and the risk of EM was greater than the best validated non-genetic risk factor, smoking.

Update of the BIPM comparison BIPM.RI(II)-K1.Cs-134 of activity measurements of the radionuclide 134Cs to include the 2008 results of the BEV (Austria), the 2009 result of the IRA (Switzerland) and the 2010 results of the NMISA (South Africa)

Purine nucleotide pyrophosphotransferase was purified to apparent homogeneity from a culture filtrate of Streptomyces morookaensis. It is a monomeric protein with a molecular weight of 24 000-25 000, and its isoelectric point is 6.9. The enzyme synthesizes purine nucleoside 5'-phosphate (mono, di, or tri) 3'-diphosphates such as pppApp, ppApp, pApp, pppGpp, ppGpp and pppIpp by transferring a pyrophosphoryl group from the 5'-position of ATP, dATP and ppApp to the 3'-position of purine nucleotides. The purified enzyme catalysed the formation of 435 mumol of pppApp and 620 mumol of pppGpp from ATP and GTP per min mg protein under the standard conditions. The enzyme requires absolutely a divalent cation for activity, and optimum pH for the enzyme activity lay above 10 for Mg2+, for Co2+ and Zn2+ from 9 to 9.5, and for Fe2+ from 7.5 to 8. The following Michaelis constants were determined: AMP, 2.78 mM; ADP, 3.23 mM; GMP, 0.89 mM; GDP, 0.46 mM and GTP, 1.54 mM, in the case of ATP donor. The enzyme is inhibited by guanine, guanosine, dGDP, dGTP, N-bromosuccinimide, iodacetate, sodium borate and mercuric acetate.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.