Screening for susceptibility genes in hereditary non-polyposis colorectal cancer

Exome-wide assessment of isolated biliary atresia: A report from the National Birth Defects Prevention Study using child-parent trios and a case-control design to identify novel rare variants

The etiology of biliary atresia (BA) is unknown, but recent studies suggest a role for rare protein-altering variants (PAVs). Exome sequencing data from the National Birth Defects Prevention Study on 54 child-parent trios, one child-mother duo, and 1513 parents of children with other birth defects were analyzed. Most (91%) cases were isolated BA. We performed (1) a trio-based analysis to identify rare de novo, homozygous, and compound heterozygous PAVs and (2) a case-control analysis using a sequence kernel-based association test to identify genes enriched with rare PAVs. While we replicated previous findings on PKD1L1, our results do not suggest that recurrent de novo PAVs play important roles in BA susceptibility. In fact, our finding in NOTCH2, a disease gene associated with Alagille syndrome, highlights the difficulty in BA diagnosis. Notably, IFRD2 has been implicated in other gastrointestinal conditions and warrants additional study. Overall, our findings strengthen the hypothesis that the etiology of BA is complex.

The role of intestinal immune cells and matrix metalloproteinases in inflammatory bowel disease

Inflammatory bowel disease (IBD) has become globally intractable. MMPs play a key role in many inflammatory diseases. However, little is known about the role of MMPs in IBD. In this study, IBD expression profiles were screened from public Gene Expression Omnibus datasets. Functional enrichment analysis revealed that IBD-related specific functions were associated with immune pathways. Five MMPS-related disease markers, namely MMP-9, CD160, PTGDS, SLC26A8, and TLR5, were selected by machine learning and the correlation between each marker and immune cells was evaluated. We then induced colitis in C57 mice using sodium dextran sulfate and validated model construction through HE staining of the mouse colon. WB and immunofluorescence experiments confirmed that the expression levels of MMP-9, PTGDS, SLC26A8, and CD160 in colitis were significantly increased, whereas that of TLR5 were decreased. Flow cytometry analysis revealed that MMPs regulate intestinal inflammation and immunity mainly through CD8 in colitis. Our findings reveal that MMPs play a crucial role in the pathogenesis of IBD and are related to the infiltration of immune cells, suggesting that MMPs may promote the development of IBD by activating immune infiltration and the immune response. This study provides insights for further studies on the occurrence and development of IBD.

Identification of Gene Coexpression Modules and Prognostic Genes Associated with Papillary Thyroid Cancer

Thyroid cancer is a great part of the endocrine tumor with an increasing incidence. Papillary thyroid carcinoma (PTC) is the most common subtype. With the enormous pace taken in the microarray technology, bioinformatics is applied in data mining more frequently. Weighted gene coexpression network analysis (WGCNA) can perform analysis combining clinic information. We performed WGCNA for prognostic genes associated with PTC. From the GEO profile, we got ten modules. We identified a key module that was closest to patients’ survival time. Then, we screened five hub genes (ATRX, BOD1L1, CEP290, DCAF16, and NEK1) from the key module based on the clinical information from TCGA. These five genes not only significantly differ between the normal and tumor groups but have prognostic value. The receiver operating characteristic (ROC) curve indicated that they had the potential to serve as prognostic genes. We performed next-generation sequencing using the PTC tissue to get more convincing evidence. Besides, we established a new signature and verified it through K-M plots and ROC. The signature could be an independent factor for the prognosis of PTC, and we built a nomogram to carry out a quantitative study. In a word, the hub genes we explored in the study deserved more basic and clinical research.

Whole exome sequencing identifies novel germline variants of SLC15A4 gene as potentially cancer predisposing in familial colorectal cancer

About 15% of colorectal cancer (CRC) patients have first-degree relatives affected by the same malignancy. However, for most families the cause of familial aggregation of CRC is unknown. To identify novel high-to-moderate-penetrance germline variants underlying CRC susceptibility, we performed whole exome sequencing (WES) on four CRC cases and two unaffected members of a Polish family without any mutation in known CRC predisposition genes. After WES, we used our in-house developed Familial Cancer Variant Prioritization Pipeline and identified two novel variants in the solute carrier family 15 member 4 (SLC15A4) gene. The heterozygous missense variant, p. Y444C, was predicted to affect the phylogenetically conserved PTR2/POT domain and to have a deleterious effect on the function of the encoded peptide/histidine transporter. The other variant was located in the upstream region of the same gene (GRCh37.p13, 12_129308531_C_T; 43 bp upstream of transcription start site, ENST00000266771.5) and it was annotated to affect the promoter region of SLC15A4 as well as binding sites of 17 different transcription factors. Our findings of two distinct variants in the same gene may indicate a synergistic up-regulation of SLC15A4 as the underlying genetic cause and implicate this gene for the first time in genetic inheritance of familial CRC.

RNA-seq profiling reveals PBMC RNA as a potential biomarker for hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors and has extremely high morbidity and mortality. Although many existing studies have focused on the identification of biomarkers, little information has been uncovered regarding the PBMC RNA profile of HCC. We attempted to create a profile throughout using expression of peripheral blood mononuclear cell (PBMC) RNA using RNA-seq technology and compared the transcriptome between HCC patients and healthy controls. Seventeen patients and 17 matched healthy controls were included in this study, and PBMC RNA was sequenced from all samples. Sequencing data were analyzed using bioinformatics tools, and quantitative reverse transcription PCR (qRT-PCR) was used for selected validation of DEGs. A total of 1,578 dysregulated genes were found in the PBMC samples, including 1,334 upregulated genes and 244 downregulated genes. GO enrichment and KEGG studies revealed that HCC is closely linked to differentially expressed genes (DEGs) implicated in the immune response. Expression of 6 selected genes (SELENBP1, SLC4A1, SLC26A8, HSPA8P4, CALM1, and RPL7p24) was confirmed by qRT-PCR, and higher sensitivity and specificity were obtained by ROC analysis of the 6 genes. CALM1 was found to gradually decrease as tumors enlarged. Nearly the opposite expression modes were obtained when compared to tumor sequencing data. Immune cell populations exhibited significant differences between HCC and controls. These findings suggest a potential biomarker for the diagnosis of HCC. This study provides new perspectives for liver cancer development and possible future successful clinical diagnosis.

MutS protein-based fiber optic particle plasmon resonance biosensor for detecting single nucleotide polymorphisms

A new biosensing method is presented to detect gene mutation by integrating the MutS protein from bacteria with a fiber optic particle plasmon resonance (FOPPR) sensing system. In this method, the MutS protein is conjugated with gold nanoparticles (AuNPs) deposited on an optical fiber core surface. The target double-stranded DNA containing an A and C mismatched base pair in a sample can be captured by the MutS protein, causing increased absorption of green light launching into the fiber and hence a decrease in transmitted light intensity through the fiber. As the signal change is enhanced through consecutive total internal reflections along the fiber, the limit of detection for an AC mismatch heteroduplex DNA can be as low as 0.49 nM. Because a microfluidic chip is used to contain the optical fiber, the narrow channel width allows an analysis time as short as 15 min. Furthermore, the label-free and real-time nature of the FOPPR sensing system enables determination of binding affinity and kinetics between MutS and single-base mismatched DNA. The method has been validated using a heterozygous PCR sample from a patient to determine the allelic fraction. The obtained allelic fraction of 0.474 reasonably agrees with the expected allelic fraction of 0.5. Therefore, the MutS-functionalized FOPPR sensor may potentially provide a convenient quantitative tool to detect single nucleotide polymorphisms in biological samples with a short analysis time at the point-of-care sites.

Candidate Gene Discovery in Hereditary Colorectal Cancer and Polyposis Syndromes-Considerations for Future Studies

To discover novel high-penetrant risk loci for hereditary colorectal cancer (hCRC) and polyposis syndromes many whole-exome and whole-genome sequencing (WES/WGS) studies have been performed. Remarkably, these studies resulted in only a few novel high-penetrant risk genes. Given this observation, the possibility and strategy to identify high-penetrant risk genes for hCRC and polyposis needs reconsideration. Therefore, we reviewed the study design of WES/WGS-based hCRC and polyposis gene discovery studies (n = 37) and provide recommendations to optimize discovery and validation strategies. The group of genetically unresolved patients is phenotypically heterogeneous, and likely composed of distinct molecular subtypes. This knowledge advocates for the screening of a homogeneous, stringently preselected discovery cohort and obtaining multi-level evidence for variant pathogenicity. This evidence can be collected by characterizing the molecular landscape of tumors from individuals with the same affected gene or by functional validation in cell-based models. Together, the combined approach of a phenotype-driven, tumor-based candidate gene search might elucidate the potential contribution of novel genetic predispositions in genetically unresolved hCRC and polyposis.

Expression of Viral microRNAs in Serum and White Blood Cells of Cows Exposed to Bovine Leukemia Virus

Bovine leukemia virus (BLV) affects the health and productivity of cattle. The virus causes abnormal immune function and immunosuppression. MicroRNAs (miRNAs) are involved in gene expression, having been associated with stress and immune response, tumor growth, and viral infection. The objective of this study was to determine the expression of circulating miRNAs produced by BLV in animals exposed to the virus. Sera from 14 animals were collected to establish IgG reactivity to BLV by ELISA, where seven animals were seropositive and seven were seronegative for BLV exposure. White blood cells (WBC) from each animal were also collected and miRNAs were identified by sequencing from sera and WBC. The seropositive group had higher counts of BLV miRNAs when compared to seronegative group in sera and WBC. Blv-miR-1-3p, blv-miR-B2-5p, blv-miR-B4-3p, and blv-miR-B5-5p were statistically significant (P < 0.00001) in serum with an average of 7 log2 fold difference between seropositive and seronegative groups. Blv-miR-B1-3p, blv-miR-B1-5p, blv-miR-B3, blv-miR-B4-3p, blv-miR-B4-5p, blv-miR-B5-5p were statistically significant (P < 1.08e⁻⁹) in WBC with an average of 7 log2 fold difference between the seropositive and the seronegative groups. Blv-miR-B2-3p and blv-miR-B2-5p were also statistically significant in WBC (P < 2.79e-17), with an average of 27 log2 fold difference between the seropositive and the seronegative groups. There were 18 genes identified as being potential targets for blv-miR-B1-5p, and 3 genes for blv-miR-B4-5p. Gene ontology analysis indicated that the target genes are mainly involved in the response to stress and in the immune system process. Several of the identified genes have been associated with leukemia development in humans and cattle. Differential expression of genes targeted by BLV miRNAs should be evaluated to determine their effect in BLV replication.

A Systematic Literature Review of Whole Exome and Genome Sequencing Population Studies of Genetic Susceptibility to Cancer

The application of next-generation sequencing (NGS) technologies in cancer research has accelerated the discovery of somatic mutations; however, progress in the identification of germline variation associated with cancer risk is less clear. We conducted a systematic literature review of cancer genetic susceptibility studies that used NGS technologies at an exome/genome-wide scale to obtain a fuller understanding of the research landscape to date and to inform future studies. The variability across studies on methodologies and reporting was considerable. Most studies sequenced few high-risk (mainly European) families, used a candidate analysis approach, and identified potential cancer-related germline variants or genes in a small fraction of the sequenced cancer cases. This review highlights the importance of establishing consensus on standards for the application and reporting of variants filtering strategies. It also describes the progress in the identification of cancer-related germline variation to date. These findings point to the untapped potential in conducting studies with appropriately sized and racially diverse families and populations, combining results across studies and expanding beyond a candidate analysis approach to advance the discovery of genetic variation that accounts for the unexplained cancer heritability.

Case report: recurrent pituitary adenoma has increased load of somatic variants

BACKGROUND

Pituitary adenomas (PA) have an increased potential for relapse in one to 5 years after resection. In this study, we investigated the genetic differences in genomic DNA of primary and rapidly recurrent tumours in the same patient to explain the causality mechanisms of PA recurrence.

CASE PRESENTATION

The patient was a 69-year-old female with non-functional pituitary macroadenoma with extension into the left cavernous sinus (Knosp grade 2) who underwent craniotomy and partial resection in August 2010. Two years later, the patient had prolonged tumour growth with an essential suprasellar extension (Knosp grade 2), and a second craniotomy with partial tumour resection was performed in September 2012. In both tumours, the KI-67 level was below 1.5%. Exome sequencing via semiconductor sequencing of patient germline DNA and somatic DNA from both tumours was performed. Tmap alignment and Platypus variant calling were performed followed by variant filtering and manual review with IGV software. We observed an increased load of missense variants in the recurrent PA tumour when compared to the original tumour. The number of detected variants increased from ten to 26 and potential clonal expansion of four variants was observed. Additionally, targeted SNP analysis revealed five rare missense SNPs with a potential impact on the function of the encoded proteins.

CONCLUSIONS

In this case study, an SNP located in HRAS is the most likely candidate inducing rapid PA progression. The relapsed PA tumour had a higher variation load and fast tumour recurrence in this patient could be caused by clonal expansion of the leftover tumour tissue.

OMICfpp: a fuzzy approach for paired RNA-Seq counts

BACKGROUND

RNA sequencing is a widely used technology for differential expression analysis. However, the RNA-Seq do not provide accurate absolute measurements and the results can be different for each pipeline used. The major problem in statistical analysis of RNA-Seq and in the omics data in general, is the small sample size with respect to the large number of variables. In addition, experimental design must be taken into account and few tools consider it.

RESULTS

We propose OMICfpp, a method for the statistical analysis of RNA-Seq paired design data. First, we obtain a p-value for each case-control pair using a binomial test. These p-values are aggregated using an ordered weighted average (OWA) with a given orness previously chosen. The aggregated p-value from the original data is compared with the aggregated p-value obtained using the same method applied to random pairs. These new pairs are generated using between-pairs and complete randomization distributions. This randomization p-value is used as a raw p-value to test the differential expression of each gene. The OMICfpp method is evaluated using public data sets of 68 sample pairs from patients with colorectal cancer. We validate our results through bibliographic search of the reported genes and using simulated data set. Furthermore, we compared our results with those obtained by the methods edgeR and DESeq2 for paired samples. Finally, we propose new target genes to validate these as gene expression signatures in colorectal cancer. OMICfpp is available at http://www.uv.es/ayala/software/OMICfpp_0.2.tar.gz .

CONCLUSIONS

Our study shows that OMICfpp is an accurate method for differential expression analysis in RNA-Seq data with paired design. In addition, we propose the use of randomized p-values pattern graphic as a powerful and robust method to select the target genes for experimental validation.