Utilizing extended pedigree information for discovery and confirmation of copy number variable regions among Mexican Americans

Characterization of Large Copy Number Variation in Mexican Type 2 Diabetes subjects

The effect of Copy Number Variants (CNVs) on Type 2 Diabetes (T2D) remains little explored. The present study characterized large rare CNVs in 686 T2D and 194 non-T2D subjects of Mexican ancestry genotyped using the Affymetrix Genome-Wide Human SNP array 5.0. Rare CNVs with ≥ 100 kb length were identified using a stringent strategy based on merging CNVs calls generated using Birdsuit, iPattern and PennCNV algorithms. We applied three different strategies to evaluate the distribution of CNVs in the T2D and non-T2D samples: 1) Burden analysis, 2) Identification of CNVs in loci previously associated to T2D, and 3) Identification of CNVs observed only in the T2D group. In the CNV burden analysis, the T2D group showed a higher proportion of CNVs, and also a higher proportion of CNVs overlapping at least one gene than the non T2D group. Five of the six loci previously associated with T2D had duplications or deletions in the T2D sample, but not the non-T2D sample. A gene-set analysis including genes with CNVs observed only in the T2D group highlighted gene-sets related with sensory perception (olfactory receptors, OR) and phenylpyruvate tautomerase/dopachrome isomerase activity (MIF and DDT genes).

The Arg59Trp variant in ANGPTL8 (betatrophin) is associated with total and HDL-cholesterol in American Indians and Mexican Americans and differentially affects cleavage of ANGPTL3

We previously identified a locus linked to total cholesterol (TC) concentration in Pima Indians on chromosome 19p. To characterize this locus, we genotyped >2000 SNPs in 1838 Pimas and assessed association with log(TC). We observed evidence for association with log(TC) with rs2278426 (3.5% decrease/copy of the T allele; P=5.045×10(-6)) in the ANGPTL8 (angiopoietin-like 8) gene. We replicated this association in 2413 participants of the San Antonio Mexican American Family Study (SAMAFS: 2.0% decrease per copy of the T allele; P=0.005842). In a meta-analysis of the combined data, we found the strongest estimated effect with rs2278426 (P=2.563×10(-7)). The variant T allele at rs2278426 predicts an Arg59Trp substitution and has previously been associated with LDL-C and HDL-C. In Pimas and SAMAFS participants, the T allele of rs2278426 was associated with reduced HDL-C levels (P=0.000741 and 0.00002, respectively), and the combined estimated effect for the two cohorts was -3.8% (P=8.526×10(-8)). ANGPTL8 transcript and protein levels increased in response to both glucose and insulin. The variant allele was associated with increased levels of cleaved ANGPTL3. We conclude that individuals with the variant allele may have lower TC and HDL-C levels due to increased activation of ANGPTL3 by ANGPTL8.

Effects of copy number variable regions on local gene expression in white blood cells of Mexican Americans

Only few systematic studies on the contribution of copy number variation to gene expression variation have been published to date. Here we identify effects of copy number variable regions (CNVRs) on nearby gene expression by investigating 909 CNVRs and expression levels of 12059 nearby genes in white blood cells from Mexican-American participants of the San Antonio Family Heart Study. We empirically evaluate our ability to detect the contribution of CNVs to proximal gene expression (presumably in cis) at various window sizes (up to a 10 Mb distance) between the gene and CNV. We found a ~1-Mb window size to be optimal for capturing cis effects of CNVs. Up to 10% of the CNVs in this study were found to be significantly associated with the expression of at least one gene within their vicinity. As expected, we find that CNVs that directly overlap gene sequences have the largest effects on gene expression (compared with non-overlapping CNVRs located nearby), with positive correlation (except for a few exceptions) between estimated genomic dosage and expression level. We find that genes whose expression level is significantly influenced by nearby CNVRs are enriched for immunity and autoimmunity related genes. These findings add to the currently limited catalog of CNVRs that are recognized as expression quantitative trait loci, and have implications for future study designs as well as for prioritizing candidate causal variants in genomic regions associated with disease.

Imputation in families using a heuristic phasing approach

Whole genome sequencing (WGS) remains prohibitively expensive, which has encouraged the development of methods to impute WGS data into nonsequenced individuals using a framework of single nucleotide polymorphisms genotyped for genome-wide association studies (GWAS). Although successful methods have been developed for cohorts of unrelated individuals, current imputation methods in related individuals are limited by pedigree size, by the distance of relationships, or by computation time. In this article, we describe a method for imputation in arbitrarily shaped multigenerational pedigrees that can impute genotypes across distantly related individuals based on identity by descent. We evaluate this approach using GWAS data and apply this approach to WGS data distributed for Genetic Analysis Workshop 18.