Associations between human disease genes and overlapping gene groups and multiple amino acid runs

Guidance for estimating penetrance of monogenic disease-causing variants in population cohorts

Penetrance is the probability that an individual with a pathogenic genetic variant develops a specific disease. Knowing the penetrance of variants for monogenic disorders is important for counseling of individuals. Until recently, estimates of penetrance have largely relied on affected individuals and their at-risk family members being clinically referred for genetic testing, a 'phenotype-first' approach. This approach substantially overestimates the penetrance of variants because of ascertainment bias. The recent availability of whole-genome sequencing data in individuals from very-large-scale population-based cohorts now allows 'genotype-first' estimates of penetrance for many conditions. Although this type of population-based study can underestimate penetrance owing to recruitment biases, it provides more accurate estimates of penetrance for secondary or incidental findings. Here, we provide guidance for the conduct of penetrance studies to ensure that robust genotypes and phenotypes are used to accurately estimate penetrance of variants and groups of similarly annotated variants from population-based studies.

Discovery of long non-coding RNAs in the liver fluke, Fasciola hepatica

Long non-coding (lnc)RNAs are a class of eukaryotic RNA that do not code for protein and are linked with transcriptional regulation, amongst a myriad of other functions. Using a custom in silico pipeline we have identified 6,436 putative lncRNA transcripts in the liver fluke parasite, Fasciola hepatica, none of which are conserved with those previously described from Schistosoma mansoni. F. hepatica lncRNAs were distinct from F. hepatica mRNAs in transcript length, coding probability, exon/intron composition, expression patterns, and genome distribution. RNA-Seq and digital droplet PCR measurements demonstrated developmentally regulated expression of lncRNAs between intra-mammalian life stages; a similar proportion of lncRNAs (14.2%) and mRNAs (12.8%) were differentially expressed (p<0.001), supporting a functional role for lncRNAs in F. hepatica life stages. While most lncRNAs (81%) were intergenic, we identified some that overlapped protein coding loci in antisense (13%) or intronic (6%) configurations. We found no unequivocal evidence for correlated developmental expression within positionally correlated lncRNA:mRNA pairs, but global co-expression analysis identified five lncRNA that were inversely co-regulated with 89 mRNAs, including a large number of functionally essential proteases. The presence of micro (mi)RNA binding sites in 3135 lncRNAs indicates the potential for miRNA-based post-transcriptional regulation of lncRNA, and/or their function as competing endogenous (ce)RNAs. The same annotation pipeline identified 24,141 putative lncRNAs in F. gigantica. This first description of lncRNAs in F. hepatica provides an avenue to future functional and comparative genomics studies that will provide a new perspective on a poorly understood aspect of parasite biology.

Amino acid homorepeats in proteins

Amino acid homorepeats, or homorepeats, are polypeptide segments found in proteins that contain stretches of identical amino acid residues. Although abnormal homorepeat expansions are linked to pathologies such as neurodegenerative diseases, homorepeats are prevalent in eukaryotic proteomes, suggesting that they are important for normal physiology. In this Review, we discuss recent advances in our understanding of the biological functions of homorepeats, which range from facilitating subcellular protein localization to mediating interactions between proteins across diverse cellular pathways. We explore how the functional diversity of homorepeat-containing proteins could be linked to the ability of homorepeats to adopt different structural conformations, an ability influenced by repeat composition, repeat length and the nature of flanking sequences. We conclude by highlighting how an understanding of homorepeats will help us better characterize and develop therapeutics against the human diseases to which they contribute.

The morbid cutaneous anatomy of the human genome revealed by a bioinformatic approach

Computational approaches have been developed to prioritize candidate genes in disease gene identification. They are based on different pieces of evidences associating each gene with the given disease. In this study, 648 genes underlying genodermatoses have been compared to 1808 genes involved in other genetic diseases using a bioinformatic approach. These genes were studied at the structural, evolutionary and functional levels. Results show that genes underlying genodermatoses present longer CDS and have more exons. Significant differences were observed in nucleotide motif and amino-acid compositions. Evolutionary conservation analysis revealed that genodermatoses genes have less paralogs, more orthologs in Mouse and Dog and are less conserved. Functional analysis revealed that genodermatosis genes seem to be involved in immune system and skin layers. The Bayesian network model returned a rate of good classification of around 80%. This computational approach could help investigators working in the field of dermatology by prioritizing positional candidate genes for mutation screening.

Overlapping protein-coding genes in human genome and their coincidental expression in tissues

The completion of human genome sequences and the advancement of next-generation sequencing technologies have engendered a clear understanding of all human genes. Overlapping genes are usually observed in compact genomes, such as those of bacteria and viruses. Notably, overlapping protein-coding genes do exist in human genome sequences. Accordingly, we used the current Ensembl gene annotations to identify overlapping human protein-coding genes. We analysed 19,200 well-annotated protein-coding genes and determined that 4,951 protein-coding genes overlapped with their adjacent genes. Approximately a quarter of all human protein-coding genes were overlapping genes. We observed different clusters of overlapping protein-coding genes, ranging from two genes (paired overlapping genes) to 22 genes. We also divided the paired overlapping protein-coding gene groups into four subtypes. We found that the divergent overlapping gene subtype had a stronger expression association than did the subtypes of 5'-tandem overlapping and 3'-tandem overlapping genes. The majority of paired overlapping genes exhibited comparable coincidental tissue expression profiles; however, a few overlapping gene pairs displayed distinctive tissue expression association patterns. In summary, we have carefully examined the genomic features and distributions about human overlapping protein-coding genes and found coincidental expression in tissues for most overlapping protein-coding genes.

The flexibility of a homeodomain transcription factor heterodimer and its allosteric regulation by DNA binding

Transcription factors are known to modify the DNA that they bind. However, DNA can also serve as an allosteric ligand whose binding modifies the conformation of transcriptional regulators. Here, we describe how heterodimer PBX1:PREP1, formed by proteins playing major roles in embryonic development and tumorigenesis, undergoes an allosteric transition upon DNA binding. We demonstrate through a number of biochemical and biophysical methods that PBX1:PREP1 exhibits a structural change upon DNA binding. Small-angle X-ray scattering (SAXS), circular dichroism (CD), isothermal titration calorimetry (ITC), and limited proteolysis demonstrate a different shape, α-helical content, thermodynamic behavior, and solution environment of the holo-complex (with DNA) compared to the apo-complex (without DNA). Given that PBX1 as such does not have a defined DNA selectivity, structural changes upon DNA binding become major factors in the function of the PBX1:PREP1 complex. The observed changes are mapped at both the amino- and carboxy-terminal regions of the two proteins thereby providing important insights to determine how PBX1:PREP1 dimer functions.

DATABASE

Small-angle scattering data are available in SASBDB under accession numbers SASDAP7, SASDAQ7, and SASDAR7.

Association of polyalanine and polyglutamine coiled coils mediates expansion disease-related protein aggregation and dysfunction

The expansion of homopolymeric glutamine (polyQ) or alanine (polyA) repeats in certain proteins owing to genetic mutations induces protein aggregation and toxicity, causing at least 18 human diseases. PolyQ and polyA repeats can also associate in the same proteins, but the general extent of their association in proteomes is unknown. Furthermore, the structural mechanisms by which their expansion causes disease are not well understood, and these repeats are generally thought to misfold upon expansion into aggregation-prone β-sheet structures like amyloids. However, recent evidence indicates a critical role for coiled-coil (CC) structures in triggering aggregation and toxicity of polyQ-expanded proteins, raising the possibility that polyA repeats may as well form these structures, by themselves or in association with polyQ. We found through bioinformatics screenings that polyA, polyQ and polyQA repeats have a phylogenetically graded association in human and non-human proteomes and associate/overlap with CC domains. Circular dichroism and cross-linking experiments revealed that polyA repeats can form—alone or with polyQ and polyQA—CC structures that increase in stability with polyA length, forming higher-order multimers and polymers in vitro. Using structure-guided mutagenesis, we studied the relevance of polyA CCs to the in vivo aggregation and toxicity of RUNX2—a polyQ/polyA protein associated with cleidocranial dysplasia upon polyA expansion—and found that the stability of its polyQ/polyA CC controls its aggregation, localization and toxicity. These findings indicate that, like polyQ, polyA repeats form CC structures that can trigger protein aggregation and toxicity upon expansion in human genetic diseases.

Genetic variants and evolutionary analyses of heparin cofactor II

Heparin cofactor II (HCII) belongs to serpin superfamily and it acts as a thrombin inhibitor in the coagulation cascade, in a glycosaminoglycan-dependent pathway using the release of a sequestered hirudin-like N-terminal tail for interaction with thrombin. This serpin belongs to multiple member group V2 of vertebrate serpin classification. However, there is no comprehensive study illustrating the exact phylogenetic history of HCII, to date. Herein, we explored phylogenetic traits of HCII genes. Structures of HCII gene from selected ray-finned fishes and lamprey varied in exon I and II with insertions of novel introns of which one in core domain for ray-finned fishes in exon II at the position 241c. We found HCII remain nested in the largest intron of phosphatidylinositol (PI) 4-kinase (PIK4CA) gene (genetic variants of this gene cause schizophrenia) at the origin of vertebrates, dated about 500MY old. We found that sequence features such as two acidic repeats (AR1-II), GAG-binding helix-D, three serpin motifs and inhibitory reactive center loop (RCL) of HCII protein are highly conserved in 55 vertebrates analyzed. We identified 985 HCII variants by analysis of 1092 human genomes with top three variation classes belongs to SNPs (84.3%), insertion (7.1%) and deletion (5.0%). We identified 37 deleterious mutations in the human HCII protein and we have described these mutations in relation to HCII sequence-structure-function relationships. These understandings may have clinical and medical importance as well.

Spatial localization of genes determined by intranuclear DNA fragmentation with the fusion proteins lamin KRED and histone KRED und visible light

The highly organized DNA architecture inside of the nuclei of cells is accepted in the scientific world. In the human genome about 3 billion nucleotides are organized as chromatin in the cell nucleus. In general, they are involved in gene regulation and transcription by histone modification. Small chromosomes are localized in a central nuclear position whereas the large chromosomes are peripherally positioned. In our experiments we inserted fusion proteins consisting of a component of the nuclear lamina (lamin B1) and also histone H2A, both combined with the light inducible fluorescence protein KillerRed (KRED). After activation, KRED generates reactive oxygen species (ROS) producing toxic effects and may cause cell death. We analyzed the spatial damage distribution in the chromatin after illumination of the cells with visible light. The extent of DNA damage was strongly dependent on its localization inside of nuclei. The ROS activity allowed to gain information about the location of genes and their functions via sequencing and data base analysis of the double strand breaks of the isolated DNA. A connection between the damaged gene sequences and some diseases was found.

Overlapping genes of Aedes aegypti: evolutionary implications from comparison with orthologs of Anopheles gambiae and other insects

Background

Although gene overlapping is a common feature of prokaryote and mitochondria genomes, such genes have also been identified in many eukaryotes. The overlapping genes in eukaryotes are extensively rearranged even between closely related species. In this study, we investigated retention and rearrangement of positionally overlapping genes between the mosquitoes Aedes aegypti (dengue virus vector) and Anopheles gambiae (malaria vector). The overlapping gene pairs of A. aegypti were further compared with orthologs of other selected insects to conduct several hypothesis driven investigations relating to the evolution and rearrangement of overlapping genes.

Results

The results show that as much as ~10% of the predicted genes of A. aegypti and A. gambiae are localized in positional overlapping manner. Furthermore, the study shows that differential abundance of introns and simple sequence repeats have significant association with positional rearrangement of overlapping genes between the two species. Gene expression analysis further suggests that antisense transcripts generated from the oppositely oriented overlapping genes are differentially regulated and may have important regulatory functions in these mosquitoes. Our data further shows that synonymous and non-synonymous mutations have differential but non-significant effect on overlapping localization of orthologous genes in other insect genomes.

Conclusion

Gene overlapping in insects may be a species-specific evolutionary process as evident from non-dependency of gene overlapping with species phylogeny. Based on the results, our study suggests that overlapping genes may have played an important role in genome evolution of insects.

Molecular pathology of polyalanine expansion disorders: new perspectives from mouse models

Disease-causing polyalanine (PA) expansion mutations have been identified in nine genes, eight of which encode transcription factors (TFs) with important roles in development. In vitro and cell overexpression studies have shown that expanded PA tracts result in protein misfolding and the formation of aggregates. This feature of PA proteins is reminiscent of the related polyglutamine (PQ) disease proteins, which have been shown to cause disease via a gain-of-function (GOF) mechanism. However, in sharp contrast to PQ disorders, the disease phenotypes associated with PA mutations are more consistent with a LOF and/or mild GOF mechanism, suggesting that their molecular pathology is inherently different to PQ disorders. Elucidating the cellular impact of PA mutations in vivo has been difficult to address as, unlike the late-onset polyglutamine disorders, all PA disorders associated with TF gene mutations are congenital. However, in recent years, significant advances have been made through the analysis of engineered (knock-in) and spontaneous PA mouse models. Here we review these recent findings and propose an updated model of the molecular and cellular mechanism of PA disorders that incorporates both LOF and GOF features.

Protein homorepeats sequences, structures, evolution, and functions

The vast majority of protein sequences are aperiodic; they do not have any strong bias in the amino acid composition, and they use a subtle mixture of all or most of the 20 amino acid residues to code a great number of various structures and functions. In this context, homorepeats, runs of a single amino acid residue, represent unusual, eye-catching motifs in proteins. Despite the sequence simplicity and relatively small size, the homorepeat runs have a strong potential for molecular interactions due to the excessively high local concentration of a certain physico-chemical property. Appearance of such runs within proteins may give them new structural and functional features. An increasing number of studies demonstrate the abundance of these motifs in proteins, their important roles in biological processes, and their link to a number of hereditary and age-related diseases. In this chapter, we summarize data on the distribution of homorepeats in proteomes and on their structural properties, evolution, and functions.

A novel protein, sperm head and tail associated protein (SHTAP), interacts with cysteine-rich secretory protein 2 (CRISP2) during spermatogenesis in the mouse

BACKGROUND INFORMATION

CRISP2 (cysteine-rich secretory protein 2) is a sperm acrosome and tail protein with the ability to regulate Ca2+ flow through ryanodine receptors. Based on these properties, CRISP2 has a potential role in fertilization through the regulation of ion signalling in the acrosome reaction and sperm motility. The purpose of the present study was to determine the expression, subcellular localization and the role in spermatogenesis of a novel CRISP2-binding partner, which we have designated SHTAP (sperm head and tail associated protein).

RESULTS

Using yeast two-hybrid screens of an adult testis expression library, we identified SHTAP as a novel mouse CRISP2-binding partner. Sequence analysis of all Shtap cDNA clones revealed that the mouse Shtap gene is embedded within a gene encoding the unrelated protein NSUN4 (NOL1/NOP2/Sun domain family member 4). Five orthologues of the Shtap gene have been annotated in public databases. SHTAP and its orthologues showed no significant sequence similarity to any known protein or functional motifs, including NSUN4. Using an SHTAP antiserum, multiple SHTAP isoforms (approximately 20-87 kDa) were detected in the testis, sperm, and various somatic tissues. Interestingly, only the approximately 26 kDa isoform of SHTAP was able to interact with CRISP2. Furthermore, yeast two-hybrid assays showed that both the CAP (CRISP/antigen 5/pathogenesis related-1) and CRISP domains of CRISP2 were required for maximal binding to SHTAP. SHTAP protein was localized to the peri-acrosomal region of round spermatids, and the head and tail of the elongated spermatids and sperm tail where it co-localized with CRISP2. During sperm capacitation, SHTAP and the SHTAP-CRISP2 complex appeared to be redistributed within the head.

CONCLUSIONS

The present study is the first report of the identification, annotation and expression analysis of the mouse Shtap gene. The redistribution observed during sperm capacitation raises the possibility that SHTAP and the SHTAP-CRISP2 complex play a role in the attainment of sperm functional competence.

No association between bipolar disorder and syngr1 or synapsin II polymorphisms in the Han Chinese population

Polymorphisms of the SYNAPTOGYRIN1 (SYNGR1) and SYNASINII (SYNII) genes have been shown to be a risk factor for bipolar disorder or schizophrenia. A case-control study with these two genes was conducted in 506 bipolar disorder patients and 507 healthy individuals from the Han Chinese population. No association was found in this study.

Performance and scalability of discriminative metrics for comparative gene identification in 12 Drosophila genomes

Comparative genomics of multiple related species is a powerful methodology for the discovery of functional genomic elements, and its power should increase with the number of species compared. Here, we use 12 Drosophila genomes to study the power of comparative genomics metrics to distinguish between protein-coding and non-coding regions. First, we study the relative power of different comparative metrics and their relationship to single-species metrics. We find that even relatively simple multi-species metrics robustly outperform advanced single-species metrics, especially for shorter exons (< or =240 nt), which are common in animal genomes. Moreover, the two capture largely independent features of protein-coding genes, with different sensitivity/specificity trade-offs, such that their combinations lead to even greater discriminatory power. In addition, we study how discovery power scales with the number and phylogenetic distance of the genomes compared. We find that species at a broad range of distances are comparably effective informants for pairwise comparative gene identification, but that these are surpassed by multi-species comparisons at similar evolutionary divergence. In particular, while pairwise discovery power plateaued at larger distances and never outperformed the most advanced single-species metrics, multi-species comparisons continued to benefit even from the most distant species with no apparent saturation. Last, we find that genes in functional categories typically considered fast-evolving can nonetheless be recovered at very high rates using comparative methods. Our results have implications for comparative genomics analyses in any species, including the human.

Interactions between homopolymeric amino acids (HPAAs)

Many human proteins contain consecutive amino acid repeats, known as homopolymeric amino acid (HPAA) tracts. Some inherited diseases are caused by proteins in which HPAAs are expanded to an excessive length. To this day, nine polyglutamine-related diseases and nine polyalanine-related diseases have been reported, including Huntington's disease and oculopharyngeal muscular dystrophy. In this study, potential HPAA-HPAA interactions were examined by yeast two-hybrid assays using HPAAs of approximately 30 residues in length. The results indicate that hydrophobic HPAAs interact with themselves and with other hydrophobic HPAAs. Previously, we reported that hydrophobic HPAAs formed large aggregates in COS-7 cells. Here, those HPAAs were shown to have significant interactions with each other, suggesting that hydrophobicity plays an important role in aggregation. Among the observed HPAA-HPAA interactions, the Ala28-Ala29 interaction was notable because polyalanine tracts of these lengths have been established to be pathogenic in several polyalanine-related diseases. By testing several constructs of different lengths, we clarified that polyalanine self-interacts at longer lengths (>23 residues) but not at shorter lengths (six to approximately 23 residues) in a yeast two-hybrid assay and a GST pulldown assay. This self-interaction was found to be SDS sensitive in SDS-PAGE and native-PAGE assays. Moreover, the intracellular localization of these long polyalanine tracts was also observed to be disturbed. Our results suggest that long tracts of polyalanine acquire SDS-sensitive self-association properties, which may be a prerequisite event for their abnormal folding. The misfolding of these tracts is thought to be a common molecular aspect underlying the pathogenesis of polyalanine-related diseases.

Oculopharyngeal muscular dystrophy: Recent advances in the understanding of the molecular pathogenic mechanisms and treatment strategies

Oculopharyngeal muscular dystrophy (OPMD) is an adult-onset disorder characterized by progressive eyelid drooping, swallowing difficulties and proximal limb weakness. OPMD is caused by a small expansion of a short polyalanine tract in the poly (A) binding protein nuclear 1 protein (PABPN1). The mechanism by which the polyalanine expansion mutation in PABPN1 causes disease is unclear. PABPN1 is a nuclear multi-functional protein which is involved in pre-mRNA polyadenylation, transcription regulation, and mRNA nucleocytoplasmic transport. The distinct pathological hallmark of OPMD is the presence of filamentous intranuclear inclusions (INIs) in patient's skeletal muscle cells. The exact relationship between mutant PABPN1 intranuclear aggregates and pathology is not clear. OPMD is a unique disease sharing common pathogenic features with other polyalanine disorders, as well as with polyglutamine and dystrophic disorders. This chapter aims to review the rapidly growing body of knowledge concerning OPMD. First, we outline the background of OPMD. Second, we compare OPMD with other trinucleotide repeat disorders. Third, we discuss the recent advances in the understanding of the molecular mechanisms underlying OPMD pathogenesis. Finally, we review recent therapeutic strategies for OPMD.

Fast sequence evolution of Hox and Hox-derived genes in the genus Drosophila

Background

It is expected that genes that are expressed early in development and have a complex expression pattern are under strong purifying selection and thus evolve slowly. Hox genes fulfill these criteria and thus, should have a low evolutionary rate. However, some observations point to a completely different scenario. Hox genes are usually highly conserved inside the homeobox, but very variable outside it.

Results

We have measured the rates of nucleotide divergence and indel fixation of three Hox genes, labial (lab), proboscipedia (pb) and abdominal-A (abd-A), and compared them with those of three genes derived by duplication from Hox3, bicoid (bcd), zerknüllt (zen) and zerknüllt-related (zen2), and 15 non-Hox genes in sets of orthologous sequences of three species of the genus Drosophila. These rates were compared to test the hypothesis that Hox genes evolve slowly. Our results show that the evolutionary rate of Hox genes is higher than that of non-Hox genes when both amino acid differences and indels are taken into account: 43.39% of the amino acid sequence is altered in Hox genes, versus 30.97% in non-Hox genes and 64.73% in Hox-derived genes. Microsatellites scattered along the coding sequence of Hox genes explain partially, but not fully, their fast sequence evolution.

Conclusion

These results show that Hox genes have a higher evolutionary dynamics than other developmental genes, and emphasize the need to take into account indels in addition to nucleotide substitutions in order to accurately estimate evolutionary rates.

The human ortholog of the rodent testis-specific ABC transporter Abca17 is a ubiquitously expressed pseudogene (ABCA17P) and shares a common 5' end with ABCA3

BACKGROUND

During the past years, we and others discovered a series of human ATP-binding cassette (ABC) transporters, now referred to as ABC A-subfamily transporters. Recently, a novel testis-specific ABC A transporter, Abca17, has been cloned in rodent. In this study, we report the identification and characterization of the human ortholog of rodent Abca17.

RESULTS

The novel human ABC A-transporter gene on chromosome 16p13.3 is ubiquitously expressed with highest expression in glandular tissues and the heart. The new ABC transporter gene exhibits striking nucleotide sequence homology with the recently cloned mouse (58%) and rat Abca17 (51%), respectively, and is located in the syntenic region of mouse Abca17 indicating that it represents the human ortholog of rodent Abca17. However, unlike in the mouse, the full-length ABCA17 transcript (4.3 kb) contains numerous mutations that preclude its translation into a bona fide ABC transporter protein strongly suggesting that the human ABCA17 gene is a transcribed pseudogene (ABCA17P). We identified numerous alternative ABCA17P splice variants which are transcribed from two distinct transcription initiation sites. Genomic analysis revealed that ABCA17P borders on another ABC A-subfamily transporter - the lung surfactant deficiency gene ABCA3. Surprisingly, we found that both genes overlap at their first exons and are transcribed from opposite strands. This genomic colocalization and the observation that the ABCA17P and ABCA3 genes share significant homologies in several exons (up to 98%) suggest that both genes have evolved by gene duplication.

CONCLUSION

Our results demonstrate that ABCA17P and ABCA3 form a complex of overlapping genes in the human genome from which both non-coding and protein-coding ABC A-transporter RNAs are expressed. The fact that both genes overlap at their 5' ends suggests interdependencies in their regulation and may have important implications for the functional analysis of the disease gene ABCA3. Moreover, this is the first demonstration of the expression of a pseudogene and its parent gene from a common overlapping DNA region in the human genome.

Foxl2 function in ovarian development

Foxl2 is a forkhead transcription factor essential for proper reproductive function in females. Human patients carrying mutations in the FOXL2 gene display blepharophimosis/ptosis/epicanthus inversus syndrome (BPES), an autosomal dominant disease associated with eyelid defects and premature ovarian failure in females. Recently, animal models for BPES have been developed that in combination with a catalogue of human FOXL2 mutations provide further insight into its molecular function. Mice homozygous mutant for Foxl2 display craniofacial malformations and female infertility. The analysis of the murine phenotype has revealed that Foxl2 is required for granulosa cell function. These ovarian somatic cells surround and nourish the oocyte and play an important role in follicle formation and activation. Mutations upstream of FOXL2 in humans, not affecting the coding sequence itself, have also been shown to cause BPES, which points to the existence of a distant regulatory element necessary for proper gene expression. The same regulatory sequences may be deleted in the goat polled intersex syndrome (PIS), in which FoxL2 expression is severely reduced. Sequence comparison of FoxL2 from several vertebrate species has shown that it is a highly conserved gene involved in ovary development. Thus, the detailed understanding of Foxl2 function and regulation and the identification of its transcriptional targets may open new avenues for the treatment of female infertility in the future.

ELMOD2 is a candidate gene for familial idiopathic pulmonary fibrosis

We performed a genomewide scan in six multiplex families with familial idiopathic pulmonary fibrosis (IPF) who originated from southeastern Finland. The majority of the Finnish multiplex families were clustered in the region, and the population history suggested that the clustering might be explained by an ancestor shared among the patients. The genomewide scan identified five loci of interest. The hierarchical fine mapping in an extended data set with 24 families originating from the same geographic region revealed a shared 110 kb to 13 Mb haplotype on chromosome 4q31, which was significantly more frequent among the patients than in population-based controls (odds ratio 6.3; 95% CI 2.5-15.9; P = .0001). The shared haplotype harbored two functionally uncharacterized genes, ELMOD2 and LOC152586, of which only ELMOD2 was expressed in lung and showed significantly decreased messenger-RNA expression in IPF lung (n = 6) when compared with that of healthy lung (n = 7; P = .05). Our results suggest ELMOD2 as a novel candidate gene for susceptibility in familial IPF.

The silence RNA keeps: cis mechanisms of RNA mediated epigenetic silencing in mammals

One of the fundamental questions of modern biology is to unravel how genes are switched on and off at the right time and in the correct tissues. It is well recognized that gene regulation depends on a dynamic balance between activating and repressing forces, and multiple mechanisms are involved in both gene silencing and activation. Work over the last decade has revealed that in some cases transcriptional silencing of specific genes is mediated by RNAs that specifically recruit repressing complexes to homologous DNA sequences. Examples of both cis and trans RNA driven transcriptional silencing have been reported. This review focuses on those examples of transcriptional gene silencing in which the RNA component seems to act uniquely in cis. Speculative models of how such cis acting transcripts may trigger transcriptional silencing are proposed. Future experimental testing of these and other mechanisms is important to gain a fuller understanding of how genes are regulated and to identify instances in which such mechanisms are defective, leading to disease. Understanding the basic molecular basis of these phenomena will provide us with invaluable tools for the future development of targeted therapies and drugs for those diseases in which they are faulty.

Nested genes in the human genome

Here we studied one special type of gene, i.e., the nested gene, in the human genome. We collected 373 reliably annotated nested genes. Two-thirds of them were on the strand opposite that of their host gene. About 58% coding nested gene pairs were conserved in mouse and some were even maintained in chicken and fish, while nested pseudogenes were poorly conserved. Ka/Ks analysis revealed that nested genes were under strong selection, although they did not demonstrate greater conservation than other genes. With microarray data we observed that two partners of one nested pair seemed to be expressed reciprocally. A significant proportion of nested genes were tissue-specifically expressed. Gene ontology analysis demonstrated that quite a number of nested genes participated in cellular signal transduction. Based on these observations, we think that nested genes are a group of genes with important physiological functions.

Molecular consequences of PHOX2B missense, frameshift and alanine expansion mutations leading to autonomic dysfunction

Heterozygous mutations of the PHOX2B gene account for a broad variety of disorders of the autonomic nervous system, either isolated or combined, including congenital central hypoventilation syndrome (CCHS), tumours of the sympathetic nervous system and Hirschsprung disease. In CCHS, the prevalent mutation is an expansion of a 20-alanine stretch ranging from +5 to +13 alanines, whereas frameshift and missense mutations are found occasionally. To determine the molecular basis of impaired PHOX2B function, we assayed the transactivation and DNA binding properties of wild-type and mutant PHOX2B proteins. Furthermore, we investigated aggregate formation by proteins with polyalanine tract expansions ranging from +5 to +13 alanines using immunofluorescence of transfected cells and gel filtration of in vitro translated proteins. We found that transactivation of the dopamine beta-hydroxylase promoter by PHOX2B proteins with frameshift and missense mutations was abolished or severely curtailed, as was in vitro DNA binding although the proteins localized to the nucleus. The transactivation potential of proteins with polyalanine tract expansions declined with increasing length of the polyalanine stretch, and DNA binding was affected for an expansion of +9 alanines and above. Cytoplasmic aggregation in transfected cells was only observed for the longest expansions, whereas even the short expansion mutants were prone to form multimers in vitro. Such a tendency to protein misfolding could explain loss of transactivation for alanine expansion mutations. However, additional mechanisms such as toxic gain-of-function may play a role in the pathogenic process.

Association study of polymorphisms in synaptic vesicle-associated genes, SYN2 and CPLX2, with schizophrenia

BACKGROUND

The occurrence of aberrant functional connectivity in the neuronal circuit is one of the integrative theories of the etiology of schizophrenia. Previous studies have reported that the protein and mRNA levels of the synapsin 2 (SYN2) and complexin 2 (CPLX2) genes were decreased in patients with schizophrenia. Synapsin 2 and complexin 2 are involved in synaptogenesis and the modulation of neurotransmitter release. This report presents a study of the association of polymorphisms of SYN2 and CPLX2 with schizophrenia in the Korean population.

METHODS

Six single nucleotide polymorphisms (SNPs) and one 5-bp insertion/deletion in SYN2 and five SNPs in CPLX2 were genotyped in 154 Korean patients with schizophrenia and 133 control patients using direct sequencing or restriction fragment length polymorphism analysis. An intermarker linkage disequilibrium map was constructed for each gene.

RESULTS

Although there was no significant difference in the genotypic distributions and allelic frequencies of either SYN2 or CPLX2 polymorphisms between the schizophrenia and control groups, the two-way haplotype analyses revealed significant associations with the disease (P < 0.05 after Bonferroni correction). The three-way haplotype analyses also revealed a significant association of SYN2 with schizophrenia (P < 0.001 after Bonferroni correction).

CONCLUSION

These results suggest that both SYN2 and CPLX2 may confer susceptibility to schizophrenia in the Korean population.

Overlapping genes in vertebrate genomes

Overlapping genes in mammalian genomes are unexpected phenomena even though hundreds of pairs of protein coding overlapping genes have been reported so far. Overlapping genes can be divided into different categories based on direction of transcription as well as on sequence segments being shared between overlapping coding regions. The biologic functions of natural antisense transcripts, their involvement in physiological processes and gene regulation in living organisms are not fully understood. Number of documented examples indicates that they may exert control at various levels of gene expression, such as transcription, mRNA processing, splicing, stability, transport, and translation. Similarly, evolutionary origin of such genes is not known, existing hypotheses can explain only selected cases of mammalian gene overlaps which could originate as result of rearrangements, overprinting and/or adoption of signals in the neighboring gene locus.

Hyper-expansion of asparagines correlates with an abundance of proteins with prion-like domains in Plasmodium falciparum

Plasmodium falciparum encodes approximately 5300 proteins of which approximately 35% have repeats of amino acids, significantly higher than in other fully sequenced eukaryotes. The proportion of proteins with amino acid homorepeats varies from 4 to 54% amongst different functional classes of proteins. These homorepeats are dominated by asparagines, which are selected over lysines despite equivalent AT codon content. Surprisingly, asparagine repeats are absent from the variant surface antigen protein families of PfEMP1s, Stevors and Rifins. The PfEMP1 protein family is instead rich in recurrences of glutamates, similar to human cell surface proteins. Structural mapping of homorepeats suggests that these segments are likely to form surface exposed structures that protrude from the main protein cores. We also found an abundance of asparagine-rich prion-like domains in P. falciparum, significantly larger than in any other eukaryote. Domains rich in glutamines and asparagines have an innate predisposition to form self-propagating amyloid fibers, which are involved both in prion-based inheritance and in human neurodegenerative disorders. Nearly 24% (1302 polypeptides) of P. falciparum proteins contain prion-forming or prion-inducing domains, in comparison to Drosophila (approximately 3.4%) which to date showed the highest number of prion-like proteins. The unexpected properties of P. falciparum revealed here open new avenues for investigating parasite biology.

Speeding disease gene discovery by sequence based candidate prioritization

Background

Regions of interest identified through genetic linkage studies regularly exceed 30 centimorgans in size and can contain hundreds of genes. Traditionally this number is reduced by matching functional annotation to knowledge of the disease or phenotype in question. However, here we show that disease genes share patterns of sequence-based features that can provide a good basis for automatic prioritization of candidates by machine learning.

Results

We examined a variety of sequence-based features and found that for many of them there are significant differences between the sets of genes known to be involved in human hereditary disease and those not known to be involved in disease. We have created an automatic classifier called PROSPECTR based on those features using the alternating decision tree algorithm which ranks genes in the order of likelihood of involvement in disease. On average, PROSPECTR enriches lists for disease genes two-fold 77% of the time, five-fold 37% of the time and twenty-fold 11% of the time.

Conclusion

PROSPECTR is a simple and effective way to identify genes involved in Mendelian and oligogenic disorders. It performs markedly better than the single existing sequence-based classifier on novel data. PROSPECTR could save investigators looking at large regions of interest time and effort by prioritizing positional candidate genes for mutation detection and case-control association studies.

Family-based association study of synapsin II and schizophrenia

Synapsin II has been proposed as a candidate gene for vulnerability to schizophrenia on the basis of its function and its location in a region of the genome implicated by linkage studies in families with schizophrenia. We recently reported positive association of synapsin II with schizophrenia in a case-control study (Chen et al. 2004). However, since case-control analyses can generate false-positive results in the presence of minor degrees of population stratification, we have performed a replication study in 366 additional Han Chinese probands and their parents by use of analyses of transmission/disequilibrium for three in/del markers and three single-nucleotide polymorphisms. Positive association was observed for rs2307981 (P =.02), rs2308169 (P =.005), rs308963 (P =.002), rs795009 (P =.02), and rs2307973 (P =.02). For transmission of six-marker haplotypes, the global P value was.0000016 (5 degrees of freedom), principally because of overtransmission of the most common haplotype, CAA/-/G/T/C/- (frequency 53.6%; chi (2) = 20.8; P =.0000051). This confirms our previous study and provides further support for the role of synapsin II variants in susceptibility to schizophrenia.

Positive association between synapsin II and schizophrenia

BACKGROUND

Synapsin II encodes a neuron-specific phosphoprotein that selectively binds to small synaptic vesicles in the presynaptic nerve terminal. The expressions of messenger ribonucleic acid and protein of synapsin II have been reported to be significantly reduced in the brains of schizophrenia patients. The synapsin II gene is located on 3p25, a region that has been implicated to be associated with schizophrenia by genetic linkage. All these findings suggest synapsin II as a candidate gene for schizophrenia.

METHODS

In this work, we studied four markers (two single nucleotide polymorphisms (SNPs): rs308963 and rs795009; and two insertion/deletion polymorphisms: rs2307981 and rs2308169) covering 144.2 kilobase pairs (kb) with an average interval of 38 kb in synapsin II in a sample of 654 schizophrenic patients and 628 normal control subjects to explore the mechanism underlying schizophrenia.

RESULTS

We found significant differences in allele frequency distribution of SNP rs795009 (p =.000018, odds ratio 1.405, 95% confidence interval 1.202-1.641) between patients and control subjects. The T allele was significantly higher in patients than in control subjects. Moreover, the overall frequency of haplotype showed significant differences between patients and control subjects (p <.000001).

CONCLUSIONS

This study suggests a positive association between synapsin II and schizophrenia, implying that synapsin II is involved in the etiology of schizophrenia.

Proteomic signatures: amino acid and oligopeptide compositions differentiate among phyla

Availability of complete genome sequences allows in-depth comparison of single-residue and oligopeptide compositions of the corresponding proteomes. We have used principal component analysis (PCA) to study the landscape of compositional motifs across more than 70 genera from all three superkingdoms. Unexpectedly, the first two principal components clearly differentiate archaea, eubacteria, and eukaryota from each other. In particular, we contrast compositional patterns typical of the three superkingdoms and characterize differences between species and phyla, as well as among patterns shared by all compositional proteomic signatures. These species-specific patterns may even extend to subsets of the entire proteome, such as proteins pertaining to individual yeast chromosomes. We identify factors that affect compositional signatures, such as living habitat, and detect strong eukaryotic preference for homopeptides and palindromic tripeptides. We further detect oligopeptides that are either universally over- or underabundant across the whole proteomic landscape, as well as oligopeptides whose over- or underabundance is phylum- or species-specific. Finally, we report that species composition signatures preserve evolutionary memory, providing a new method to compare phylogenetic relationships among species that avoids problems of sequence alignment and ortholog detection.

Human disease genes: patterns and predictions

We compared genes at which mutations are known to cause human disease (disease genes) with other human genes (nondisease genes) using a large set of human-rodent alignments to infer evolutionary patterns. Such comparisons may be of use both in predicting disease genes and in understanding the general evolution of human genes. Four features were found to differ significantly between disease and nondisease genes, with disease genes (i) evolving with higher nonsynonymous/synonymous substitution rate ratios (Ka/Ks), (ii) evolving at higher synonymous substitution rates, (iii) with longer protein-coding sequences, and (iv) expressed in a narrower range of tissues. Discriminant analysis showed that these differences may help to predict human disease genes. We also investigated other factors affecting the mode of evolution in the disease genes: Ka/Ks is significantly affected by protein function, mode of inheritance, and the reduction of life expectancy caused by disease.

Genome sequence and splice site analysis of low-fidelity DNA polymerases H and I involved in replication of damaged DNA

POLH and POLI are paralogs encoding low-fidelity, class Y, DNA polymerases involved in replication of damaged DNA in the human disease xeroderma pigmentosum variant. Analysis of genomic regions for human and mouse homologs, employing the analytic tool Genome Cryptographer, detected low-repetitive or unique regions at exons and other potential control regions, especially within intron I of human POLH. The human and mouse homologs are structurally similar, but the paralogs have undergone evolutionary divergence. The information content of splice sites for human POLH, the probability that a base would contribute to splicing, was low only for the acceptor site of exon II, which is preceded by a region of high information content that could contain sequences controlling splicing. This analysis explains previous observations of tissue-specific skipping during mRNA processing, resulting in the loss of the transcription start site in exon II, in human tissues.