Strategies for genetic model specification in the screening of genome-wide meta-analysis signals for further replication

Background Meta-analysis is increasingly being employed as a screening procedure in large-scale association studies to select promising variants for follow-up studies. However, standard methods for meta-analysis require the assumption of an underlying genetic model, which is typically unknown a priori. This drawback can introduce model misspecifications, causing power to be suboptimal, or the evaluation of multiple genetic models, which augments the number of false-positive associations, ultimately leading to waste of resources with fruitless replication studies. We used simulated meta-analyses of large genetic association studies to investigate naive strategies of genetic model specification to optimize screenings of genome-wide meta-analysis signals for further replication. Methods Different methods, meta-analytical models and strategies were compared in terms of power and type-I error. Simulations were carried out for a binary trait in a wide range of true genetic models, genome-wide thresholds, minor allele frequencies (MAFs), odds ratios and between-study heterogeneity (tau(2)). Results Among the investigated strategies, a simple Bonferroni-corrected approach that fits both multiplicative and recessive models was found to be optimal in most examined scenarios, reducing the likelihood of false discoveries and enhancing power in scenarios with small MAFs either in the presence or in absence of heterogeneity. Nonetheless, this strategy is sensitive to tau(2) whenever the susceptibility allele is common (MAF epsilon 30%), resulting in an increased number of false-positive associations compared with an analysis that considers only the multiplicative model. Conclusion Invoking a simple Bonferroni adjustment and testing for both multiplicative and recessive models is fast and an optimal strategy in large meta-analysis-based screenings. However, care must be taken when examined variants are common, where specification of a multiplicative model alone may be preferable.

A high-resolution genetic map of the familial Mediterranean fever candidate region allows identification of haplotype-sharing among ethnic groups

Familial Mediterranean fever (FMF) is a recessive disorder of inflammation caused by mutations in a gene (designated MEFV) on chromosome 16p13.3, We have recently constructed a 1-Mb cosmid contig that includes the FMF critical region. Here we show genotype data for 12 markers from our physical map, including 5 newly identified microsatellites, in FMF families. Intrafamilial recombinations placed MEFV in the similar to 285 kb between D16S468/D16S3070 and D16S3376. We observed significant linkage disequilibrium in the North African Jewish population, and historical recombinants in the founder haplotype placed MEFV between D16S3082 and D16S3373 (similar to 200 kb). In smaller panels of Iraqi Jewish, Arab, and Armenian families, there were significant allelic associations only for D16S3370 and D16S2617 among the Armenians. A sizable minority of Iraqi Jewish and Armenian carrier chromosomes appeared to be derived from the North African Jewish ancestral haplotype. We observed a unique FMF haplotype common to Iraqi Jews, Arabs, and Armenians and two other haplotypes restricted to either the Iraqi Jewish or the Armenian population. These data support the view that a few major mutations account for a large percentage of the cases of FMF and suggest that same of these mutations arose before the affected Middle Eastern populations diverged from one another. (C) 1997 Academic Press.

Recessive COL6A2 C-globular Missense Mutations in Ullrich Congenital Muscular Dystrophy ROLE OF THE C2a SPLICE VARIANT

Ullrich congenital muscular dystrophy (UCMD) is a disabling and life-threatening disorder resulting from either recessive or dominant mutations in genes encoding collagen VI. Although the majority of the recessive UCMD cases have frameshift or nonsense mutations in COL6A1, COL6A2, or COL6A3, recessive structural mutations in the COL6A2 C-globular region are emerging also. However, the underlying molecular mechanisms have remained elusive. Here we identified a homozygous COL6A2 E624K mutation (C1 subdomain) and a homozygous COL6A2 R876S mutation (C2 subdomain) in two UCMD patients. The consequences of the mutations were investigated using fibroblasts from patients and cells stably transfected with the mutant constructs. In contrast to expectations based on the clinical severity of these two patients, secretion and assembly of collagen VI were moderately affected by the E624K mutation but severely impaired by the R876S substitution. The E624K substitution altered the electrostatic potential of the region surrounding the metal ion-dependent adhesion site, resulting in a collagen VI network containing thick fibrils and spots with densely packed microfibrils. The R876S mutation prevented the chain from assembling into triple-helical collagen VI molecules. The minute amount of collagen VI secreted by the R876S fibroblasts was solely composed of a faster migrating chain corresponding to the C2a splice variant with an alternative C2 subdomain. In transfected cells, the C2a splice variant was able to assemble into short microfibrils. Together, the results suggest that the C2a splice variant may functionally compensate for the loss of the normal COL6A2 chain when mutations occur in the C2 subdomain.

Hepatitis Delta virus genotype 8 infection in Northeast Brazil: Inheritance from African slaves?

Hepatitis Delta virus (HDV) is endemic worldwide, but its prevalence varies in different geographical areas. While in the Brazilian Amazon, HDV is known to be endemic and to represent a significant public health problem, few studies have assessed its prevalence in other regions in the country. This study evaluated the seroprevalence of HDV among HBsAg chronic carriers from Maranhao state, a region located in the Northeast of Brazil. Among 133 patients, 5 had anti-HD, of whom 3 had HDV RNA. HDV genotypes were characterized by Bayesian phylogenetic analysis of nucleotide sequences from the HDAg coding region. HDV-3 was identified in one patient who lives in Maranhao, but was born in Amazonas state (Western Amazon basin). Phylogenetic analysis shows that this HDV-3 sequence grouped with other HDV-3 sequences isolated in this state, which suggests that the patient probably contracted HDV infection there. Surprisingly, the other two patients were infected with HDV-8, an African genotype. These patients were born and have always lived in Urbano Santos, a rural county of Maranhao state, moreover they had never been to Africa and denied any contact with people from that continent. This is the first description of the HDV-8 in non-native African populations. This genotype may have been introduced to Brazil through the slaves brought to the country from the West Africa regions during the 16-18th centuries. Our results indicate that the need of clinical and epidemiological studies to investigate the presence of this infection in other areas in Brazil. (C) 2011 Elsevier B.V. All rights reserved.

Genetic association and cellular function of MC1R variant alleles in human pigmentation

We have examined MC1R variant allele frequencies in the general population of South East Queensland and in a collection of adolescent dizygotic and monozygotic twins and family members to define statistical associations with hair and skin color, freckling, and mole count. Results of these studies are consistent with a linear recessive allelic model with multiplicative penetrance in the inheritance of red hair. Four alleles, D84E, R151C, R160W, and D294H, are strongly associated with red hair and fair skin with multinomial regression analysis showing odds ratios of 63, 118, 50, and 94, respectively. An additional three low-penetrance alleles V60L, V92M, and R163Q have odds ratios 6, 5, and 2 relative to the wild-type allele. To address the cellular effects of MC1R variant alleles in signal transduction, we expressed these receptors in permanently transfected HEK293 cells. Measurement of receptor activity via induction of a cAMP-responsive luciferase reporter gene found that the R151C and R160W receptors were active in the presence of NDP-MSH ligand, but at much reduced levels compared with that seen with the wild-type receptor. The ability to stimulate phosphorylation of the cAMP response element binding protein (CREB) transcription factor was also apparent in all stimulated MC1R variant allele-expressing HEK293 cell extracts as assessed by immunoblotting. In contrast, human melanoma cell lines showed wide variation in the their ability to undergo cAMP-mediated CREB phosphorylation. Culture of human melanocytes of known MC1R genotype may provide the best experimental approach to examine the functional consequences for each MC1R variant allele. With this objective, we have established more than 300 melanocyte cell strains of defined MC1R genotype.

New SMS mutation leads to a striking reduction in spermine synthase protein function and a severe form of Snyder-Robinson X-linked recessive mental retardation syndrome

We report the identification of a novel mutation at a highly conserved residue within the N-terminal region of spermine synthase (SMS) in a second family with Snyder-Robinson X-linked mental retardation syndrome ( OMIM 309583). This missense mutation, p.G56S, greatly reduces SMS activity and leads to severe epilepsy and cognitive impairment. Our findings contribute to a better delineation and expansion of the clinical spectrum of Snyder-Robinson syndrome, support the important role of the N-terminus in the function of the SMS protein, and provide further evidence for the importance of SMS activity in the development of intellectual processing and other aspects of human development.

Screening of autosomal gene deletions in patients with hypogonadotropic hypogonadism using multiplex ligation-dependent probe amplification: detection of a hemizygosis for the fibroblast growth factor receptor 1

P>Objective Congenital hypogonadotropic hypogonadism with anosmia (Kallmann syndrome) or with normal sense of smell is a heterogeneous genetic disorder caused by defects in the synthesis, secretion and action of gonadotrophin-releasing hormone (GnRH). Mutations involving autosomal genes have been identified in approximately 30% of all cases of hypogonadotropic hypogonadism. However, most studies that screened patients with hypogonadotropic hypogonadism for gene mutations did not include gene dosage methodologies. Therefore, it remains to be determined whether patients without detected point mutation carried a heterozygous deletion of one or more exons. Measurements We used the multiplex ligation-dependent probe amplification (MLPA) assay to evaluate the potential contribution of heterozygous deletions of FGFR1, GnRH1, GnRHR, GPR54 and NELF genes in the aetiology of GnRH deficiency. Patients We studied a mutation-negative cohort of 135 patients, 80 with Kallmann syndrome and 55 with normosmic hypogonadotropic hypogonadism. Results One large heterozygous deletion involving all FGFR1 exons was identified in a female patient with sporadic normosmic hypogonadotropic hypogonadism and mild dimorphisms as ogival palate and cavus foot. FGFR1 hemizygosity was confirmed by gene dosage with comparative multiplex and real-time PCRs. Conclusions FGFR1 or other autosomal gene deletion is a possible but very rare event and does not account for a significant number of sporadic or inherited cases of isolated GnRH deficiency.

Linkage to chromosome 2q36.1 in autosomal dominant Dandy-Walker malformation with occipital cephalocele and evidence for genetic heterogeneity

We previously reported a Vietnamese-American family with isolated autosomal dominant occipital cephalocele. Upon further neuroimaging studies, we have recharacterized this condition as autosomal dominant Dandy-Walker with occipital cephalocele (ADDWOC). A similar ADDWOC family from Brazil was also recently described. To determine the genetic etiology of ADDWOC, we performed genome-wide linkage analysis on members of the Vietnamese-American and Brazilian pedigrees. Linkage analysis of the Vietnamese-American family identified the ADDWOC causative locus on chromosome 2q36.1 with a multipoint parametric LOD score of 3.3, while haplotype analysis refined the locus to 1.1 Mb. Sequencing of the five known genes in this locus did not identify any protein-altering mutations. However, a terminal deletion of chromosome 2 in a patient with an isolated case of Dandy-Walker malformation also encompassed the 2q36.1 chromosomal region. The Brazilian pedigree did not show linkage to this 2q36.1 region. Taken together, these results demonstrate a locus for ADDWOC on 2q36.1 and also suggest locus heterogeneity for ADDWOC.

Association of Dominant and Recessive Genes Confers Anthracnose Resistance in Stem and Leaves of Common Bean

Different genes might be involved in Colletotrichum lindemuthianum resistance in leaves and stem of common bean. This work aimed to study the genetic mechanisms of the resistance in the leaf and stem in segregating populations from backcrosses involving resistant cultivar AN 910408 and susceptible cultivar Ruda inoculated with spore suspensions of C. lindemuthianum race 83. Our results indicate that two genes which interact epistatically, one dominant and one recessive, are involved in the genetic control of leaf anthracnose resistance. As for stem anthracnose resistance, two genes also epistatic, one dominant and one recessive, explain the resistance to C. lindemuthianum race 83. The recessive gene is the same for leaf and stem resistance; however, the dominant genes are distinct and independent from each other. The three independent resistance genes of AN 910408 observed in this work could be derived from Guanajuato 31.

Identification of an apolipoprotein(e) variant associated with type III hyperlipoproteinaemia in an indigenous Australian

As a result of testing for lipid and apolipoprotein(e) (apo E) phenotype status of an indigenous Australian community, an apo E variant associated with type III hyperlipoproteinaemia has been identified. Apo E phenotype was determined by analysis of VLDL by isoelectric focusing, and genotype on DNA amplified by polymerase chain reaction, using two different restriction enzyme isotyping assays. Phenotypes and genotypes were discordant in samples from two subjects and an abnormal-sized restriction fragment was also observed in their genotyping gel patterns. DNA sequencing studies revealed this was due to a single nucleotide deletion. 3817delC, at amino acid 136 on apo E. This resulted in a new reading frame and the premature termination of the apo E protein due to a stop codon (TGA) at nucleotide 4105. The variant apo E null allele showed a recessive mode of inheritance and, in combination with the E2 allele, resulted in the type III hyperlipoproteinaemic phenotype but when inherited with the E4 allele had no marked effect on plasma lipids.

CADASIL (cerebral autosomal dominant arteriopathy with subcortical infarcts and leucoencephalopathy): an Australian perspective

Cerebral Autosomal Dominant Arteriopathy with Subcortical Infarcts and Leucoencephalopathy (CADASIL) is a recently described cause of stroke or stroke-like episodes. It is caused by mutations in the Notch3 gene on chromosome 19p. We sought to demonstrate mutations of the Notch3 gene in Australian patients suspected of having CADASIL. Patients from several families were referred to the study. A diagnosis was determined clinically and by neuroimaging. Those suspected of having CADASIL had sequencing of exons 3 and 4 of the Notch3 gene. Eight patients, two of whom were siblings, were suspected of having CADASIL. Five patients (including the siblings) had mutations. Because of strong clustering of Notch3 mutations in CADASIL, this has potential as a reliable test for the disease in Australian patients. (C) 2001 Harcourt Publishers Ltd.

Power of the joint segregation analysis method for testing mixed major-gene and polygene inheritance models of quantitative traits

Understanding the genetic architecture of quantitative traits can greatly assist the design of strategies for their manipulation in plant-breeding programs. For a number of traits, genetic variation can be the result of segregation of a few major genes and many polygenes (minor genes). The joint segregation analysis (JSA) is a maximum-likelihood approach for fitting segregation models through the simultaneous use of phenotypic information from multiple generations. Our objective in this paper was to use computer simulation to quantify the power of the JSA method for testing the mixed-inheritance model for quantitative traits when it was applied to the six basic generations: both parents (P-1 and P-2), F-1, F-2, and both backcross generations (B-1 and B-2) derived from crossing the F-1 to each parent. A total of 1968 genetic model-experiment scenarios were considered in the simulation study to quantify the power of the method. Factors that interacted to influence the power of the JSA method to correctly detect genetic models were: (1) whether there were one or two major genes in combination with polygenes, (2) the heritability of the major genes and polygenes, (3) the level of dispersion of the major genes and polygenes between the two parents, and (4) the number of individuals examined in each generation (population size). The greatest levels of power were observed for the genetic models defined with simple inheritance; e.g., the power was greater than 90% for the one major gene model, regardless of the population size and major-gene heritability. Lower levels of power were observed for the genetic models with complex inheritance (major genes and polygenes), low heritability, small population sizes and a large dispersion of favourable genes among the two parents; e.g., the power was less than 5% for the two major-gene model with a heritability value of 0.3 and population sizes of 100 individuals. The JSA methodology was then applied to a previously studied sorghum data-set to investigate the genetic control of the putative drought resistance-trait osmotic adjustment in three crosses. The previous study concluded that there were two major genes segregating for osmotic adjustment in the three crosses. Application of the JSA method resulted in a change in the proposed genetic model. The presence of the two major genes was confirmed with the addition of an unspecified number of polygenes.

Inheritance and mapping of 11 avirulence genes in Phytophthora sojae

Two new crosses involving four races (races 7, 16, 17, and 25) of the soybean root and stem rot pathogen Phytophthora sojae were established (7/16 cross; 17/25 cross). An F-2 Population derived from each cross was used to determine the genetic basis of avirulence towards 11 different resistance genes in soybean. Avirulence was found to be dominant and determined by a single locus for Avr1b, 1d, 1k, 3b, 4, and 6, as expected for a simple gene-for-gene model. We also observed several cases of segregation, inconsistent with a single dominant gene being solely responsible for avirulence, which suggests that the genetic background of the different crosses can affect avirulence. Avr4 and 6 cosegregated in both the 7/16 and 17/25 crosses and, in the 7/16 cross, Avr1b and 1k were closely linked. Information from segregating RAPD, RFLP, and AFLP markers screened on F-2 progeny from the two new crosses and two crosses described previously (a total of 212 F-2 individuals, 53 from each cross) were used to construct an integrated genetic linkage map of P. sojae. This revised genetic linkage map consists of 386 markers comprising 35 RFLP, 236 RAPD, and 105 AFLP markers, as well as 10 avirulence genes. The map is composed of 21 major linkage groups and seven minor linkage groups covering a total map distance of 1640.4 cM. (C) 2002 Elsevier Science (USA). All rights reserved.