Reconstruction of a functional human gene network, with an application for prioritizing positional candidate genes

Most common genetic disorders have a complex inheritance and may result from variants in many genes, each contributing only weak effects to the disease. Pinpointing these disease genes within the myriad of susceptibility loci identified in linkage studies is difficult because these loci may contain hundreds of genes. However, in any disorder, most of the disease genes will be involved in only a few different molecular pathways. If we know something about the relationships between the genes, we can assess whether some genes (which may reside in different loci) functionally interact with each other, indicating a joint basis for the disease etiology. There are various repositories of information on pathway relationships. To consolidate this information, we developed a functional human gene network that integrates information on genes and the functional relationships between genes, based on data from the Kyoto Encyclopedia of Genes and Genomes, the Biomolecular Interaction Network Database, Reactome, the Human Protein Reference Database, the Gene Ontology database, predicted protein-protein interactions, human yeast two-hybrid interactions, and microarray co-expressions. We applied this network to interrelate positional candidate genes from different disease loci and then tested 96 heritable disorders for which the Online Mendelian Inheritance in Man database reported at least three disease genes. Artificial susceptibility loci, each containing 100 genes, were constructed around each disease gene, and we used the network to rank these genes on the basis of their functional interactions. By following up the top five genes per artificial locus, we were able to detect at least one known disease gene in 54% of the loci studied, representing a 2.8-fold increase over random selection. This suggests that our method can significantly reduce the cost and effort of pinpointing true disease genes in analyses of disorders for which numerous loci have been reported but for which most of the genes are unknown.     

The stationary distribution of allele frequencies when selection acts at unlinked loci

We consider population genetics models where selection acts at a set of unlinked loci. It is known that if the fitness of an individual is multiplicative across loci, then these loci are independent. We consider general selection models, but assume parent-independent mutation at each locus. For such a model, the joint stationary distribution of allele frequencies is proportional to the stationary distribution under neutrality multiplied by a known function of the mean fitness of the population. We further show how knowledge of this stationary distribution enables direct simulation of the genealogy of a sample at a single-locus. For a specific selection model appropriate for complex disease genes, we use simulation to determine what features of the genealogy differ between our general selection model and a multiplicative model.     

DIFFERENTIATION OF A MULTIGENE FAMILY BETWEEN POPULATIONS

Differentiation of a multigene family between populations is investigated under the forces of intrachromosomal unbiased gene conversion, neutral mutation, recombination, and random drift. The infinite-site model without recombination within the gene is employed, and the average numbers of the different sites between pairs of genes at the same locus and at different loci from two populations that divided some time ago are calculated as functions of time. Since differentiation is a relative concept in which differences between populations are expressed in comparison with those within a population and since there is most likely a reduction in population size after the division of populations, the time-dependent behavior of the corresponding quantities within a population is also investigated using approximate and numerical methods. It is found that intrachromosomal unbiased gene conversion promotes differentiation at the same locus but promotes or retards that between different loci, depending on the parameters in a recently divided population.     

Genetic identity between populations when mutation rates vary within and across loci

A formula is derived for the probability that two genes taken at random from the same locus in two populations isolated at time t ago are of the same allelic type. The model assumed is a neutral one where there are possibly different mutation rates between different alleles. Inequalities are derived for this probability. A particular result is that for a fixed overall mutation rate, the probability is least for the infinite alleles model. Inequalities and approximations are found for Nei's genetic identity at one locus when mutation rates vary, and also for the identity across loci when the overall mutation rates per locus vary. Genetic identity at the molecular level is considered and a probability generating function found for the number of segregating sites between two randomly chosen gametes from two divergent populations, under various models.     

Utility and efficiency of linked marker genes for genetic counseling. III. Proportion of informative families under linkage disequilibrium.

A marker locus closely linked to a disease locus is often useful for genetic counseling provided that a counselee is heterozygous at both disease and marker loci. Furthermore, the linkage phase of these genes in the counselee must be known. When the linkage between the disease and marker loci is very close, one often finds linkage disequilibrium between the loci. To evaluate the effect of such nonrandom associations on the utility of linked marker genes for genetic counseling, the proportion of informative families is studied for X-linked recessive and autosomal dominant diseases. This proportion is higher for X-linked genes than for autosomal genes, if other factors are the same. In general, codominant markers are more useful than dominant markers. Also, under appropriate conditions, the proportion of informative families is higher when linkage disequilibrium is present. The results obtained in this paper are useful for evaluating the utility of polymorphic restriction endonuclease cleavage sites as markers in genetic counseling.     

A Novel Locus for Ectodermal Dysplasia of Hair,Nail and Skin Pigmentation Anomalies Maps to Chromosome 18p11.32-p11.31

Ectodermal dysplasias (EDs) are a large heterogeneous group of inherited disorders exhibiting abnormalities in ectodermally derived appendages such as hair, nails, teeth and sweat glands. EDs associated with reticulated pigmentation phenotype are rare entities for which the genetic basis and pathophysiology are not well characterized. The present study describes a five generation consanguineous Pakistani family segregating an autosomal recessive form of a novel type of ectodermal dysplasia. The affected members present with sparse and woolly hair, severe nail dystrophy and reticulate skin pigmentation. After exclusion of known gene loci related with other skin disorders, genome-wide linkage analysis was performed using Illumina HumanOmniExpress beadchip SNP arrays. We linked this form of ED to human chromosome 18p11.32-p11.31 flanked by the SNPs rs9284390 (0.113Mb) and rs4797100 (3.14 Mb). A maximum two-point LOD score of 3.3 was obtained with several markers along the disease interval. The linkage interval of 3.03 Mb encompassed seventeen functional genes. However, sequence analysis of all these genes did not discover any potentially disease causing-variants. The identification of this novel locus provides additional information regarding the mapping of a rare form of ED. Further research, such as the use of whole-genome sequencing, would be expected to reveal any pathogenic mutation within the disease locus.     

Gene organization of the mating type regions in the ectomycorrhizal fungus Laccaria bicolor reveals distinct evolution between the two mating type loci

In natural conditions, basidiomycete ectomycorrhizal fungi such as Laccaria bicolor are typically in the dikaryotic state when forming symbioses with trees, meaning that two genetically different individuals have to fuse or 'mate'. Nevertheless, nothing is known about the molecular mechanisms of mating in these ecologically important fungi. Here, advantage was taken of the first sequenced genome of the ectomycorrhizal fungus, Laccaria bicolor, to determine the genes that govern the establishment of cell-type identity and orchestrate mating. The L. bicolor mating type loci were identified through genomic screening. The evolutionary history of the genomic regions that contained them was determined by genome-wide comparison of L. bicolor sequences with those of known tetrapolar and bipolar basidiomycete species, and by phylogenetic reconstruction of gene family history. It is shown that the genes of the two mating type loci, A and B, are conserved across the Agaricales, but they are contained in regions of the genome with different evolutionary histories. The A locus is in a region where the gene order is under strong selection across the Agaricales. By contrast, the B locus is in a region where the gene order is likely under a low selection pressure but where gene duplication, translocation and transposon insertion are frequent.     

Mutations in STIL, Encoding a Pericentriolar and Centrosomal Protein, Cause Primary Microcephaly

Primary microcephaly (MCPH) is an autosomal-recessive congenital disorder characterized by smaller-than-normal brain size and mental retardation. MCPH is genetically heterogeneous with six known loci: MCPH1–MCPH6. We report mapping of a novel locus, MCPH7, to chromosome 1p32.3–p33 between markers D1S2797 and D1S417, corresponding to a physical distance of 8.39 Mb. Heterogeneity analysis of 24 families previously excluded from linkage to the six known MCPH loci suggested linkage of five families (20.83%) to the MCPH7 locus. In addition, four families were excluded from linkage to the MCPH7 locus as well as all of the six previously known loci, whereas the remaining 15 families could not be conclusively excluded or included. The combined maximum two-point LOD score for the linked families was 5.96 at marker D1S386 at θ = 0.0. The combined multipoint LOD score was 6.97 between markers D1S2797 and D1S417. Previously, mutations in four genes, MCPH1, CDK5RAP2, ASPM, and CENPJ, that code for centrosomal proteins have been shown to cause this disorder. Three different homozygous mutations in STIL, which codes for a pericentriolar and centrosomal protein, were identified in patients from three of the five families linked to the MCPH7 locus; all are predicted to truncate the STIL protein. Further, another recently ascertained family was homozygous for the same mutation as one of the original families. There was no evidence for a common haplotype. These results suggest that the centrosome and its associated structures are important in the control of neurogenesis in the developing human brain.     

Identification of a novel SCA locus ( SCA19) in a Dutch autosomal dominant cerebellar ataxia family on chromosome region 1p21-q21

We present a linkage study in a four-generation autosomal dominant cerebellar ataxia (ADCA) family of Dutch ancestry. The family shows a clinically and genetically distinct form of ADCA. This neurodegenerative disorder manifests in the family as a relatively mild ataxia syndrome with some additional characteristic symptoms. We have identified a SCA19 locus, approved by the Human Genome Nomenclature Committee that can be assigned to the chromosome region 1p21-q21. Our mutation analysis failed to identify any mutations in the known spinocerebellar ataxia ( SCA) genes and linkage analysis excluded the remaining SCA loci. We therefore performed a genome-wide scan with 350 microsatellite markers to identify the location of the disease-causing gene in this family. Multi-point analysis was performed and exclusion maps were generated. Linkage and haplotype analysis revealed linkage to an interval located on chromosome 1. The estimated minimal prevalence of ADCA in the Netherlands is about 3:100,000. To date, sixteen different SCA loci have been identified in ADCA ( SCA1-8 and SCA10-17). However, mutation analysis has been commercially available only for the SCA1, 2, 3, 6 and 7 genes. So far, a molecular analysis in these SCA genes cannot be made in about one-third of the ADCA families. Thus, the identification of this new, additional SCA19 locus will contribute to expanding the DNA diagnostic possibilities.     

Genotyping microsatellite DNA markers at putative disease loci in inbred/multiplex families with respiratory chain complex I deficiency allows rapid identification of a novel nonsense mutation (IVS1nt -1) in the NDUFS4 gene in Leigh syndrome

Complex I deficiency, the most common cause of mitochondrial disorders, accounts for a variety of clinical symptoms and its genetic heterogeneity makes identification of the disease genes particularly tedious. Indeed, most of the 43 complex I subunits are encoded by nuclear genes, only seven of them being mitochondrially encoded. In order to offer urgent prenatal diagnosis, we have studied an inbred/multiplex family with complex I deficiency by using microsatellite DNA markers flanking the putative disease loci. Microsatellite DNA markers have allowed us to exclude the NDUFS7, NDUFS8, NDUFV1 and NDUFS1 genes and to find homozygosity at the NDUFS4 locus. Direct sequencing has led to identification of a homozygous splice acceptor site mutation in intron 1 of the NDUFS4 gene (IVS1nt -1, G-->A); this was not found in chorion villi of the ongoing pregnancy. We suggest that genotyping microsatellite DNA markers at putative disease loci in inbred/multiplex families helps to identify the disease-causing mutation. More generally, we suggest giving consideration to a more systematic microsatellite analysis of putative disease loci for identification of disease genes in inbred/multiplex families affected with genetically heterogeneous conditions.     

Multipoint linkage analysis and heterogeneity testing in 20 X-linked retinitis pigmentosa families

Using multipoint linkage analysis in 20 families segregating for X-linked retinitis pigmentosa (XLRP), the lod scores on a map of eight RFLP loci were obtained. Our results indicate that under the hypothesis of homogeneity the maximal multipoint lod score supports one disease locus located slightly distal to OTC at Xp21.1. Heterogeneity testing for two XLRP loci suggested that a second XLRP locus may be located 8.5 cM proximal to DXS28 at Xp21.3. Further heterogeneity testing for three disease loci failed to detect a third XLRP locus proximal to DXS7 in any of our 20 XLRP families.     

Neural networks for modeling gene-gene interactions in association studies

Background

The monogenic disease osteogenesis imperfecta (OI) is due to single mutations in either of the collagen genes ColA1 or ColA2, but within the same family a given mutation is accompanied by a wide range of disease severity. Although this phenotypic variability implies the existence of modifier gene variants, genome wide scanning of DNA from OI patients has not been reported. Promising genome wide marker-independent physical methods for identifying disease-related loci have lacked robustness for widespread applicability. Therefore we sought to improve these methods and demonstrate their performance to identify known and novel loci relevant to OI.

Results

We have improved methods for enriching regions of identity-by-descent (IBD) shared between related, afflicted individuals. The extent of enrichment exceeds 10- to 50-fold for some loci. The efficiency of the new process is shown by confirmation of the identification of the Col1A2 locus in osteogenesis imperfecta patients from Amish families. Moreover the analysis revealed additional candidate linkage loci that may harbour modifier genes for OI; a locus on chromosome 1q includes COX-2, a gene implicated in osteogenesis.

Conclusion

Technology for physical enrichment of IBD loci is now robust and applicable for finding genes for monogenic diseases and genes for complex diseases. The data support the further investigation of genetic loci other than collagen gene loci to identify genes affecting the clinical expression of osteogenesis imperfecta. The discrimination of IBD mapping will be enhanced when the IBD enrichment procedure is coupled with deep resequencing.     

Using Linkage Analysis to Detect Gene-Gene Interaction by Stratifying Family Data on Known Disease,or Disease-Associated,Alleles

Detecting gene-gene interaction in complex diseases is a major challenge for common disease genetics. Most interaction detection approaches use disease-marker associations and such methods have low power and unknown reliability in real data. We developed and tested a powerful linkage-analysis-based gene-gene interaction detection strategy based on conditioning the family data on a known disease-causing allele or disease-associated marker allele. We computer-generated multipoint linkage data for a disease caused by two epistatically interacting loci (A and B). We examined several two-locus epistatic inheritance models: dominant-dominant, dominant-recessive, recessive-dominant, recessive-recessive. At one of the loci (A), there was a known disease-related allele. We stratified the family data on the presence of this allele, eliminating family members who were without it. This elimination step has the effect of raising the “penetrance” at the second locus (B). We then calculated the lod score at the second locus (B) and compared the pre- and post-stratification lod scores at B. A positive difference indicated interaction. We also examined if it was possible to detect interaction with locus B based on a disease-marker association (instead of an identified disease allele) at locus A. We also tested whether the presence of genetic heterogeneity would generate false positive evidence of interaction. The power to detect interaction for a known disease allele was 60–90%. The probability of false positives, based on heterogeneity, was low. Decreasing linkage disequilibrium between the disease and marker at locus A decreased the likelihood of detecting interaction. The allele frequency of the associated marker made little difference to the power.     

林木全同胞群体的多位点连锁分析

多位点连锁分析是构建人类以及动植物的遗传连锁图谱的关键步骤之一。但是由于林木遗传背景的复杂性,多位点连锁分析在林木的全同胞群体中还没有得到应用．本文将多位点连锁分析应用到林木的F1代全同胞群体中．对于全同胞群体的任意分离比的两个位点,给出了在不同的连锁相下从一个位点到另一个位点的转移概率矩阵．对于给定的一列标记位点,考虑了不同分离比位点以及两位点间的连锁相信息,采用隐马尔可夫链模型计算极大似然函数和相邻位点间的重组率．本文的方法有助于构建完整的高密度的林木遗传连锁图谱．     

Utility and efficiency of linked marker genes for genetic counseling. II. Identification of linkage phase by offspring phenotypes

For a linked marker locus to be useful for genetic counseling, the counselee must be heterozygous for both disease and marker loci and his or her linkage phase must be known. It is shown that when the phenotypes of the counselee's previous children for the disease and marker loci are known, the linkage phase can often be inferred with a high probability, and thus it is possible to conduct genetic counseling. To evaluate the utility of linked marker genes for genetic counseling, the accuracy of prediction of the risk for a prospective child with a given marker gene to develop the genetic disease and the proportion of families in which a particular marker locus can be used for genetic counseling are studied for X-linked recessive, autosomal dominant, and autosomal recessive diseases. In the case of X-linked genetic diseases, information from children is very useful for determining the linkage phase of the counselee and predicting the genetic disease. In the case of autosomal dominant diseases, not all children are informative, but if the number of children is large, the phenotypes of children are often more informative than the information from grandparents. In the case of autosomal recessive diseases, information from grandparents is usually useless, since they show a normal phenotype for the disease locus. If we use information on the phenotypes of children, however, the linkage phase of the counselee and the risk of a prospective child can be inferred with a high probability. The proportion of informative families depends on the dominance relationship and frequencies of marker alleles, and the number of children. In general, codominant markers are more useful than are dominant markers, and a locus with high heterozygosity is more useful than is a locus with low heterozygosity.     

CtIP Mutations Cause Seckel and Jawad Syndromes

Seckel syndrome is a recessively inherited dwarfism disorder characterized by microcephaly and a unique head profile. Genetically, it constitutes a heterogeneous condition, with several loci mapped (SCKL1-5) but only three disease genes identified: the ATR, CENPJ, and CEP152 genes that control cellular responses to DNA damage. We previously mapped a Seckel syndrome locus to chromosome 18p11.31-q11.2 (SCKL2). Here, we report two mutations in the CtIP (RBBP8) gene within this locus that result in expression of C-terminally truncated forms of CtIP. We propose that these mutations are the molecular cause of the disease observed in the previously described SCKL2 family and in an additional unrelated family diagnosed with a similar form of congenital microcephaly termed Jawad syndrome. While an exonic frameshift mutation was found in the Jawad family, the SCKL2 family carries a splicing mutation that yields a dominant-negative form of CtIP. Further characterization of cell lines derived from the SCKL2 family revealed defective DNA damage induced formation of single-stranded DNA, a critical co-factor for ATR activation. Accordingly, SCKL2 cells present a lowered apoptopic threshold and hypersensitivity to DNA damage. Notably, over-expression of a comparable truncated CtIP variant in non-Seckel cells recapitulates SCKL2 cellular phenotypes in a dose-dependent manner. This work thus identifies CtIP as a disease gene for Seckel and Jawad syndromes and defines a new type of genetic disease mechanism in which a dominant negative mutation yields a recessively inherited disorder.     

Screening for homozygosity by descent in families with autosomal recessive retinitis pigmentosa

Retinitis pigmentosa (RP) is a genetically heterogeneous disease and an important cause of blindness in the state of Andhra Pradesh in India. In an attempt to identify the disease locus in families with the recessive form of the disease, we used the approach of screening for homozygosity by descent in offspring of consanguineous and nonconsanguineous families with RP. Microsatellite markers closely flanking 21 known candidate genes for RP were genotyped in parents and affected offspring to determine whether there was homozygosity at these loci that was shared by affected individuals of a family. This screening approach may be a rapid preliminary method to test known loci for possible cosegregation with disease.     

Identification, by homozygosity mapping, of a novel locus for autosomal recessive congenital ichthyosis on chromosome 17p, and evidence for further genetic heterogeneity

Autosomal recessive congenital ichthyosis (ARCI) comprises a group of severe disorders of keratinization, characterized by variable erythema and skin scaling. It is known for its high degree of genetic and clinical heterogeneity. Mutations in the gene for keratinocyte transglutaminase (TGM1) on chromosome 14q11 were shown in patients with ARCI, and a second locus was described, on chromosome 2q, in families from northern Africa. Three other loci for ARCI, on chromosomes 3p and 19p, were identified recently. We have embarked on a whole-genome scan for further loci for ARCI in four families from Germany, Turkey, and the United Arab Emirates. A novel ARCI locus was identified on chromosome 17p, between the markers at D17S938 and D17S1856, with a maximum LOD score of 3.38, at maximum recombination fraction 0.00, at D17S945, under heterogeneity. This locus is linked to the disease in the Turkish family and in the German family. Extensive genealogical studies revealed that the parents of the German patients with ARCI were eighth cousins. By homozygosity mapping, the localization of the gene could then be refined to the 8.4-cM interval between D17S938 and D17S1879. It could be shown, however, that ARCI in the two Arab families is linked neither to the new locus on chromosome 17p nor to one of the five loci known previously. Our findings give evidence of further genetic heterogeneity that is not linked to distinctive phenotypes.     

Mismatch Repair Genes on Chromosomes 2p and 3p Account for a Major Share of Hereditary Nonpolyposis Colorectal Cancer Families Evaluable by Linkage

Two susceptibility loci for hereditary nonpolyposis colo-rectal cancer (HNPCC) have been identified, and each contains a mismatch repair gene: MSH2 on chromosome 2p and MLH1 on chromosome 3p. We studied the involvement of these loci in 13 large HNPCC kindreds originating from three different continents. Six families showed close linkage to the 2p locus, and a heritable mutation of the MSH2 gene was subsequently found in four. The 2p-linked kindreds included a family characterized by the lack of extracolonic manifestations (Lynch I syndrome), as well as two families with cutaneous manifestations typical of the Muir-Torre syndrome. Four families showed evidence for linkage to the 3p locus, and a heritable mutation of the MLH1 gene was later detected in three. One 3p-linked kindred was of Amerindian origin. Of the remaining three families studied for linkage, one showed lod scores compatible with exclusion of both MSH2 and MLH1, while lod scores obtained in the other two families suggested exclusion of one HNPCC locus (MSH2 or MLH1) but were uninformative for markers flanking the other locus. Our results suggest that mismatch repair genes on 2p and 3p account for a major share of HNPCC in kindreds that can be evaluated by linkage analysis.     

Genetic linkage studies of autosomal dominant polycystic kidney disease: search for the second gene in a large Sicilian family

Summary Polycystic kidney disease (PKD) is a common autosomal dominant genetic disorder caused by mutation in at least two different gene loci. The PKD1 gene has been localized on the short arm of chromosome 16. The location of a second genetic locus in the human genome is not yet known. A large PKD kindred, unlinked to chromosome 16, with over 250 members was studied using both DNA and classical markers. In total, 29 informative marker loci on 11 autosomes have been analyzed for linkage with PKD. The data significantly exclude the linkage with disease locus from 17 marker loci and show no evidence of close linkage with the other loci.