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.

Autosomal Recessive Posterior Microphthalmos/Nanophthalmos is Caused by Loss-of-function Mutations in LOC646960, a Novel Gene Encoding a Trypsin-like Serine Protease

Purpose: Posterior microphthalmos (MCOP)/nanophthalmos (NNO) is a developmental anomaly characterized by extreme hyperopia due to short axial length. The population of the Faroe Islands shows a high prevalence of an autosomal recessive form (arMCOP). The gene mutated in arMCOP is not yet known.Methods: Genetic mapping by linkage analysis using microsatellite and single nucleotide polymorphisms, mutation analysis by PCR and sequencing, molecular modellingResults: Having refined the position of the disease locus (MCOP6) in an interval of 250 kb in chromosome 2q37.1 in Faroese families, we detected 3 mutations in a novel gene, LOC646960: Patients of 10 different Faroese families were either homozygous (n=22) for c.926G>C (p.Trp309Ser) or compound heterozygous (n=6) for c.926G>C and c.526C>G (p.Arg176Gly), whereas a homozygous 1 bp duplication (c.1066dupC) was identified in patients with arNNO from a Tunisian family. In two unrelated patients with MCOP, no LOC646960 mutation was found. LOC646960 is expressed in the human adult retina and RPE. The expression of the mouse homologue in the eye can be first detected at E17 and is highest in adults. The predicted protein is a 603 amino acid long secreted trypsin-like serine peptidase. c.1066dupC should result in a functional null allele. Molecular modelling of the p.Trp309Ser mutant suggests that both affinity and reactivity of the enzyme towards in vivo substrates are substantially reduced.Conclusions: Postnatal growth of the eye is important for proper development of the refractive components (emmetropization), and is mainly due to elongation of the posterior segment from 10-11 mm at birth to 15-16 mm at the age of 13 years. Optical defocus leads to changes in axial length by moving the retina towards the image plane. arMCOP may theoretically be explained, in line with the expression pattern of LOC646960, by a postnatal growth retardation of the posterior segment.

Developmental and Environmental Regulation of Aquaporin Gene Expression across Populus Species: Divergence or Redundancy?

Aquaporins (AQPs) are membrane channels belonging to the major intrinsic proteins family and are known for their ability to facilitate water movement. While in Populus trichocarpa, AQP proteins form a large family encompassing fifty-five genes, most of the experimental work focused on a few genes or subfamilies. The current work was undertaken to develop a comprehensive picture of the whole AQP gene family in Populus species by delineating gene expression domain and distinguishing responsiveness to developmental and environmental cues. Since duplication events amplified the poplar AQP family, we addressed the question of expression redundancy between gene duplicates. On these purposes, we carried a meta-analysis of all publicly available Affymetrix experiments. Our in-silico strategy controlled for previously identified biases in cross-species transcriptomics, a necessary step for any comparative transcriptomics based on multispecies design chips. Three poplar AQPs were not supported by any expression data, even in a large collection of situations (abiotic and biotic constraints, temporal oscillations and mutants). The expression of 11 AQPs was never or poorly regulated whatever the wideness of their expression domain and their expression level. Our work highlighted that PtTIP1;4 was the most responsive gene of the AQP family. A high functional divergence between gene duplicates was detected across species and in response to tested cues, except for the root-expressed PtTIP2;3/PtTIP2;4 pair exhibiting 80% convergent responses. Our meta-analysis assessed key features of aquaporin expression which had remained hidden in single experiments, such as expression wideness, response specificity and genotype and environment interactions. By consolidating expression profiles using independent experimental series, we showed that the large expansion of AQP family in poplar was accompanied with a strong divergence of gene expression, even if some cases of functional redundancy could be suspected.

Mirror extreme BMI phenotypes associated with gene dosage at the chromosome 16p11.2 locus.

Both obesity and being underweight have been associated with increased mortality. Underweight, defined as a body mass index (BMI) ≤ 18.5 kg per m(2) in adults and ≤ -2 standard deviations from the mean in children, is the main sign of a series of heterogeneous clinical conditions including failure to thrive, feeding and eating disorder and/or anorexia nervosa. In contrast to obesity, few genetic variants underlying these clinical conditions have been reported. We previously showed that hemizygosity of a ∼600-kilobase (kb) region on the short arm of chromosome 16 causes a highly penetrant form of obesity that is often associated with hyperphagia and intellectual disabilities. Here we show that the corresponding reciprocal duplication is associated with being underweight. We identified 138 duplication carriers (including 132 novel cases and 108 unrelated carriers) from individuals clinically referred for developmental or intellectual disabilities (DD/ID) or psychiatric disorders, or recruited from population-based cohorts. These carriers show significantly reduced postnatal weight and BMI. Half of the boys younger than five years are underweight with a probable diagnosis of failure to thrive, whereas adult duplication carriers have an 8.3-fold increased risk of being clinically underweight. We observe a trend towards increased severity in males, as well as a depletion of male carriers among non-medically ascertained cases. These features are associated with an unusually high frequency of selective and restrictive eating behaviours and a significant reduction in head circumference. Each of the observed phenotypes is the converse of one reported in carriers of deletions at this locus. The phenotypes correlate with changes in transcript levels for genes mapping within the duplication but not in flanking regions. The reciprocal impact of these 16p11.2 copy-number variants indicates that severe obesity and being underweight could have mirror aetiologies, possibly through contrasting effects on energy balance.

Vitellogenesis and the vitellogenin gene family.

Vitellogenin is synthesized under estrogen control in the liver, extensively modified, transported to the ovary, and there processed to the yolk proteins lipovitellin and phosvitin. In the frog Xenopus laevis there are at least four distinct but related vitellogenin genes. The two genes A1 and A2 have a 95 percent sequence homology in their messenger RNA coding regions, and contain 33 introns that interrupt the coding region (exons) at homologous positions. Sequences and lengths of analogous introns differ, and many introns contain repetitive DNA elements. The introns in these two genes that have apparently arisen by duplication have diverged extensively by events that include deletions, insertions, and probably duplications. Rapid evolutionary change involving rearrangements and the presence of repeated DNA suggests that the bulk of the sequences within introns may not have any specific function.

RPS4Y gene family evolution in primates

Background: The RPS4 gene codifies for ribosomal protein S4, a very well-conserved protein present in all kingdoms. In primates, RPS4 is codified by two functional genes located on both sex chromosomes: the RPS4X and RPS4Y genes. In humans, RPS4Y is duplicated and the Y chromosome therefore carries a third functional paralog: RPS4Y2, which presents a testis-specific expression pattern. Results: DNA sequence analysis of the intronic and cDNA regions of RPS4Y genes from species covering the entire primate phylogeny showed that the duplication event leading to the second Y-linked copy occurred after the divergence of New World monkeys, about 35 million years ago. Maximum likelihood analyses of the synonymous and non-synonymous substitutions revealed that positive selection was acting on RPS4Y2 gene in the human lineage, which represents the first evidence of positive selection on a ribosomal protein gene. Putative positive amino acid replacements affected the three domains of the protein: one of these changes is located in the KOW protein domain and affects the unique invariable position of this motif, and might thus have a dramatic effect on the protein function.Conclusion: Here, we shed new light on the evolutionary history of RPS4Y gene family, especially on that of RPS4Y2. The results point that the RPS4Y1 gene might be maintained to compensate gene dosage between sexes, while RPS4Y2 might have acquired a new function, at least in the lineage leading to humans.

Genesis and evolution of the Evx and Mox genes and the extended Hox and ParaHox gene clusters

Background: Hox and ParaHox gene clusters are thought to have resulted from the duplication of a ProtoHox gene cluster early in metazoan evolution. However, the origin and evolution of the other genes belonging to the extended Hox group of homeobox-containing genes, that is, Mox and Evx, remains obscure. We constructed phylogenetic trees with mouse, amphioxus and Drosophila extended Hox and other related Antennapedia-type homeobox gene sequences and analyzed the linkage data available for such genes.Results: We claim that neither Mox nor Evx is a Hox or ParaHox gene. We propose a scenariothat reconciles phylogeny with linkage data, in which an Evx/Mox ancestor gene linked to aProtoHox cluster was involved in a segmental tandem duplication event that generated an arrayof all Hox-like genes, referred to as the `coupled¿ cluster. A chromosomal breakage within thiscluster explains the current composition of the extended Hox cluster (with Evx, Hox and Moxgenes) and the ParaHox cluster.Conclusions: Most studies dealing with the origin and evolution of Hox and ParaHox clustershave not included the Hox-related genes Mox and Evx. Our phylogenetic analyses and theavailable linkage data in mammalian genomes support an evolutionary scenario in which anancestor of Evx and Mox was linked to the ProtoHox cluster, and that a tandem duplication of alarge genomic region early in metazoan evolution generated the Hox and ParaHox clusters, plusthe cluster-neighbors Evx and Mox. The large `coupled¿ Hox-like cluster EvxHox/MoxParaHox wassubsequently broken, thus grouping the Mox and Evx genes to the Hox clusters, and isolating theParaHox cluster.

Repeated Whole-Genome Duplication, Karyotype Reshuffling, and Biased Retention of Stress-Responding Genes in Buckler Mustard.

Whole-genome duplication (WGD) is usually followed by gene loss and karyotype repatterning. Despite evidence of new adaptive traits associated with WGD, the underpinnings and evolutionary significance of such genome fractionation remain elusive. Here, we use Buckler mustard (Biscutella laevigata) to infer processes that have driven the retention of duplicated genes after recurrent WGDs. In addition to the β- and α-WGD events shared by all Brassicaceae, cytogenetic and transcriptome analyses revealed two younger WGD events that occurred at times of environmental changes in the clade of Buckler mustard (Biscutelleae): a mesopolyploidy event from the late Miocene that was followed by considerable karyotype reshuffling and chromosome number reduction and a neopolyploidy event during the Pleistocene. Although a considerable number of the older duplicates presented signatures of retention under positive selection, the majority of retained duplicates arising from the younger mesopolyploidy WGD event matched predictions of the gene balance hypothesis and showed evidence of strong purifying selection as well as enrichment in gene categories responding to abiotic stressors. Retention of large stretches of chromosomes for both genomic copies supported the hypothesis that cycles of WGD and biased fractionation shaped the genome of this stress-tolerant polypolyloid, promoting the adaptive recruitment of stress-responding genes in the face of environmental challenges.

Functional and evolutionary analysis of DXL1, a non-essential gene encoding a 1-deoxy-D-xylulose 5-phosphate synthase like protein in Arabidopsis thaliana

The synthesis of 1-deoxy-D-xylulose 5-phosphate (DXP), catalyzed by the enzyme DXP synthase (DXS), represents a key regulatory step of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway for isoprenoid biosynthesis. In plants DXS is encoded by small multigene families that can be classified into, at least, three specialized subfamilies. Arabidopsis thaliana contains three genes encoding proteins with similarity to DXS, including the well-known DXS1/CLA1 gene, which clusters within subfamily I. The remaining proteins, initially named DXS2 and DXS3, have not yet been characterized. Here we report the expression and functional analysis of A. thaliana DXS2. Unexpectedly, the expression of DXS2 failed to rescue Escherichia coli and A. thaliana mutants defective in DXS activity. Coherently, we found that DXS activity was negligible in vitro, being renamed as DXL1 following recent nomenclature recommendation. DXL1 is targeted to plastids as DXS1, but shows a distinct expression pattern. The phenotypic analysis of a DXL1 defective mutant revealed that the function of the encoded protein is not essential for growth and development. Evolutionary analyses indicated that DXL1 emerged from DXS1 through a recent duplication apparently specific of the Brassicaceae lineage. Divergent selective constraints would have affected a significant fraction of sites after diversification of the paralogues. Furthermore, amino acids subjected to divergent selection and likely critical for functional divergence through the acquisition of a novel, although not yet known, biochemical function, were identified. Our results provide with the first evidences of functional specialization at both the regulatory and biochemical level within the plant DXS family.

Consequences of lineage-specific gene loss on functional evolution of surviving paralogs: ALDH1A and retinoic acid signaling in vertebrate genomes

Genome duplications increase genetic diversity and may facilitate the evolution of gene subfunctions. Little attention, however, has focused on the evolutionary impact of lineage-specific gene loss. Here, we show that identifying lineage-specific gene loss after genome duplication is important for understanding the evolution of gene subfunctions in surviving paralogs and for improving functional connectivity among human and model organism genomes. We examine the general principles of gene loss following duplication, coupled with expression analysis of the retinaldehyde dehydrogenase Aldh1a gene family during retinoic acid signaling in eye development as a case study. Humans have three ALDH1A genes, but teleosts have just one or two. We used comparative genomics and conserved syntenies to identify loss of ohnologs (paralogs derived from genome duplication) and to clarify uncertain phylogenies. Analysis showed that Aldh1a1 and Aldh1a2 form a clade that is sister to Aldh1a3-related genes. Genome comparisons showed secondarily loss of aldh1a1 in teleosts, revealing that Aldh1a1 is not a tetrapod innovation and that aldh1a3 was recently lost in medaka, making it the first known vertebrate with a single aldh1a gene. Interestingly, results revealed asymmetric distribution of surviving ohnologs between co-orthologous teleost chromosome segments, suggesting that local genome architecture can influence ohnolog survival. We propose a model that reconstructs the chromosomal history of the Aldh1a family in the ancestral vertebrate genome, coupled with the evolution of gene functions in surviving Aldh1a ohnologs after R1, R2, and R3 genome duplications. Results provide evidence for early subfunctionalization and late subfunction-partitioning and suggest a mechanistic model based on altered regulation leading to heterochronic gene expression to explain the acquisition or modification of subfunctions by surviving ohnologs that preserve unaltered ancestral developmental programs in the face of gene loss.

Lateral gene exchanges shape the genomes of amoeba-resisting microorganisms.

Based on Darwin's concept of the tree of life, vertical inheritance was thought to be dominant, and mutations, deletions, and duplication were streaming the genomes of living organisms. In the current genomic era, increasing data indicated that both vertical and lateral gene inheritance interact in space and time to trigger genome evolution, particularly among microorganisms sharing a given ecological niche. As a paradigm to their diversity and their survival in a variety of cell types, intracellular microorganisms, and notably intracellular bacteria, were considered as less prone to lateral genetic exchanges. Such specialized microorganisms generally have a smaller gene repertoire because they do rely on their host's factors for some basic regulatory and metabolic functions. Here we review events of lateral gene transfer (LGT) that illustrate the genetic exchanges among intra-amoebal microorganisms or between the microorganism and its amoebal host. We tentatively investigate the functions of laterally transferred genes in the light of the interaction with their host as they should confer a selective advantage and success to the amoeba-resisting microorganisms (ARMs).

Thr(118)Met amino acid substitution in the peripheral myelin protein 22 does not influence the clinical phenotype of Charcot-Marie-Tooth disease type 1A due to the 17p11.2-p12 duplication

The Thr(118)Met substitution in the peripheral myelin protein 22 (PMP22) gene has been detected in a number of families with demyelinating Charcot-Marie-Tooth (CMT1) neuropathy or with the hereditary neuropathy with liability to pressure palsy, but in none of them has it consistently segregated with the peripheral neuropathy. We describe here a CMT1 family (a 63-year-old man, his brother and his niece) in which two mutations on different chromosomes were found in the PMP22 gene, the 17p duplication, detected by fluorescent semiquantitative polymerase chain reaction (PCR) of microsatellite markers localized within the duplicated region on chromosome 17p11.2-p12, and the Thr(118)Met substitution, detected by direct sequencing the four coding exons of the PMP22 gene. A genotype/phenotype correlation study showed that the neuropathy segregates with the duplication and that the amino acid substitution does not seem to modify the clinical characteristics or the severity of the peripheral neuropathy. We did not find any evidence to characterize this substitution as a polymorphism in the population studied and we propose that the high frequency reported for this point mutation in the literature suggests that the Thr(118)Met substitution may be a hotspot for mutations in the PMP22 gene.

Approches algorithmiques pour l’inférence d’histoires de duplication en tandem avec inversions et délétions pour des familles multigéniques

[Français] Une fraction importante des génomes eucaryotes est constituée de Gènes Répétés en Tandem (GRT). Un mécanisme fondamental dans l’évolution des GRT est la recombinaison inégale durant la méiose, entrainant la duplication locale (en tandem) de segments chromosomiques contenant un ou plusieurs gènes adjacents. Différents algorithmes ont été proposés pour inférer une histoire de duplication en tandem pour un cluster de GRT. Cependant, leur utilisation est limitée dans la pratique, car ils ne tiennent pas compte d’autres événements évolutifs pourtant fréquents, comme les inversions, les duplications inversées et les délétions. Cette thèse propose différentes approches algorithmiques permettant d’intégrer ces événements dans le modèle de duplication en tandem classique. Nos contributions sont les suivantes: • Intégrer les inversions dans un modèle de duplication en tandem simple (duplication d’un gène à la fois) et proposer un algorithme exact permettant de calculer le nombre minimal d’inversions s’étant produites dans l’évolution d’un cluster de GRT. • Généraliser ce modèle pour l’étude d’un ensemble de clusters orthologues dans plusieurs espèces. • Proposer un algorithme permettant d’inférer l’histoire évolutive d’un cluster de GRT en tenant compte des duplications en tandem, duplications inversées, inversions et délétions de segments chromosomiques contenant un ou plusieurs gènes adjacents.

On the function of the Dictyostelium Argonaute A protein (AgnA) in epigenetic gene regulation

With molecular biology methods and bioinformatics, the Argonaute proteins in Dictyostelium discoideum were characterized, and the function of the AgnA protein in RNAi and DNA methylation was investigated, as well as cellular features. Also interaction partners of the PAZ-Piwi domain of AgnA (PAZ-PiwiAgnA) were discovered. The Dictyostelium genome encodes five Argonaute proteins, termed AgnA/B/C/D/E. The expression level of Argonaute proteins was AgnB/D/E > AgnA > AgnC. All these proteins contain the characteristic conserved of PAZ and Piwi domains. Fluorescence microscopy revealed that the overexpressed C-terminal GFP-fusion of PAZ-PiwiAgnA (PPWa-GFP) localized to the cytoplasm. Overexpression of PPWa-GFP leaded to an increased gene silencing efficiency mediated by RNAi but not by antisense RNA. This indicated that PAZ-PiwiAgnA is involved in the RNAi pathway, but not in the antisense pathway. An analysis of protein-protein interactions by a yeast-two-hybrid screen on a cDNA library from vegetatively grown Dictyostelium revealed that several proteins, such as EF2, EF1-I, IfdA, SahA, SamS, RANBP1, UAE1, CapA, and GpdA could interact with PAZ-PiwiAgnA. There was no interaction between PAZ-PiwiAgnA and HP1, HelF and DnmA detected by direct yeast-two-hybrid analysis. The fluorescence microscopy images showed that the overexpressed GFP-SahA or IfdA fusion proteins localized to both cytoplasm and nuclei, while the overexpressed GFP-SamS localized to the cytoplasm. The expression of SamS in AgnA knock down mutants was strongly down regulated on cDNA and mRNA level in, while the expression of SahA was only slightly down regulated. AgnA knock down mutants displayed defects in growth and phagocytosis, which suggested that AgnA affects also cell biological features. The inhibition of DNA methylation on DIRS-1 and Skipper retroelements, as well as the endogenous mvpB and telA gene, observed for the same strains, revealed that AgnA is involved in the DNA methylation pathway. Northern blot analysis showed that Skipper and DIRS-1 were rarely expressed in Ax2, but the expression of Skipper was upregulated in AgnA knock down mutants, while the expression of DIRS-1 was not changed. A knock out of the agnA gene failed even though the homologous recombination of the disruption construct occurred at the correct site, which indicated that there was a duplication of the agnA gene in the genome. The same phenomenon was also observed in ifdA knock out experiments.

Genome dynamics explain the evolution of flowering time CCT domain gene families in the Poaceae

Numerous CCT domain genes are known to control flowering in plants. They belong to the CONSTANS-like (COL) and PREUDORESPONSE REGULATOR (PRR) gene families, which in addition to a CCT domain possess B-box or response-regulator domains, respectively. Ghd7 is the most recently identified COL gene to have a proven role in the control of flowering time in the Poaceae. However, as it lacks B-box domains, its inclusion within the COL gene family, technically, is incorrect. Here, we show Ghd7 belongs to a larger family of previously uncharacterized Poaceae genes which possess just a single CCT domain, termed here CCT MOTIF FAMILY (CMF) genes. We molecularly describe the CMF (and related COL and PRR) gene families in four sequenced Poaceae species, as well as in the draft genome assembly of barley (Hordeum vulgare). Genetic mapping of the ten barley CMF genes identified, as well as twelve previously unmapped HvCOL and HvPRR genes, finds the majority map to colinear positions relative to their Poaceae orthologues. Combined inter-/intra-species comparative and phylogenetic analysis of CMF, COL and PRR gene families indicates they evolved prior to the monocot/dicot divergence ~200 mya, with Poaceae CMF evolution described as the interplay between whole genome duplication in the ancestral cereal, and subsequent clade-specific mutation, deletion and duplication events. Given the proven role of CMF genes in the modulation of cereals flowering, the molecular, phylogenetic and comparative analysis of the Poaceae CMF, COL and PRR gene families presented here provides the foundation from which functional investigation can be undertaken.