Information content of the Centre d'Etude du Polymorphisme Humain (CEPH) family structures for linkage studies

Summary This paper derives theoretical values for joint polymorphism information content for two markers from a family structure consisting of four grandparents, two parents, and many offspring. These data determine the efficiency of linkage map construction.     

Mean and variance of FST in a finite number of incompletely isolated populations

In the presence of migration F_ST in a finite number of incompletely isolated populations first increases, but after reaching a certain maximum value, it starts to decline and eventually becomes 0. The mean and variance of F_ST in this process are studied by using the recurrence formulas for the moments of gene frequencies in the island model of finite size as well as by using Monte Carlo simulation. The mean and variance in the early generations can be predicted by the approximate formulas developed. On the other hand, if we exclude the cases of an allele being fixed in all subpopulations, the mean of F_ST eventually reaches a steady-state value. This value is given by 1 − 2N_T(1 − λ) approximately, where N_T is the total population size and λ is the rate of decay of heterozygosity at steady state. It is shown that the mean and variance of F_ST depend on the initial gene frequency and when this is close to 0 or 1, Lewontin and Krakauer's test of the neutrality of polymorphic genes is not valid.     

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.     

Segregation analysis of blood pressure and body mass index in a rural US community

To assess evidence for a gene with large effect on systolic blood pressure (SBP), diastolic blood pressure (DBP), and body mass index (BMI), we conducted segregation analyses on 261 nuclear families collected from a rural Caucasian community in Michigan. The families were ascertained through a hypertensive proband. Each phenotype was adjusted for significant covariate effects (e.g., gender and age). We used class D regressive models to conduct the segregation analyses. Our analysis results support the segregation of a major gene for BMI, but not for SBP or DBP. A recessive locus effect provided the best explanation for BMI where approximately 43% of the variance of BMI was due to this gene.     

An evaluation of the assembly of an approximately 15-Mb region on human chromosome 13q32-q33 linked to bipolar disorder and schizophrenia

The human 13q32-q33 region has been linked to both bipolar disorder and schizophrenia. Before completion of the draft sequences, we developed an approximately 15-Mb comprehensive map for the region extending from D13S1300 to ATA35H12. This map was assembled using publicly available mapping data and sequence-tagged site (STS)-based PCR confirmation. We then compared this map with the NCBI, Celera Genomics, and UCSC Golden Path data in February, June, and September 2001. All data sets showed gaps, misassignment of STSs, and errors in orientation and marker order. Surprisingly, the completed sequences of many bacterial artificial chromosomes (BACs) had been truncated. Of 21 gaps that were detected, 4 were present in both the NCBI and Celera databases. All gaps could be filled using 1-2 BAC clones. A total of 39 loci mapped to additional sites within the human genome, providing evidence of segmental duplications. Additionally, 61 unique cDNA clones were sequenced to increase available transcribed sequence, and 11,353 reference single-nucleotide polymorphisms (SNPs) with an average density of 1 SNP/3720 bases were identified. Overall, integration of the data from multiple sources is still needed for complete assembly of the 13q32-q33 region. (c)     

2013 William Allan Award: My Multifactorial Journey

Gonadotropin-releasing hormone (GnRH) neurons originate outside the CNS in the olfactory placode and migrate into the CNS, where they become integral components of the hypothalamic-pituitary-gonadal (HPG) axis. Disruption of this migration results in Kallmann syndrome (KS), which is characterized by anosmia and pubertal failure due to hypogonadotropic hypogonadism. Using candidate-gene screening, autozygosity mapping, and whole-exome sequencing in a cohort of 30 individuals with KS, we searched for genes newly associated with KS. We identified homozygous loss-of-function mutations in FEZF1 in two independent consanguineous families each with two affected siblings. The FEZF1 product is known to enable axons of olfactory receptor neurons (ORNs) to penetrate the CNS basal lamina in mice. Because a subset of axons in these tracks is the migratory pathway for GnRH neurons, in FEZF1 deficiency, GnRH neurons also fail to enter the brain. These results indicate that FEZF1 is required for establishment of the central component of the HPG axis in humans.     

Meta-analysis of Correlated Traits via Summary Statistics from GWASs with an Application in Hypertension

Genome-wide association studies (GWASs) have identified many genetic variants underlying complex traits. Many detected genetic loci harbor variants that associate with multiple—even distinct—traits. Most current analysis approaches focus on single traits, even though the final results from multiple traits are evaluated together. Such approaches miss the opportunity to systemically integrate the phenome-wide data available for genetic association analysis. In this study, we propose a general approach that can integrate association evidence from summary statistics of multiple traits, either correlated, independent, continuous, or binary traits, which might come from the same or different studies. We allow for trait heterogeneity effects. Population structure and cryptic relatedness can also be controlled. Our simulations suggest that the proposed method has improved statistical power over single-trait analysis in most of the cases we studied. We applied our method to the Continental Origins and Genetic Epidemiology Network (COGENT) African ancestry samples for three blood pressure traits and identified four loci (CHIC2, HOXA-EVX1, IGFBP1/IGFBP3, and CDH17; p < 5.0 × 10⁻⁸) associated with hypertension-related traits that were missed by a single-trait analysis in the original report. Six additional loci with suggestive association evidence (p < 5.0 × 10⁻⁷) were also observed, including CACNA1D and WNT3. Our study strongly suggests that analyzing multiple phenotypes can improve statistical power and that such analysis can be executed with the summary statistics from GWASs. Our method also provides a way to study a cross phenotype (CP) association by using summary statistics from GWASs of multiple phenotypes.     

Functional Loss of Semaphorin 3C and/or Semaphorin 3D and Their Epistatic Interaction with Ret Are Critical to Hirschsprung Disease Liability

Innervation of the gut is segmentally lost in Hirschsprung disease (HSCR), a consequence of cell-autonomous and non-autonomous defects in enteric neuronal cell differentiation, proliferation, migration, or survival. Rare, high-penetrance coding variants and common, low-penetrance non-coding variants in 13 genes are known to underlie HSCR risk, with the most frequent variants in the ret proto-oncogene (RET). We used a genome-wide association (220 trios) and replication (429 trios) study to reveal a second non-coding variant distal to RET and a non-coding allele on chromosome 7 within the class 3 Semaphorin gene cluster. Analysis in Ret wild-type and Ret-null mice demonstrates specific expression of Sema3a, Sema3c, and Sema3d in the enteric nervous system (ENS). In zebrafish embryos, sema3 knockdowns show reduction of migratory ENS precursors with complete ablation under conjoint ret loss of function. Seven candidate receptors of Sema3 proteins are also expressed within the mouse ENS and their expression is also lost in the ENS of Ret-null embryos. Sequencing of SEMA3A, SEMA3C, and SEMA3D in 254 HSCR-affected subjects followed by in silico protein structure modeling and functional analyses identified five disease-associated alleles with loss-of-function defects in semaphorin dimerization and binding to their cognate neuropilin and plexin receptors. Thus, semaphorin 3C/3D signaling is an evolutionarily conserved regulator of ENS development whose dys-regulation is a cause of enteric aganglionosis.     

Quantifying and modeling birth order effects in autism

Autism is a complex genetic disorder with multiple etiologies whose molecular genetic basis is not fully understood. Although a number of rare mutations and dosage abnormalities are specific to autism, these explain no more than 10% of all cases. The high heritability of autism and low recurrence risk suggests multifactorial inheritance from numerous loci but other factors also intervene to modulate risk. In this study, we examine the effect of birth rank on disease risk which is not expected for purely hereditary genetic models. We analyzed the data from three publicly available autism family collections in the USA for potential birth order effects and studied the statistical properties of three tests to show that adequate power to detect these effects exist. We detect statistically significant, yet varying, patterns of birth order effects across these collections. In multiplex families, we identify V-shaped effects where middle births are at high risk; in simplex families, we demonstrate linear effects where risk increases with each additional birth. Moreover, the birth order effect is gender-dependent in the simplex collection. It is currently unknown whether these patterns arise from ascertainment biases or biological factors. Nevertheless, further investigation of parental age-dependent risks yields patterns similar to those observed and could potentially explain part of the increased risk. A search for genes considering these patterns is likely to increase statistical power and uncover novel molecular etiologies.     

Structure and assembly of designed beta-hairpin peptides in crystals as models for beta-sheet aggregation

De novo designed beta-hairpin peptides have generally been recalcitrant to crystallization. The crystal structures of four synthetic peptide beta-hairpins, Boc-Leu-Val-Val-DPro-Gly-Leu-Phe-Val-OMe (1), Boc-Leu-Phe-Val-DPro-Ala-Leu-Phe-Val-OMe (2), Boc-Leu-Val-Val-DPro-Aib-Leu-Val-Val-OMe (3), and Boc-Met-Leu-Phe-Val-DPro-Ala-Leu-Val-Val-Phe-OMe (4), are described. The centrally positioned DPro-Xxx segment promotes prime beta-turn formation, thereby nucleating beta-hairpin structures. In all four peptides well-defined beta-hairpins nucleated by central type II' DPro-Xxx beta-turns have been characterized by X-ray diffraction, providing a view of eight crystallographically independent hairpins. In peptides 1-3 three intramolecular cross-strand hydrogen bonds stabilized the observed beta-hairpin, with some fraying of the structures at the termini. In peptide 4, four intramolecular cross-strand hydrogen bonds stabilized the hairpin. Peptides 1-4 reveal common features of packing of beta-hairpins into crystals. Two-dimensional sheet formation mediated by intermolecular hydrogen bonds formed between antiparallel strands of adjacent molecule is a recurrent theme. The packing of two-dimensional sheets into the crystals is mediated in the third dimension by bridging solvents and interactions of projecting side chains, which are oriented on either face of the sheet. In all cases, solvation of the central DPro-Xxx peptide unit beta-turn is observed. The hairpins formed in the octapeptides are significantly buckled as compared to the larger hairpin in peptide 4, which is much flatter. The crystal structures provide insights into the possible modes of beta-sheet packing in regular crystalline arrays, which may provide a starting point for understanding beta-sandwich and cross-beta-structures in amyloid fibrils.