Reliability of genomic predictions of complex human phenotypes

Genome-wide association studies have helped us identify a wealth of genetic variants associated with complex human phenotypes. Because most variants explain a small portion of the total phenotypic variation, however, marker-based studies remain limited in their ability to predict such phenotypes. Here, we show how modern statistical genetic techniques borrowed from animal breeding can be employed to increase the accuracy of genomic prediction of complex phenotypes and the power of genetic mapping studies. Specifically, using the triglyceride data of the GAW20 data set, we apply genomic-best linear unbiased prediction (G-BLUP) methods to obtain empirical genetic values (EGVs) for each triglyceride phenotype and each individual. We then study 2 different factors that influence the prediction accuracy of G-BLUP for the analysis of human data: (a) the choice of kinship matrix, and (b) the overall level of relatedness. The resulting genetic values represent the total genetic component for the phenotype of interest and can be used to represent a trait without its environmental component. Finally, using empirical data, we demonstrate how this method can be used to increase the power of genetic mapping studies. In sum, our results show that dense genome-wide data can be used in a wider scope than previously anticipated

Heritability and genetic associations of triglyceride and HDL-C levels using pedigree-based and empirical kinships

The heritability of a phenotype is an estimation of the percent of variance in that phenotype that is attributable to additive genetic factors. Heritability is optimally estimated in family-based sample populations. Traditionally, this involves use of a pedigree-based kinship coefficient generated from the collected genealogical relationships between family members. An alternative, when dense genotype data are available, is to directly measure the empirical kinship between samples. This study compares the use of pedigree and empirical kinships in the GAW20 data set. Two phenotypes were assessed: triglyceride levels and high-density lipoprotein cholesterol (HDL-C) levels pre- and postintervention with the cholesterol-reducing drug fenofibrate. Using SOLAR (Sequential Oligogenic Linkage Analysis Routines), pedigree-based kinships and empirically calculated kinships (using IBDLD and LDAK) were used to calculate phenotype heritability. In addition, a genome-wide association study was conducted using each kinship model for each phenotype to identify genetic variants significantly associated with phenotypic variation. The variant rs247617 was significantly associated with HDL-C levels both pre- and post-fenofibrate intervention. Overall, the phenotype heritabilities calculated using pedigree based kinships or either of the empirical kinships generated using IBDLD or LDAK were comparable. Phenotype heritabilities estimated from empirical kinships generated using IBDLD were closest to the pedigree-based estimations. Given that there was not an appreciable amount of unknown relatedness between the pedigrees in this data set, a large increase in heritability in using empirical kinship was not expected, and our calculations support this. Importantly, these results demonstrate that when sufficient genotypic data are available, empirical kinship estimation is a practical alternative to using pedigree-based kinships

Genome-wide linkage scan for loci influencing plasma triglyceride levels

We conducted a genome-wide linkage scan to detect loci that influence the levels of fasting triglycerides in plasma. Fasting triglyceride levels were available at 4 time points (visits), 2 pre- and 2 post-fenofibrate intervention. Multipoint identity-by-descent (MIBD) matrices were derived from genotypes using IBDLD. Variance-component linkage analyses were then conducted using SOLAR (Sequential Oligogenic Linkage Analysis Routines). We found evidence of linkage (logarithm of odds [LOD] ≥3) at 5 chromosomal regions with triglyceride levels in plasma. The highest LOD scores were observed for linkage to the estimated genetic value (additive genetic component) of the log-normalized triglyceride levels in plasma. Our results suggest that a chromosome 10 locus at 37 cM (LODpre = 3.01, LODpost = 3.72) influences fasting triglyceride levels in plasma regardless of the fenofibrate intervention, and that loci in chromosomes 1 at 170 cM and 4 at 24 cM ceases to affect the triglyceride levels when fenofibrate is present, while the regions in chromosomes 6 at 136 to 162 cM and 11 at 39 to 40 cM appear to influence triglyceride levels in response to fenofibrate

Comparative Analysis of Immune Checkpoint Molecules and Their Potential Role in the Transmissible Tasmanian Devil Facial Tumor Disease

Immune checkpoint molecules function as a system of checks and balances that enhance or inhibit immune responses to infectious agents, foreign tissues, and cancerous cells. Immunotherapies that target immune checkpoint molecules, particularly the inhibitory molecules programmed cell death 1 and cytotoxic T-lymphocyte-associated protein 4 (CTLA-4), have revolutionized human oncology in recent years, yet little is known about these key immune signaling molecules in species other than primates and rodents. The Tasmanian devil facial tumor disease is caused by transmissible cancers that have resulted in a massive decline in the wild Tasmanian devil population. We have recently demonstrated that the inhibitory checkpoint molecule PD-L1 is upregulated on Tasmanian devil (Sarcophilus harrisii) facial tumor cells in response to the interferon-gamma cytokine. As this could play a role in immune evasion by tumor cells, we performed a thorough comparative analysis of checkpoint molecule protein sequences among Tasmanian devils and eight other species. We report that many of the key signaling motifs and ligand-binding sites in the checkpoint molecules are highly conserved across the estimated 162 million years of evolution since the last common ancestor of placental and non-placental mammals. Specifically, we discovered that the CTLA-4 (MYPPPY) ligand-binding motif and the CTLA-4 (GVYVKM) inhibitory domain are completely conserved across all nine species used in our comparative analysis, suggesting that the function of CTLA-4 is likely conserved in these species. We also found that cysteine residues for intra- and intermolecular disulfide bonds were also highly conserved. For instance, all 20 cysteine residues involved in disulfide bonds in the human 4-1BB molecule were also present in devil 4-1BB. Although many key sequences were conserved, we have also identified immunoreceptor tyrosine-based inhibitory motifs (ITIMs) and immunoreceptor tyrosine-based switch motifs (ITSMs) in genes and protein domains that have not been previously reported in any species. This checkpoint molecule analysis and review of salient features for each of the molecules presented here can serve as road map for the development of a Tasmanian devil facial tumor disease immunotherapy. Finally, the strategies can be used as a guide for veterinarians, ecologists, and other researchers willing to venture into the nascent field of wild immunology

Dose response of the 16p11.2 distal copy number variant on intracranial volume and basal ganglia

Carriers of large recurrent copy number variants (CNVs) have a higher risk of developing neurodevelopmental disorders. The 16p11.2 distal CNV predisposes carriers to e.g., autism spectrum disorder and schizophrenia. We compared subcortical brain volumes of 12 16p11.2 distal deletion and 12 duplication carriers to 6882 non-carriers from the large-scale brain Magnetic Resonance Imaging collaboration, ENIGMA-CNV. After stringent CNV calling procedures, and standardized FreeSurfer image analysis, we found negative dose-response associations with copy number on intracranial volume and on regional caudate, pallidum and putamen volumes (β = −0.71 to −1.37; P \u3c 0.0005). In an independent sample, consistent results were obtained, with significant effects in the pallidum (β = −0.95, P = 0.0042). The two data sets combined showed significant negative dose-response for the accumbens, caudate, pallidum, putamen and ICV (P=0.0032, 8.9 × 10−6, 1.7 × 10−9, 3.5 × 10−12 and 1.0 × 10−4, respectively). Full scale IQ was lower in both deletion and duplication carriers compared to noncarriers. This is the first brain MRI study of the impact of the 16p11.2 distal CNV, and we demonstrate a specific effect on subcortical brain structures, suggesting a neuropathological pattern underlying the neurodevelopmental syndromes

Molecular characterization of a marine turtle tumor epizootic, profiling external, internal and postsurgical regrowth tumors

Sea turtle populations are under threat from an epizootic tumor disease (animal epidemic) known as fibropapillomatosis. Fibropapillomatosis continues to spread geographically, with prevalence of the disease also growing at many longer-affected sites globally. However, we do not yet understand the precise environmental, mutational and viral events driving fibropapillomatosis tumor formation and progression. Here we perform transcriptomic and immunohistochemical profiling of five fibropapillomatosis tumor types: external new, established and postsurgical regrowth tumors, and internal lung and kidney tumors. We reveal that internal tumors are molecularly distinct from the more common external tumors. However, they have a small number of conserved potentially therapeutically targetable molecular vulnerabilities in common, such as the MAPK, Wnt, TGFβ and TNF oncogenic signaling pathways. These conserved oncogenic drivers recapitulate remarkably well the core pan-cancer drivers responsible for human cancers. Fibropapillomatosis has been considered benign, but metastatic-related transcriptional signatures are strongly activated in kidney and established external tumors. Tumors in turtles with poor outcomes (died/euthanized) have genes associated with apoptosis and immune function suppressed, with these genes providing putative predictive biomarkers. Together, these results offer an improved understanding of fibropapillomatosis tumorigenesis and provide insights into the origins, inter-tumor relationships, and therapeutic treatment for this wildlife epizootic

Direct observation of phase coherence in 3-<b>k</b> magnetic configurations

International audienceWe report the observation by neutron diffraction of phase coherent Bragg reflections in a multi-k magnetic configuration with a spatial periodicity outside the conventional scattering cross-section. The peaks, which exist in the 3-k state of UAs0.8Se0.2, display long-range order with a wavevector dependence characteristic of a magnetic interaction. The results confirm the long-range order and temperature dependence reported in an earlier study of similar peaks in this material using x-ray resonant scattering by (a) the non-trivial extension to the technique of neutron diffraction, and (b) the observation of similar 3-k phase-coherent reflections in other samples by x-ray resonant scattering. The importance of the neutron diffraction results lies primarily in the fact that magnetic neutron diffraction is well established as a weak probe operating on thermodynamic time scales. This alleviates concern that the rapid (10-15 - 10-14 s), strong interaction, characteristic of the resonant x-ray technique, is imaging a transient or non-equilibrium configuration. Likewise, the extension of the x-ray resonant scattering results to other samples establishes the generality of this effect. The enigma of how to understand the observed diffraction, which appears to lie strictly outside both the conventional neutron and x-ray scattering cross sections, remains

Rare, Potentially Pathogenic Variants in ZNF469 Are Not Enriched in Keratoconus in a Large Australian Cohort of European Descent

Purpose: The Zinc Finger Protein 469 (ZNF469) gene has been proposed as a candidate gene for keratoconus due to the association of an upstream polymorphism (rs9938149) with the disease in two independent studies, and the role of the gene in the autosomal recessive disease Brittle Cornea Syndrome. Coding variants in ZNF469 have been assessed for association with keratoconus in several small studies, with conflicting results. We assessed rare, potentially pathogenic variants in ZNF469 for enrichment in keratoconus patients in a cohort larger than all previous studies combined. Methods: ZNF469 was sequenced in 385 Australian keratoconus patients of European descent, 346 population controls, and 230 ethnically matched screened controls by either whole exome sequencing or targeted gene sequencing. The frequency of rare and very rare potentially pathogenic variants was compared between cases and controls using χ2 or Fisher\u27s exact tests and further explored using a gene based test (Sequence Kernel Association Test [SKAT]), weighting on the rarity of variants. Results: A total of 49 rare, including 33 very rare, potentially pathogenic variants were identified across all groups. No enrichment of rare or very rare potentially pathogenic variants in ZNF469 was observed in our cases compared to the control groups following analysis using χ2 or Fisher\u27s exact tests. This finding was further supported by the SKAT results, which found no significant difference in the frequency of variants predicted to be damaging between cases and either control group (P = 0.06). Conclusions: Rare variants in ZNF469 do not contribute to keratoconus susceptibility and do not account for the association at rs9938149

Genetic analysis localizes a novel locus on chromosome 4q for the glaucoma endophenotype, cup-to-disc ratio: The Jiri Eye Study

Purpose: Glaucoma is a heterogeneous disease influenced by genetic risk factors. However, not all genetic risk factors have been identified. The aim of this project is to localize genetic factors influencing known glaucoma endophenotypes: intraocular pressure (IOP), central corneal thickness (CCT), and vertical cup-to-disc ratio (VCDR). Methods: This family-based study design utilizes phenotypic and genomic data from a single well-characterized pedigree residing in the Jiri region of Nepal. Measures of IOP, CCT and VCDR were obtained by Goldmann applanation tonometry, OCT, and slit lamp biomicroscopy, respectively. Using a genome-wide genotype data set (~550,000 SNPs), we performed a genome-wide linkage scan for IOP, CCT, and VCDR adopting a quantitative approach in SOLAR. For localized quantitative trait locus (QTL) signals, we screened all SNPs within the 1-LOD (95% confidence) interval using the classical measured genotype approach to association analysis and allowing for non-independence amongst the pedigree members. Results: For this study, phenotypic and genotype data from 1,163 (55% female) members of the Jirel population were available. The mean age of the sample is 43.8 (SD=15.7) years. IOP (h2=19%, p=6.1×10-5), CCT (h2=57%, p=1.6×10-26), and VCDR (h2=48%, p=9.7×10-22) were significantly heritable. We localized a significant QTL for VCDR on chromosome 4 (LOD=3.05 at 86.83 Mb). The top association signal within this QTL was for an intronic SNP (rs4148155; p=2.01×10-6, b=0.24) in the ABCG2 (ATP binding cassette subfamily G member 2) gene, which satisfied our QTL-specific Bonferroni-corrected significance criterion (p\u3c6.59×10-5). ABCG2 is a known stem cell marker, which is positively expressed in clonal human trabecular meshwork stem cells. Another positional candidate gene of note is SCD5 (Stearoyl-CoA desaturase 5), which is shown to suppress neurite outgrowth, a marker of neuronal differentiation. SCD5 is of significant interest given that expression of myocilin (MYOC) also inhibits neurite outgrowth. Conclusions: To our knowledge, the VCDR QTL on chromosome 4 is a novel locus and does not overlap with other glaucoma endophenotypes or glaucoma disease status. These results highlight the importance of continued evaluation of genetic factors influencing glaucoma endophenotypes in under-studied populations, such as the Jirels, as new information may be elucidated

1q21.1 distal copy number variants are associated with cerebral and cognitive alterations in humans

Low-frequency 1q21.1 distal deletion and duplication copy number variant (CNV) carriers are predisposed to multiple neurodevelopmental disorders, including schizophrenia, autism and intellectual disability. Human carriers display a high prevalence of micro- and macrocephaly in deletion and duplication carriers, respectively. The underlying brain structural diversity remains largely unknown. We systematically called CNVs in 38 cohorts from the large-scale ENIGMACNV collaboration and the UK Biobank and identified 28 1q21.1 distal deletion and 22 duplication carriers and 37,088 non-carriers (48% male) derived from 15 distinct magnetic resonance imaging scanner sites. With standardized methods, we compared subcortical and cortical brain measures (all) and cognitive performance (UK Biobank only) between carrier groups also testing for mediation of brain structure on cognition. We identified positive dosage effects of copy number on intracranial volume (ICV) and total cortical surface area, with the largest effects in frontal and cingulate cortices, and negative dosage effects on caudate and hippocampal volumes. The carriers displayed distinct cognitive deficit profiles in cognitive tasks from the UK Biobank with intermediate decreases in duplication carriers and somewhat larger in deletion carriers—the latter potentially mediated by ICV or cortical surface area. These results shed light on pathobiological mechanisms of neurodevelopmental disorders, by demonstrating gene dose effect on specific brain structures and effect on cognitive function