Genometrics as an essential tool for the assembly of whole genome sequences: the example of the chromosome of Bifidobacterium longum NCC2705.

BACKGROUND: Analysis of the first reported complete genome sequence of Bifidobacterium longum NCC2705, an actinobacterium colonizing the gastrointestinal tract, uncovered its proteomic relatedness to Streptomyces coelicolor and Mycobacterium tuberculosis. However, a rapid scrutiny by genometric methods revealed a genome organization totally different from all so far sequenced high-GC Gram-positive chromosomes. RESULTS: Generally, the cumulative GC- and ORF orientation skew curves of prokaryotic genomes consist of two linear segments of opposite slope: the minimum and the maximum of the curves correspond to the origin and the terminus of chromosome replication, respectively. However, analyses of the B. longum NCC2705 chromosome yielded six, instead of two, linear segments, while its dnaA locus, usually associated with the origin of replication, was not located at the minimum of the curves. Furthermore, the coorientation of gene transcription with replication was very low. Comparison with closely related actinobacteria strongly suggested that the chromosome of B. longum was misassembled, and the identification of two pairs of relatively long homologous DNA sequences offers the possibility for an alternative genome assembly proposed here below. By genometric criteria, this configuration displays all of the characters common to bacteria, in particular to related high-GC Gram-positives. In addition, it is compatible with the partially sequenced genome of DJO10A B. longum strain. Recently, a corrected sequence of B. longum NCC2705, with a configuration similar to the one proposed here below, has been deposited in GenBank, confirming our predictions. CONCLUSION: Genometric analyses, in conjunction with standard bioinformatic tools and knowledge of bacterial chromosome architecture, represent fast and straightforward methods for the evaluation of chromosome assembly.

Meta-analysis of gene-level associations for rare variants based on single-variant statistics.

Meta-analysis of genome-wide association studies (GWASs) has led to the discoveries of many common variants associated with complex human diseases. There is a growing recognition that identifying "causal" rare variants also requires large-scale meta-analysis. The fact that association tests with rare variants are performed at the gene level rather than at the variant level poses unprecedented challenges in the meta-analysis. First, different studies may adopt different gene-level tests, so the results are not compatible. Second, gene-level tests require multivariate statistics (i.e., components of the test statistic and their covariance matrix), which are difficult to obtain. To overcome these challenges, we propose to perform gene-level tests for rare variants by combining the results of single-variant analysis (i.e., p values of association tests and effect estimates) from participating studies. This simple strategy is possible because of an insight that multivariate statistics can be recovered from single-variant statistics, together with the correlation matrix of the single-variant test statistics, which can be estimated from one of the participating studies or from a publicly available database. We show both theoretically and numerically that the proposed meta-analysis approach provides accurate control of the type I error and is as powerful as joint analysis of individual participant data. This approach accommodates any disease phenotype and any study design and produces all commonly used gene-level tests. An application to the GWAS summary results of the Genetic Investigation of ANthropometric Traits (GIANT) consortium reveals rare and low-frequency variants associated with human height. The relevant software is freely available.

Geographic variation in sex-chromosome differentiation in the common frog (Rana temporaria).

In sharp contrast with birds and mammals, sex-determination systems in ectothermic vertebrates are often highly dynamic and sometimes multifactorial. Both environmental and genetic effects have been documented in common frogs (Rana temporaria). One genetic linkage group, mapping to the largest pair of chromosomes and harbouring the candidate sex-determining gene Dmrt1, associates with sex in several populations throughout Europe, but association varies both within and among populations. Here, we show that sex association at this linkage group differs among populations along a 1500-km transect across Sweden. Genetic differentiation between sexes is strongest (FST  = 0.152) in a northern-boreal population, where male-specific alleles and heterozygote excesses (FIS  = -0.418 in males, +0.025 in females) testify to a male-heterogametic system and lack of X-Y recombination. In the southernmost population (nemoral climate), in contrast, sexes share the same alleles at the same frequencies (FST  = 0.007 between sexes), suggesting unrestricted recombination. Other populations show intermediate levels of sex differentiation, with males falling in two categories: some cluster with females, while others display male-specific Y haplotypes. This polymorphism may result from differences between populations in the patterns of X-Y recombination, co-option of an alternative sex-chromosome pair, or a mixed sex-determination system where maleness is controlled either by genes or by environment depending on populations or families. We propose approaches to test among these alternative models, to disentangle the effects of climate and phylogeography on the latitudinal trend, and to sort out how this polymorphism relates to the 'sexual races' described in common frogs in the 1930s.

Post-axial polydactyly type A2, overgrowth and autistic traits associated with a chromosome 13q31.3 microduplication encompassing miR-17-92 and GPC5.

Genomic rearrangements at chromosome 13q31.3q32.1 have been associated with digital anomalies, dysmorphic features, and variable degree of mental disability. Microdeletions leading to haploinsufficiency of miR17∼92, a cluster of micro RNA genes closely linked to GPC5 in both mouse and human genomes, has recently been associated with digital anomalies in the Feingold like syndrome. Here, we report on a boy with familial dominant post-axial polydactyly (PAP) type A, overgrowth, significant facial dysmorphisms and autistic traits who carries the smallest germline microduplication known so far in that region. The microduplication encompasses the whole miR17∼92 cluster and the first 5 exons of GPC5. This report supports the newly recognized role of miR17∼92 gene dosage in digital developmental anomalies, and suggests a possible role of GPC5 in growth regulation and in cognitive development.

ParA2, a Vibrio cholerae chromosome partitioning protein, forms left-handed helical filaments on DNA.

Most bacterial chromosomes contain homologs of plasmid partitioning (par) loci. These loci encode ATPases called ParA that are thought to contribute to the mechanical force required for chromosome and plasmid segregation. In Vibrio cholerae, the chromosome II (chrII) par locus is essential for chrII segregation. Here, we found that purified ParA2 had ATPase activities comparable to other ParA homologs, but, unlike many other ParA homologs, did not form high molecular weight complexes in the presence of ATP alone. Instead, formation of high molecular weight ParA2 polymers required DNA. Electron microscopy and three-dimensional reconstruction revealed that ParA2 formed bipolar helical filaments on double-stranded DNA in a sequence-independent manner. These filaments had a distinct change in pitch when ParA2 was polymerized in the presence of ATP versus in the absence of a nucleotide cofactor. Fitting a crystal structure of a ParA protein into our filament reconstruction showed how a dimer of ParA2 binds the DNA. The filaments formed with ATP are left-handed, but surprisingly these filaments exert no topological changes on the right-handed B-DNA to which they are bound. The stoichiometry of binding is one dimer for every eight base pairs, and this determines the geometry of the ParA2 filaments with 4.4 dimers per 120 A pitch left-handed turn. Our findings will be critical for understanding how ParA proteins function in plasmid and chromosome segregation.

Limited clinical benefit of minority K103N and Y181C-variant detection in addition to routine genotypic resistance testing in antiretroviral therapy-naive patients.

OBJECTIVE: The presence of minority nonnucleoside reverse transcriptase inhibitor (NNRTI)-resistant HIV-1 variants prior to antiretroviral therapy (ART) has been linked to virologic failure in treatment-naive patients. DESIGN: We performed a large retrospective study to determine the number of treatment failures that could have been prevented by implementing minority drug-resistant HIV-1 variant analyses in ART-naïve patients in whom no NNRTI resistance mutations were detected by routine resistance testing. METHODS: Of 1608 patients in the Swiss HIV Cohort Study, who have initiated first-line ART with two nucleoside reverse transcriptase inhibitors (NRTIs) and one NNRTI before July 2008, 519 patients were eligible by means of HIV-1 subtype, viral load and sample availability. Key NNRTI drug resistance mutations K103N and Y181C were measured by allele-specific PCR in 208 of 519 randomly chosen patients. RESULTS: Minority K103N and Y181C drug resistance mutations were detected in five out of 190 (2.6%) and 10 out of 201 (5%) patients, respectively. Focusing on 183 patients for whom virologic success or failure could be examined, virologic failure occurred in seven out of 183 (3.8%) patients; minority K103N and/or Y181C variants were present prior to ART initiation in only two of those patients. The NNRTI-containing, first-line ART was effective in 10 patients with preexisting minority NNRTI-resistant HIV-1 variant. CONCLUSION: As revealed in settings of case-control studies, minority NNRTI-resistant HIV-1 variants can have an impact on ART. However, the implementation of minority NNRTI-resistant HIV-1 variant analysis in addition to genotypic resistance testing (GRT) cannot be recommended in routine clinical settings. Additional associated risk factors need to be discovered.

Structural variation-associated expression changes are paralleled by chromatin architecture modifications.

Copy number variants (CNVs) influence the expression of genes that map not only within the rearrangement, but also to its flanks. To assess the possible mechanism(s) underlying this "neighboring effect", we compared intrachromosomal interactions and histone modifications in cell lines of patients affected by genomic disorders and control individuals. Using chromosome conformation capture (4C-seq), we observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together. The newly identified interacting genes include AUTS2, mutations of which are associated with autism and intellectual disabilities. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We also pinpointed concomitant changes in histone modifications between samples. We conclude that large genomic rearrangements can lead to chromatin conformation changes that extend far away from the structural variant, thereby possibly modulating expression globally and modifying the phenotype. GEO SERIES ACCESSION NUMBER: GSE33784, GSE33867.

Phylogeographical structure, postglacial recolonization and barriers to gene flow in the distinctive Valais chromosome race of the common shrew (Sorex araneus).

Using one male-inherited and eight biparentally inherited microsatellite markers, we investigate the population genetic structure of the Valais chromosome race of the common shrew (Sorex araneus) in the Central Alps of Europe. Unexpectedly, the Y-chromosome microsatellite suggests nearly complete absence of male gene flow among populations from the St-Bernard and Simplon regions (Switzerland). Autosomal markers also show significant genetic structuring among these two geographical areas. Isolation by distance is significant and possible barriers to gene flow exist in the study area. Two different approaches are used to better understand the geographical patterns and the causes of this structuring. Using a principal component analysis for which testing procedure exists, and partial Mantel tests, we show that the St-Bernard pass does not represent a significant barrier to gene flow although it culminates at 2469 m, close to the highest altitudinal record for this species. Similar results are found for the Simplon pass, indicating that both passes represented potential postglacial recolonization routes into Switzerland from Italian refugia after the last Pleistocene glaciations. In contrast with the weak effect of these mountain passes, the Rhône valley lowlands significantly reduce gene flow in this species. Natural obstacles (the large Rhône river) and unsuitable habitats (dry slopes) are both present in the valley. Moreover, anthropogenic changes to landscape structures are likely to have strongly reduced available habitats for this shrew in the lowlands, thereby promoting genetic differentiation of populations found on opposite sides of the Rhône valley.

Identification of a new cancer/testis gene family, CT47, among expressed multicopy genes on the human X chromosome.

Cancer/testis (CT) genes are normally expressed in germ cells only, yet are reactivated and expressed in some tumors. Of the approximately 40 CT genes or gene families identified to date, 20 are on the X chromosome and are present as multigene families, many with highly conserved members. This indicates that novel CT gene families may be identified by detecting duplicated expressed genes on chromosome X. By searching for transcript clusters that map to multiple locations on the chromosome, followed by in silico analysis of their gene expression profiles, we identified five novel gene families with testis-specific expression and >98% sequence identity among family members. The expression of these genes in normal tissues and various tumor cell lines and specimens was evaluated by qualitative and quantitative RT-PCR, and a novel CT gene family with at least 13 copies was identified on Xq24, designated as CT47. mRNA expression of CT47 was found mainly in the testes, with weak expression in the placenta. Brain tissue was the only positive somatic tissue tested, with an estimated CT47 transcript level 0.09% of that found in testis. Among the tumor specimens tested, CT47 expression was found in approximately 15% of lung cancer and esophageal cancer specimens, but not in colorectal cancer or breast cancer. The putative CT47 protein consists of 288 amino acid residues, with a C-terminus rich in alanine and glutamic acid. The only species other than human in which a gene homologous to CT47 has been detected is the chimpanzee, with the predicted protein showing approximately 80% identity in its carboxy terminal region.

Infantile spasms is associated with deletion of the MAGI2 gene on chromosome 7q11.23-q21.11.

Infantile spasms (IS) is the most severe and common form of epilepsy occurring in the first year of life. At least half of IS cases are idiopathic in origin, with others presumed to arise because of brain insult or malformation. Here, we identify a locus for IS by high-resolution mapping of 7q11.23-q21.1 interstitial deletions in patients. The breakpoints delineate a 500 kb interval within the MAGI2 gene (1.4 Mb in size) that is hemizygously disrupted in 15 of 16 participants with IS or childhood epilepsy, but remains intact in 11 of 12 participants with no seizure history. MAGI2 encodes the synaptic scaffolding protein membrane-associated guanylate kinase inverted-2 that interacts with Stargazin, a protein also associated with epilepsy in the stargazer mouse.

Mutations leading to X-linked hypohidrotic ectodermal dysplasia affect three major functional domains in the tumor necrosis factor family member ectodysplasin-A.

Mutations in the epithelial morphogen ectodysplasin-A (EDA), a member of the tumor necrosis factor (TNF) family, are responsible for the human disorder X-linked hypohidrotic ectodermal dysplasia (XLHED) characterized by impaired development of hair, eccrine sweat glands, and teeth. EDA-A1 and EDA-A2 are two splice variants of EDA, which bind distinct EDA-A1 and X-linked EDA-A2 receptors. We identified a series of novel EDA mutations in families with XLHED, allowing the identification of the following three functionally important regions in EDA: a C-terminal TNF homology domain, a collagen domain, and a furin protease recognition sequence. Mutations in the TNF homology domain impair binding of both splice variants to their receptors. Mutations in the collagen domain can inhibit multimerization of the TNF homology region, whereas those in the consensus furin recognition sequence prevent proteolytic cleavage of EDA. Finally, a mutation affecting an intron splice donor site is predicted to eliminate specifically the EDA-A1 but not the EDA-A2 splice variant. Thus a proteolytically processed, oligomeric form of EDA-A1 is required in vivo for proper morphogenesis.

Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila.

Ionotropic glutamate receptors (iGluRs) mediate neuronal communication at synapses throughout vertebrate and invertebrate nervous systems. We have characterized a family of iGluR-related genes in Drosophila, which we name ionotropic receptors (IRs). These receptors do not belong to the well-described kainate, AMPA, or NMDA classes of iGluRs, and they have divergent ligand-binding domains that lack their characteristic glutamate-interacting residues. IRs are expressed in a combinatorial fashion in sensory neurons that respond to many distinct odors but do not express either insect odorant receptors (ORs) or gustatory receptors (GRs). IR proteins accumulate in sensory dendrites and not at synapses. Misexpression of IRs in different olfactory neurons is sufficient to confer ectopic odor responsiveness. Together, these results lead us to propose that the IRs comprise a novel family of chemosensory receptors. Conservation of IR/iGluR-related proteins in bacteria, plants, and animals suggests that this receptor family represents an evolutionarily ancient mechanism for sensing both internal and external chemical cues.