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.     

1,2,3-Triazole Fused Quinoline-Peptidomimetics: Studies on Synthesis, DNA Binding and Photonuclease Activity

Abstract  

We present, simple approach for the accession of 1,2,3-triazole fused quinoline peptide analogues from 3-(azidomethyl)-2-chloroquinoline in a three-step mechanistic pathway. The UV–Visible absorbance plot shows dynamic interaction of parent triazole derivative with CT DNA as efficient DNA intercalator (K _b = 4.6 × 10⁻⁴ M⁻¹). Finally, the efficient DNA damage was observed on photo-irradiation at 360 nm in the presence of 2-(9H-Fluoren-9-ylmethoxycarbonylamino)-propionic acid 1-(2-chloro-quinolin-3-ylmethyl)-1H-[1,2,3]triazole-4-ylmethyl ester (6a).     

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.     

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.     

Crystallographic characterization of the α,γ C12 helix in hybrid peptide sequences

The solid‐state conformations of two αγ hybrid peptides Boc‐[Aib‐γ⁴(R)Ile]₄‐OMe 1 and Boc‐[Aib‐γ⁴(R)Ile]₅‐OMe 2 are described. Peptides 1 and 2 adopt C₁₂‐helical conformations in crystals. The structure of octapeptide 1 is stabilized by six intramolecular 4 → 1 hydrogen bonds, forming 12 atom C₁₂ motifs. The structure of peptide 2 reveals the formation of eight successive C₁₂ hydrogen‐bonded turns. Average backbone dihedral angles for αγ C₁₂ helices are peptide 1 , Aib; φ (°) = ?57.2 ± 0.8, ψ (°) = ?44.5 ± 4.7; γ⁴(R)Ile; φ (°) = ?127.3 ± 7.3, θ₁ (°) = 58.5 ± 12.1, θ₂ (°) = 67.6 ± 10.1, ψ (°) = ?126.2 ± 16.1; peptide 2 , Aib; φ (°) = ?58.8 ± 5.1, ψ (°) = ?40.3 ± 5.5; ψ⁴(R)Ile; φ (°) = ?123.9 ± 2.7, θ₁ (°) = 53.3 θ 4.9, θ ₂ (°) = 61.2 ± 1.6, ψ (°) = ?121.8 ± 5.1. The tendency of γ⁴‐substituted residues to adopt gauche–gauche conformations about the C^α–C^β and C^β–C^γ bonds facilitates helical folding. The αγ C₁₂ helix is a backbone expanded analog of α peptide 3₁₀ helix. The hydrogen bond parameters for α peptide 3₁₀ and α‐helices are compared with those for αγ hybrid C₁₂ helix. Copyright © 2016 European Peptide Society and John Wiley & Sons.     

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.     

A common genetic variant in the neurexin superfamily member CNTNAP2 increases familial risk of autism

Autism is a childhood neuropsychiatric disorder that, despite exhibiting high heritability, has largely eluded efforts to identify specific genetic variants underlying its etiology. We performed a two-stage genetic study in which genome-wide linkage and family-based association mapping was followed up by association and replication studies in an independent sample. We identified a common polymorphism in contactin-associated protein-like 2 (CNTNAP2), a member of the neurexin superfamily, that is significantly associated with autism susceptibility. Importantly, the genetic variant displays a parent-of-origin and gender effect recapitulating the inheritance of autism.     

Development and characterization of a reverse genetic system for studying dengue virus serotype 3 strain variation and neutralization

Dengue viruses (DENV) are enveloped single-stranded positive-sense RNA viruses transmitted by Aedes spp. mosquitoes. There are four genetically distinct serotypes designated DENV-1 through DENV-4, each further subdivided into distinct genotypes. The dengue scientific community has long contended that infection with one serotype confers lifelong protection against subsequent infection with the same serotype, irrespective of virus genotype. However this hypothesis is under increased scrutiny and the role of DENV genotypic variation in protection from repeated infection is less certain. As dengue vaccine trials move increasingly into field-testing, there is an urgent need to develop tools to better define the role of genotypic variation in DENV infection and immunity. To better understand genotypic variation in DENV-3 neutralization and protection, we designed and constructed a panel of isogenic, recombinant DENV-3 infectious clones, each expressing an envelope glycoprotein from a different DENV-3 genotype; Philippines 1982 (genotype I), Thailand 1995 (genotype II), Sri Lanka 1989 and Cuba 2002 (genotype III) and Puerto Rico 1977 (genotype IV). We used the panel to explore how natural envelope variation influences DENV-polyclonal serum interactions. When the recombinant viruses were tested in neutralization assays using immune sera from primary DENV infections, neutralization titers varied by as much as ～19-fold, depending on the expressed envelope glycoprotein. The observed variability in neutralization titers suggests that relatively few residue changes in the E glycoprotein may have significant effects on DENV specific humoral immunity and influence antibody mediated protection or disease enhancement in the setting of both natural infection and vaccination. These genotypic differences are also likely to be important in temporal and spatial microevolution of DENV-3 in the background of heterotypic neutralization. The recombinant and synthetic tools described here are valuable for testing hypotheses on genetic determinants of DENV-3 immunopathogenesis.     

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.