Resequencing and fine-mapping of the chromosome 12q13-14 locus associated with multiple sclerosis refines the number of implicated genes

Multiple sclerosis (MS) is a common chronic inflammatory disease of the central nervous system. Susceptibility to the disease is affected by both environmental and genetic factors. Genetic factors include haplotypes in the histocompatibility complex (MHC) and over 50 non-MHC loci reported by genome-wide association studies. Amongst these, we previously reported polymorphisms in chromosome 12q13-14 with a protective effect in individuals of European descent. This locus spans 288 kb and contains 17 genes, including several candidate genes which have potentially significant pathogenic and therapeutic implications. In this study, we aimed to fine-map this locus. We have implemented a two-phase study: a variant discovery phase where we have used next-generation sequencing and two target-enrichment strategies [long-range polymerase chain reaction (PCR) and Nimblegen's solution phase hybridization capture] in pools of 25 samples; and a genotyping phase where we genotyped 712 variants in 3577 healthy controls and 3269 MS patients. This study confirmed the association (rs2069502, P = 9.9 × 10−11, OR = 0.787) and narrowed down the locus of association to an 86.5 kb region. Although the study was unable to pinpoint the key-associated variant, we have identified a 42 (genotyped and imputed) single-nucleotide polymorphism haplotype block likely to harbour the causal variant. No evidence of association at previously reported low-frequency variants in CYP27B1 was observed. As part of the study we compared variant discovery performance using two target-enrichment strategies. We concluded that our pools enriched with Nimblegen's solution phase hybridization capture had better sensitivity to detect true variants than the pools enriched with long-range PCR, whilst specificity was better in the long-range PCR-enriched pools compared with solution phase hybridization capture enriched pools; this result has important implications for the design of future fine-mapping studies.

Resequencing of the CCL5 and CCR5 genes and investigation of variants for association with diabetic nephropathy

Chemokine (C-C motif) ligand 5 (CCL5) and chemokine (C-C motif) receptor 5 are implicated in the pathogenesis of diabetic nephropathy (DN). We hypothesize that variants in these genes may be associated with DN. The CCL5 and chemokine receptor type 5 (CCR5) genes were resequenced, variants identified (n=58), allele frequencies determined in 46 individuals (92 chromosomes) and efficient haplotype tag single-nucleotide polymorphisms (htSNPs) selected to effectively evaluate the common variation in these genes. One reportedly functional gene variant and eight htSNPs were genotyped in a case-control association study involving Caucasian individuals with type 1 diabetes (267 cases with DN and 442 non-nephropathic diabetic controls). Genotyping was performed using MassARRAY iPLEX, TaqMan, gel electrophoresis and direct capillary sequencing. After correction for multiple testing, there were no statistically significant associations between variants in the CCL5 and CCR5 genes and DN. Journal of Human Genetics (2010) 55, 248-251; doi:10.1038/jhg.2010.15; published online 5 March 2010

Comprehensive Investigation of the Caveolin 2 Gene: Resequencing and Association for Kidney Transplant Outcomes

Caveolae are plasma membrane structures formed from a complex of the proteins caveolin-1 and caveolin-2. Caveolae interact with pro-inflammatory cytokines and are dysregulated in fibrotic disease. Although caveolae are present infrequently in healthy kidneys, they are abundant during kidney injury. An association has been identified between a CAV1 gene variant and long term kidney transplant survival. Chronic, gradual decline in transplant function is a persistent problem in kidney transplantation. The aetiology of this is diverse but fibrosis within the transplanted organ is the common end point. This study is the first to investigate the association of CAV2 gene variants with kidney transplant outcomes. Genomic DNA from donors and recipients of 575 kidney transplants performed in Belfast was investigated for common variation in CAV2 using a tag SNP approach. The CAV2 SNP rs13221869 was nominally significant for kidney transplant failure. Validation was sought in an independent group of kidney transplant donors and recipients from Dublin, Ireland using a second genotyping technology. Due to the unexpected absence of rs13221869 from this cohort, the CAV2 gene was resequenced. One novel SNP and a novel insertion/deletion in CAV2 were identified; rs13221869 is located in a repetitive region and was not a true variant in resequenced populations. CAV2 is a plausible candidate gene for association with kidney transplant outcomes given its proximity to CAV1 and its role in attenuating fibrosis. This study does not support an association between CAV2 variation and kidney transplant survival. Further analysis of CAV2 should be undertaken with an awareness of the sequence complexities and genetic variants highlighted by this study. 

Sequence capture and next-generation resequencing of multiple tagged nucleic acid samples for mutation screening of urea cycle disorders

Molecular genetic testing is commonly used to confirm clinical diagnoses of inherited urea cycle disorders (UCDs); however, conventional mutation screenings encompassing only the coding regions of genes may not detect disease-causing mutations occurring in regulatory elements and introns. Microarray-based target enrichment and next-generation sequencing now allow more-comprehensive genetic screening. We applied this approach to UCDs and combined it with the use of DNA bar codes for more cost-effective, parallel analyses of multiple samples.

Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease

More than 1,000 susceptibility loci have been identified through genome-wide association studies (GWAS) of common variants; however, the specific genes and full allelic spectrum of causal variants underlying these findings have not yet been defined. Here we used pooled next-generation sequencing to study 56 genes from regions associated with Crohn's disease in 350 cases and 350 controls. Through follow-up genotyping of 70 rare and low-frequency protein-altering variants in nine independent case-control series (16,054 Crohn's disease cases, 12,153 ulcerative colitis cases and 17,575 healthy controls), we identified four additional independent risk factors in NOD2, two additional protective variants in IL23R, a highly significant association with a protective splice variant in CARD9 (P < 1 × 10(-16), odds ratio ≈ 0.29) and additional associations with coding variants in IL18RAP, CUL2, C1orf106, PTPN22 and MUC19. We extend the results of successful GWAS by identifying new, rare and probably functional variants that could aid functional experiments and predictive models.

The Car Resequencing Problem with Pull-Off Tables

The car sequencing problem determines sequences of different car models launched down a mixed-model assembly line. To avoid work overloads of workforce, car sequencing restricts the maximum occurrence of labor-intensive options, e.g., a sunroof, by applying sequencing rules. We consider this problem in a resequencing context, where a given number of buffers (denoted as pull-off tables) is available for rearranging a stirred sequence. The problem is formalized and suited solution procedures are developed. A lower bound and a dominance rule are introduced which both reduce the running time of our graph approach. Finally, a real-world resequencing setting is investigated.

Genetic Heterogeneity in Autism Spectrum Disorders in a Population Isolate

Positional cloning has enabled hypothesis-free, genome-wide scans for genetic factors contributing to disorders or traits. Traditionally linkage analysis has been used to identify regions of interest, followed by meticulous fine mapping and candidate gene screening using association methods and finally sequencing of regions of interest. More recently, genome-wide association analysis has enabled a more direct approach to identify specific genetic variants explaining a part of the variance of the phenotype of interest. Autism spectrum disorders (ASDs) are a group of childhood onset neuropsychiatric disorders with shared core symptoms but varying severity. Although a strong genetic component has been established in ASDs, genetic susceptibility factors have largely eluded characterization. Here, we have utilized modern molecular genetic methods combined with the advantages provided by the special population structure in Finland to identify genetic risk factors for ASDs. The results of this study show that numerous genetic risk factors exist for ASDs even within a population isolate. Stratification based on clinical phenotype resulted in encouraging results, as previously identified linkage to 3p14-p24 was replicated in an independent family set of families with Asperger syndrome, but no other ASDs. Fine-mapping of the previously identified linkage peak for ASDs at 3q25-q27 revealed association between autism and a subunit of the 5-hydroxytryptamine receptor 3C (HTR3C). We also used dense, genome-wide single nucleotide polymorphism (SNP) data to characterize the population structure of Finns. We observed significant population substructure which correlates with the known history of multiple consecutive bottle-necks experienced by the Finnish population. We used this information to ascertain a genetically homogenous subset of autism families to identify possible rare, enriched risk variants using genome-wide SNP data. No rare enriched genetic risk factors were identified in this dataset, although a subset of families could be genealogically linked to form two extended pedigrees. The lack of founder mutations in this isolated population suggests that the majority of genetic risk factors are rare, de novo mutations unique to individual nuclear families. The results of this study are consistent with others in the field. The underlying genetic architecture for this group of disorders appears highly heterogeneous, with common variants accounting for only a subset of genetic risk. The majority of identified risk factors have turned out to be exceedingly rare, and only explain a subset of the genetic risk in the general population in spite of their high penetrance within individual families. The results of this study, together with other results obtained in this field, indicate that family specific linkage, homozygosity mapping and resequencing efforts are needed to identify these rare genetic risk factors.

Analysis of the barley bract suppression gene Trd1

A typical barley (Hordeum vulgare) floret consists of reproductive organs three stamens and a pistil, and non-reproductive organs-lodicules and two floral bracts, abaxial called 'lemma' and adaxial 'palea'. The floret is subtended by two additional bracts called outer or empty glumes. Together these organs form the basic structural unit of the grass inflorescence, a spikelet. There are commonly three spikelets at each rachis (floral stem of the barley spike) node, one central and two lateral spikelets. Rare naturally occurring or induced phenotypic variants that contain a third bract subtending the central spikelets have been described in barley. The gene responsible for this phenotype was called the THIRD OUTER GLUME1 (Trd1). The Trd1 mutants fail to suppress bract growth and as a result produce leaf-like structures that subtend each rachis node in the basal portion of the spike. Also, floral development at the collar is not always suppressed. In rice and maize, recessive mutations in NECK LEAF1 (Nl1) and TASSEL SHEATH1 (Tsh1) genes, respectively, have been shown to be responsible for orthologous phenotypes. Fine mapping of the trd1 phenotype in an F-3 recombinant population enabled us to position on the long arm of chromosome 1H to a 10 cM region. We anchored this to a conserved syntenic region on rice chromosome Os05 and selected a set of candidate genes for validation by resequencing PCR amplicons from a series of independent mutant alleles. This analysis revealed that a GATA transcription factor, recently proposed to be Trd1, contained mutations in 10 out of 14 independent trd1 mutant alleles that would generate non-functional TRD1 proteins. Together with genetic linkage data, we confirm the identity of Trd1 as the GATA transcription factor ortholog of rice Nl1 and maize Tsh1 genes.