Accelerating Gene Discovery by Phenotyping Whole-Genome Sequenced Multi-mutation Strains and Using the Sequence Kernel Association Test (SKAT)

The Caenorhabditis elegans Shugoshin regulates TAC-1 in cilia

The conserved Shugoshin (SGO) protein family is essential for mediating proper chromosome segregation from yeast to humans but has also been implicated in diverse roles outside of the nucleus. SGO's roles include inhibiting incorrect spindle attachment in the kinetochore, regulating the spindle assembly checkpoint (SAC), and ensuring centriole cohesion in the centrosome, all functions that involve different microtubule scaffolding structures in the cell. In Caenorhabditis elegans, a species with holocentric chromosomes, SGO-1 is not required for cohesin protection or spindle attachment but appears important for licensing meiotic recombination. Here we provide the first functional evidence that in C. elegans, Shugoshin functions in another extranuclear, microtubule-based structure, the primary cilium. We identify the centrosomal and microtubule-regulating transforming acidic coiled-coil protein, TACC/TAC-1, which also localizes to the basal body, as an SGO-1 binding protein. Genetic analyses indicate that TAC-1 activity must be maintained below a threshold at the ciliary base for correct cilia function, and that SGO-1 likely participates in constraining TAC-1 to the basal body by influencing the function of the transition zone 'ciliary gate'. This research expands our understanding of cellular functions of Shugoshin proteins and contributes to the growing examples of overlap between kinetochore, centrosome and cilia proteomes.

Cell-type-specific promoters for C. elegans glia

Glia shape the development and function of the C. elegans nervous system, especially its sense organs and central neuropil (nerve ring). Cell-type-specific promoters allow investigators to label or manipulate individual glial cell types, and therefore provide a key tool for deciphering glial function. In this technical resource, we compare the specificity, brightness, and consistency of cell-type-specific promoters for C. elegans glia. We identify a set of promoters for the study of seven glial cell types (F16F9.3, amphid and phasmid sheath glia; F11C7.2, amphid sheath glia only; grl-2, amphid and phasmid socket glia; hlh-17, cephalic (CEP) sheath glia; and grl-18, inner labial (IL) socket glia) as well as a pan-glial promoter (mir-228). We compare these promoters to promoters that are expressed more variably in combinations of glial cell types (delm-1 and itx-1). We note that the expression of some promoters depends on external conditions or the internal state of the organism, such as developmental stage, suggesting glial plasticity. Finally, we demonstrate an approach for prospectively identifying cell-type-specific glial promoters using existing single-cell sequencing data, and we use this approach to identify two novel promoters specific to IL socket glia (col-53 and col-177).

Association Study of Genetic Variants in Autophagy Pathway and Risk of Non-syndromic Cleft Lip With or Without Cleft Palate

Although genetic variants in autophagy pathway genes were associated with the risk of oral cancers and early development in embryos, their associations with non-syndromic cleft lip with or without cleft palate (NSCL/P) risk remained unclear. A two-stage case-control study (2,027 NSCL/P cases and 1,843 controls) was performed to investigate the associations between single nucleotide polymorphisms (SNPs) in 23 autophagy pathway genes and NSCL/P susceptibility. The logistic regression model was used to calculate effects of SNPs on NSCL/P susceptibility. Gene-based analysis was performed via the sequence kernel association test (SKAT) and multi-marker analysis of genomic annotation (MAGMA) methods. Expression quantitative trait loci (eQTL) analysis was conducted using NSCL/P lip tissue samples. Gene expression during embryonic development was evaluated using RNA-Seq. Functional roles were explored by luciferase activity assay, cell apoptosis, proliferation, and cycle in vitro. Rs2301104 in HIF1A was significantly associated with NSCL/P susceptibility in the combined analysis (OR: 1.29, 95% CI: 1.09-1.29, P = 3.39 × 10^-03), and showed strong evidence of association heterogeneity (P = 9.06 × 10^-03) with obvious association in the female (OR: 1.80; 95% CI: 1.32-2.45; P = 1.79 × 10^-04). The G allele of rs2301104 was associated with enhanced transcription activity and high expression of HIF1A compared with that of C allele. Moreover, rs2301104 exhibited an eQTL effect for HIF1A with its GC/CC genotypes associated with decreased HIF1A expression compared with those with GG genotypes (P = 3.1 × 10^-2). Knockdown of HIF1A induced cell apoptosis and inhibited cell proliferation in human embryonic palate mesenchyme (HEPM) and human oral epithelium cells (HOEC). This study demonstrated that rs2301104 in autophagy pathway gene HIF1A was associated with susceptibility of NSCL/P.

Systematic analyses of regulatory variants in DNase I hypersensitive sites identified two novel lung cancer susceptibility loci

DNase I hypersensitive sites (DHS) are abundant in regulatory elements, such as promoter, enhancer and transcription factor binding sites. Many studies have revealed that disease-associated variants were concentrated in DHS-related regions. However, limited studies are available on the roles of DHS-related variants in lung cancer. In this study, we performed a large-scale case-control study with 20 871 lung cancer cases and 15 971 controls to evaluate the associations between regulatory genetic variants in DHS and lung cancer susceptibility. The expression quantitative trait loci (eQTL) analysis and pathway-enrichment analysis were performed to identify the possible target genes and pathways. In addition, we performed motif-based analysis to explore the lung-cancer-related motifs using sequence kernel association test. Two novel variants, rs186332 in 20q13.3 (C>T, odds ratio [OR] = 1.17, 95% confidence interval [95% CI]: 1.10-1.24, P = 8.45 × 10-7) and rs4839323 in 1p13.2 (T>C, OR = 0.92, 95% CI: 0.89-0.95, P = 1.02 × 10-6) showed significant association with lung cancer risk. The eQTL analysis suggested that these two SNPs might regulate the expression of MRGBP and SLC16A1, respectively. What's more, the expression of both MRGBP and SLC16A1 was aberrantly elevated in lung tumor tissues. The motif-based analysis identified 10 motifs related to the risk of lung cancer (P < 1.71 × 10-4). Our findings suggested that variants in DHS might modify lung cancer susceptibility through regulating the expression of surrounding genes. This study provided us a deeper insight into the roles of DHS-related genetic variants for lung cancer.

Exploring Impact of Rare Variation in Systemic Lupus Erythematosus by a Genome Wide Imputation Approach

The importance of low frequency and rare variation in complex disease genetics is difficult to estimate in patient populations. Genome-wide association studies are therefore, underpowered to detect rare variation. We have used a combined approach of genome-wide-based imputation with a highly stringent sequence kernel association (SKAT) test and a case-control burden test. We identified 98 candidate genes containing rare variation that in aggregate show association with SLE many of which have recognized immunological function, but also function and expression related to relevant tissues such as the joints, skin, blood or central nervous system. In addition we also find that there is a significant enrichment of genes annotated for disease-causing mutations in the OMIM database, suggesting that in complex diseases such as SLE, such mutations may be involved in subtle or combined phenotypes or could accelerate specific organ abnormalities found in the disease. We here provide an important resource of candidate genes for SLE.

Multi-marker analysis of genomic annotation on gastric cancer GWAS data from Chinese populations

BACKGROUND

Gastric cancer (GC) is one of the high-incidence and high-mortality cancers all over the world. Though genome-wide association studies (GWASs) have found some genetic loci related to GC, they could only explain a small fraction of the potential pathogenesis for GC.

METHODS

We used multi-marker analysis of genomic annotation (MAGMA) to analyze pathways from four public pathway databases based on Chinese GWAS data including 2631 GC cases and 4373 controls. The differential expressions of selected genes in certain pathways were assessed on the basis of The Cancer Genome Atlas database. Immunohistochemistry was also conducted on 55 GC and paired normal tissues of Chinese patients to localize the expression of genes and further validate the differential expression.

RESULTS

We identified three pathways including chemokine signaling pathway, potassium ion import pathway, and interleukin-7 (IL7) pathway, all of which were associated with GC risk. NMI in IL7 pathway and RAC1 in chemokine signaling pathway might be two new candidate genes involved in GC pathogenesis. Additionally, NMI and RAC1 were overexpressed in GC tissues than normal tissues.

CONCLUSION

Immune and inflammatory associated processes and potassium transporting might participate in the development of GC. Besides, NMI and RAC1 might represent two new key genes related to GC. Our findings might give new insight into the biological mechanism and immunotherapy for GC.

Nanopore sequencing: An enrichment-free alternative to mitochondrial DNA sequencing

Mitochondrial DNA sequence data are often utilized in disease studies, conservation genetics and forensic identification. The current approaches for sequencing the full mtGenome typically require several rounds of PCR enrichment during Sanger or MPS protocols followed by fairly tedious assembly and analysis. Here we describe an efficient approach to sequencing directly from genomic DNA samples without prior enrichment or extensive library preparation steps. A comparison is made between libraries sequenced directly from native DNA and the same samples sequenced from libraries generated with nine overlapping mtDNA amplicons on the Oxford Nanopore MinION™ device. The native and amplicon library preparation methods and alternative base calling strategies were assessed to establish error rates and identify trends of discordance between the two library preparation approaches. For the complete mtGenome, 16 569 nucleotides, an overall error rate of approximately 1.00% was observed. As expected with mtDNA, the majority of error was detected in homopolymeric regions. The use of a modified basecaller that corrects for ambiguous signal in homopolymeric stretches reduced the error rate for both library preparation methods to approximately 0.30%. Our study indicates that direct mtDNA sequencing from native DNA on the MinION™ device provides comparable results to those obtained from common mtDNA sequencing methods and is a reliable alternative to approaches using PCR-enriched libraries.

Whole Exome Sequencing of Patients from Multicase Families with Systemic Lupus Erythematosus Identifies Multiple Rare Variants

In an effort to identify rare alleles associated with SLE, we have performed whole exome sequencing of the most distantly related affected individuals from two large Icelandic multicase SLE families followed by Ta targeted genotyping of additional relatives. We identified multiple rare likely pathogenic variants in nineteen genes co-segregating with the disease through multiple generations. Gene co-expression and protein-protein interaction analysis identified a network of highly connected genes comprising several loci previously implicated in autoimmune diseases. These genes were significantly enriched for immune system development, lymphocyte activation, DNA repair, and V(D)J gene recombination GO-categories. Furthermore, we found evidence of aggregate association and enrichment of rare variants at the FAM71E1/EMC10 locus in an independent set of 4,254 European SLE-cases and 4,349 controls. Our study presents evidence supporting that multiple rare likely pathogenic variants, in newly identified genes involved in known disease pathogenic pathways, segregate with SLE at the familial and population level.

Different mutational function of low- and high-linear energy transfer heavy-ion irradiation demonstrated by whole-genome resequencing of Arabidopsis mutants

Heavy-ion irradiation is a powerful mutagen that possesses high linear energy transfer (LET). Several studies have indicated that the value of LET affects DNA lesion formation in several ways, including the efficiency and the density of double-stranded break induction along the particle path. We assumed that the mutation type can be altered by selecting an appropriate LET value. Here, we quantitatively demonstrate differences in the mutation type induced by irradiation with two representative ions, Ar ions (LET: 290 keV μm^-1 ) and C ions (LET: 30.0 keV μm^-1 ), by whole-genome resequencing of the Arabidopsis mutants produced by these irradiations. Ar ions caused chromosomal rearrangements or large deletions (≥100 bp) more frequently than C ions, with 10.2 and 2.3 per mutant genome under Ar- and C-ion irradiation, respectively. Conversely, C ions induced more single-base substitutions and small indels (<100 bp) than Ar ions, with 28.1 and 56.9 per mutant genome under Ar- and C-ion irradiation, respectively. Moreover, the rearrangements induced by Ar-ion irradiation were more complex than those induced by C-ion irradiation, and tended to accompany single base substitutions or small indels located close by. In conjunction with the detection of causative genes through high-throughput sequencing, selective irradiation by beams with different effects will be a powerful tool for forward genetics as well as studies on chromosomal rearrangements.

Systematical analyses of variants in DNase I hypersensitive sites to identify hepatocellular carcinoma susceptibility loci in a Chinese population

BACKGROUND AND AIM

Although several variants located at coding and non-coding regions were evaluated by previous studies, the evidence for associations between variants located in DNase I-hypersensitive sites (DHSs) and hepatocellular carcinoma (HCC) risk was still limited. Recent advances using ENCODE data indicated that genetic variants in DHSs played an important role in carcinogenesis. Therefore, systematically investigating the associations between regulatory variants in DHSs and HCC risk should be put on the agenda.

METHODS

We conducted a case-control design (1538 HCC cases vs 1465 normal controls) to evaluate the effects on HCC for the variants located at the uniform DNase I hypersensitive sites sequencing peaks in a Chinese population.

RESULTS

We found two novel single nucleotide polymorphisms rs12309362 (odds ratio = 0.64, P = 5.61 × 10^-6 ) and rs9970827 (odds ratio = 0.73, P = 7.23 × 10^-6 ) significantly associated with decreased risk of HCC. Conditional analysis proved that both of them independently contributed to the susceptibility of HCC. Expression quantitative trait loci analysis found that A allele of rs12309362 was significantly associated with an elevated expression of phosphatase phosphoglycerate mutase 5 in liver tissues. In addition, gene-based analysis indicated that CEBPB (P = 1 × 10^-4 ) was associated with the risk of HCC, and the expression of CEBPB was significantly lower in 50 The Cancer Genome Atlas HCC tumor tissues compared with matched normal tissues.

CONCLUSIONS

Our results indicated that rs12309362 (G > A), rs9970827 (A > G) in DHSs, and elevated expression of CEBPB were associated with a decreased risk of HCC. These results may contribute us to understand the function of regulatory DNA sequences in HCC development.

Genes and molecular pathways underpinning ciliopathies

Motile and non-motile (primary) cilia are nearly ubiquitous cellular organelles. The dysfunction of cilia causes diseases known as ciliopathies. The number of reported ciliopathies (currently 35) is increasing, as is the number of established (187) and candidate (241) ciliopathy-associated genes. The characterization of ciliopathy-associated proteins and phenotypes has improved our knowledge of ciliary functions. In particular, investigating ciliopathies has helped us to understand the molecular mechanisms by which the cilium-associated basal body functions in early ciliogenesis, as well as how the transition zone functions in ciliary gating, and how intraflagellar transport enables cargo trafficking and signalling. Both basic biological and clinical studies are uncovering novel ciliopathies and the ciliary proteins involved. The assignment of these proteins to different ciliary structures, processes and ciliopathy subclasses (first order and second order) provides insights into how this versatile organelle is built, compartmentalized and functions in diverse ways that are essential for human health.