Single Gene Level Research Articles

Genetics of autosomal recessive intellectual disability.

In the last few years, next-generation sequencing has led to enormous progress in deciphering monogenic forms of intellectual disability. Autosomal dominant intellectual disability (ADID) and X chromosomal intellectual disability (XLID) have been the focus of research. Apart from metabolic disorders, autosomal recessive intellectual disability (ARID) is still behind, probably because it is more heterogeneous and less prevalent in industrial populations. The prevalence of ARID in a cohort of affected children of an outbred population is estimated to be about 10%, with an upward tendency in still unclarified cases. The risk for ARID in children of first cousins or closer is a magnitude higher than for children of unrelated parents. Taken together, it seems that children of related parents are at a 2 to 3 times higher risk for ID. There are no prevalent ARID genes, pathways, or protein complexes and the functions of the affected proteins are very diverse and limited not only to neurological aspects. Thus, in a regular case, there is no reasoning for picking a few genes for a first diagnostic step, and a genetic diagnosis of ID in general, and ARID specifically, is better made using large panels or exome sequencing. In addition, in the last few months, evidence has been growing that many ARID genes are pleiotropic and that the resulting phenotypes may have a broad spectrum. For an exhaustive deciphering of the genetics of ARID, we suggest research at the level of single genes rather than large meta-analyses.

Gene Expression Order Attributed to Genome Reduction and the Steady Cellular State in Escherichia coli.

Transcriptomes not only reflect the growth status but also link to the genome in bacteria. To investigate if and how genome or cellular state changes contribute to the gene expression order, the growth profile-associated transcriptomes of an assortment of genetically differentiated Escherichia coli either exponentially growing under varied conditions or in response to environmental disturbance were analyzed. A total of 168 microarray data sets representing 56 transcriptome variations, were categorized by genome size (full length or reduced) and cellular state (steady or unsteady). At the genome-wide level, the power-law distribution of gene expression was found to be significantly disturbed by the genome size but not the cellular state. At the regulatory network level, more networks with improved coordination of growth rates were observed in genome reduction than at the steady state. At the single-gene level, both genome reduction and steady state increased the correlation of gene expression to growth rate, but the enriched gene categories with improved correlations were different. These findings not only illustrate the order of gene expression attributed to genome reduction and steady cellular state but also indicate that the accessory sequences acquired during genome evolution largely participated in the coordination of transcriptomes to growth fitness.

Genomic Prediction of Complex Phenotypes Using Genic Similarity Based Relatedness Matrix.

In the last years, a series of methods for genomic prediction (GP) have been established, and the advantages of GP over pedigree best linear unbiased prediction (BLUP) have been reported. However, the majority of previously proposed GP models are purely based on mathematical considerations while seldom take the abundant biological knowledge into account. Prediction ability of those models largely depends on the consistency between the statistical assumptions and the underlying genetic architectures of traits of interest. In this study, gene annotation information was incorporated into GP models by constructing haplotypes with SNPs mapped to genic regions. Haplotype allele similarity between pairs of individuals was measured through different approaches at single gene level and then converted into whole genome level, which was then treated as a special kernel and used in kernel based GP models. Results shown that the gene annotation guided methods gave higher or at least comparable predictive ability in some traits, especially in the Arabidopsis dataset and the rice breeding population. Compared to SNP models and haplotype models without gene annotation, the gene annotation based models improved the predictive ability by 0.56~26.67% in the Arabidopsis and 1.62~16.53% in the rice breeding population, respectively. However, incorporating gene annotation slightly improved the predictive ability for several traits but did not show any extra gain for the rest traits in a chicken population. In conclusion, integrating gene annotation into GP models could be beneficial for some traits, species, and populations compared to SNP models and haplotype models without gene annotation. However, more studies are yet to be conducted to implicitly investigate the characteristics of these gene annotation guided models.

Rainbow-Seq: Combining Cell Lineage Tracing with Single-Cell RNA Sequencing in Preimplantation Embryos

SummaryWe developed the Rainbow-seq technology to trace cell division history and reveal single-cell transcriptomes. With distinct fluorescent protein genes as lineage markers, Rainbow-seq enables each single-cell RNA sequencing (RNA-seq) experiment to simultaneously decode the lineage marker genes and read single-cell transcriptomes. We triggered lineage tracking in each blastomere at the 2-cell stage, observed microscopically inequivalent contributions of the progeny to the two embryonic poles at the blastocyst stage, and analyzed every single cell at either 4- or 8-cell stage with deep paired-end sequencing of full-length transcripts. Although lineage difference was not marked unequivocally at a single-gene level, it became clear when the transcriptome was analyzed as a whole. Moreover, several groups of novel transcript isoforms with embedded repeat sequences exhibited lineage difference, suggesting a possible link between DNA demethylation and cell fate decision. Rainbow-seq bridged a critical gap between division history and single-cell RNA-seq assays.

Genomic Characterization of Testicular Germ Cell Tumors Relapsing After Chemotherapy

BackgroundAlthough both seminomatous and nonseminomatous testicular germ cell tumors (TGCTs) have favorable outcomes with chemotherapy, a subset is chemorefractory, and novel therapeutic options are needed. ObjectiveTo molecularly characterize chemotherapy-refractory TGCTs. Design, setting, and participantsArchival tissues from 107 chemotherapy-treated and relapsed TGCT patients (23 seminomas; 84 nonseminomas) underwent hybrid-capture–based genomic profiling to evaluate four classes of genomic alterations (GAs). Tumor mutational burden (TMB) and microsatellite instability (MSI) were also measured. InterventionGenomic profiling on tumor samples from chemotherapy-refractory TGCTs. Outcome measurements and statistical analysisDescriptive analyses and differences between seminoma and nonseminoma subgroups were reported. Results and limitationsThe mean GA/tumor was 2.9 for seminomas and 4.0 for nonseminomas (p=0.04). KRAS alterations (mainly amplifications) were the most common GAs at the single-gene level (47.8% of seminomas and 51.2% of nonseminomas). RAS-RAF pathway (56.5% vs 52.3%) and cell-cycle pathway (52.2% vs 56.0%) were the most common GA classes in seminomas and nonseminomas, respectively. Receptor tyrosine kinase pathway and PI3K pathway GAs were more frequent in seminomas (p=0.02). Median TMB was 1.8 mutations/Mb for seminomas and 2.7 mutations/Mb for nonseminomas (p=0.098), and MSI-high status was found in one nonseminoma only (1.2%). A lack of clinical outcome correlation is a limitation of the present analyses. ConclusionsIn chemotherapy-refractory TGCTs, trials with agents targeting the KRAS pathway may be pursued due to the high frequency of KRAS GAs. Overall, the GAs found in refractory seminomas and nonseminomas differ significantly. Considering the frequency of high TMB or MSI-high status, immunotherapy may benefit a small subset of nonseminomas. Patient summaryTesticular cancers that are resistant to or relapse after standard chemotherapy may harbor genomic alterations that are potentially druggable, particularly in the clinical trial setting, and genomic profiling can guide clinical research and disclose therapeutic opportunities for these patients.

Comprehensive analysis of mutation and expression based pathways in head and neck squamous cell carcinoma.

e18015Background: Over- or under expression of mRNA result from genetic alteration. mRNA expression based pathway are as important as single gene level mutations analyzed by next generation sequenc...

Discovery and Analysis of Invertebrate IgVJ-C2 Structure from Amphioxus Provides Insight into the Evolution of the Ig Superfamily.

The emergence of adaptive immunity in jawed vertebrates depended on the appearance of variable immune receptors, BCRs and TCRs, which exhibit variable-J-constant (VJ-C)-type Ig superfamily folds. Hitherto, however, the structures of IgV-J-IgC-type molecules had never been characterized in invertebrates, leaving the origin of BCR/TCR-type molecules unknown. Using x-ray crystallography, the structure of a VJ-C2 molecule, named AmpIgVJ-C2, was determined in amphioxus (Branchiostoma floridae). The first domain shows typical V folding, including the hydrophobic core, CDR analogs, and eight conserved residues. The second domain is a C2-type Ig superfamily domain, as defined by its short length and the absence of β-strand D- and C1-typical motifs. AmpIgVJ-C2 molecules form homodimers, using "three-layer packing dimerization," as described for TCRs and BCRs. The AmpIgVJ-C2 V domain harbors a diglycine motif in β-strand G and forms a β-bulge structure participating in V-V intermolecular interaction. By immunohistochemistry, AmpIgVJ-C2 molecules were primarily found in mucosal tissues, whereas PCR and sequence analysis indicated considerable genetic variation at the single-gene level; these findings would be consistent with an immune function and a basic ability to adapt to binding different immune targets. Our results show a BCR/TCR-ancestral like molecule in amphioxus and help us to understand the evolution of the adaptive immune system.

Into the Fourth Dimension: Dysregulation of Genome Architecture in Aging and Alzheimer's Disease.

Alzheimer’s disease (AD) is a progressive neurodegenerative disease characterized by synapse dysfunction and cognitive impairment. Understanding the development and progression of AD is challenging, as the disease is highly complex and multifactorial. Both environmental and genetic factors play a role in AD pathogenesis, highlighted by observations of complex DNA modifications at the single gene level, and by new evidence that also implicates changes in genome architecture in AD patients. The four-dimensional structure of chromatin in space and time is essential for context-dependent regulation of gene expression in post-mitotic neurons. Dysregulation of epigenetic processes have been observed in the aging brain and in patients with AD, though there is not yet agreement on the impact of these changes on transcription. New evidence shows that proteins involved in genome organization have altered expression and localization in the AD brain, suggesting that the genomic landscape may play a critical role in the development of AD. This review discusses the role of the chromatin organizers and epigenetic modifiers in post-mitotic cells, the aging brain, and in the development and progression of AD. How these new insights can be used to help determine disease risk and inform treatment strategies will also be discussed.

Large-scale genomic analysis shows association between homoplastic genetic variation in Mycobacterium tuberculosis genes and meningeal or pulmonary tuberculosis

BackgroundMeningitis is the most severe manifestation of tuberculosis. It is largely unknown why some people develop pulmonary TB (PTB) and others TB meningitis (TBM); we examined if the genetic background of infecting M. tuberculosis strains may be relevant.MethodsWe whole-genome sequenced M. tuberculosis strains isolated from 322 HIV-negative tuberculosis patients from Indonesia and compared isolates from patients with TBM (n = 106) and PTB (n = 216). Using a phylogeny-adjusted genome-wide association method to count homoplasy events we examined phenotype-related changes at specific loci or genes in parallel branches of the phylogenetic tree. Enrichment scores for the TB phenotype were calculated on single nucleotide polymorphism (SNP), gene, and pathway level. Genetic associations were validated in an independent set of isolates.ResultsStrains belonged to the East-Asian lineage (36.0%), Euro-American lineage (61.5%), and Indo-Oceanic lineage (2.5%). We found no association between lineage and phenotype (Chi-square = 4.556; p = 0.207). Large genomic differences were observed between isolates; the minimum pairwise genetic distance varied from 17 to 689 SNPs. Using the phylogenetic tree, based on 28,544 common variable positions, we selected 54 TBM and 54 PTB isolates in terminal branch sets with distinct phenotypes. Genetic variation in Rv0218, and absence of Rv3343c, and nanK were significantly associated with disease phenotype in these terminal branch sets, and confirmed in the validation set of 214 unpaired isolates.ConclusionsUsing homoplasy counting we identified genetic variation in three separate genes to be associated with the TB phenotype, including one (Rv0218) which encodes a secreted protein that could play a role in host-pathogen interaction by altering pathogen recognition or acting as virulence effector.

DNA methylation in the cancerogenesis process and methods of its detection

Epigenetic modifications, apart from affecting gene expression, play an important role in the chromatin structure stabilization, embryonic development and the genomic imprinting. Recent studies have shown that they also play a vital role in other biological processes, including silencing of the expression and mobility of transposable elements and resistance to viral infections by blocking the expression of viral genes. The stability of the genome and the expression of genes in normal cells are strongly dependent on the DNA methylation pattern, which is visibly disturbed in tumor cells. These alterations may be a consequence of the attachment of methyl groups to cytosines in unmethylated DNA sequences, resulting in an increase in the degree of methylation or can be a result of demethylation, i.e. a reduction in the level of DNA methylation. Currently, many techniques are available to determine the level of methylcytosine in DNA, both at the level of single genes and the whole genome. However, each method has its advantages and disadvantages, not being universal in relation to the type of research material and the purpose of planned analyses..

An Exome-Wide Association Study Identifies New Susceptibility Loci for Age of Smoking Initiation in African- and European-American Populations.

Cigarette smoking is one of the largest causes of preventable death worldwide. This study aimed to identify susceptibility loci for age at smoking initiation (ASI) by performing an exome-wide association analysis. A total of 2510 smokers of either African-American (AA) or European-American (EA) origin were genotyped and analyzed at both the single nucleotide polymorphism (SNP) and gene levels. After removal of those SNPs with a minor allele frequency (<0.01), 48091 and 34933 SNPs for AAs and EAs, respectively, were used to conduct a SNP-based association analysis. Gene-based analyses were then performed for all SNPs examined within each gene. Further, we estimated the proportion of variance explained by all common SNPs included in the analysis. The strongest signals were detected for SNPs rs17849904 in the pitrilysin metallopeptidase 1 gene (PITRM1) in the AA sample (p = 9.02 × 10-7) and rs34722354 in the discoidin domain of the receptor tyrosine kinase 2 gene (DDR2) in the EA sample (p = 9.74 × 10-7). Both SNPs remained significant after Bonferroni correction for the number of SNPs tested. Subsequently, the gene-based association analysis revealed a significantly associated gene, DHRS7, in the AA sample (p = 5.00 × 10-6), a gene previously implicated in nicotine metabolism. Our study revealed two susceptibility loci for age of smoking initiation in the two ethnic samples, with the first being PITRM1 for AA smokers and the second DDR2 for EA smokers. In addition, we found DHRS7 to be a plausible candidate for ASI in the AA sample from our gene-based association analysis. PITRM1 and DHRS7 for African-American smokers and DDR2 for European-American smokers are new candidate genes for smoking initiation. These genes represent new additions to smoking initiation, an important but less studied phenotype in nicotine dependence research.

Reflection on teaching of boinformatics course for medical students in the big data era

Bioinformatics is an intercross field of biology and computer science. With the develop-ment of computer network information technology andprecision medical care, biomedical field ushered in the era of big data. The demand for bioinformatics knowledge to medical students has shifted from single gene level to the analysis on massive medical data. This paper analyzed the characteristics of the big data era, as well as the problems of the teaching idea, curriculum coverage and teaching platform could not meet the requirements of the new era in teaching of bioinformatics course for medical students. Finally, we pro-posed several measures to strengthen the reform of bioinformatics teaching to meet the challenge of the big data era, including optimizing courses, changing educational philosophy, introducing new teaching methods, and setting up teaching platform with times. Key words: Bioinformatics; Big data era; Teaching methods

A putative role for genome-wide epigenetic regulatory mechanisms in Huntington’s disease: A computational assessment

Background: Huntington's Disease (HD) is an incurable disease of the adult brain. Massive changes in gene expression are a prominent feature. Epigenetic effects have been reported to be implicated in HD, but the role of chromatin is not well understood. We tested if the chromatin state of dysregulated genes in HD is affected at a genome-wide scale and examined how epigenetic processes are associated with CpG-island-mediated gene expression. Methods: Our general approach incorporates computational and functional analysis of public data before embarking on expensive wet-lab experiments. We compared the location in the genome of the genes that were deregulated in HD human brain, obtained from public gene expression data, to the location of particular chromatin marks in reference tissues using data from the ENCODE project. Results: We found that differentially expressed genes were enriched in the active chromatin state, but not enriched in the silent state. In the caudate nucleus, the most highly affected brain region in HD, genes in the active state were associated with transcription, cell cycle, protein transport and modification, RNA splicing, histone post-translational modifications and RNA processing. Genes in the repressed state were linked with developmental processes and responses related to zinc and cadmium stimulus. We confirmed that genes within CpG-islands are enriched among HD dysregulated genes in human and mouse in HD. Epigenetic processes were associated more with genes that overlap with CpG-islands than genes that do not. Conclusion: Our results suggest that massive transcriptional dysregulation in HD is not matched by large-scale relocation of gene activity, i.e. inactive chromatin regions are altered into actively expressed chromatin regions and vice versa. We expect that changes in epigenetic chromatin state might occur at the level of single genes (e.g. promoters, gene body) and scattered genomic sites (e.g. CTCF sites, enhancer regions) instead of large-scale genomic regions.

Genetic and Transcriptomic Bases of Intestinal Epithelial Barrier Dysfunction in Inflammatory Bowel Disease.

Intestinal barrier defects are common in patients with inflammatory bowel disease (IBD). To identify which components could underlie these changes, we performed an in-depth analysis of epithelial barrier genes in IBD. A set of 128 intestinal barrier genes was selected. Polygenic risk scores were generated based on selected barrier gene variants that were associated with Crohn's disease (CD) or ulcerative colitis (UC) in our study. Gene expression was analyzed using microarray and quantitative reverse transcription polymerase chain reaction. Influence of barrier gene variants on expression was studied by cis-expression quantitative trait loci mapping and comparing patients with low- and high-risk scores. Barrier risk scores were significantly higher in patients with IBD than controls. At single-gene level, the associated barrier single-nucleotide polymorphisms were most significantly enriched in PTGER4 for CD and HNF4A for UC. As a group, the regulating proteins were most enriched for CD and UC. Expression analysis showed that many epithelial barrier genes were significantly dysregulated in active CD and UC, with overrepresentation of mucus layer genes. In uninflamed CD ileum and IBD colon, most barrier gene levels restored to normal, except for MUC1 and MUC4 that remained persistently increased compared with controls. Expression levels did not depend on cis-regulatory variants nor combined genetic risk. We found genetic and transcriptomic dysregulations of key epithelial barrier genes and components in IBD. Of these, we believe that mucus genes, in particular MUC1 and MUC4, play an essential role in the pathogenesis of IBD and could represent interesting targets for treatment.

Maltose Utilization as a Novel Selection Strategy for Continuous Evolution of Microbes with Enhanced Metabolite Production.

We have developed a novel selection circuit based on carbon source utilization that establishes and sustains growth-production coupling over several generations in a medium with maltose as the sole carbon source. In contrast to traditional antibiotic resistance-based circuits, we first proved that coupling of cell fitness to metabolite production by our circuit was more robust with a much lower escape risk even after many rounds of selection. We then applied the selection circuit to the optimization of L-tryptophan (l-Trp) production. We demonstrated that it enriched for specific mutants with increased l-Trp productivity regardless of whether it was applied to a small and defined mutational library or a relatively large and undefined one. From the latter, we identified four novel mutations with enhanced l-Trp output. Finally, we used it to select for several high l-Trp producers with randomly generated genome-wide mutations and obtained strains with up to 65% increased l-Trp production. This selection circuit provides new perspectives for the optimization of microbial cell factories for diverse metabolite production and the discovery of novel genotype-phenotype associations at the single-gene and whole-genome levels.

Intracoronary Imaging, Cholesterol Efflux, and Transcriptomics after Intensive Statin Treatment in Diabetes

Residual atherothrombotic risk remains higher in patients with versus without diabetes mellitus (DM) despite statin therapy. The underlying mechanisms are unclear. This is a retrospective post-hoc analysis of the YELLOW II trial, comparing patients with and without DM (non-DM) who received rosuvastatin 40 mg for 8–12 weeks and underwent intracoronary multimodality imaging of an obstructive nonculprit lesion, before and after therapy. In addition, blood samples were drawn to assess cholesterol efflux capacity (CEC) and changes in gene expression in peripheral blood mononuclear cells (PBMC). There was a significant reduction in low density lipoprotein-cholesterol (LDL-C), an increase in CEC and beneficial changes in plaque morphology including increase in fibrous cap thickness and decrease in the prevalence of thin cap fibro-atheroma by optical coherence tomography in DM and non-DM patients. While differential gene expression analysis did not demonstrate differences in PBMC transcriptome between the two groups on the single-gene level, weighted gene coexpression network analysis revealed two modules of coexpressed genes associated with DM, Collagen Module and Platelet Module, related to collagen catabolism and platelet function respectively. Bayesian network analysis revealed key driver genes within these modules. These transcriptomic findings might provide potential mechanisms responsible for the higher cardiovascular risk in DM patients.

Association of Intron Microsatellite Instability and Exon Mutational Profile of TP53 in Human Gastric Cancers.

Microsatellite instability (MSI) is a hallmark of genomic instability. In gastric cancer (GC), MSI phenotype is an important molecular subgroup that is often closely correlated with elevated mutation rates on the whole genome level. However, on a single gene level, it is still unknown whether the MSI status of a gene is correlated with the mutational profile of the gene itself. We analyzed intron MSI status and exon mutational profile of TP53 through short tandem repeat (STR) scanning and direct sequencing respectively in gastric cancers and their matched normal tissues. MSI status of the TP53 intron was significantly associated with the mutational profile of seven TP53 exon positions (p=0.0416), male patients (p=0.0095), and drinking (p=0.0474), and showed a mild correlation with longer survival time (p=0.0584) and increasing age (p=0.0611). TP53 exons tended to mutation in the status of TP53 intron MSI.

Detecting the genetic link between Alzheimer's disease and obesity using bioinformatics analysis of GWAS data.

Alzheimer's disease (AD) represents the major form of dementia in the elderly. In recent years, accumulating evidence indicate that obesity may act as a risk factor for AD, while the genetic link between the two conditions remains unclear. This bioinformatics analysis aimed to detect the genetic link between AD and obesity on single nucleotide polymorphisms (SNPs), gene, and pathway levels based on genome-wide association studies data. A total of 31 SNPs were found to be shared by AD and obesity, which were linked to 7 genes. These genes included PSMC3, CELF1, MYBPC3, SPI1, APOE, MTCH2 and RAPSN. Further functional enrichment analysis of these genes revealed the following biological pathways, including proteasome, osteoclast differentiation, hypertrophic cardiomyopathy, dilated cardiomyopathy, Epstein-Barr virus and TLV-I infection, as well as several cancer associated pathways, to be common among AD and obesity. The findings deepened our understanding on the genetic basis linking obesity and AD and may help shape possible prevention and treatment strategies.

Immune Gene Expression Is Associated with Genomic Aberrations in Breast Cancer.

The presence of tumor-infiltrating lymphocytes (TIL) is a favorable prognostic factor in breast cancer, but what drives immune infiltration remains unknown. Here we examine if clonal heterogeneity, total mutation load, neoantigen load, copy number variations (CNV), gene- or pathway-level somatic mutations, or germline polymorphisms (SNP) are associated with immune metagene expression in breast cancer subtypes. Thirteen published immune metagenes correlated separately with genomic metrics in the three major breast cancer subtypes. We analyzed RNA-Seq, DNA copy number, mutation and germline SNP data of 627 ER+, 207 HER2+, and 191 triple-negative (TNBC) cancers from The Cancer Genome Atlas. P-values were adjusted for multiple comparisons, and permutation testing was used to assess false discovery rates. Increased immune metagene expression associated significantly with lower clonal heterogeneity estimated by MATH score in all subtypes and with a trend for lower overall mutation, neoantigen, and CNV loads in TNBC and HER2+ cancers. In ER+ cancers, mutation load, neoantigen load, and CNV load weakly but positively associated with immune infiltration, which reached significance for overall mutation load only. No highly recurrent single gene or pathway level mutations associated with immune infiltration. High immune gene expression and lower clonal heterogeneity in TNBC and HER2+ cancers suggest an immune pruning effect and equilibrium between immune surveillance and clonal expansion. Thus, immune checkpoint inhibitors may tip the balance in favor of immune surveillance in these cancers. Cancer Res; 77(12); 3317-24. ©2017 AACR.

Orthogonal Genetic Regulation in Human Cells Using Chemically Induced CRISPR/Cas9 Activators.

The concerted action of multiple genes in a time-dependent manner controls complex cellular phenotypes, yet the temporal regulation of gene expressions is restricted on a single-gene level, which limits our ability to control higher-order gene networks and understand the consequences of multiplex genetic perturbations. Here we developed a system for temporal regulation of multiple genes. This system combines the simplicity of CRISPR/Cas9 activators for orthogonal targeting of multiple genes and the orthogonality of chemically induced dimerizing (CID) proteins for temporal control of CRISPR/Cas9 activator function. In human cells, these transcription activators exerted simultaneous activation of multiple genes and orthogonal regulation of different genes in a ligand-dependent manner with minimal background. We envision that our system will enable the perturbation of higher-order gene networks with high temporal resolution and accelerate our understanding of gene-gene interactions in a complex biological setting.