Visualization and quality assessment of de novo genome assemblies.

Recent technological progress has greatly facilitated de novo genome sequencing. However, de novo assemblies consist in many pieces of contiguous sequence (contigs) arranged in thousands of scaffolds instead of small numbers of chromosomes. Confirming and improving the quality of such assemblies is critical for subsequent analysis. We present a method to evaluate genome scaffolding by aligning independently obtained transcriptome sequences to the genome and visually summarizing the alignments using the Cytoscape software. Applying this method to the genome of the red fire ant Solenopsis invicta allowed us to identify inconsistencies in 7%, confirm contig order in 20% and extend 16% of scaffolds.Scripts that generate tables for visualization in Cytoscape from FASTA sequence and scaffolding information files are publicly available at https://github.com/ksanao/TGNet.

Proximal 11p deletion syndrome (P11pDS): additional evaluation of the clinical and molecular aspects.

The combination of multiple exostoses (EXT) and enlarged parietal foramina (foramina parietalia permagna, FPP) represent the main features of the proximal 11p deletion syndrome (P11pDS), a contiguous gene syndrome (MIM 601224) caused by an interstitial deletion on the short arm of chromosome 11. Here we present clinical aspects of two new P11pDS patients and the clinical follow-up of one patient reported in the original paper describing this syndrome. Recognised clinical signs include EXT, FPP, mental retardation, facial asymmetry, asymmetric calcification of coronary sutures, defective vision (severe myopia, nystagmus, strabismus), skeletal anomalies (small hands and feet, tapering fingers), heart defect, and anal stenosis. In addition fluorescence in situ hybridisation and molecular analysis were performed to gain further insight in potential candidate genes involved in P11pDS.

Contribution à la cytotaxonomie des Soricidés (Mammalia, Insectivora) de l'Afrique occidentale.

Chromosomes analysis of 28 shrews belonging to the genus Crocidura from West Africa allowed the detection of five new karyotypes: C. wimmeri (2n = 50, NF = 84), C. cf. gracilipes (56/86), C. cf. nimbae (46/68), and C. cf. planiceps (44/72). Among the individuals identified morphologically as C. poensis (52/70) a distinct species was recognized on the basis of its chromosome complement, C. nigeriae (50/76). The karyotype of C. occidentalis was confirmed with the study of another subspecies, C. occidentalis manni (50/66). An identical chromosome set characterizes C. odorata giffardi and provides evidence for its relationship with C. occidentalis.

Characterization of a new clinical yeast species, Candida tunisiensis sp. nov., isolated from a strain collection from Tunisian hospitals.

From a collection of yeast isolates isolated from patients in Tunisian hospitals between September 2006 and July 2010, the yeast strain JEY63 (CBS 12513), isolated from a 50-year-old male that suffered from oral thrush, could not be identified to the species level using conventional methods used in clinical laboratories. These methods include matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS), germ tube formation, and the use of CHROMagar Candida and metabolic galleries. Sequence analysis of the nuclear rRNA (18S rRNA, 5.8S rRNA, and 26S rRNA) and internal transcribed spacer regions (ITS1 and ITS2) indicated that the ribosomal DNA sequences of this species were not yet reported. Multiple gene phylogenic analyses suggested that this isolate clustered at the base of the Dipodascaceae (Saccharomycetales, Saccharomycetes, and Ascomycota). JEY63 was named Candida tunisiensis sp. nov. according to several phenotypic criteria and its geographical origin. C. tunisiensis was able to grow at 42°C and does not form chlamydospores and hyphae but could grow as yeast and pseudohyphal forms. C. tunisiensis exhibited most probably a haploid genome with an estimated size of 10 Mb on at least three chromosomes. Using European Committee for Antimicrobial Susceptibility Testing (EUCAST) and Clinical and Laboratory Standards Institute (CLSI) Candida albicans susceptibility breakpoints as a reference, C. tunisiensis was resistant to fluconazole (MIC = 8 μg/ml), voriconazole (MIC = 0.5 μg/ml), itraconazole (MIC = 16 μg/ml), and amphotericin B (MIC = 4 μg/ml) but still susceptible to posaconazole (MIC = 0.008 μg/ml) and caspofungin (MIC = 0.5 μg/ml). In conclusion, MALDI-TOF MS permitted the early selection of an unusual isolate, which was still unreported in molecular databases but could not be unambiguously classified based on phylogenetic approaches.

Amyloid and non-amyloid forms of 5q31-linked corneal dystrophy resulting from kerato-epithelin mutations at Arg-124 are associated with abnormal turnover of the protein.

Mutations in kerato-epithelin are responsible for a group of hereditary cornea-specific deposition diseases, 5q31-linked corneal dystrophies. These conditions are characterized by progressive accumulation of protein deposits of different ultrastructure. Herein, we studied the corneas with mutations at kerato-epithelin residue Arg-124 resulting in amyloid (R124C), non-amyloid (R124L), and a mixed pattern of deposition (R124H). We found that aggregated kerato-epithelin comprised all types of pathological deposits. Each mutation was associated with characteristic changes of protein turnover in corneal tissue. Amyloidogenesis in R124C corneas was accompanied by the accumulation of N-terminal kerato-epithelin fragments, whereby species of 44 kDa were the major constituents of amyloid fibrils. R124H corneas with prevailing non-amyloid inclusions showed accumulation of a new 66-kDa species altogether with the full-size 68-kDa form. Finally, in R124L cornea with non amyloid deposits, we found only the accumulation of the 68-kDa form. Two-dimensional gels revealed mutation-specific changes in the processing of the full-size protein in all affected corneas. It appears that substitutions at the same residue (Arg-124) result in cornea-specific deposition of kerato-epithelin via distinct aggregation pathways each involving altered turnover of the protein in corneal tissue.

Functional analysis of human POT1 and RPA : insights into the role of single stranded DNA binding proteins at telomeres

Abstract Telomeres, the natural ends of chromosomes, need to be protected from chromosome end fusions, aberrant homologous recombination and degradation. In humans, chromosome ends are specified through arrays of tandemly repeated 5'-TTAGGG-3' hexamers, ending in a 3' overhang. A complex formed by the six proteins TRF1, TRF2, hRap1, TIN2, TPP1 and POT1 specifically assocìates with and protects telomeres. Telomeres are maintained by semiconservative DNA replication and by a specialized reverse transcriptase, telomerase, that carries an RNA subunit which templates new telomeric repeat synthesis. The telomeric single stranded (ss) DNA binding protein POT1 protects the telomeric 3' overhang and modulates telomerase-mediated telomere elongation. It is possible that POT1 also influences DNA synthesis during semiconservative DNA replication, which is initiated by the DNA polymerase alpha-primase complex. The heterotrimeric ss DNA-binding protein RPA plays essential roles during DNA replication. RPA binds to ss DNA with high affinity in order to stabilize ss DNA and facilitate nascent strand synthesis at the replication fork. Here we investigate how the two proteins RPA and POT1 contribute to telomere maintenance by regulating semi-conservative DNA replication and telomerase. Using chromatin immunoprecipitation experiments, we show that RPA associates with telomeres during S-phase. Analysis of telomere structure in cells shRNA-depleted for RPA and POT1 reveals that loss of RPA and POT1 causes exposure of single-stranded DNA at telomeres, suggestive of incomplete DNA replication. Biochemical experiments using purified recombinant POT1 and RPA show that saturating telomeric oligonucleotides with POT1 or RPA reduces the primase activity of the DNA polymerase alpha-primase complex and the overall activity of telomerase. POT1 and RPA also increase the primer extension by DNA polymerase alpha-primase complex and the processivity of telomerase under certain conditions, although POT1 increases the activities to a greater extent than RPA. We propose that POT1 is required for proper replication of the lagging strand of telomeres and that some phenotypes observed in POT1-depleted cells may stern from incomplete DNA replication rather than de-protection of the single-stranded overhang. Résumé Les télomères, les extrémités normales des chromosomes linéaires, doivent être protégés des fusions chromosomiques, d'événements de recombinaison homologue aberrants et de phénomènes de dégradation. Chez l'Homme, les extrémités des chromosomes sont constitués d'ADN double brin répétitif de séquence 5'-TTAGGG-3', d'une extension simple brin 3' sortante et d'un complexe protéique formé des six facteurs TRF1, TRF2, hRap1, TIN2, TPP1 et POT1 qui, s'associant à cette séquence, protègent l'ADN télomèrique. Les télomères sont maintenus par la télomérase, une transcriptase inverse capable d'allonger l'extension 3' sortante télomérique. POT1 lie l'ADN simple brin télomérique et module l'élongation des télomères par la télomérase. POT1 pourrait en théorie également influencer la réplication semi-conservative de l'ADN. L'ADN-polymérase Pal alpha-primase amorce et initie la synthèse d'ADN. Pendant la réplication, l'ADN simple brin est stabilisé par RPA, un complexe hétérotrimèrique qui lie l'ADN simple brin. RPA facilite la synthèse du brin naissant à la fourche de réplication. Ici nous avons étudié comment ces deux protéines qui lient l'ADN simple brin, RPA et POT1, régulent la réplication des télomères par la télomérase et la machinerie classique de réplication de l'ADN. Par immunoprécipitation de chromatine (ChIP), nous montrons que RPA est localisé aux télomères lors de la phase S du cycle cellulaire. De plus, l'analyse de la structure des télomeres indique que !a perte de RPA ou de POT1 conduit à l'apparition d'ADN simple brin télomérique, suggérant une réplication incomplète de l'ADN télomérique in vivo. Par une approche complémentaire biochimique utilisant les protéines POT1 et RPA recombinantes purifiées, nous montrons également que la liaison de POT1 ou de RPA à des oligonucléotides télomériques bloque l'activité primase du complexe polymérase alpha/primase et réduit l'activité télomérase sur ces substrats. En revanche, leur liaison augmente l'activité ADN-polymérase du complexe polymérase alpha/primase, ainsi que fa processivité de la télomérase dans certaines conditions, POT1 étant le plus efficace des deux facteurs. Nous proposons que POT1 est nécessaire à la réplication du brin retardé au niveau des télomères, ce qui suggère que certains phénotypes des cellules déplétés en POT1 puissent résulter d'une réplication incomplète de l'ADN télémétrique plutôt que d'une déprotection de l'extrémité sortante des télomères.

Comparison of human chromosome 21 conserved nongenic sequences (CNGs) with the mouse and dog genomes shows that their selective constraint is independent of their genic environment.

The analysis of conservation between the human and mouse genomes resulted in the identification of a large number of conserved nongenic sequences (CNGs). The functional significance of this nongenic conservation remains unknown, however. The availability of the sequence of a third mammalian genome, the dog, allows for a large-scale analysis of evolutionary attributes of CNGs in mammals. We have aligned 1638 previously identified CNGs and 976 conserved exons (CODs) from human chromosome 21 (Hsa21) with their orthologous sequences in mouse and dog. Attributes of selective constraint, such as sequence conservation, clustering, and direction of substitutions were compared between CNGs and CODs, showing a clear distinction between the two classes. We subsequently performed a chromosome-wide analysis of CNGs by correlating selective constraint metrics with their position on the chromosome and relative to their distance from genes. We found that CNGs appear to be randomly arranged in intergenic regions, with no bias to be closer or farther from genes. Moreover, conservation and clustering of substitutions of CNGs appear to be completely independent of their distance from genes. These results suggest that the majority of CNGs are not typical of previously described regulatory elements in terms of their location. We propose models for a global role of CNGs in genome function and regulation, through long-distance cis or trans chromosomal interactions.

Genetic aspects of normal and disturbed sleep.

Sleep disorders commonly involve genetic susceptibility, environmental effects, and interactions between these factors. The heritability of sleep patterns has been shown in studies of monozygotic twins, and sleep electroencephalogram patterns offer a unique genetic fingerprint which may assist in the identification of genes involved in the regulation of sleep. Genetic factors are also thought to play a role in sleep disorders; narcolepsy is a disabling sleep condition and research has revealed the complexity of underlying genetic and environmental influences in the development of this disorder. An understanding of sleep regulation at the molecular level is essential in the identification of new targets for the treatment of sleep disorders, and genome-wide association studies for both normal sleep and sleep disorders may shed new light on the molecular architecture of mechanisms regulating these behaviours.

Extensive gene traffic on the mammalian X chromosome.

Mammalian sex chromosomes have undergone profound changes since evolving from ancestral autosomes. By examining retroposed genes in the human and mouse genomes, we demonstrate that, during evolution, the mammalian X chromosome has generated and recruited a disproportionately high number of functional retroposed genes, whereas the autosomes experienced lower gene turnover. Most autosomal copies originating from X-linked genes exhibited testis-biased expression. Such export is incompatible with mutational bias and is likely driven by natural selection to attain male germline function. However, the excess recruitment is consistent with a combination of both natural selection and mutational bias.

A burst of segmental duplications in the genome of the African great ape ancestor

It is generally accepted that the extent of phenotypic change between human and great apes is dissonant with the rate of molecular change. Between these two groups, proteins are virtually identical, cytogenetically there are few rearrangements that distinguish ape-human chromosomes, and rates of single-base-pair change and retrotransposon activity have slowed particularly within hominid lineages when compared to rodents or monkeys. Studies of gene family evolution indicate that gene loss and gain are enriched within the primate lineage. Here, we perform a systematic analysis of duplication content of four primate genomes (macaque, orang-utan, chimpanzee and human) in an effort to understand the pattern and rates of genomic duplication during hominid evolution. We find that the ancestral branch leading to human and African great apes shows the most significant increase in duplication activity both in terms of base pairs and in terms of events. This duplication acceleration within the ancestral species is significant when compared to lineage-specific rate estimates even after accounting for copy-number polymorphism and homoplasy. We discover striking examples of recurrent and independent gene-containing duplications within the gorilla and chimpanzee that are absent in the human lineage. Our results suggest that the evolutionary properties of copy-number mutation differ significantly from other forms of genetic mutation and, in contrast to the hominid slowdown of single-base-pair mutations, there has been a genomic burst of duplication activity at this period during human evolution.

Recurrent deletions and reciprocal duplications of 10q11.21q11.23 including CHAT and SLC18A3 are likely mediated by complex low-copy repeats.

We report 24 unrelated individuals with deletions and 17 additional cases with duplications at 10q11.21q21.1 identified by chromosomal microarray analysis. The rearrangements range in size from 0.3 to 12 Mb. Nineteen of the deletions and eight duplications are flanked by large, directly oriented segmental duplications of >98% sequence identity, suggesting that nonallelic homologous recombination (NAHR) caused these genomic rearrangements. Nine individuals with deletions and five with duplications have additional copy number changes. Detailed clinical evaluation of 20 patients with deletions revealed variable clinical features, with developmental delay (DD) and/or intellectual disability (ID) as the only features common to a majority of individuals. We suggest that some of the other features present in more than one patient with deletion, including hypotonia, sleep apnea, chronic constipation, gastroesophageal and vesicoureteral refluxes, epilepsy, ataxia, dysphagia, nystagmus, and ptosis may result from deletion of the CHAT gene, encoding choline acetyltransferase, and the SLC18A3 gene, mapping in the first intron of CHAT and encoding vesicular acetylcholine transporter. The phenotypic diversity and presence of the deletion in apparently normal carrier parents suggest that subjects carrying 10q11.21q11.23 deletions may exhibit variable phenotypic expressivity and incomplete penetrance influenced by additional genetic and nongenetic modifiers.

The evolution of gene expression levels in mammalian organs.

Changes in gene expression are thought to underlie many of the phenotypic differences between species. However, large-scale analyses of gene expression evolution were until recently prevented by technological limitations. Here we report the sequencing of polyadenylated RNA from six organs across ten species that represent all major mammalian lineages (placentals, marsupials and monotremes) and birds (the evolutionary outgroup), with the goal of understanding the dynamics of mammalian transcriptome evolution. We show that the rate of gene expression evolution varies among organs, lineages and chromosomes, owing to differences in selective pressures: transcriptome change was slow in nervous tissues and rapid in testes, slower in rodents than in apes and monotremes, and rapid for the X chromosome right after its formation. Although gene expression evolution in mammals was strongly shaped by purifying selection, we identify numerous potentially selectively driven expression switches, which occurred at different rates across lineages and tissues and which probably contributed to the specific organ biology of various mammals.

Solving ambiguities in contig assembly of Idiomarina loihiensis L2TR chromosome by in silico analyses.

Nucleotide composition analyses of bacterial genomes such as cumulative GC skew highlight the atypical, strongly asymmetric architecture of the recently published chromosome of Idiomarina loihiensis L2TR, suggesting that an inversion of a 600-kb chromosomal segment occurred. The presence of 3.4-kb inverted repeated sequences at the borders of the putative rearrangement supports this hypothesis. Reverting in silico this segment restores (1) a symmetric chromosome architecture; (2) the co-orientation of transcription of all rRNA operons with DNA replication; and (3) a better conservation of gene order between this chromosome and other gamma-proteobacterial ones. Finally, long-range PCRs encompassing the ends of the 600-kb segment reveal the existence of the reverted configuration but not of the published one. This demonstrates how cumulative nucleotide-skew analyses can validate genome assemblies.

Population consequences of environmental sex reversal.

When sex determination in a species is predominantly genetic but environmentally reversible, exposure to (anthropogenic) changes in the environment can lead to shifts in a population's sex ratio. Such scenarios may be common in many fishes and amphibians, yet their ramifications remain largely unexplored. We used a simple model to study the (short-term) population consequences of environmental sex reversal (ESR). We examined the effects on sex ratios, sex chromosome frequencies, and population growth and persistence after exposure to environmental forces with feminizing or masculinizing tendencies. When environmental feminization was strong, X chromosomes were driven to extinction. Analogously, extinction of normally male-linked genetic factors (e.g., Y chromosomes) was caused by continuous environmental masculinization. Although moderate feminization was beneficial for population growth in the absence of large viability effects, our results suggest that the consequences of ESR are generally negative in terms of population size and the persistence of sex chromosomes. Extreme sex ratios resulting from high rates of ESR also reduced effective population sizes considerably. This may limit any evolutionary response to the deleterious effects of ESR. Our findings suggest that ESR changes population growth and sex ratios in some counter-intuitive ways and can change the predominant factor in sex determination from genetic to fully environmental, often within only a few tens of generations. Populations that lose genetic sex determination may quickly go extinct if the environmental forces that cause sex reversal cease.