Search

Patients with inherited retinal dystrophies (IRDs) were recruited from two understudied populations: Mexico and Pakistan as well as a third well-studied population of European Americans to define the genetic architecture of IRD by performing whole-genome sequencing (WGS). Whole-genome analysis was performed on 409 individuals from 108 unrelated pedigrees with IRDs. All patients underwent an ophthalmic evaluation to establish the retinal phenotype. Although the 108 pedigrees in this study had previously been examined for mutations in known IRD genes using a wide range of methodologies including targeted gene(s) or mutation(s) screening, linkage analysis and exome sequencing, the gene mutations responsible for IRD in these 108 pedigrees were not determined. WGS was performed on these pedigrees using Illumina X10 at a minimum of 30X depth. The sequence reads were mapped against hg19 followed by variant calling using GATK. The genome variants were annotated using SnpEff, PolyPhen2, and CADD score; the structural variants (SVs) were called using GenomeSTRiP and LUMPY. We identified potential causative sequence alterations in 61 pedigrees (57%), including 39 novel and 54 reported variants in IRD genes. For 57 of these pedigrees the observed genotype was consistent with the initial clinical diagnosis, the remaining 4 had the clinical diagnosis reclassified based on our findings. In seven pedigrees (12%) we observed atypical causal variants, i.e. unexpected genotype(s), including 4 pedigrees with causal variants in more than one IRD gene within all affected family members, one pedigree with intrafamilial genetic heterogeneity (different affected family members carrying causal variants in different IRD genes), one pedigree carrying a dominant causative variant present in pseudo-recessive form due to consanguinity and one pedigree with a de-novo variant in the affected family member. Combined atypical and large structural variants contributed to about 20% of cases. Among the novel mutations, 75% were detected in Mexican and 50% found in European American pedigrees and have not been reported in any other population while only 20% were detected in Pakistani pedigrees and were not previously reported. The remaining novel IRD causative variants were listed in gnomAD but were found to be very rare and population specific. Mutations in known IRD associated genes contributed to pathology in 63% Mexican, 60% Pakistani and 45% European American pedigrees analyzed. Overall, contribution of known IRD gene variants to disease pathology in these three populations was similar to that observed in other populations worldwide. This study revealed a spectrum of mutations contributing to IRD in three populations, identified a large proportion of novel potentially causative variants that are specific to the corresponding population or not reported in gnomAD and shed light on the genetic architecture of IRD in these diverse global populations.

To investigate the genetic diversity, spatiotemporal dynamics, and transmission networks of HIV-1 CRF55_01B epidemic in China. A total of 209 partial pol gene sequences of HIV-1 CRF55_01B were sampled during 2007-2015 from 7 provinces of China. Phylogenetic analyses and trait diffusion process of these sequences were performed using Bayesian methods. Distance-based molecular network analyses were performed to infer putative relationships. Characteristics of genetically linked individuals were analyzed. Our study identified that HIV-1 CRF55_01B likely originated among men who have sex with men (MSM) in Guangdong province in January 2003 (April 2000-April 2005), and that Guangdong province and MSM are major hubs for the spread of the HIV-1 CRF55_01B epidemic in China. A Bayesian Skygrid plot revealed that the effective population size of HIV-1 CRF55_01B experienced increased phase followed by a plateau. All sequences from persons of unknown risk clustered within groups who reported MSM risk. This could be because Chinese MSM may not report such risk due to HIV/AIDS-related stigmatization and discrimination. This study inferred the transmission dynamics of the HIV-1 CRF55_01B epidemic in China at high resolution. The methods developed in this study may be critical for designing effective HIV prevention strategies in China and beyond.

PurposePrimary congenital glaucoma (PCG) is a genetically heterogeneous disorder caused by developmental defects in the anterior chamber and trabecular meshwork. This disease is an important cause of childhood blindness. In this study, we aim to identify the genetic determinants of PCG in three consanguineous families of Pakistani descent.MethodsAffected members of all three families underwent detailed ophthalmological examination including slit-lamp biomicroscopy. Blood samples were collected from affected and healthy members of all three families, and genomic DNA was extracted. Linkage analysis was performed for the known or reported loci of PCG to localize the disease interval, and logarithm of odds (LOD) scores were calculated. All protein-coding exons of the candidate gene, latent transforming growth factor-beta binding protein 2 (LTBP2), were bidirectionally sequenced to identify the disease-causing mutation.ResultsShort tandem repeat (STR) marker-based linkage analysis localized the critical interval to chromosome 14q with a maximum two-point LOD score of 2.86 (PKGL076), 2.8 (PKGL015), and 2.92 (PKGL042). Bidirectional Sanger sequencing of LTBP2 revealed three novel pathogenic variants, i.e., c.3028G>A (p.Asp1010Asn), c.3427delC (p.Gln1143Argfs*35), and c.5270G>A (p.Cys1757Tyr) in PKGL076, PKGL015, and PKGL042, respectively. All three mutations segregated with the disease phenotype in their respective families and were absent in 200 ethnically matched normal chromosomes.ConclusionsWe identified three novel mutations, p.D1010N, p.Q1143Rfs*35, and p.C1757Y, in LTBP2 responsible for PCG.

Keratoconus (KC) is the most common corneal ectatic disorder affecting >300,000 people in the US. KC normally has its onset in adolescence, progressively worsening through the third to fourth decades of life. KC patients report significant impaired vision-related quality of life. Genetic factors play an important role in KC pathogenesis. To identify novel genes in familial KC patients, we performed whole exome and genome sequencing in a four-generation family. We identified potential variants in the PPIP5K2 and PCSK1 genes. Using in vitro cellular model and in vivo gene-trap mouse model, we found critical evidence to support the role of PPIP5K2 in normal corneal function and KC pathogenesis. The gene-trap mouse showed irregular corneal surfaces and pathological corneal thinning resembling KC. For the first time, we have integrated corneal tomography and pachymetry mapping into characterization of mouse corneal phenotypes which could be widely implemented in basic and translational research for KC diagnosis and therapy in the future.

BackgroundThis study addresses the strength of associations between trichotillomania (TTM) and other DSM-IV Axis I conditions in a large sample (n = 2606) enriched for familial obsessive-compulsive disorder (OCD), to inform TTM classification.MethodsWe identified participants with TTM in the Johns Hopkins OCD Family Study (153 families) and the OCD Collaborative Genetics Study, a six-site genetic linkage study of OCD (487 families). We used logistic regression (with generalized estimating equations) to assess the strength of associations between TTM and other DSM-IV disorders.ResultsTTM had excess comorbidity with a number of conditions from different DSM-IV chapters, including tic disorders, alcohol dependence, mood disorders, anxiety disorders, impulse-control disorders, and bulimia nervosa. However, association strengths (odds ratios) were highest for kleptomania (6.6), pyromania (5.8), OCD (5.6), skin picking disorder (4.4), bulimia nervosa (3.5), and pathological nail biting (3.4).ConclusionsTTM is comorbid with a number of psychiatric conditions besides OCD, and it is strongly associated with other conditions involving impaired impulse control. Though DSM-5 includes TTM as an OCD-related disorder, its comorbidity pattern also emphasizes the impulsive, appetitive aspects of this condition that may be relevant to classification.

ObjectiveAutosomal dominant optic atrophy (ADOA) starts in early childhood with loss of visual acuity and color vision deficits. OPA1 mutations are responsible for the majority of cases, but in a portion of patients with a clinical diagnosis of ADOA, the cause remains unknown. This study aimed to identify novel ADOA-associated genes and explore their causality.MethodsLinkage analysis and sequencing were performed in multigeneration families and unrelated patients to identify disease-causing variants. Functional consequences were investigated in silico and confirmed experimentally using the zebrafish model.ResultsWe defined a new ADOA locus on 7q33-q35 and identified 3 different missense variants in SSBP1 (NM_001256510.1; c.113G>A [p.(Arg38Gln)], c.320G>A [p.(Arg107Gln)] and c.422G>A [p.(Ser141Asn)]) in affected individuals from 2 families and 2 singletons with ADOA and variable retinal degeneration. The mutated arginine residues are part of a basic patch that is essential for single-strand DNA binding. The loss of a positive charge at these positions is very likely to lower the affinity of SSBP1 for single-strand DNA. Antisense-mediated knockdown of endogenous ssbp1 messenger RNA (mRNA) in zebrafish resulted in compromised differentiation of retinal ganglion cells. A similar effect was achieved when mutated mRNAs were administered. These findings point toward an essential role of ssbp1 in retinal development and the dominant-negative nature of the identified human variants, which is consistent with the segregation pattern observed in 2 multigeneration families studied.InterpretationSSBP1 is an essential protein for mitochondrial DNA replication and maintenance. Our data have established pathogenic variants in SSBP1 as a cause of ADOA and variable retinal degeneration. ANN NEUROL 2019;86:368-383.

Despite intensive study, most of the specific genetic factors that contribute to variation in human height remain undiscovered. We conducted a family-based linkage study of height in a unique cohort of very large nuclear families from a founder (Jewish) population. This design allowed for increased power to detect linkage, compared to previous family-based studies. Loci we identified in discovery families could explain an estimated lower bound of 6% of the variance in height in validation families. We showed that these loci are not tagging known common variants associated with height. Rather, we suggest that the observed signals arise from variants with large effects that are rare globally but elevated in frequency in the Jewish population.

BackgroundGeographical allocation of interventions focusing on hotspots of human immunodeficiency virus (HIV) transmission has the potential to improve efficiency. We used phylogeographic analyses to identify hotspots of the HIV transmission in Cologne-Bonn, Germany.MethodsWe included 714 HIV-1 infected individuals, followed up at the University Hospitals Cologne and Bonn. Distance-based molecular network analyses were performed to infer putative relationships. Characteristics of genetically linked individuals and assortativity (shared characteristics) were analyzed. Geospatial diffusion (ie, viral gene flow) was evaluated using a Slatkin-Maddison approach. Geospatial dispersal was determined by calculating the average distance between the residences of linked individuals (centroids of 3-digit zip code).ResultsIn sum, 217/714 (30.4%) sequences had a putative genetic linkage, forming 77 clusters (size range: 2-8). Linked individuals were more likely to live in areas surrounding the city center (P = .043),

Large structural variants (SVs) in the human genome are difficult to detect and study by conventional sequencing technologies. With long-range genome analysis platforms, such as optical mapping, one can identify large SVs (>2 kb) across the genome in one experiment. Analyzing optical genome maps of 154 individuals from the 26 populations sequenced in the 1000 Genomes Project, we find that phylogenetic population patterns of large SVs are similar to those of single nucleotide variations in 86% of the human genome, while ~2% of the genome has high structural complexity. We are able to characterize SVs in many intractable regions of the genome, including segmental duplications and subtelomeric, pericentromeric, and acrocentric areas. In addition, we discover ~60 Mb of non-redundant genome content missing in the reference genome sequence assembly. Our results highlight the need for a comprehensive set of alternate haplotypes from different populations to represent SV patterns in the genome.

In this study, we investigated low-frequency and rare variants associated with blood pressure (BP) by focusing on a linkage region on chromosome 16p13. We used whole genome sequencing (WGS) data obtained through the NHLBI Trans-Omics for Precision Medicine (TOPMed) program on 395 Cleveland Family Study (CFS) European Americans (CFS-EA). By analyzing functional coding variants and non-coding rare variants with CADD score > 10 residing within the chromosomal region in families with linkage evidence, we observed 25 genes with nominal statistical evidence (burden or SKAT p

Primary ciliary dyskinesia (PCD) is a hereditary defect of motile cilia in humans and several domestic animal species. Typical clinical findings are chronic recurrent infections of the respiratory tract and fertility problems. We analyzed an Alaskan Malamute family, in which two out of six puppies were affected by PCD. The parents were unaffected suggesting autosomal recessive inheritance. Linkage and homozygosity mapping defined critical intervals comprising ~118 Mb. Whole genome sequencing of one case and comparison to 601 control genomes identified a disease associated frameshift variant, c.43delA, in the NME5 gene encoding a sparsely characterized protein associated with ciliary function. Nme5-/- knockout mice exhibit doming of the skull, hydrocephalus and sperm flagellar defects. The genotypes at NME5:c.43delA showed the expected co-segregation with the phenotype in the Alaskan Malamute family. An additional unrelated Alaskan Malamute with PCD and hydrocephalus that became available later in the study was also homozygous mutant at the NME5:c.43delA variant. The mutant allele was not present in more than 1000 control dogs from different breeds. Immunohistochemistry demonstrated absence of the NME5 protein from nasal epithelia of an affected dog. We therefore propose NME5:c.43delA as the most likely candidate causative variant for PCD in Alaskan Malamutes. These findings enable genetic testing to avoid the unintentional breeding of affected dogs in the future. Furthermore, the results of this study identify NME5 as a novel candidate gene for unsolved human PCD and/or hydrocephalus cases.

Genomic datasets continue to increase in number due to the ease of production for a wider selection of species including non-model organisms. For many of these species, especially those with large or polyploid genomes, highly contiguous and well-annotated genomes are still rare due to the complexity and cost involved in their assembly. As a result, a common starting point for genomic work in non-model species is the production of a linkage map. Dense linkage maps facilitate the analysis of genomic data in a variety of ways, from broad scale observations regarding genome structure e.g., chromosome number and type or sex-related structural differences, to fine scale patterns e.g., recombination rate variation and co-localization of differentiated regions. Here we present both sex-averaged and sex-specific linkage maps for Coregonus sp. Albock, a member of the European whitefish lineage (C. lavaretus spp. complex), containing 5395 single nucleotide polymorphism (SNP) loci across 40 linkage groups to facilitate future investigation into the genomic basis of whitefish adaptation and speciation. The map was produced using restriction-site associated digestion (RAD) sequencing data from two wild-caught parents and 156 F1 offspring. We discuss the differences between our sex-averaged and sex-specific maps and identify genome-wide synteny between C. sp. Albock and Atlantic Salmon (Salmo salar), which have diverged following the salmonid-specific whole genome duplication. Our analysis confirms that many patterns of synteny observed between Atlantic Salmon and Oncorhynchus and Salvelinus species are also shared by members of the Coregoninae subfamily. We also show that regions known for their species-specific rediploidization history can pose challenges for synteny identification since these regions have diverged independently in each salmonid species following the salmonid-specific whole genome duplication. The European whitefish map provided here will enable future studies to understand the distribution of loci of interest, e.g., FST outliers, along the whitefish genome as well as assisting with the de novo assembly of a whitefish reference genome.

PurposeCongenital cataract, opacification of the ocular lens, is clinically and genetically a heterogeneous childhood disease. In this study we aimed to identify the underlying genetic cause of isolated autosomal-dominant lamellar cataract in a multi-generation English family.MethodsWhole-genome sequencing (WGS) was undertaken in two affected subjects and one unaffected individual. Segregation analysis was performed and a known cataract-causing mutation was identified. Segregation was further validated by sanger sequencing in the entire pedigree.ResultsA heterozygous mutation c.7 G > T; p.D3Y was identified in an NH2-terminal region of the gap junction protein GJA3 and found to co-segregate with disease.ConclusionWe have identified a recurrent mutation in GJA3 in a large British pedigree causing the novel phenotype of autosomal-dominant congenital lamellar cataract. Previously, p.D3Y was found in a Hispanic family causing pulverulent cataract. WGS proved an efficient method to find the underlying molecular cause in this large family, which could not be mapped due to uninformative markers.

The recent biochemical distinction between antibodies against leucine-rich, glioma-inactivated-1 (LGI1), contactin-associated protein-2 (CASPR2) and intracellular epitopes of voltage-gated potassium-channels (VGKCs) demands aetiological explanations. Given established associations between human leucocyte antigen (HLA) alleles and adverse drug reactions, and our clinical observation of frequent adverse drugs reactions in patients with LGI1 antibodies, we compared HLA alleles between healthy controls (n = 5553) and 111 Caucasian patients with VGKC-complex autoantibodies. In patients with LGI1 antibodies (n = 68), HLA-DRB1*07:01 was strongly represented [odds ratio = 27.6 (95% confidence interval 12.9-72.2), P = 4.1 × 10-26]. In contrast, patients with CASPR2 antibodies (n = 31) showed over-representation of HLA-DRB1*11:01 [odds ratio = 9.4 (95% confidence interval 4.6-19.3), P = 5.7 × 10-6]. Other allelic associations for patients with LGI1 antibodies reflected linkage, and significant haplotypic associations included HLA-DRB1*07:01-DQA1*02:01-DQB1*02:02, by comparison to DRB1*11:01-DQA1*05:01-DQB1*03:01 in CASPR2-antibody patients. Conditional analysis in LGI1-antibody patients resolved further independent class I and II associations. By comparison, patients with both LGI1 and CASPR2 antibodies (n = 3) carried yet another complement of HLA variants, and patients with intracellular VGKC antibodies (n = 9) lacked significant HLA associations. Within LGI1- or CASPR2-antibody patients, HLA associations did not correlate with clinical features. In silico predictions identified unique CASPR2- and LGI1-derived peptides potentially presented by the respective over-represented HLA molecules. These highly significant HLA associations dichotomize the underlying immunology in patients with LGI1 or CASPR2 antibodies, and inform T cell specificities and cellular interactions at disease initiation.10.1093/brain/awy109_video1awy109media15796480660001.

Two obsessive-compulsive disorder (OCD) genome-wide association studies (GWASs) have been published by independent OCD consortia, the International Obsessive-Compulsive Disorder Foundation Genetics Collaborative (IOCDF-GC) and the OCD Collaborative Genetics Association Study (OCGAS), but many of the top-ranked signals were supported in only one study. We therefore conducted a meta-analysis from the two consortia, investigating a total of 2688 individuals of European ancestry with OCD and 7037 genomically matched controls. No single-nucleotide polymorphisms (SNPs) reached genome-wide significance. However, in comparison with the two individual GWASs, the distribution of P-values shifted toward significance. The top haplotypic blocks were tagged with rs4733767 (P=7.1 × 10-7; odds ratio (OR)=1.21; confidence interval (CI): 1.12-1.31, CASC8/CASC11), rs1030757 (P=1.1 × 10-6; OR=1.18; CI: 1.10-1.26, GRID2) and rs12504244 (P=1.6 × 10-6; OR=1.18; CI: 1.11-1.27, KIT). Variants located in or near the genes ASB13, RSPO4, DLGAP1, PTPRD, GRIK2, FAIM2 and CDH20, identified in linkage peaks and the original GWASs, were among the top signals. Polygenic risk scores for each individual study predicted case-control status in the other by explaining 0.9% (P=0.003) and 0.3% (P=0.0009) of the phenotypic variance in OCGAS and the European IOCDF-GC target samples, respectively. The common SNP heritability in the combined OCGAS and IOCDF-GC sample was estimated to be 0.28 (s.e.=0.04). Strikingly, ∼65% of the SNP-based heritability in the OCGAS sample was accounted for by SNPs with minor allele frequencies of ⩾40%. This joint analysis constituting the largest single OCD genome-wide study to date represents a major integrative step in elucidating the genetic causes of OCD.

A genome-wide association study (GWAS) of 94,674 ancestrally diverse Kaiser Permanente members using 478,866 longitudinal electronic health record (EHR)-derived measurements for untreated serum lipid levels empowered multiple new findings: 121 new SNP associations (46 primary, 15 conditional, and 60 in meta-analysis with Global Lipids Genetic Consortium data); an increase of 33-42% in variance explained with multiple measurements; sex differences in genetic impact (greater impact in females for LDL, HDL, and total cholesterol and the opposite for triglycerides); differences in variance explained among non-Hispanic whites, Latinos, African Americans, and East Asians; genetic dominance and epistatic interaction, with strong evidence for both at the ABO and FUT2 genes for LDL; and tissue-specific enrichment of GWAS-associated SNPs among liver, adipose, and pancreas eQTLs. Using EHR pharmacy data, both LDL and triglyceride genetic risk scores (477 SNPs) were strongly predictive of age at initiation of lipid-lowering treatment. These findings highlight the value of longitudinal EHRs for identifying new genetic features of cholesterol and lipoprotein metabolism with implications for lipid treatment and risk of coronary heart disease.

Background:The Caribbean vervet monkey (Chlorocebus aethiops sabaeus) is a potentially valuable animal model of neurodegenerative disease. However, the trajectory of aging in vervets and its relationship to human disease is incompletely understood. Methods:To characterize biomarkers associated with neurodegeneration, we measured cerebrospinal fluid (CSF) concentrations of Aβ1-40, Aβ1-42, total tau, and p-tau181 in 329 members of a multigenerational pedigree. Linkage and genome-wide association were used to elucidate a genetic contribution to these traits. Results:Aβ1-40 concentrations were significantly correlated with age, brain total surface area, and gray matter thickness. Levels of p-tau181 were associated with cerebral volume and brain total surface area. Among the measured analytes, only CSF Aβ1-40 was heritable. No significant linkage (LOD > 3.3) was found, though suggestive linkage was highlighted on chromosomes 4 and 12. Genome-wide association identified a suggestive locus near the chromosome 4 linkage peak. Conclusions:Overall, these results support the vervet as a non-human primate model of amyloid-related neurodegeneration, such as Alzheimer's disease and cerebral amyloid angiopathy, and highlight Aβ1-40 and p-tau181 as potentially valuable biomarkers of these processes.

Genome-wide association studies have identified numerous variants associated with lipid levels; yet, the majority are located in non-coding regions with unclear mechanisms. In the Insulin Resistance Atherosclerosis Family Study (IRASFS), heritability estimates suggest a strong genetic basis: low-density lipoprotein (LDL, h2 = 0.50), high-density lipoprotein (HDL, h2 = 0.57), total cholesterol (TC, h2 = 0.53), and triglyceride (TG, h2 = 0.42) levels. Exome sequencing of 1,205 Mexican Americans (90 pedigrees) from the IRASFS identified 548,889 variants and association and linkage analyses with lipid levels were performed. One genome-wide significant signal was detected in APOA5 with TG (rs651821, PTG = 3.67 × 10-10, LODTG = 2.36, MAF = 14.2%). In addition, two correlated SNPs (r2 = 1.0) rs189547099 (PTG = 6.31 × 10-08, LODTG = 3.13, MAF = 0.50%) and chr4:157997598 (PTG = 6.31 × 10-08, LODTG = 3.13, MAF = 0.50%) reached exome-wide significance (P  3 with the strongest signal at rs1141070 (LODLDL = 4.30, PLDL = 0.33, MAF = 21.6%) in DFFB. A total of 53 nominally associated variants (P

Vernal pool clam shrimp (Eulimnadia texana) are a promising model system due to their ease of lab culture, short generation time, modest sized genome, a somewhat rare stable androdioecious sex determination system, and a requirement to reproduce via desiccated diapaused eggs. We generated a highly contiguous genome assembly using 46× of PacBio long read data and 216× of Illumina short reads, and annotated using Illumina RNAseq obtained from adult males or hermaphrodites. Of the 120 Mb genome 85% is contained in the largest eight contigs, the smallest of which is 4.6 Mb. The assembly contains 98% of transcripts predicted via RNAseq. This assembly is qualitatively different from scaffolded Illumina assemblies: It is produced from long reads that contain sequence data along their entire length, and is thus gap free. The contiguity of the assembly allows us to order the HOX genes within the genome, identifying two loci that contain HOX gene orthologs, and which approximately maintain the order observed in other arthropods. We identified a partial duplication of the Antennapedia complex adjacent to the few genes homologous to the Bithorax locus. Because the sex chromosome of an androdioecious species is of special interest, we used existing allozyme and microsatellite markers to identify the E. texana sex chromosome, and find that it comprises nearly half of the genome of this species. Linkage patterns indicate that recombination is extremely rare and perhaps absent in hermaphrodites, and as a result the location of the sex determining locus will be difficult to refine using recombination mapping.

PurposeThe Pakistan Punjab population has been a rich source for identifying genes causing or contributing to autosomal recessive retinal degenerations (arRD). This study was carried out to delineate the genetic architecture of arRD in the Pakistani population.MethodsThe genetic origin of arRD in a total of 144 families selected only for having consanguineous marriages and multiple members affected with arRD was examined. Of these, causative mutations had been identified in 62 families while only the locus had been identified for an additional 15. The remaining 67 families were subjected to homozygosity exclusion mapping by screening of closely flanking microsatellite markers at 180 known candidate genes/loci followed by sequencing of the candidate gene for pathogenic changes.ResultsOf these 67 families subjected to homozygosity mapping, 38 showed homozygosity for at least one of the 180 regions, and sequencing of the corresponding genes showed homozygous cosegregating mutations in 27 families. Overall, mutations were detected in approximately 61.8 % (89/144) of arRD families tested, with another 10.4% (15/144) being mapped to a locus but without a gene identified.ConclusionsThese results suggest the involvement of unmapped novel genes in the remaining 27.8% (40/144) of families. In addition, this study demonstrates that homozygosity mapping remains a powerful tool for identifying the genetic defect underlying genetically heterogeneous arRD disorders in consanguineous marriages for both research and clinical applications.

Insertions and deletions (INDELs) represent a significant fraction of interindividual variation in the human genome yet their contribution to phenotypes is poorly understood. To confirm the quality of imputed INDELs and investigate their roles in mediating cardiometabolic phenotypes, genome-wide association and linkage analyses were performed for 15 phenotypes with 1,273,952 imputed INDELs in 1,024 Mexican-origin Americans. Imputation quality was validated using whole exome sequencing with an average kappa of 0.93 in common INDELs (minor allele frequencies [MAFs] ≥ 5%). Association analysis revealed one genome-wide significant association signal for the cholesterylester transfer protein gene (CETP) with high-density lipoprotein levels (rs36229491, P = 3.06 × 10-12 ); linkage analysis identified two peaks with logarithm of the odds (LOD) > 5 (rs60560566, LOD = 5.36 with insulin sensitivity (SI ) and rs5825825, LOD = 5.11 with adiponectin levels). Suggestive overlapping signals between linkage and association were observed: rs59849892 in the WSC domain containing 2 gene (WSCD2) was associated and nominally linked with SI (P = 1.17 × 10-7 , LOD = 1.99). This gene has been implicated in glucose metabolism in human islet cell expression studies. In addition, rs201606363 was linked and nominally associated with low-density lipoprotein (P = 4.73 × 10-4 , LOD = 3.67), apolipoprotein B (P = 1.39 × 10-3 , LOD = 4.64), and total cholesterol (P = 1.35 × 10-2 , LOD = 3.80) levels. rs201606363 is an intronic variant of the UBE2F-SCLY (where UBE2F is ubiquitin-conjugating enzyme E2F and SCLY is selenocysteine lyase) fusion gene that may regulate cholesterol through selenium metabolism. In conclusion, these results confirm the feasibility of imputing INDELs from array-based single nucleotide polymorphism (SNP) genotypes. Analysis of these variants using association and linkage replicated previously identified SNP signals and identified multiple novel INDEL signals. These results support the inclusion of INDELs into genetic studies to more fully interrogate the spectrum of genetic variation.

Sex influences risk for opioid dependence (OD). We hypothesized that sex might interact with genetic loci that influence the risk for OD. Therefore we performed an analysis to identify sex-specific genomic susceptibility regions for OD using linkage. Over 6,000 single nucleotide polymorphism (SNP) markers were genotyped for 1,758 African- and European-American (AA and EA) individuals from 739 families, ascertained via affected sib-pairs with OD and/or cocaine dependence. Autosomewide non-parametric linkage scans, stratified by sex and population, were performed. We identified one significant linkage region, segregating with OD in EA men, at 71.1 cM on chromosome 4 (LOD = 3.29; point-wise P = 0.00005; empirical autosome-wide P = 0.042), which significantly differed from the linkage signal at the same location in EA women (empirical P = 0.002). Three suggestive linkage signals were identified at 181.3 cM on chromosome 7 (LOD = 2.18), 104 cM on chromosome 11 (LOD = 1.85), and 60.9 cM on chromosome 16 (LOD = 1.93) in EA women. In AA men, four suggestive linkage signals were detected at 201.1 cM on chromosome 3 (LOD = 2.32), 152.9 cM on chromosome 6 (LOD = 1.86), 16.8 cM on chromosome 7 (LOD = 1.95), and 36.1 cM on chromosome 17 (LOD = 1.99). The significant region, mapping to 4q12-4q13.1, harbors several OD candidate genes with interconnected functionality, including VEGFR, CLOCK, PDCL2, NMU, NRSF, and IGFBP7. In conclusion, these results provide an evidence for the existence of sex-specific and population-specific differences in OD. Furthermore, these results provide positional information that will facilitate the use of targeted next-generation sequencing to search for genes that contribute to sex-specific differences in OD. © 2016 Wiley Periodicals, Inc.

Linkage studies of complex genetic diseases have been largely replaced by genome-wide association studies, due in part to limited success in complex trait discovery. However, recent interest in rare and low-frequency variants motivates re-examination of family-based methods. In this study, we investigated the performance of two-point linkage analysis for over 1.6 million single-nucleotide polymorphisms (SNPs) combined with single variant association analysis to identify high impact variants, which are both strongly linked and associated with cardiometabolic traits in up to 1414 Hispanics from the Insulin Resistance Atherosclerosis Family Study (IRASFS). Evaluation of all 50 phenotypes yielded 83 557 000 LOD (logarithm of the odds) scores, with 9214 LOD scores ⩾3.0, 845 ⩾4.0 and 89 ⩾5.0, with a maximal LOD score of 6.49 (rs12956744 in the LAMA1 gene for tumor necrosis factor-α (TNFα) receptor 2). Twenty-seven variants were associated with P4, including variants in the NFIB gene under a linkage peak with TNFα receptor 2 levels on chromosome 9. Linkage regions of interest included a broad peak (31 Mb) on chromosome 1q with acute insulin response (max LOD=5.37). This region was previously documented with type 2 diabetes in family-based studies, providing support for the validity of these results. Overall, we have demonstrated the utility of two-point linkage and association in comprehensive genome-wide array-based SNP genotypes.

Genetic determinants of sleep-disordered breathing (SDB), a common set of disorders that contribute to significant cardiovascular and neuropsychiatric morbidity, are not clear. Overnight nocturnal oxygen saturation (SaO2) is a clinically relevant and easily measured indicator of SDB severity but its genetic contribution has never been studied. Our recent study suggests nocturnal SaO2 is heritable. We performed linkage analysis, association analysis and haplotype analysis of average nocturnal oxyhaemoglobin saturation in participants in the Cleveland Family Study (CFS), followed by gene-based association and additional tests in four independent samples. Linkage analysis identified a peak (LOD = 4.29) on chromosome 8p23. Follow-up association analysis identified two haplotypes in angiopoietin-2 (ANGPT2) that significantly contributed to the variation of SaO2 (P = 8 × 10-5) and accounted for a portion of the linkage evidence. Gene-based association analysis replicated the association of ANGPT2 and nocturnal SaO2. A rare missense SNP rs200291021 in ANGPT2 was associated with serum angiopoietin-2 level (P = 1.29 × 10-4), which was associated with SaO2 (P = 0.002). Our study provides the first evidence for the association of ANGPT2, a gene previously implicated in acute lung injury syndromes, with nocturnal SaO2, suggesting that this gene has a broad range of effects on gas exchange, including influencing oxygenation during sleep.

In many taxa, sex chromosomes are heteromorphic and largely non-recombining. Evolutionary models predict that spread of recombination suppression on the Y chromosome is fueled by the accumulation of sexually antagonistic alleles in close linkage to the sex determination region. However, empirical evidence for the existence of sexually antagonistic alleles is scarce. In the mosquito Aedes aegypti, the sex-determining chromosomes are homomorphic. The region of suppressed recombination, which surrounds the male-specific sex-determining gene, remains very small, despite ancient origin of the sex chromosomes in the Aedes lineage. We conducted a genetic analysis of the A. aegypti chromosome region tightly linked to the sex locus. We used a strain with an enhanced green fluorescent protein (EGFP)-tagged transgene inserted near the male-determining gene to monitor crossing-over events close to the boundary of the sex-determining region (SDR), and to trace the inheritance pattern of the transgene in relation to sex. In a series of crossing experiments involving individuals with a recombinant sex chromosome we found developmental abnormalities leading to 1:2 sex biases, caused by lethality of half of the male or female progeny. Our results suggest that various factors causing sex-specific lethal effects are clustered within the neighborhood of the SDR, which in the affected sex are likely lost or gained through recombination, leading to death. These may include genes that are recessive lethal, vital for development and/or sexually antagonistic. The sex chromosome fragment in question represents a fascinating test case for the analysis of processes that shape stable boundaries of a non-recombining region.

The patterns of sex determination and sex differentiation have been shown to differ among geographic populations of common frogs. Notably, the association between phenotypic sex and linkage group 2 (LG2) has been found to be perfect in a northern Swedish population, but weak and variable among families in a southern one. By analyzing these populations with markers from other linkage groups, we bring two new insights: (1) the variance in phenotypic sex not accounted for by LG2 in the southern population could not be assigned to genetic factors on other linkage groups, suggesting an epigenetic component to sex determination; (2) a second linkage group (LG7) was found to co-segregate with sex and LG2 in the northern population. Given the very short timeframe since post-glacial colonization (in the order of 1000 generations) and its seemingly localized distribution, this neo-sex chromosome system might be the youngest one described so far. It does not result from a fusion, but more likely from a reciprocal translocation between the original Y chromosome (LG2) and an autosome (LG7), causing their co-segregation during male meiosis. By generating a strict linkage between several important genes from the sex-determination cascade (Dmrt1, Amh and Amhr2), this neo-sex chromosome possibly contributes to the 'differentiated sex race' syndrome (strictly genetic sex determination and early gonadal development) that characterizes this northern population.

High-resolution genetic and physical maps are invaluable tools for building accurate genome assemblies, and interpreting results of genome-wide association studies (GWAS). Previous genetic and physical maps anchored good quality draft assemblies of the domestic cat genome, enabling the discovery of numerous genes underlying hereditary disease and phenotypes of interest to the biomedical science and breeding communities. However, these maps lacked sufficient marker density to order thousands of shorter scaffolds in earlier assemblies, which instead relied heavily on comparative mapping with related species. A high-resolution map would aid in validating and ordering chromosome scaffolds from existing and new genome assemblies. Here, we describe a high-resolution genetic linkage map of the domestic cat genome based on genotyping 453 domestic cats from several multi-generational pedigrees on the Illumina 63K SNP array. The final maps include 58,055 SNP markers placed relative to 6637 markers with unique positions, distributed across all autosomes and the X chromosome. Our final sex-averaged maps span a total autosomal length of 4464 cM, the longest described linkage map for any mammal, confirming length estimates from a previous microsatellite-based map. The linkage map was used to order and orient the scaffolds from a substantially more contiguous domestic cat genome assembly (Felis catus v8.0), which incorporated ∼20 × coverage of Illumina fragment reads. The new genome assembly shows substantial improvements in contiguity, with a nearly fourfold increase in N50 scaffold size to 18 Mb. We use this map to report probable structural errors in previous maps and assemblies, and to describe features of the recombination landscape, including a massive (∼50 Mb) recombination desert (of virtually zero recombination) on the X chromosome that parallels a similar desert on the porcine X chromosome in both size and physical location.

Identifying homology between sex chromosomes of different species is essential to understanding the evolution of sex determination. Here, we show that the identity of a homomorphic sex chromosome pair can be established using a linkage map, without information on offspring sex. By comparing sex-specific maps of the European tree frog Hyla arborea, we find that the sex chromosome (linkage group 1) shows a threefold difference in marker number between the male and female maps. In contrast, the number of markers on each autosome is similar between the two maps. We also find strongly conserved synteny between H. arborea and Xenopus tropicalis across 200 million years of evolution, suggesting that the rate of chromosomal rearrangement in anurans is low. Finally, we show that recombination in males is greatly reduced at the centers of large chromosomes, consistent with previous cytogenetic findings. Our research shows the importance of high-density linkage maps for studies of recombination, chromosomal rearrangement and the genetic architecture of ecologically or economically important traits.

Three strikingly different alternative male mating morphs (aggressive 'independents', semicooperative 'satellites' and female-mimic 'faeders') coexist as a balanced polymorphism in the ruff, Philomachus pugnax, a lek-breeding wading bird. Major differences in body size, ornamentation, and aggressive and mating behaviors are inherited as an autosomal polymorphism. We show that development into satellites and faeders is determined by a supergene consisting of divergent alternative, dominant and non-recombining haplotypes of an inversion on chromosome 11, which contains 125 predicted genes. Independents are homozygous for the ancestral sequence. One breakpoint of the inversion disrupts the essential CENP-N gene (encoding centromere protein N), and pedigree analysis confirms the lethality of homozygosity for the inversion. We describe new differences in behavior, testis size and steroid metabolism among morphs and identify polymorphic genes within the inversion that are likely to contribute to the differences among morphs in reproductive traits.

PurposeThis study was undertaken to identify causal mutations responsible for autosomal recessive retinitis pigmentosa (arRP) in consanguineous families.MethodsLarge consanguineous families were ascertained from the Punjab province of Pakistan. An ophthalmic examination consisting of a fundus evaluation and electroretinography (ERG) was completed, and small aliquots of blood were collected from all participating individuals. Genomic DNA was extracted from white blood cells, and a genome-wide linkage or a locus-specific exclusion analysis was completed with polymorphic short tandem repeats (STRs). Two-point logarithm of odds (LOD) scores were calculated, and all coding exons and exon-intron boundaries of RP1 were sequenced to identify the causal mutation.ResultsThe ophthalmic examination showed that affected individuals in all families manifest cardinal symptoms of RP. Genome-wide scans localized the disease phenotype to chromosome 8q, a region harboring RP1, a gene previously implicated in the pathogenesis of RP. Sanger sequencing identified a homozygous single base deletion in exon 4: c.3697delT (p.S1233Pfs22*), a single base substitution in intron 3: c.787+1G>A (p.I263Nfs8*), a 2 bp duplication in exon 2: c.551_552dupTA (p.Q185Yfs4*) and an 11,117 bp deletion that removes all three coding exons of RP1. These variations segregated with the disease phenotype within the respective families and were not present in ethnically matched control samples.ConclusionsThese results strongly suggest that these mutations in RP1 are responsible for the retinal phenotype in affected individuals of all four consanguineous families.

BackgroundMouse chromosome 2 is linked to growth and body fat phenotypes in many mouse crosses. With the goal to identify the underlying genes regulating growth and body fat on mouse chromosome 2, we developed five overlapping subcongenic strains that contained CAST/EiJ donor regions in a C57BL/6J (hg/hg) background (hg is a spontaneous deletion of 500 Kb on mouse chromosome 10). To fine map QTL on distal mouse chromosome 2 a total of 1,712 F2 mice from the five subcongenic strains, plus 278 F2 mice from the HG2D founder congenic strain were phenotyped and analyzed. Interval mapping (IM) and composite IM (CIM) were performed on body weight and body fat traits on a combination of SNP and microsatellite markers, which generated a high-density genotyping panel.ResultsPhenotypic analysis and interval mapping of total fat mass identified two QTL on distal mouse chromosome 2. One QTL between 150 and 161 Mb, Fatq2a, and the second between 173.3 and 175.6 Mb, Fatq2b. The two QTL reside in different congenic strains with significant total fat differences between homozygous cast/cast and b6/b6 littermates. Both of these QTL were previously identified only as a single QTL affecting body fat, Fatq2. Furthermore, through a novel approach referred here as replicated CIM, Fatq2b was mapped to the Gnas imprinted locus.ConclusionsThe integration of subcongenic strains, high-density genotyping, and CIM succesfully partitioned two previously linked QTL 20 Mb apart, and the strongest QTL, Fatq2b, was fine mapped to a ~2.3 Mb region interval encompassing the Gnas imprinted locus.

BackgroundIn order to study the genetics of diseases more accurately and effectively, one often collects large families. Different members of a large family may provide differing information about the structure and make-up of their pedigree. Thus, software is needed to facilitate reconciliation of pedigrees collected independently from multiple informants from a single large family to create a unified pedigree that is based on the best composite information available.FindingsOur implementation demonstrates that Genetic ME performs merging in terms of adding, replacing and combining information from two pedigrees. Through a tracking process, all of the changes made to the data set for the individuals can be traced back to their original source material. A new pedigree structure can be easily visualized while reconciling disparate information from multiple pedigrees.MethodsWe developed the Genetic Merging & Editing (Genetic ME) program, an open source Java application built on top of CraneFoot and Ghostscript, to support comparing, editing and merging of pedigrees collected from multiple sources in a visually-oriented manner.ConclusionsGenetic ME constitutes an ideal addition to software packages for reconciling pedigree information from multiple sources. Genetic ME provides a friendly graphical user interface, traces the changes made by users, and produces viewable merged pedigree structures able to be further used by other popular analysis programs.

Arsenic toxicokinetics are important for disease risks in exposed populations, but genetic determinants are not fully understood. We examined urine arsenic species patterns measured by HPLC-ICPMS among 2189 Strong Heart Study participants 18 years of age and older with data on ~400 genome-wide microsatellite markers spaced ~10 cM and arsenic speciation (683 participants from Arizona, 684 from Oklahoma, and 822 from North and South Dakota). We logit-transformed % arsenic species (% inorganic arsenic, %MMA, and %DMA) and also conducted principal component analyses of the logit % arsenic species. We used inverse-normalized residuals from multivariable-adjusted polygenic heritability analysis for multipoint variance components linkage analysis. We also examined the contribution of polymorphisms in the arsenic metabolism gene AS3MT via conditional linkage analysis. We localized a quantitative trait locus (QTL) on chromosome 10 (LOD 4.12 for %MMA, 4.65 for %DMA, and 4.84 for the first principal component of logit % arsenic species). This peak was partially but not fully explained by measured AS3MT variants. We also localized a QTL for the second principal component of logit % arsenic species on chromosome 5 (LOD 4.21) that was not evident from considering % arsenic species individually. Some other loci were suggestive or significant for 1 geographical area but not overall across all areas, indicating possible locus heterogeneity. This genome-wide linkage scan suggests genetic determinants of arsenic toxicokinetics to be identified by future fine-mapping, and illustrates the utility of principal component analysis as a novel approach that considers % arsenic species jointly.

Individuals with the short variant of the serotonin transporter linked polymorphic region gene are more susceptible than individuals homozygous for the long allele to the effects of stressful life events on risk for internalizing and externalizing problems. We tested whether individual differences in coping style explained this increased risk for problem behavior among youth who were at both genetic and environmental risk. Participants included 279 children, ages 8-11, from the Children's Experiences and Development Study. Caregivers and teachers reported on children's internalizing and externalizing symptoms, and caregivers and children on children's exposure to harsh parenting and parental warmth in middle childhood, and traumatic events. Children reported how frequently they used various coping strategies. Results revealed that short/short homozygotes had higher levels of internalizing problems compared with long allele carriers and that short allele carriers had higher levels of externalizing problems compared with long/long homozygotes under conditions of high cumulative risk. Moreover, among children who were homozygous for the short allele, those who had more cumulative risk indicators less frequently used distraction coping strategies, which partly explained why they had higher levels of internalizing problems. Coping strategies did not significantly mediate Gene × Environment effects on externalizing symptoms.

Cerebral dominance of language function and hand preference are suggested to be heritable traits with possible shared genetic background. However, joined genetic studies of these traits have never been conducted. We performed a genetic linkage study in 37 multigenerational human pedigrees of both sexes (consisting of 355 subjects) enriched with left-handedness in which we also measured language lateralization. Hand preference was measured with the Edinburgh Handedness Inventory, and language lateralization was measured with functional transcranial Doppler during language production. The estimated heritability of left-handedness and language lateralization in these pedigrees is 0.24 and 0.31, respectively. A parametric major gene model was tested for left-handedness. Nonparametric analyses were performed for left-handedness, atypical lateralization, and degree of language lateralization. We did not observe genome-wide evidence for linkage in the parametric or nonparametric analyses for any of the phenotypes tested. However, multiple regions showed suggestive evidence of linkage. The parametric model showed suggestive linkage for left-handedness in the 22q13 region [heterogeneity logarithm of odds (HLOD) = 2.18]. Nonparametric multipoint analysis of left-handedness showed suggestive linkage in the same region [logarithm of odds (LOD) = 2.80]. Atypical language lateralization showed suggestive linkage in the 7q34 region (LODMax = 2.35). For strength of language lateralization, we observed suggestive linkage in the 6p22 (LODMax = 2.54), 7q32 (LODMax = 1.93), and 9q33 (LODMax = 2.10) regions. We did not observe any overlap of suggestive genetic signal between handedness and the extent of language lateralization. The absence of significant linkage argues against the presence of a major gene coding for both traits; rather, our results are suggestive of these traits being two independent polygenic complex traits.

ObjectiveTo identify novel genetic markers of obesity-related traits and to identify gene-diet interactions with n-3 polyunsaturated fatty acid (n-3 PUFA) intake in Yup'ik people.Material and methodsWe measured body composition, plasma adipokines and ghrelin in 982 participants enrolled in the Center for Alaska Native Health Research (CANHR) Study. We conducted a genome-wide SNP linkage scan and targeted association analysis, fitting additional models to investigate putative gene-diet interactions. Finally, we performed bioinformatic analysis to uncover likely candidate genes within the identified linkage peaks.ResultsWe observed evidence of linkage for all obesity-related traits, replicating previous results and identifying novel regions of interest for adiponectin (10q26.13-2) and thigh circumference (8q21.11-13). Bioinformatic analysis revealed DOCK1, PTPRE (10q26.13-2) and FABP4 (8q21.11-13) as putative candidate genes in the newly identified regions. Targeted SNP analysis under the linkage peaks identified associations between three SNPs and obesity-related traits: rs1007750 on chromosome 8 and thigh circumference (P=0.0005), rs878953 on chromosome 5 and thigh skinfold (P=0.0004), and rs1596854 on chromosome 11 for waist circumference (P=0.0003). Finally, we showed that n-3 PUFA modified the association between obesity related traits and two additional variants (rs2048417 on chromosome 3 for adiponectin, P for interaction=0.0006 and rs730414 on chromosome 11 for percentage body fat, P for interaction=0.0004).ConclusionsThis study presents evidence of novel genomic regions and gene-diet interactions that may contribute to the pathophysiology of obesity-related traits among Yup'ik people.

Tourette syndrome (TS) is a neuropsychiatric disorder characterized by recurrent motor and vocal tics, often accompanied by obsessive-compulsive disorder and/or attention-deficit/hyperactivity disorder. While the evidence for a genetic contribution is strong, its exact nature has yet to be clarified fully. There is now mounting evidence that the genetic risks for TS include both common and rare variants and may involve complex multigenic inheritance or, in rare cases, a single major gene. Based on recent progress in many other common disorders with apparently similar genetic architectures, it is clear that large patient cohorts and open-access repositories will be essential to further advance the field. To that end, the large multicenter Tourette International Collaborative Genetics (TIC Genetics) study was established. The goal of the TIC Genetics study is to undertake a comprehensive gene discovery effort, focusing both on familial genetic variants with large effects within multiply affected pedigrees and on de novo mutations ascertained through the analysis of apparently simplex parent-child trios with non-familial tics. The clinical data and biomaterials (DNA, transformed cell lines, RNA) are part of a sharing repository located within the National Institute for Mental Health Center for Collaborative Genomics Research on Mental Disorders, USA, and will be made available to the broad scientific community. This resource will ultimately facilitate better understanding of the pathophysiology of TS and related disorders and the development of novel therapies. Here, we describe the objectives and methods of the TIC Genetics study as a reference for future studies from our group and to facilitate collaboration between genetics consortia in the field of TS.

We previously identified a low-frequency (1.1 %) coding variant (G45R; rs200573126) in the adiponectin gene (ADIPOQ) which was the basis for a multipoint microsatellite linkage signal (LOD = 8.2) for plasma adiponectin levels in Hispanic families. We have empirically evaluated the ability of data from targeted common variants, exome chip genotyping, and genome-wide association study data to detect linkage and association to adiponectin protein levels at this locus. Simple two-point linkage and association analyses were performed in 88 Hispanic families (1,150 individuals) using 10,958 SNPs on chromosome 3. Approaches were compared for their ability to map the functional variant, G45R, which was strongly linked (two-point LOD = 20.98) and powerfully associated (p value = 8.1 × 10(-50)). Over 450 SNPs within a broad 61 Mb interval around rs200573126 showed nominal evidence of linkage (LOD > 3) but only four other SNPs in this region were associated with p values

Cave animals, compared to surface-dwelling relatives, tend to have reduced eyes and pigment, longer appendages, and enhanced mechanosensory structures. Pressing questions include how certain cave-related traits are gained and lost, and if they originate through the same or different genetic programs in independent lineages. An excellent system for exploring these questions is the isopod, Asellus aquaticus. This species includes multiple cave and surface populations that have numerous morphological differences between them. A key feature is that hybrids between cave and surface individuals are viable, which enables genetic crosses and linkage analyses. Here, we advance this system by analyzing single animal transcriptomes of Asellus aquaticus. We use high throughput sequencing of non-normalized cDNA derived from the head of a surface-dwelling male, the head of a cave-dwelling male, the head of a hybrid male (produced by crossing a surface individual with a cave individual), and a pooled sample of surface embryos and hatchlings. Assembling reads from surface and cave head RNA pools yielded an integrated transcriptome comprised of 23,984 contigs. Using this integrated assembly as a reference transcriptome, we aligned reads from surface-, cave- and hybrid- head tissue and pooled surface embryos and hatchlings. Our approach identified 742 SNPs and placed four new candidate genes to an existing linkage map for A. aquaticus. In addition, we examined SNPs for allele-specific expression differences in the hybrid individual. All of these resources will facilitate identification of genes and associated changes responsible for cave adaptation in A. aquaticus and, in concert with analyses of other species, will inform our understanding of the evolutionary processes accompanying adaptation to the subterranean environment.

Here, we provide revised gene models for D. ananassae, D. yakuba, and D. simulans, which include untranslated regions and empirically verified intron-exon boundaries, as well as ortholog groups identified using a fuzzy reciprocal-best-hit blast comparison. Using these revised annotations, we perform differential expression testing using the cufflinks suite to provide a broad overview of differential expression between reproductive tissues and the carcass. We identify thousands of genes that are differentially expressed across tissues in D. yakuba and D. simulans, with roughly 60% agreement in expression patterns of orthologs in D. yakuba and D. simulans. We identify several cases of putative polycistronic transcripts, pointing to a combination of transcriptional read-through in the genome as well as putative gene fusion and fission events across taxa. We furthermore identify hundreds of lineage specific genes in each species with no blast hits among transcripts of any other Drosophila species, which are candidates for neofunctionalized proteins and a potential source of genetic novelty.

Hereditary eye diseases of animals serve as excellent models of human ocular disorders and assist in the development of gene and drug therapies for inherited forms of blindness. Several primary hereditary eye conditions affecting various ocular tissues and having different rates of progression have been documented in domestic cats. Gene therapy for canine retinopathies has been successful, thus the cat could be a gene therapy candidate for other forms of retinal degenerations. The current study investigates a hereditary, autosomal recessive, retinal degeneration specific to Persian cats. A multi-generational pedigree segregating for this progressive retinal atrophy was genotyped using a 63 K SNP array and analyzed via genome-wide linkage and association methods. A multi-point parametric linkage analysis localized the blindness phenotype to a ~1.75 Mb region with significant LOD scores (Z ≈ 14, θ = 0.00) on cat chromosome E1. Genome-wide TDT, sib-TDT, and case-control analyses also consistently supported significant association within the same region on chromosome E1, which is homologous to human chromosome 17. Using haplotype analysis, a ~1.3 Mb region was identified as highly associated for progressive retinal atrophy in Persian cats. Several candidate genes within the region are reasonable candidates as a potential causative gene and should be considered for molecular analyses.

BackgroundRecurrent bacterial and fungal infections, eczema, and increased serum IgE levels characterize patients with the hyper-IgE syndrome (HIES). Known genetic causes for HIES are mutations in signal transducer and activator of transcription 3 (STAT3) and dedicator of cytokinesis 8 (DOCK8), which are involved in signal transduction pathways. However, glycosylation defects have not been described in patients with HIES. One crucial enzyme in the glycosylation pathway is phosphoglucomutase 3 (PGM3), which catalyzes a key step in the synthesis of uridine diphosphate N-acetylglucosamine, which is required for the biosynthesis of N-glycans.ObjectiveWe sought to elucidate the genetic cause in patients with HIES who do not carry mutations in STAT3 or DOCK8.MethodsAfter establishing a linkage interval by means of SNPchip genotyping and homozygosity mapping in 2 families with HIES from Tunisia, mutational analysis was performed with selector-based, high-throughput sequencing. Protein expression was analyzed by means of Western blotting, and glycosylation was profiled by using mass spectrometry.ResultsMutational analysis of candidate genes in an 11.9-Mb linkage region on chromosome 6 shared by 2 multiplex families identified 2 homozygous mutations in PGM3 that segregated with disease status and followed recessive inheritance. The mutations predict amino acid changes in PGM3 (p.Glu340del and p.Leu83Ser). A third homozygous mutation (p.Asp502Tyr) and the p.Leu83Ser variant were identified in 2 other affected families, respectively. These hypomorphic mutations have an effect on the biosynthetic reactions involving uridine diphosphate N-acetylglucosamine. Glycomic analysis revealed an aberrant glycosylation pattern in leukocytes demonstrated by a reduced level of tri-antennary and tetra-antennary N-glycans. T-cell proliferation and differentiation were impaired in patients. Most patients had developmental delay, and many had psychomotor retardation.ConclusionImpairment of PGM3 function leads to a novel primary (inborn) error of development and immunity because biallelic hypomorphic mutations are associated with impaired glycosylation and a hyper-IgE-like phenotype.

Linkage analysis of complex traits has had limited success in identifying trait-influencing loci. Recently, coding variants have been implicated as the basis for some biomedical associations. We tested whether coding variants are the basis for linkage peaks of complex traits in 42 African-American (n = 596) and 90 Hispanic (n = 1,414) families in the Insulin Resistance Atherosclerosis Family Study (IRASFS) using Illumina HumanExome Beadchips. A total of 92,157 variants in African Americans (34%) and 81,559 (31%) in Hispanics were polymorphic and tested using two-point linkage and association analyses with 37 cardiometabolic phenotypes. In African Americans 77 LOD scores greater than 3 were observed. The highest LOD score was 4.91 with the APOE SNP rs7412 (MAF = 0.13) with plasma apolipoprotein B (ApoB). This SNP was associated with ApoB (P-value = 4 × 10(-19)) and accounted for 16.2% of the variance in African Americans. In Hispanic families, 104 LOD scores were greater than 3. The strongest evidence of linkage (LOD = 4.29) was with rs5882 (MAF = 0.46) in CETP with HDL. CETP variants were strongly associated with HDL (0.00049 < P-value

ImportanceGenetic factors contribute to risk for bipolar disorder (BP), but its pathogenesis remains poorly understood. A focus on measuring multisystem quantitative traits that may be components of BP psychopathology may enable genetic dissection of this complex disorder, and investigation of extended pedigrees from genetically isolated populations may facilitate the detection of specific genetic variants that affect BP as well as its component phenotypes.ObjectiveTo identify quantitative neurocognitive, temperament-related, and neuroanatomical phenotypes that appear heritable and associated with severe BP (bipolar I disorder [BP-I]) and therefore suitable for genetic linkage and association studies aimed at identifying variants contributing to BP-I risk.Design, setting, and participantsMultigenerational pedigree study in 2 closely related, genetically isolated populations: the Central Valley of Costa Rica and Antioquia, Colombia. A total of 738 individuals, all from Central Valley of Costa Rica and Antioquia pedigrees, participated; among them, 181 have BP-I.Main outcomes and measuresFamilial aggregation (heritability) and association with BP-I of 169 quantitative neurocognitive, temperament, magnetic resonance imaging, and diffusion tensor imaging phenotypes.ResultsOf 169 phenotypes investigated, 126 (75%) were significantly heritable and 53 (31%) were associated with BP-I. About one-quarter of the phenotypes, including measures from each phenotype domain, were both heritable and associated with BP-I. Neuroimaging phenotypes, particularly cortical thickness in prefrontal and temporal regions as well as volume and microstructural integrity of the corpus callosum, represented the most promising candidate traits for genetic mapping related to BP based on strong heritability and association with disease. Analyses of phenotypic and genetic covariation identified substantial correlations among the traits, at least some of which share a common underlying genetic architecture.Conclusions and relevanceTo our knowledge, this is the most extensive investigation of BP-relevant component phenotypes to date. Our results identify brain and behavioral quantitative traits that appear to be genetically influenced and show a pattern of BP-I association within families that is consistent with expectations from case-control studies. Together, these phenotypes provide a basis for identifying loci contributing to BP-I risk and for genetic dissection of the disorder.

Congenital cataracts are an important cause of bilateral visual impairment in infants. Through genome-wide linkage analysis in a four-generation family of Irish descent, the disease-associated gene causing autosomal-dominant congenital nuclear cataract was mapped to chromosome 4p16.1. The maximum logarithm of odds (LOD) score was 2.62 at a recombination fraction θ=0, obtained for marker D4S432 physically close to the Wolfram gene (WFS1). By sequencing the coding regions and intron-exon boundaries of WFS1, we identified a DNA substitution (c.1385A-to-G) in exon 8, causing a missense mutation at codon 462 (E462G) of the Wolframin protein. This is the first report of a mutation in this gene causing an isolated nuclear congenital cataract. These findings suggest that the membrane trafficking protein Wolframin may be important for supporting the developing lens.

ObjectiveTo identify a genetic cause for migrating partial seizures in infancy (MPSI).MethodsWe characterized a consanguineous pedigree with MPSI and obtained DNA from affected and unaffected family members. We analyzed single nucleotide polymorphism 500K data to identify regions with evidence of linkage. We performed whole exome sequencing and analyzed homozygous variants in regions of linkage to identify a candidate gene and performed functional studies of the candidate gene SLC25A22.ResultsIn a consanguineous pedigree with 2 individuals with MPSI, we identified 2 regions of linkage, chromosome 4p16.1-p16.3 and chromosome 11p15.4-pter. Using whole exome sequencing, we identified 8 novel homozygous variants in genes in these regions. Only 1 variant, SLC25A22 c.G328C, results in a change of a highly conserved amino acid (p.G110R) and was not present in control samples. SLC25A22 encodes a glutamate transporter with strong expression in the developing brain. We show that the specific G110R mutation, located in a transmembrane domain of the protein, disrupts mitochondrial glutamate transport.InterpretationWe have shown that MPSI can be inherited and have identified a novel homozygous mutation in SLC25A22 in the affected individuals. Our data strongly suggest that SLC25A22 is responsible for MPSI, a severe condition with few known etiologies. We have demonstrated that a combination of linkage analysis and whole exome sequencing can be used for disease gene discovery. Finally, as SLC25A22 had been implicated in the distinct syndrome of neonatal epilepsy with suppression bursts on electroencephalogram, we have expanded the phenotypic spectrum associated with SLC25A22.

BackgroundThe many attempts that have been made to identify genes for bipolar disorder (BD) have met with limited success, which may reflect an inadequacy of diagnosis as an informative and biologically relevant phenotype for genetic studies. Here we have explored aspects of personality as quantitative phenotypes for bipolar disorder through the use of the Temperament and Character Inventory (TCI), which assesses personality in seven dimensions. Four temperament dimensions are assessed: novelty seeking (NS), harm avoidance (HA), reward dependence (RD), and persistence (PS). Three character dimensions are also included: self-directedness (SD), cooperativeness (CO), and self-transcendence (ST).MethodsWe compared personality scores between diagnostic groups and assessed heritability in a sample of 101 families collected for genetic studies of BD. A genome-wide SNP linkage analysis was then performed in the subset of 51 families for which genetic data was available.ResultsSignificant group differences were observed between BD subjects, their first-degree relatives, and independent controls for all but RD and PS, and all but HA and RD were found to be significantly heritable in this sample. Linkage analysis of the heritable dimensions produced several suggestive linkage peaks for NS (chromosomes 7q21 and 10p15), PS (chromosomes 6q16, 12p13, and 19p13), and SD (chromosomes 4q35, 8q24, and 18q12).LimitationsThe relatively small size of our linkage sample likely limited our ability to reach genome-wide significance in this study.ConclusionsWhile not genome-wide significant, these results suggest that aspects of personality may prove useful in the identification of genes underlying BD susceptibility.

We report an analysis of allele-specific expression (ASE) and parent-of-origin expression in adult mouse liver using next generation sequencing (RNA-Seq) of reciprocal crosses of heterozygous F1 mice from the parental strains C57BL/6J and DBA/2J. We found a 60% overlap between genes exhibiting ASE and putative cis-acting expression quantitative trait loci (cis-eQTL) identified in an intercross between the same strains. We discuss the various biological and technical factors that contribute to the differences. We also identify genes exhibiting parental imprinting and complex expression patterns. Our study demonstrates the importance of biological replicates to limit the number of false positives with RNA-Seq data.

BackgroundThe many attempts to identify genes for bipolar disorder (BD) have met with limited success, which has generally been attributed to genetic heterogeneity and small gene effects. However, it is also possible that the categorical phenotypes used in genetic studies of BD are not the most informative or biologically relevant. We have explored aspects of temperament as quantitative phenotypes for BD through the use of the Temperament Evaluation of Memphis, Pisa, Paris, and San Diego Auto-questionnaire (TEMPS-A), which is designed to assess lifelong, milder aspects of bipolar symptomatology and defines five temperaments: hyperthymic, dysthymic, cyclothymic, irritable, and anxious.MethodsWe compared temperament scores between diagnostic groups and assessed heritability in a sample of 101 families collected for genetic studies of BD. A genome-wide SNP linkage study was then performed in the subset of 51 families for which genetic data was available.ResultsSignificant group differences were observed between BD subjects, their first-degree relatives, and independent controls, and all five temperaments were found to be significantly heritable, with heritabilities ranging from 21% for the hyperthymic to 52% for the irritable temperaments. Suggestive evidence for linkage was observed for the hyperthymic (chromosomes 1q44, 2p16, 6q16, and 14q23), dysthymic (chromosomes 3p21 and 13q34), and irritable (chromosome 6q24) temperaments.LimitationsThe relatively small size of our linkage sample likely limited our ability to reach genome-wide significance in this study.ConclusionsWhile not genome-wide significant, these results suggest that aspects of temperament may prove useful in the identification of genes underlying BD susceptibility.

ObjectivesTo identify genomic regions associated with fasting plasma lipid profiles, insulin, glucose, and glycosylated hemoglobin in a Yup'ik study population, and to evaluate whether the observed associations between genetic factors and metabolic traits were modified by dietary intake of marine derived omega-3 polyunsaturated acids (n-3 PUFA).MethodsA genome-wide linkage scan was conducted among 982 participants of the Center for Alaska Native Health Research study. n-3 PUFA intake was estimated using the nitrogen stable isotope ratio (δ(15) N) of erythrocytes. All genotyped SNPs located within genomic regions with LOD scores > 2 were subsequently tested for individual SNP associations with metabolic traits using linear models that account for familial correlation as well as age, sex, community group, and n-3 PUFA intake. Separate linear models were fit to evaluate interactions between the genotype of interest and n-3 PUFA intake.ResultsWe identified several chromosomal regions linked to serum apolipoprotein A2, high density lipoprotein-, low density lipoprotein-, and total cholesterol, insulin, and glycosylated hemoglobin. Genetic variants found to be associated with total cholesterol mapped to a region containing previously validated lipid loci on chromosome 19, and additional novel peaks of biological interest were identified at 11q12.2-11q13.2. We did not observe any significant interactions between n-3 PUFA intake, genotypes, and metabolic traits.ConclusionsWe have completed a whole genome linkage scan for metabolic traits in Native Alaskans, confirming previously identified loci, and offering preliminary evidence of novel loci implicated in chronic disease pathogenesis in this population.