Cytogenetic and molecular relationships between Zarudny's rock shrew (Crocidura zarudnyi; Mammalia : Soricomorpha) and Eurasian taxa

We karyotyped and sequenced 1,140 base pairs of the mitochondrial DNA cytochrome b of a specimen of Zarudny's rock shrew (Crocidura zarudnyi) from Baluchestan, southeastern Iran, to clarify its cytogenetic and molecular relationships with other Eurasian species of Crocidura. According to the karyotype (2N = 40, FN = 50), Zarudny's rock shrew belongs to the group of the lesser white-toothed shrew (C. suaveolens), which is different from other known crocidurine karyotypes, considering the combination of the diploid and fundamental number of chromosomes. Molecular results revealed that C. zarudnyi is included in a monophyletic clade with the C. suaveolens group, where it is a sister taxon to the others (mean Kimura 2-parameter distance = 9.7%).

Molecular changes underlying mammalian phenotypic innovation

Mammals are characterized by specific phenotypic traits that include lactation, hair, and relatively large brains with unique structures. Individual mammalian lineages have, in turn, evolved characteristic traits that distinguish them from others. These include obvious anatom¬ical differences but also differences related to reproduction, life span, cognitive abilities, be¬havior. and disease susceptibility. However, the molecular basis of the diverse mammalian phenotypes and the selective pressures that shaped their evolution remain largely unknown. In the first part of my thesis, I analyzed the genetic factors associated with the origin of a unique mammalian phenotype lactation and I studied the selective pressures that forged the transition from oviparity to viviparity. Using a comparative genomics approach and evolutionary simulations, I showed that the emergence of lactation, as well as the appear¬ance of the casein gene family, significantly reduced selective pressure on the major egg-yolk proteins (the vitellogenin family). This led to a progressive loss of vitellogenins, which - in oviparous species - act as storage proteins for lipids, amino acids, phosphorous and calcium in the isolated egg. The passage to internal fertilization and placentation in therian mam¬mals rendered vitellogenins completely dispensable, which ended in the loss of the whole gene family in this lineage. As illustrated by the vitellogenin study, changes in gene content are one possible underlying factor for the evolution of mammalian-specific phenotypes. However, more subtle genomic changes, such as mutations in protein-coding sequences, can also greatly affect the phenotypes. In particular, it was proposed that changes at the level of gene reg¬ulation could underlie many (or even most) phenotypic differences between species. In the second part of my thesis, I participated in a major comparative study of mammalian tissue transcriptomes, with the goal of understanding how evolutionary forces affected expression patterns in the past 200 million years of mammalian evolution. I showed that, while com¬parisons of gene expressions are in agreement with the known species phylogeny, the rate of expression evolution varies greatly among lineages. Species with low effective population size, such as monotremes and hominoids, showed significantly accelerated rates of gene expression evolution. The most likely explanation for the high rate of gene expression evolution in these lineages is the accumulation of mildly deleterious mutations in regulatory regions, due to the low efficiency of purifying selection. Thus, our observations are in agreement with the nearly neutral theory of molecular evolution. I also describe substantial differences in evolutionary rates between tissues, with brain being the most constrained (especially in primates) and testis significantly accelerated. The rate of gene expression evolution also varies significantly between chromosomes. In particular, I observed an acceleration of gene expression changes on the X chromosome, probably as a result of adaptive processes associated with the origin of therian sex chromosomes. Lastly, I identified several individual genes as well as co-regulated expression modules that have undergone lineage specific expression changes and likely under¬lie various phenotypic innovations in mammals. The methods developed during my thesis, as well as the comprehensive gene content analyses and transcriptomics datasets made available by our group, will likely prove to be useful for further exploratory analyses of the diverse mammalian phenotypes.

Mutations in the heparan-sulfate proteoglycan glypican 6 (GPC6) impair endochondral ossification and cause recessive omodysplasia.

Glypicans are a family of glycosylphosphatidylinositol (GPI)-anchored, membrane-bound heparan sulfate (HS) proteoglycans. Their biological roles are only partly understood, although it is assumed that they modulate the activity of HS-binding growth factors. The involvement of glypicans in developmental morphogenesis and growth regulation has been highlighted by Drosophila mutants and by a human overgrowth syndrome with multiple malformations caused by glypican 3 mutations (Simpson-Golabi-Behmel syndrome). We now report that autosomal-recessive omodysplasia, a genetic condition characterized by short-limbed short stature, craniofacial dysmorphism, and variable developmental delay, maps to chromosome 13 (13q31.1-q32.2) and is caused by point mutations or by larger genomic rearrangements in glypican 6 (GPC6). All mutations cause truncation of the GPC6 protein and abolish both the HS-binding site and the GPI-bearing membrane-associated domain, and thus loss of function is predicted. Expression studies in microdissected mouse growth plate revealed expression of Gpc6 in proliferative chondrocytes. Thus, GPC6 seems to have a previously unsuspected role in endochondral ossification and skeletal growth, and its functional abrogation results in a short-limb phenotype.

A general definition and nomenclature for alternative splicing events

Understanding the molecular mechanisms responsible for the regulation of the transcriptome present in eukaryotic cells isone of the most challenging tasks in the postgenomic era. In this regard, alternative splicing (AS) is a key phenomenoncontributing to the production of different mature transcripts from the same primary RNA sequence. As a plethora ofdifferent transcript forms is available in databases, a first step to uncover the biology that drives AS is to identify thedifferent types of reflected splicing variation. In this work, we present a general definition of the AS event along with anotation system that involves the relative positions of the splice sites. This nomenclature univocally and dynamically assignsa specific ‘‘AS code’’ to every possible pattern of splicing variation. On the basis of this definition and the correspondingcodes, we have developed a computational tool (AStalavista) that automatically characterizes the complete landscape of ASevents in a given transcript annotation of a genome, thus providing a platform to investigate the transcriptome diversityacross genes, chromosomes, and species. Our analysis reveals that a substantial part—in human more than a quarter—ofthe observed splicing variations are ignored in common classification pipelines. We have used AStalavista to investigate andto compare the AS landscape of different reference annotation sets in human and in other metazoan species and found thatproportions of AS events change substantially depending on the annotation protocol, species-specific attributes, andcoding constraints acting on the transcripts. The AStalavista system therefore provides a general framework to conductspecific studies investigating the occurrence, impact, and regulation of AS.

X-chromosome tiling path array detection of copy number variants in patients with chromosome X-linked mental retardation

Background: Aproximately 5–10% of cases of mental retardation in males are due to copy number variations (CNV) on the X chromosome. Novel technologies, such as array comparative genomic hybridization (aCGH), may help to uncover cryptic rearrangements in X-linked mental retardation (XLMR) patients. We have constructed an X-chromosome tiling path array using bacterial artificial chromosomes (BACs) and validated it using samples with cytogenetically defined copy number changes. We have studied 54 patients with idiopathic mental retardation and 20 controls subjects. Results: Known genomic aberrations were reliably detected on the array and eight novel submicroscopic imbalances, likely causative for the mental retardation (MR) phenotype, were detected. Putatively pathogenic rearrangements included three deletions and five duplications (ranging between 82 kb to one Mb), all but two affecting genes previously known to be responsible for XLMR. Additionally, we describe different CNV regions with significant different frequencies in XLMR and control subjects (44% vs. 20%). Conclusion:This tiling path array of the human X chromosome has proven successful for the detection and characterization of known rearrangements and novel CNVs in XLMR patients.

Effects of supercoiling on enhancer-promoter contacts.

Using Brownian dynamics simulations, we investigate here one of possible roles of supercoiling within topological domains constituting interphase chromosomes of higher eukaryotes. We analysed how supercoiling affects the interaction between enhancers and promoters that are located in the same or in neighbouring topological domains. We show here that enhancer-promoter affinity and supercoiling act synergistically in increasing the fraction of time during which enhancer and promoter stay in contact. This stabilizing effect of supercoiling only acts on enhancers and promoters located in the same topological domain. We propose that the primary role of recently observed supercoiling of topological domains in interphase chromosomes of higher eukaryotes is to assure that enhancers contact almost exclusively their cognate promoters located in the same topological domain and avoid contacts with very similar promoters but located in neighbouring topological domains.

Chromosomal rearrangements and genetic structure at different evolutionary levels of the Sorex araneus group.

Robertsonian (Rb) fusions received large theoretical support for their role in speciation, but empirical evidence is often lacking. Here, we address the role of Rb rearrangements on the genetic differentiation of the karyotypically diversified group of shrews, Sorex araneus. We compared genetic structure between 'rearranged' and 'common' chromosomes in pairwise comparisons of five karyotypic taxa of the group. Considering all possible comparisons, we found a significantly greater differentiation at rearranged chromosomes, supporting the role of chromosomal rearrangements in the general genetic diversification of this group. Intertaxa structure and distance were larger across rearranged chromosomes for most of the comparisons, although these differences were not significant. This last result could be explained by the large variance observed among microsatellite-based estimates. The differences observed among the pairs of taxa analysed support the role of both the hybrid karyotypic complexity and the level of evolutionary divergence.

Diversity of staphylococcal cassette chromosome mec elements in predominant methicillin-resistant Staphylococcus aureus clones in a small geographic area.

Recent population genetic studies suggest that staphylococcal cassette chromosome mec (SCCmec) was acquired much more frequently than previously thought. In the present study, we aimed to investigate the diversity of SCCmec elements in a local methicillin-resistant Staphylococcus aureus (MRSA) population. Each MRSA isolate (one per patient) recovered in the Vaud canton of Switzerland from January 2005 to December 2008 was analyzed by the double-locus sequence typing (DLST) method and SCCmec typing. DLST analysis indicated that 1,884/2,036 isolates (92.5%) belong to four predominant clones. As expected from the local spread of a clone, most isolates within clones harbored an identical SCCmec type. However, three to seven SCCmec types have been recovered in every predominant DLST clone, suggesting that some of these elements might have been acquired locally. This pattern could also be explained by distinct importations of related isolates into the study region. The addition of a third highly variable locus to further increase the discriminatory power of typing as well as epidemiological data suggested that most ambiguous situations were explained by the second hypothesis. In conclusion, our study showed that even if the acquisition of new SCCmec elements at a local level likely occurs, it does not explain all the diversity observed in the study region.

The TPTE gene family: cellular expression, subcellular localization and alternative splicing.

The human TPTE (Transmembrane Phosphatase with TEnsin homology) gene family encodes a PTEN-related tyrosine phosphatase with four potential transmembrane domains. Chromosomal mapping revealed multiple copies of the TPTE gene on chromosomes 13, 15, 21, 22 and Y. Human chromosomes 13 and 21 copies encode two functional proteins, TPIP (TPTE and PTEN homologous Inositol lipid Phosphatase) and TPTE, respectively, whereas only one copy of the gene exists in the mouse genome. In the present study, we show that TPTE and TPIP proteins are expressed in secondary spermatocytes and/or prespermatids. In addition, we report the existence of several novel alternatively spliced isoforms of these two proteins with variable number of transmembrane domains. The latter has no influence on the subcellular localization of these different peptides as shown by co-immunofluorescence experiments. Finally, we identify another expressed TPTE copy, mapping to human chromosome 22, whose transcription appears to be under the control of the LTR of human endogenous retrovirus RTVL-H3.

Identification of an RP1 Prevalent Founder Mutation and Related Phenotype in Spanish Patients with Early-Onset Autosomal Recessive Retinitis.

OBJECTIVE: To identify the genetic causes underlying early-onset autosomal recessive retinitis pigmentosa (arRP) in the Spanish population and describe the associated phenotype. DESIGN: Case series. PARTICIPANTS: A total of 244 unrelated families affected by early-onset arRP. METHODS: Homozygosity mapping or exome sequencing analysis was performed in 3 families segregating arRP. A mutational screening was performed in 241 additional unrelated families for the p.Ser452Stop mutation. Haplotype analysis also was conducted. Individuals who were homozygotes, double heterozygotes, or carriers of mutations in RP1 underwent an ophthalmic evaluation to establish a genotype-phenotype correlation. MAIN OUTCOME MEASURES: DNA sequence variants, homozygous regions, haplotypes, best-corrected visual acuity, visual field assessments, electroretinogram responses, and optical coherence tomography images. RESULTS: Four novel mutations in RP1 were identified. The new mutation p.Ser542Stop was present in 11 of 244 (4.5%) of the studied families. All chromosomes harboring this mutation shared the same haplotype. All patients presented a common phenotype with an early age of onset and a prompt macular degeneration, whereas the heterozygote carriers did not show any signs of retinitis pigmentosa (RP). CONCLUSIONS: p.Ser542Stop is a single founder mutation and the most prevalent described mutation in the Spanish population. It causes early-onset RP with a rapid macular degeneration and is responsible for 4.5% of all cases. Our data suggest that the implication of RP1 in arRP may be underestimated. FINANCIAL DISCLOSURE(S): The author(s) have no proprietary or commercial interest in any materials discussed in this article.

Androgenetic alopecia: identification of four genetic risk loci and evidence for the contribution of WNT signaling to its etiology.

The pathogenesis of androgenetic alopecia (AGA, male-pattern baldness) is driven by androgens, and genetic predisposition is the major prerequisite. Candidate gene and genome-wide association studies have reported that single-nucleotide polymorphisms (SNPs) at eight different genomic loci are associated with AGA development. However, a significant fraction of the overall heritable risk still awaits identification. Furthermore, the understanding of the pathophysiology of AGA is incomplete, and each newly associated locus may provide novel insights into contributing biological pathways. The aim of this study was to identify unknown AGA risk loci by replicating SNPs at the 12 genomic loci that showed suggestive association (5 × 10(-8)<P<10(-5)) with AGA in a recent meta-analysis. We analyzed a replication set comprising 2,759 cases and 2,661 controls of European descent to confirm the association with AGA at these loci. Combined analysis of the replication and the meta-analysis data identified four genome-wide significant risk loci for AGA on chromosomes 2q35, 3q25.1, 5q33.3, and 12p12.1. The strongest association signal was obtained for rs7349332 (P=3.55 × 10(-15)) on chr2q35, which is located intronically in WNT10A. Expression studies in human hair follicle tissue suggest that WNT10A has a functional role in AGA etiology. Thus, our study provides genetic evidence supporting an involvement of WNT signaling in AGA development.

Range-Wide Sex-Chromosome Sequence Similarity Supports Occasional XY Recombination in European Tree Frogs (Hyla arborea).

In contrast with mammals and birds, most poikilothermic vertebrates feature structurally undifferentiated sex chromosomes, which may result either from frequent turnovers, or from occasional events of XY recombination. The latter mechanism was recently suggested to be responsible for sex-chromosome homomorphy in European tree frogs (Hyla arborea). However, no single case of male recombination has been identified in large-scale laboratory crosses, and populations from NW Europe consistently display sex-specific allelic frequencies with male-diagnostic alleles, suggesting the absence of recombination in their recent history. To address this apparent paradox, we extended the phylogeographic scope of investigations, by analyzing the sequences of three sex-linked markers throughout the whole species distribution. Refugial populations (southern Balkans and Adriatic coast) show a mix of X and Y alleles in haplotypic networks, and no more within-individual pairwise nucleotide differences in males than in females, testifying to recurrent XY recombination. In contrast, populations of NW Europe, which originated from a recent postglacial expansion, show a clear pattern of XY differentiation; the X and Y gametologs of the sex-linked gene Med15 present different alleles, likely fixed by drift on the front wave of expansions, and kept differentiated since. Our results support the view that sex-chromosome homomorphy in H. arborea is maintained by occasional or historical events of recombination; whether the frequency of these events indeed differs between populations remains to be clarified.

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.

Whole genome sequencing in patients with retinitis pigmentosa reveals pathogenic DNA structural changes and NEK2 as a new disease gene.

We performed whole genome sequencing in 16 unrelated patients with autosomal recessive retinitis pigmentosa (ARRP), a disease characterized by progressive retinal degeneration and caused by mutations in over 50 genes, in search of pathogenic DNA variants. Eight patients were from North America, whereas eight were Japanese, a population for which ARRP seems to have different genetic drivers. Using a specific workflow, we assessed both the coding and noncoding regions of the human genome, including the evaluation of highly polymorphic SNPs, structural and copy number variations, as well as 69 control genomes sequenced by the same procedures. We detected homozygous or compound heterozygous mutations in 7 genes associated with ARRP (USH2A, RDH12, CNGB1, EYS, PDE6B, DFNB31, and CERKL) in eight patients, three Japanese and five Americans. Fourteen of the 16 mutant alleles identified were previously unknown. Among these, there was a 2.3-kb deletion in USH2A and an inverted duplication of ∼446 kb in EYS, which would have likely escaped conventional screening techniques or exome sequencing. Moreover, in another Japanese patient, we identified a homozygous frameshift (p.L206fs), absent in more than 2,500 chromosomes from ethnically matched controls, in the ciliary gene NEK2, encoding a serine/threonine-protein kinase. Inactivation of this gene in zebrafish induced retinal photoreceptor defects that were rescued by human NEK2 mRNA. In addition to identifying a previously undescribed ARRP gene, our study highlights the importance of rare structural DNA variations in Mendelian diseases and advocates the need for screening approaches that transcend the analysis of the coding sequences of the human genome.

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.