A Streamlined Method for Detecting Structural Variants in Cancer Genomes by Short Read Paired-End Sequencing

Defining the architecture of a specific cancer genome, including its structural variants, is essential for understanding tumor biology, mechanisms of oncogenesis, and for designing effective personalized therapies. Short read paired-end sequencing is currently the most sensitive method for detecting somatic mutations that arise during tumor development. However, mapping structural variants using this method leads to a large number of false positive calls, mostly due to the repetitive nature of the genome and the difficulty of assigning correct mapping positions to short reads. This study describes a method to efficiently identify large tumor-specific deletions, inversions, duplications and translocations from low coverage data using SVDetect or BreakDancer software and a set of novel filtering procedures designed to reduce false positive calls. Applying our method to a spontaneous T cell lymphoma arising in a core RAG2/p53-deficient mouse, we identified 40 validated tumor-specific structural rearrangements supported by as few as 2 independent read pairs.     

A Unified Method for Detecting Secondary Trait Associations with Rare Variants: Application to Sequence Data

Next-generation sequencing has made possible the detection of rare variant (RV) associations with quantitative traits (QT). Due to high sequencing cost, many studies can only sequence a modest number of selected samples with extreme QT. Therefore association testing in individual studies can be underpowered. Besides the primary trait, many clinically important secondary traits are often measured. It is highly beneficial if multiple studies can be jointly analyzed for detecting associations with commonly measured traits. However, analyzing secondary traits in selected samples can be biased if sample ascertainment is not properly modeled. Some methods exist for analyzing secondary traits in selected samples, where some burden tests can be implemented. However p-values can only be evaluated analytically via asymptotic approximations, which may not be accurate. Additionally, potentially more powerful sequence kernel association tests, variable selection-based methods, and burden tests that require permutations cannot be incorporated. To overcome these limitations, we developed a unified method for analyzing secondary trait associations with RVs (STAR) in selected samples, incorporating all RV tests. Statistical significance can be evaluated either through permutations or analytically. STAR makes it possible to apply more powerful RV tests to analyze secondary trait associations. It also enables jointly analyzing multiple cohorts ascertained under different study designs, which greatly boosts power. The performance of STAR and commonly used RV association tests were comprehensively evaluated using simulation studies. STAR was also implemented to analyze a dataset from the SardiNIA project where samples with extreme low-density lipoprotein levels were sequenced. A significant association between LDLR and systolic blood pressure was identified, which is supported by pharmacogenetic studies. In summary, for sequencing studies, STAR is an important tool for detecting secondary-trait RV associations.     

Application of Combined Long Amplicon Sequencing (CoLAS) for Genetic Analysis of Neurofibromatosis Type 1: A Pilot Study

Elaborate analyses of the status of gene mutations in neurofibromatosis type 1 (NF1) are still difficult nowadays due to the large gene sizes, broad mutation spectrum, and the various effects of mutations on mRNA splicing. These problems cannot be solved simply by sequencing the entire coding region using next-generation sequencing (NGS). We recently developed a new strategy, named combined long amplicon sequencing (CoLAS), which is a method for simultaneously analysing the whole genomic DNA region and, also, the full-length cDNA of the disease-causative gene with long-range PCR-based NGS. In this study, CoLAS was specifically arranged for NF1 genetic analysis, then applied to 20 patients (five previously reported and 15 newly recruited patients, including suspicious cases) for optimising the method and to verify its efficacy and benefits. Among new cases, CoLAS detected not only 10 mutations, including three unreported mutations and one mosaic mutation, but also various splicing abnormalities and allelic expression ratios quantitatively. In addition, heterozygous mapping by polymorphisms, including introns, showed copy number monitoring of the entire NF1 gene region was possible in the majority of patients tested. Moreover, it was shown that, when a chromosomal level microdeletion was suspected from heterozygous mapping, it could be detected directly by breakpoint-specific long PCR. In conclusion, CoLAS not simply detect the causative mutation but accurately elucidated the entire structure of the NF1 gene, its mRNA expression, and also the splicing status, which reinforces its high usefulness in the gene analysis of NF1.     

Systematic Analysis of Small RNAs Associated with Human Mitochondria by Deep Sequencing: Detailed Analysis of Mitochondrial Associated miRNA

Mitochondria are one of the central regulators of many cellular processes beyond its well established role in energy metabolism. The inter-organellar crosstalk is critical for the optimal function of mitochondria. Many nuclear encoded proteins and RNA are imported to mitochondria. The translocation of small RNA (sRNA) including miRNA to mitochondria and other sub-cellular organelle is still not clear. We characterized here sRNA including miRNA associated with human mitochondria by cellular fractionation and deep sequencing approach. Mitochondria were purified from HEK293 and HeLa cells for RNA isolation. The sRNA library was generated and sequenced using Illumina system. The analysis showed the presence of unique population of sRNA associated with mitochondria including miRNA. Putative novel miRNAs were characterized from unannotated sRNA sequences. The study showed the association of 428 known, 196 putative novel miRNAs to mitochondria of HEK293 and 327 known, 13 putative novel miRNAs to mitochondria of HeLa cells. The alignment of sRNA to mitochondrial genome was also studied. The targets were analyzed using DAVID to classify them in unique networks using GO and KEGG tools. Analysis of identified targets showed that miRNA associated with mitochondria regulates critical cellular processes like RNA turnover, apoptosis, cell cycle and nucleotide metabolism. The six miRNAs (counts >1000) associated with mitochondria of both HEK293 and HeLa were validated by RT-qPCR. To our knowledge, this is the first systematic study demonstrating the associations of sRNA including miRNA with mitochondria that may regulate site-specific turnover of target mRNA important for mitochondrial related functions.     

MosaicBase:A Knowledgebase of Postzygotic Mosaic Variants in Noncancer Disease-related and Healthy Human Individuals

Mosaic variants resulting from postzygotic mutations are prevalent in the human genome and play important roles in human diseases. However, except for cancer-related variants, there is no collection of postzygotic mosaic variants in noncancer disease-related and healthy individuals. Here, we present MosaicBase, a comprehensive database that includes 6698 mosaic variants related to 266 noncancer diseases and 27,991 mosaic variants identified in 422 healthy individuals. Genomic and phenotypic information of each variant was manually extracted and curated from 383 publications. MosaicBase supports the query of variants with Online Mendelian Inheritance in Man (OMIM) entries, genomic coordinates, gene symbols, or Entrez IDs. We also provide an integrated genome browser for users to easily access mosaic variants and their related annotations for any genomic region. By analyzing the variants collected in MosaicBase, we find that mosaic variants that directly contribute to disease phenotype show features distinct from those of variants in individuals with mild or no phenotypes, in terms of their genomic distribution, mutation signatures, and fraction of mutant cells. MosaicBase will not only assist clinicians in genetic counseling and diagnosis but also provide a useful resource to understand the genomic baseline of postzygotic mutations in the general human population. MosaicBase is publicly available at http://mosaicbase.com/ or http://49.4.21.8:8000.     

Population Specific Impact of Genetic Variants in KCNJ11 Gene to Type 2 Diabetes: A Case-Control and Meta-Analysis Study

Background and Objectives

Potassium inwardly rectifying channel, subfamily J, member 11 (KCNJ11) gene have a key role in insulin secretion and is of substantial interest as a candidate gene for type 2 diabetes (T2D). The current work was performed to delineate the genetic influence of KCNJ11 polymorphisms on risk of T2D in South Indian population through case-control association study along with systematic review and meta-analysis.

Methods

A case-control study of 400 T2D cases and controls of South Indian origin were performed to analyze the association of KCNJ11 polymorphisms (rs5219, rs5215, rs41282930, rs1800467) and copy number variations (CNV) on the risk of T2D. In addition a systematic review and meta-analysis for KCNJ11 rs5219 was conducted in 3,831 cases and 3,543 controls from 5 published reports from South-Asian population by searching various databases. Odds ratio with 95% confidence interval (CI) was used to assess the association strength. Cochran''s Q, I² statistics were used to study heterogeneity between the eligible studies.

Results

KCNJ11 rs5215, C-G-C-C haplotype and two loci analysis (rs5219 vs rs1800467) showed a significant association with T2D but CNV analysis did not show significant variation between T2D cases and control subjects. Lower age of disease onset (P = 0.04) and higher body mass index (BMI) (P = 0.04) were associated with rs5219 TT genotype in T2D patients. The meta-analysis of KCNJ11 rs5219 on South Asian population showed no association on susceptibility to T2D with an overall pooled OR = 0.98, 95% CI = 0.83–1.16. Stratification analysis showed East Asian population and global population were associated with T2D when compared to South Asians.

Conclusion

KCNJ11 rs5219 is not independently associated with T2D in South-Indian population and our meta-analysis suggests that KCNJ11 polymorphism (rs5219) is associated with risk of T2D in East Asian population and global population but this outcome could not be replicated in South Asian sub groups.     

Association of Genetic Variants in the Adiponectin Gene with Metabolic Syndrome: A Case-Control Study and a Systematic Meta-Analysis in the Chinese Population

Background

The prevalence of metabolic syndrome has been rising worldwide, including in China, but knowledge on specific genetic determinants of metabolic syndrome is very limited. A number of studies have reported that polymorphisms in the ADIPOQ gene are associated with metabolic syndrome in Chinese Han populations. However, data is still conflicting. The objective of this study was to examine the associations of the adiponectin genetic variants with metabolic syndrome by a case-control study and meta-analyses in Chinese.

Methods

We first investigated the association of ADIPOQ rs2241766 (+45T>G in exon 2), rs266729 (−11377C>G in promoter) and rs1501299 (+276G>T in intron 2) polymorphisms with metabolic syndrome in a Hubei Han Chinese population with 322 metabolic syndrome patients and 161 normal controls recruited from the Yichang, Hubei. Then we comprehensively reviewed the association between ADIPOQ rs2241766/rs266729/rs1501299 and metabolic syndrome in the Chinese populations via a meta-analysis. The strength of association was assessed by odds ratios (ORs) with 95% confidence intervals (CI).

Results

The G allele frequency of rs2241766 in metabolic syndrome patients was significantly higher than those of controls group (29.8% vs 23.3%, OR = 1.40, P = 0.033). The logistic regression analysis adjusted by gender and age showed a nominally significant association for rs2241766 GG+GT genotype (P = 0.065, OR = 1.55) and rs1501299 GG genotype in recessive model (OR = 1.54, P = 0.066). However, no association was observed for rs266729 in our sample. We identified thirteen studies for rs2241766 (2,684 metabolic syndrome patients and 2,864 controls), three studies for rs266729, and eleven studies for rs1501299 (2,889 metabolic syndrome patients and 3,304 controls) in Chinese. Meta-analysis indicated significant associations for the rs2241766 G allele (OR = 1.14, 95%CI = 1.05–1.24, P = 0.003), rs266729 GG+GT genotypes (OR = 0.80, 95%CI = 0.68–0.92, P = 0.003) and rs1501299 GG+TG genotypes (OR = 1.42, 95%CI 1.16–1.75, P = 0.001).

Conclusions

Our results demonstrated ADIPOQ as a pleiotropic locus for metabolic syndrome and its components in the Han Chinese population.     

VDJ-Seq: Deep Sequencing Analysis of Rearranged Immunoglobulin Heavy Chain Gene to Reveal Clonal Evolution Patterns of B Cell Lymphoma

Understanding tumor clonality is critical to understanding the mechanisms involved in tumorigenesis and disease progression. In addition, understanding the clonal composition changes that occur within a tumor in response to certain micro-environment or treatments may lead to the design of more sophisticated and effective approaches to eradicate tumor cells. However, tracking tumor clonal sub-populations has been challenging due to the lack of distinguishable markers. To address this problem, a VDJ-seq protocol was created to trace the clonal evolution patterns of diffuse large B cell lymphoma (DLBCL) relapse by exploiting VDJ recombination and somatic hypermutation (SHM), two unique features of B cell lymphomas.In this protocol, Next-Generation sequencing (NGS) libraries with indexing potential were constructed from amplified rearranged immunoglobulin heavy chain (IgH) VDJ region from pairs of primary diagnosis and relapse DLBCL samples. On average more than half million VDJ sequences per sample were obtained after sequencing, which contain both VDJ rearrangement and SHM information. In addition, customized bioinformatics pipelines were developed to fully utilize sequence information for the characterization of IgH-VDJ repertoire within these samples. Furthermore, the pipeline allows the reconstruction and comparison of the clonal architecture of individual tumors, which enables the examination of the clonal heterogeneity within the diagnosis tumors and deduction of clonal evolution patterns between diagnosis and relapse tumor pairs. When applying this analysis to several diagnosis-relapse pairs, we uncovered key evidence that multiple distinctive tumor evolutionary patterns could lead to DLBCL relapse. Additionally, this approach can be expanded into other clinical aspects, such as identification of minimal residual disease, monitoring relapse progress and treatment response, and investigation of immune repertoires in non-lymphoma contexts.     

Correction to: In-Silico analyses of Nonsynonymous Variants in the BRCA1 Gene

Biochemical Genetics -     

基于新一代测序技术的mtDNA突变检测方法学的建立

目的:大量研究证实线粒体DNA(mtDNA)突变与肿瘤发生及进展密切相关,但使用传统测序方法难以高通量、高精确度的检测mtDNA突变,为此本研究建立了基于新一代测序技术的mtDNA突变检测方法.方法:提取肝癌患者癌、癌旁组织以及外周血细胞总DNA,利用PCR技术对线粒体基因组进行富集并对PCR产物进行平末端、粘性末端连接或对PCR引物进行氨基修饰,构建mtDNA测序文库.经Illumina HiSeq 2000平台测序后利用生物信息学方法与人类mtDNA参考序列进行比对,并进行测序数据分析.结果:通过对不同质量基因组DNA进行评估后,发现三对引物法适用于大部分DNA样本的mtDNA富集.进一步我们发现PCR引物的氨基修饰可显著提高测序数据覆盖均一性,降低测序成本.结论:本研究利用新一代测序技术通过对线粒体DNA富集方法以及测序覆盖度均一性进行优化,建立了一套灵敏、特异、高通量的mtDNA突变检测策略,为mtDNA突变与疾病研究提供了新方法.     

Assessment of Trust in Physician: A Systematic Review of Measures

Over the last decades, trust in physician has gained in importance. Studies have shown that trust in physician is associated with positive health behaviors in patients. However, the validity of empirical findings fundamentally depends on the quality of the measures in use. Our aim was to provide an overview of trust in physician measures and to evaluate the methodological quality of the psychometric studies and the quality of psychometric properties of identified measures. We conducted an electronic search in three databases (Medline, EMBASE and PsycInfo). The secondary search strategy included reference and citation tracking of included full texts and consultation of experts in the field. Retrieved records were screened independently by two reviewers. Full texts that reported on testing of psychometric properties of trust in physician measures were included in the review. Study characteristics and psychometric properties were extracted. We evaluated the quality of design, methods and reporting of studies with the COnsensus based Standards for the selection of health status Measurement INstruments (COSMIN) checklist. The quality of psychometric properties was assessed with Terwee’s 2007 quality criteria. After screening 3284 records and assessing 169 full texts for eligibility, fourteen studies on seven trust in physician measures were included. Most of the studies were conducted in the USA and used English measures. All but one measure were generic. Sample sizes range from 25 to 1199 participants, recruited in very heterogeneous settings. Quality assessments revealed several flaws in the methodological quality of studies. COSMIN scores were mainly fair or poor. The overall quality of measures’ psychometric properties was intermediate. Several trust in physician measures have been developed over the last years, but further psychometric evaluation of these measures is strongly recommended. The methodological quality of psychometric property studies could be improved by adhering to quality criteria like the COSMIN checklist.     

Monitoring Pharmacologically Induced Immunosuppression by Immune Repertoire Sequencing to Detect Acute Allograft Rejection in Heart Transplant Patients: A Proof-of-Concept Diagnostic Accuracy Study

Background

It remains difficult to predict and to measure the efficacy of pharmacological immunosuppression. We hypothesized that measuring the B-cell repertoire would enable assessment of the overall level of immunosuppression after heart transplantation.

Methods and Findings

In this proof-of-concept study, we implemented a molecular-barcode-based immune repertoire sequencing assay that sensitively and accurately measures the isotype and clonal composition of the circulating B cell repertoire. We used this assay to measure the temporal response of the B cell repertoire to immunosuppression after heart transplantation. We selected a subset of 12 participants from a larger prospective cohort study (ClinicalTrials.gov NCT01985412) that is ongoing at Stanford Medical Center and for which enrollment started in March 2010. This subset of 12 participants was selected to represent post-heart-transplant events, with and without acute rejection (six participants with moderate-to-severe rejection and six without). We analyzed 130 samples from these patients, with an average follow-up period of 15 mo. Immune repertoire sequencing enables the measurement of a patient’s net state of immunosuppression (correlation with tacrolimus level, r = −0.867, 95% CI −0.968 to −0.523, p = 0.0014), as well as the diagnosis of acute allograft rejection, which is preceded by increased immune activity with a sensitivity of 71.4% (95% CI 30.3% to 94.9%) and a specificity of 82.0% (95% CI 72.1% to 89.1%) (cell-free donor-derived DNA as noninvasive gold standard). To illustrate the potential of immune repertoire sequencing to monitor atypical post-transplant trajectories, we analyzed two more patients, one with chronic infections and one with amyloidosis. A larger, prospective study will be needed to validate the power of immune repertoire sequencing to predict rejection events, as this proof-of-concept study is limited to a small number of patients who were selected based on several criteria including the availability of a large number of samples and the absence or presence of rejection events.

Conclusions

If confirmed in larger, prospective studies, the method described here has potential applications in the tailored management of post-transplant immunosuppression and, more broadly, as a method for assessing the overall activity of the immune system.     

Gene Control in Somatic Embryos of Digitalis lanata: Expression of the β-Glucuronidase Gene Fused to a Plastocyanin Promoter

Digitalis lanata was transformed by agrobacteria-mediated gene transfer with a chimeric reporter gene encoding for β-glucuronidase (CUS) from Escherichia coll under the control of the plastocyanin 3 (Pc3) promoter from Spinada oleracea (Pc3::uidA fusion gene). Transformed cell lines were regenerated to plants via somatic embryos. CUS activity was determined fluorometrically and histochemically. The Pc3::uidA fusion gene was expressed in the late globular and bipolar stages of somatic embryos. Expression started in globular embryos (stage-1-globules) in that part of the parenchymatic tissue which later on formed the cotyledons. No GUS activity was detectable in the parenchymatic tissue forming the root pole, in cells of the developing procambium or in epidermal cells. These tissues were free of GUS activity also in bipolar embryos. The parenchymatic cells of the cotyledons and the primary cortex of the hypocotyl of germinating embryos showed GUS activity, in contrast to the epidermal cells and the cells of the central cylinder.     

With Reference to Reference Genes: A Systematic Review of Endogenous Controls in Gene Expression Studies

The choice of reference genes that are stably expressed amongst treatment groups is a crucial step in real-time quantitative PCR gene expression studies. Recent guidelines have specified that a minimum of two validated reference genes should be used for normalisation. However, a quantitative review of the literature showed that the average number of reference genes used across all studies was 1.2. Thus, the vast majority of studies continue to use a single gene, with β-actin (ACTB) and/or glyceraldehyde 3-phosphate dehydrogenase (GAPDH) being commonly selected in studies of vertebrate gene expression. Few studies (15%) tested a panel of potential reference genes for stability of expression before using them to normalise data. Amongst studies specifically testing reference gene stability, few found ACTB or GAPDH to be optimal, whereby these genes were significantly less likely to be chosen when larger panels of potential reference genes were screened. Fewer reference genes were tested for stability in non-model organisms, presumably owing to a dearth of available primers in less well characterised species. Furthermore, the experimental conditions under which real-time quantitative PCR analyses were conducted had a large influence on the choice of reference genes, whereby different studies of rat brain tissue showed different reference genes to be the most stable. These results highlight the importance of validating the choice of normalising reference genes before conducting gene expression studies.     

A 3.7 Mb Deletion Encompassing ZEB2 Causes a Novel Polled and Multisystemic Syndrome in the Progeny of a Somatic Mosaic Bull

Polled and Multisystemic Syndrome (PMS) is a novel developmental disorder occurring in the progeny of a single bull. Its clinical spectrum includes polledness (complete agenesis of horns), facial dysmorphism, growth delay, chronic diarrhea, premature ovarian failure, and variable neurological and cardiac anomalies. PMS is also characterized by a deviation of the sex-ratio, suggesting male lethality during pregnancy. Using Mendelian error mapping and whole-genome sequencing, we identified a 3.7 Mb deletion on the paternal bovine chromosome 2 encompassing ARHGAP15, GTDC1 and ZEB2 genes. We then produced control and affected 90-day old fetuses to characterize this syndrome by histological and expression analyses. Compared to wild type individuals, affected animals showed a decreased expression of the three deleted genes. Based on a comparison with human Mowat-Wilson syndrome, we suggest that deletion of ZEB2, is responsible for most of the effects of the mutation. Finally sperm-FISH, embryo genotyping and analysis of reproduction records confirmed somatic mosaicism in the founder bull and male-specific lethality during the first third of gestation. In conclusion, we identified a novel locus involved in bovid horn ontogenesis and suggest that epithelial-to-mesenchymal transition plays a critical role in horn bud differentiation. We also provide new insights into the pathogenicity of ZEB2 loss of heterozygosity in bovine and humans and describe the first case of male-specific lethality associated with an autosomal locus in a non-murine mammalian species. This result sets PMS as a unique model to study sex-specific gene expression/regulation.     

RNA Deep Sequencing Analysis Reveals That Nicotine Restores Impaired Gene Expression by Viral Proteins in the Brains of HIV-1 Transgenic Rats

Persons infected with HIV-1 often develop neurologic disorders despite receiving highly active anti-retroviral therapy. Although the underlying mechanism is largely undetermined, our previous RNA-seq-based study showed that the expression of many genes was altered in the central nervous system (CNS) of HIV-1 transgenic (HIV-1Tg) rats. Because nicotine, a natural agonist of nicotinic acetylcholine receptors, exhibits a neuroprotective effect, we presently tested the hypothesis that nicotine restores the expression of altered genes in the CNS of HIV-1Tg rats. Adult male HIV-1Tg and F344 control strain rats were injected with either nicotine (0.25 mg/kg) or saline subcutaneously twice a day for 17 days. Gene expression in the prefrontal cortex (PFC), dorsal hippocampus (HIP), and dorsal striatum (STR) was evaluated using the RNA deep sequencing technique. We found that about 20% of the altered genes in the HIV-1Tg rat were affected by nicotine in each brain region, with the expression of most restored. Analysis of the restored genes showed distinct pathways corrected by nicotine in different brain regions of HIV-1Tg rats. Specifically, the two most significantly restored pathways were Wnt/β-catenin signaling and ephrin B signaling in the PFC, cAMP-responsive element-binding protein (CREB) signaling and glutathione metabolism pathway in the HIP, and tricarboxylic acid (TCA) cycle and calcium signaling in the STR. Together, our findings indicate that cholinergic modulators such as nicotine have beneficial effects on HIV-1-induced neurologic deficits.