Interpretation of array comparative genome hybridization data: a major challenge.

The advent and application of high-resolution array-based comparative genome hybridization (array CGH) has led to the detection of large numbers of copy number variants (CNVs) in patients with developmental delay and/or multiple congenital anomalies as well as in healthy individuals. The notion that CNVs are also abundantly present in the normal population challenges the interpretation of the clinical significance of detected CNVs in patients. In this review we will illustrate a general clinical workflow based on our own experience that can be used in routine diagnostics for the interpretation of CNVs.

16p11.2 600 kb Duplications confer risk for typical and atypical Rolandic epilepsy.

Rolandic epilepsy (RE) is the most common idiopathic focal childhood epilepsy. Its molecular basis is largely unknown and a complex genetic etiology is assumed in the majority of affected individuals. The present study tested whether six large recurrent copy number variants at 1q21, 15q11.2, 15q13.3, 16p11.2, 16p13.11 and 22q11.2 previously associated with neurodevelopmental disorders also increase risk of RE. Our association analyses revealed a significant excess of the 600 kb genomic duplication at the 16p11.2 locus (chr16: 29.5-30.1 Mb) in 393 unrelated patients with typical (n = 339) and atypical (ARE; n = 54) RE compared with the prevalence in 65,046 European population controls (5/393 cases versus 32/65,046 controls; Fisher's exact test P = 2.83 × 10(-6), odds ratio = 26.2, 95% confidence interval: 7.9-68.2). In contrast, the 16p11.2 duplication was not detected in 1738 European epilepsy patients with either temporal lobe epilepsy (n = 330) and genetic generalized epilepsies (n = 1408), suggesting a selective enrichment of the 16p11.2 duplication in idiopathic focal childhood epilepsies (Fisher's exact test P = 2.1 × 10(-4)). In a subsequent screen among children carrying the 16p11.2 600 kb rearrangement we identified three patients with RE-spectrum epilepsies in 117 duplication carriers (2.6%) but none in 202 carriers of the reciprocal deletion. Our results suggest that the 16p11.2 duplication represents a significant genetic risk factor for typical and atypical RE.

Marked hemiatrophy in carriers of Duchenne muscular dystrophy.

OBJECTIVE: To describe the clinical and molecular genetic findings in 2 carriers of Duchenne muscular dystrophy (DMD) who exhibited marked hemiatrophy. Duchenne muscular dystrophy is an X-linked disorder in which affected male patients harbor mutations in the dystrophin gene. Female patients with heterozygous mutations may be manifesting carriers. DESIGN: Case study. SETTING: Neurology clinic. PATIENTS: Two manifesting carriers of DMD. INTERVENTIONS: Clinical and radiologic examinations along with histologic and molecular investigations. RESULTS: Both patients had marked right-sided hemiatrophy on examination with radiologic evidence of muscle atrophy and fatty replacement on the affected side. In each case, histologic analysis revealed a reduction in dystrophin staining on the right side. Genetic analysis of the dystrophin gene revealed a tandem exonic duplication in patient 1 and a multiexonic deletion in patient 2 with no further point mutations identified on the other chromosome. CONCLUSIONS: Marked hemiatrophy can occur in DMD manifesting carriers. This is likely to result from a combination of skewed X-inactivation and somatic mosaicism.

The serine repeat antigen (SERA) gene family phylogeny in Plasmodium: the impact of GC content and reconciliation of gene and species trees.

Plasmodium falciparum is the parasite responsible for the most acute form of malaria in humans. Recently, the serine repeat antigen (SERA) in P. falciparum has attracted attention as a potential vaccine and drug target, and it has been shown to be a member of a large gene family. To clarify the relationships among the numerous P. falciparum SERAs and to identify orthologs to SERA5 and SERA6 in Plasmodium species affecting rodents, gene trees were inferred from nucleotide and amino acid sequence data for 33 putative SERA homologs in seven different species. (A distance method for nucleotide sequences that is specifically designed to accommodate differing GC content yielded results that were largely compatible with the amino acid tree. Standard-distance and maximum-likelihood methods for nucleotide sequences, on the other hand, yielded gene trees that differed in important respects.) To infer the pattern of duplication, speciation, and gene loss events in the SERA gene family history, the resulting gene trees were then "reconciled" with two competing Plasmodium species tree topologies that have been identified by previous phylogenetic studies. Parsimony of reconciliation was used as a criterion for selecting a gene tree/species tree pair and provided (1) support for one of the two species trees and for the core topology of the amino acid-derived gene tree, (2) a basis for critiquing fine detail in a poorly resolved region of the gene tree, (3) a set of predicted "missing genes" in some species, (4) clarification of the relationship among the P. falciparum SERA, and (5) some information about SERA5 and SERA6 orthologs in the rodent malaria parasites. Parsimony of reconciliation and a second criterion--implied mutational pattern at two key active sites in the SERA proteins-were also seen to be useful supplements to standard "bootstrap" analysis for inferred topologies.

The amphioxus genome and the evolution of the chordate karyotype.

Lancelets ('amphioxus') are the modern survivors of an ancient chordate lineage, with a fossil record dating back to the Cambrian period. Here we describe the structure and gene content of the highly polymorphic approximately 520-megabase genome of the Florida lancelet Branchiostoma floridae, and analyse it in the context of chordate evolution. Whole-genome comparisons illuminate the murky relationships among the three chordate groups (tunicates, lancelets and vertebrates), and allow not only reconstruction of the gene complement of the last common chordate ancestor but also partial reconstruction of its genomic organization, as well as a description of two genome-wide duplications and subsequent reorganizations in the vertebrate lineage. These genome-scale events shaped the vertebrate genome and provided additional genetic variation for exploitation during vertebrate evolution.

Transcriptional regulation of RNA polymerase III in mammalian cells

L'ARN Polymérase III (Pol III) transcrit un ensemble de petits ARN non traduits impliqués dans des processus cellulaires tels que la biosynthèse des protéines, la maturation des ARNs ou le contrôle transcriptionnel. De ce fait, la Pol III joue un rôle important dans la régulation de la croissance et la prolifération cellulaire. L'initiation de la transcription par la Pol III nécessite l'interaction entre des facteurs de transcription et le complexe de la Pol III lui-même. Un sous- complexe de la Pol III, composé de 3 sous-unités, HsRPC3, HsRPC6 et HsRPC7 sert d'intermédiaire dans cette interaction. Dans cette étude, nous avons caractérisé une nouvelle sous-unité de la Pol III, HsRPC7-Like, homologue à HsRPC7. Nous avons montré que ces deux homologues se trouvent spécifiquement chez les vertébrés. Ils proviennent d'un ancêtre commun qui, après duplication il y a 600 millions d'années, a donné naissance à ces deux paralogues. Dans les cellules humaines, deux formes de Pol III coexistent : l'une contientt HsRPC7, l'autre HsRPC7-Like. Nous avons localisé, à l'échelle du génome entier, la présence de ces deux formes de Pol III dans des cellules humaines et dans le foie de souris. Les deux sous-unités ont démontré des caractéristiques identiques, suggérant qu'elles possèdent des fonctions similaires. Cependant, nous avons analysé les motifs d'expression des gènes codant pour RPC7 et RPC7-Like dans des lignées cellulaires dans des conditions variées telles que la concentration de sérum et la densité cellulaire, ainsi que les motifs d'expression dans le foie de souris et des cellules d'hépatocarcinome de souris. Nos résultats suggèrent que l'expression de ces deux sous-untiés varie en fonction de l'activité de prolifération de la cellule. - RNA polymerase III (Pol III) transcribes a set of genes coding for short untranslated RNAs involved in essential cellular processes as for example protein biosynthesis, RNA maturation, and transcriptional control. Thereby Pol III plays an important role in regulating cell growth and proliferation. Initiation of Pol III transcription requires interactions between transcription factors and the Pol III core complex. A Pol III sub-complex composed of three subunits, HsRPC3, HsRPC6, and HsRPC7 mediates this interaction. In this study, we have characterized a new Pol III subunit, HsRPC7-Like, an homologue of HsRPC7. We have shown that these two homologues are specific to vertebrates and originate from an ancestor gene that duplicated 600 mio years ago to give birth to two paralogues. In human cells, two forms of Pol III coexist, one containing HsRPC7 and the other HsRPC7-Like. We have localized, genome-wide, these two Pol III forms in human cells and mouse liver. Both subunits were found on all types of Pol III genes, suggesting that they share similar function. However, we analysed the expression patterns of the RPC7 and RPC7-Like coding genes under various conditions of serum concentration and cell density in different cell lines, as well as expression patterns in mouse liver and mouse hepatocarcinoma cells. Our results suggest that the expression of these two subunits varies with the proliferation rate of the cell.

Mutations in CNNM4 cause recessive cone-rod dystrophy with amelogenesis imperfecta.

Cone-rod dystrophies are inherited dystrophies of the retina characterized by the accumulation of deposits mainly localized to the cone-rich macular region of the eye. Dystrophy can be limited to the retina or be part of a syndrome. Unlike nonsyndromic cone-rod dystrophies, syndromic cone-rod dystrophies are genetically heterogeneous with mutations in genes encoding structural, cell-adhesion, and transporter proteins. Using a genome-wide single-nucleotide polymorphism (SNP) haplotype analysis to fine map the locus and a gene-candidate approach, we identified homozygous mutations in the ancient conserved domain protein 4 gene (CNNM4) that either generate a truncated protein or occur in highly conserved regions of the protein. Given that CNNM4 is implicated in metal ion transport, cone-rod dystrophy and amelogenesis imperfecta may originate from abnormal ion homeostasis.

The Intra-Alpine terrain: A Paleotethyan remnant in the Alpine Variscides

A Cordilleran type evolution is proposed for the Variscan orogen of middle Europe. This orogenesis is regarded as mainly evolving through terrain accretion and subsequent collapse of the overthickened crust. A major terrain accretion took place between late Devonian and early Carboniferous when the Intra-Alpine terrain collided with the Ligerian-Moldanubian active margin. This terrain is regarded as being a segment of the northern margin of Paleotethys. Oblique subduction of Paleotethys under the newly accreted terrain is responsible for the voluminous calc-alkaline magmatism in late Carboniferous. The Paleotethys subduction has generated a lateral displacement of the eastern part of the Intra-Alpine terrain inducing a duplication of its western end. The late Carboniferous closure of Paleotethys in middle Europe is not found eastward where this closure happened only in early-Triassic times, following the simultaneous opening of the Neotethys ocean and the Meliata back-arc. Palinspastic models of the western Tethyan realm are proposed from the Carboniferous to early Jurassic.

Copy number variation and chromatin structure

The functional consequences of structural variation in the human genome range from adaptation, to phenotypic variation, to predisposition to diseases. Copy number variation (CNV) was shown to influence the phenotype by modifying, in a somewhat dose-dependent manner, the expression of genes that map within them, as well as that of genes located on their flanks. To assess the possible mechanism(s) behind this neighboring effect, we compared histone modification status of cell lines from patients affected by Williams-Beuren, Williams-Beuren region duplication, Smith-Magenis or DiGeorge Syndrome and control individuals using a high-throughput version of chromatin immuno-precipitation method (ChIP), called ChlP-seq. We monitored monomethylation of lysine K20 on histone H4 and trimethylation of lysine K27 on histone H3, as proxies for open and condensed chromatin, respectively. Consistent with the changes in expression levels observed for multiple genes mapping on the entire length of chromosomes affected by structural variants, we also detected regions with modified histone status between samples, up- and downstream from the critical regions, up to the end of the rearranged chromosome. We also gauged the intrachromosomal interactions of these cell lines utilizing chromosome conformation capture (4C-seq) technique. We observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together, possibly forming an interacting cluster with each other and the WBSCR. 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 conclude, that large genomic rearrangements can lead to changes in the state of the chromatin spreading far away from the critical region, thus possibly affecting expression globally and as a result modifying the phenotype of the patients. - Les conséquences fonctionnelles des variations structurelles dans le génome humain sont vastes, allant de l'adaptation, en passant par les variations phénotypiques, aux prédispositions à certaines maladies. Il a été démontré que les variations du nombre de copies (CNV) influencent le phénotype en modifiant, d'une manière plus ou moins dose-dépendante, l'expression des gènes se situant à l'intérieur de ces régions, mais également celle des gènes se trouvant dans les régions flanquantes. Afin d'étudier les mécanismes possibles sous-jacents à cet effet de voisinage, nous avons comparé les états de modification des histones dans des lignées cellulaires dérivées de patients atteints du syndrome de Williams-Beuren, de la duplication de la région Williams-Beuren, du syndrome de Smith-Magenis ou du syndrome de Di- George et d'individus contrôles en utilisant une version haut-débit de la méthode d'immunoprécipitation de la chromatine (ChIP), appelée ChIP-seq. Nous avons suivi la mono-méthylation de la lysine K20 sur l'histone H4 et la tri-méthylation de la lysine K27 sur l'histone H3, marqueurs respectifs de la chromatine ouverte et fermée. En accord avec les changements de niveaux d'expression observés pour de multiples gènes tout le long des chromosomes affectés par les CNVs, nous avons aussi détecté des régions présentant des modifications d'histones entre les échantillons, situées de part et d'autre des régions critiques, jusqu'aux extrémités du chromosome réarrangé. Nous avons aussi évalué les interactions intra-chromosomiques ayant lieu dans ces cellules par l'utilisation de la technique de capture de conformation des chromosomes (4C-seq). Nous avons observé qu'un groupe de gènes flanquants la région critique du syndrome de Williams-Beuren (WBSCR) forment souvent une boucle, constituant un groupe d'interactions privilégiées entre ces gènes et la WBSCR. La délétion de la WBSCR perturbe l'expression de ce groupe de gènes flanquants, mais également les interactions à grande échelle entre eux et la région réarrangée. Nous en concluons que les larges réarrangements génomiques peuvent aboutir à des changements de l'état de la chromatine pouvant s'étendre bien plus loin que la région critique, affectant donc potentiellement l'expression de manière globale et ainsi modifiant le phénotype des patients.

How confident can we be that orthologs are similar, but paralogs differ?

Homologous genes are classified into orthologs and paralogs, depending on whether they arose by speciation or duplication. It is widely assumed that orthologs share similar functions, whereas paralogs are expected to diverge more from each other. But does this assumption hold up on further examination? We present evidence that orthologs and paralogs are not so different in either their evolutionary rates or their mechanisms of divergence. We emphasize the importance of appropriately designed studies to test models of gene evolution between orthologs and between paralogs. Thus, functional change between orthologs might be as common as between paralogs, and future studies should be designed to test the impact of duplication against this alternative model.

Birth and expression evolution of mammalian microRNA genes.

MicroRNAs (miRNAs) are major post-transcriptional regulators of gene expression, yet their origins and functional evolution in mammals remain little understood due to the lack of appropriate comparative data. Using RNA sequencing, we have generated extensive and comparable miRNA data for five organs in six species that represent all main mammalian lineages and birds (the evolutionary outgroup) with the aim to unravel the evolution of mammalian miRNAs. Our analyses reveal an overall expansion of miRNA repertoires in mammals, with threefold accelerated birth rates of miRNA families in placentals and marsupials, facilitated by the de novo emergence of miRNAs in host gene introns. Generally, our analyses suggest a high rate of miRNA family turnover in mammals with many newly emerged miRNA families being lost soon after their formation. Selectively preserved mammalian miRNA families gradually evolved higher expression levels, as well as altered mature sequences and target gene repertoires, and were apparently mainly recruited to exert regulatory functions in nervous tissues. However, miRNAs that originated on the X chromosome evolved high expression levels and potentially diverse functions during spermatogenesis, including meiosis, through selectively driven duplication-divergence processes. Overall, our study thus provides detailed insights into the birth and evolution of mammalian miRNA genes and the associated selective forces.

Comparative genomics of ParaHox clusters of teleost fishes: gene cluster breakup and the retention of gene sets following whole genome duplications.

BACKGROUND: The evolutionary lineage leading to the teleost fish underwent a whole genome duplication termed FSGD or 3R in addition to two prior genome duplications that took place earlier during vertebrate evolution (termed 1R and 2R). Resulting from the FSGD, additional copies of genes are present in fish, compared to tetrapods whose lineage did not experience the 3R genome duplication. Interestingly, we find that ParaHox genes do not differ in number in extant teleost fishes despite their additional genome duplication from the genomic situation in mammals, but they are distributed over twice as many paralogous regions in fish genomes. RESULTS: We determined the DNA sequence of the entire ParaHox C1 paralogon in the East African cichlid fish Astatotilapia burtoni, and compared it to orthologous regions in other vertebrate genomes as well as to the paralogous vertebrate ParaHox D paralogons. Evolutionary relationships among genes from these four chromosomal regions were studied with several phylogenetic algorithms. We provide evidence that the genes of the ParaHox C paralogous cluster are duplicated in teleosts, just as it had been shown previously for the D paralogon genes. Overall, however, synteny and cluster integrity seems to be less conserved in ParaHox gene clusters than in Hox gene clusters. Comparative analyses of non-coding sequences uncovered conserved, possibly co-regulatory elements, which are likely to contain promoter motives of the genes belonging to the ParaHox paralogons. CONCLUSION: There seems to be strong stabilizing selection for gene order as well as gene orientation in the ParaHox C paralogon, since with a few exceptions, only the lengths of the introns and intergenic regions differ between the distantly related species examined. The high degree of evolutionary conservation of this gene cluster's architecture in particular - but possibly clusters of genes more generally - might be linked to the presence of promoter, enhancer or inhibitor motifs that serve to regulate more than just one gene. Therefore, deletions, inversions or relocations of individual genes could destroy the regulation of the clustered genes in this region. The existence of such a regulation network might explain the evolutionary conservation of gene order and orientation over the course of hundreds of millions of years of vertebrate evolution. Another possible explanation for the highly conserved gene order might be the existence of a regulator not located immediately next to its corresponding gene but further away since a relocation or inversion would possibly interrupt this interaction. Different ParaHox clusters were found to have experienced differential gene loss in teleosts. Yet the complete set of these homeobox genes was maintained, albeit distributed over almost twice the number of chromosomes. Selection due to dosage effects and/or stoichiometric disturbance might act more strongly to maintain a modal number of homeobox genes (and possibly transcription factors more generally) per genome, yet permit the accumulation of other (non regulatory) genes associated with these homeobox gene clusters.

Evolutionary histories of expanded peptidase families in Schistosoma mansoni

Schistosoma mansoni is one of the three main causative agents of human schistosomiasis, a major health problem with a vast socio-economic impact. Recent advances in the proteomic analysis of schistosomes have revealed that peptidases are the main virulence factors involved in the pathogenesis of this disease. In this context, evolutionary studies can be applied to identify peptidase families that have been expanded in genomes over time in response to different selection pressures. Using a phylogenomic approach, we searched for expanded endopeptidase families in the S. mansoni predicted proteome with the aim of contributing to the knowledge of such enzymes as potential therapeutic targets. We found three endopeptidase families that comprise leishmanolysins (metallopeptidase M8 family), cercarial elastases (serine peptidase S1 family) and cathepsin D proteins (aspartic peptidase A1 family). Our results suggest that the Schistosoma members of these families originated from successive gene duplication events in the parasite lineage after its diversification from other metazoans. Overall, critical residues are conserved among the duplicated genes/proteins. Furthermore, each protein family displays a distinct evolutionary history. Altogether, this work provides an evolutionary view of three S. mansoni peptidase families, which allows for a deeper understanding of the genomic complexity and lineage-specific adaptations potentially related to the parasitic lifestyle.

Role of the HCF-1 basic region in sustaining cell proliferation.

BACKGROUND: The human herpes simplex virus-associated host cell factor 1 (HCF-1) is a conserved human transcriptional co-regulator that links positive and negative histone modifying activities with sequence-specific DNA-binding transcription factors. It is synthesized as a 2035 amino acid precursor that is cleaved to generate an amino- (HCF-1(N)) terminal subunit, which promotes G1-to-S phase progression, and a carboxy- (HCF-1(C)) terminal subunit, which controls multiple aspects of cell division during M phase. The HCF-1(N) subunit contains a Kelch domain that tethers HCF-1 to sequence-specific DNA-binding transcription factors, and a poorly characterized so called "Basic" region (owing to a high ratio of basic vs. acidic amino acids) that is required for cell proliferation and has been shown to associate with the Sin3 histone deacetylase (HDAC) component. Here we studied the role of the Basic region in cell proliferation and G1-to-S phase transition assays. METHODOLOGY/PRINCIPAL FINDINGS: Surprisingly, much like the transcriptional activation domains of sequence-specific DNA-binding transcription factors, there is no unique sequence within the Basic region required for promoting cell proliferation or G1-to-S phase transition. Indeed, the ability to promote these activities is size dependent such that the shorter the Basic region segment the less activity observed. We find, however, that the Basic region requirements for promoting cell proliferation in a temperature-sensitive tsBN67 cell assay are more stringent than for G1-to-S phase progression in an HCF-1 siRNA-depletion HeLa-cell assay. Thus, either half of the Basic region alone can support G1-to-S phase progression but not cell proliferation effectively in these assays. Nevertheless, the Basic region displays considerable structural plasticity because each half is able to promote cell proliferation when duplicated in tandem. Consistent with a potential role in promoting cell-cycle progression, the Sin3a HDAC component can associate independently with either half of the Basic region fused to the HCF-1 Kelch domain. CONCLUSIONS/SIGNIFICANCE: While conserved, the HCF-1 Basic region displays striking structural flexibility for controlling cell proliferation.