Human genetic variation recognizes functional elements in noncoding sequence

Genetic variants in ATP2B2 as risk factors for mortality in patients unrelated but not associated with families with severe COVID-19

Introduction

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of Coronavirus Disease 2019 (COVID-19). The disease has a wide range of clinical manifestations, from asymptomatic to severe. Ancestral contribution, sex, immune response, and genetic factors influence the presentation of the disease. The objective of the present study was to validate these genetic variants in patients with severe COVID-19 who died and in survivor patients. Methods: Single nucleotide variants (SNVs) in six genes: ATPase plasma membrane Ca2+ transporting 2 (ATP2B2), transmembrane serine protease 2 (TMPRSS2), dedicator of cytokinesis 2 (DOCK2), (interferon alpha and beta receptor subunit 2) IFNAR2, tumor necrosis factor receptor superfamily, member 1A (TNFRSF1A), and tumor necrosis factor receptor superfamily, member 1B (TNFRSF1B), were explored in two groups: the first consisted of severe COVID-19-related patients (familial cases from 58 families, n = 130), and the second group of unrelated severe COVID-19 patients (n = 1045). In each study group, death was evaluated as the outcome.

Results

In non-related patients with severe COVID-19, carriers of GG genotype (rs2289274) in the ATP2B2 gene showed a high-risk probability of non-surviving (OR = 1.43). Survival analysis to 75 days indicates that carriers of GG have a higher risk than GA or AA genotypes (p = 0.0059). The haplotype GG (rs2289273-rs2289274) in ATP2B2 was found to be associated with a high risk of death in severe non-related COVID-19 patients. No significant associations were found between severe COVID-19-related patients and SNVs in ATP2B2, TMPRSS2, DOCK2, IFNAR2, TNFRSF1A, or TNFRSF1B.

Conclusions

Unrelated patients with severe COVID-19 that carry the GG genotype (rs2289274) in ATP2B2 showed a high death risk. Survival analysis to 75 days indicates that carriers of GG have a higher risk of non-survival compared to GA or AA genotypes.

Enhancer-promoter activation by the Kaposi sarcoma-associated herpesvirus episome maintenance protein LANA

Higher-order genome structure influences the transcriptional regulation of cellular genes through the juxtaposition of regulatory elements, such as enhancers, close to promoters of target genes. While enhancer activation has emerged as an important facet of Kaposi sarcoma-associated herpesvirus (KSHV) biology, the mechanisms controlling enhancer-target gene expression remain obscure. Here, we discover that the KSHV genome tethering protein latency-associated nuclear antigen (LANA) potentiates enhancer-target gene expression in primary effusion lymphoma (PEL), a highly aggressive B cell lymphoma causally associated with KSHV. Genome-wide analyses demonstrate increased levels of enhancer RNA transcription as well as activating chromatin marks at LANA-bound enhancers. 3D genome conformation analyses identified genes critical for latency and tumorigenesis as targets of LANA-occupied enhancers, and LANA depletion results in their downregulation. These findings reveal a mechanism in enhancer-gene coordination and describe a role through which the main KSHV tethering protein regulates essential gene expression in PEL.

Fyn and Lyn gene polymorphisms impact the risk of thyroid cancer

Thyroid cancer is the most common malignancy of the endocrine glands, and during last couple of decades, its incidence has risen alarmingly, across the globe. Etiology of thyroid cancer is still debatable. There are a few worth mentioning risk factors which contribute to initiation of abnormalities in thyroid gland leading to cancer. Genetic instability is major risk factors in thyroid carcinogenesis. Among the genetic factors, the Src family of genes (Src, Yes1, Fyn and Lyn) have been implicated in many cancers but there is little data regarding the association of these (Src, Yes1, Fyn and Lyn) genes with thyroid carcinogenesis. Fyn and Lyn genes of Src family found engaged in proliferation, migration, invasion, angiogenesis, and metastasis in different cancers. This study was planned to examine the effect of Fyn and Lyn SNPs on thyroid cancer risk in Pakistani population in 500 patients and 500 controls. Three polymorphisms of Fyn gene (rs6916861, rs2182644 and rs12910) and three polymorphisms of Lyn gene (rs2668011, rs45587541 and rs45489500) were analyzed using Tetra-primer ARMS-PCR followed by DNA sequencing. SNP rs6916861 of Fyn gene mutant genotype (CC) showed statistically significant threefold increased risk of thyroid cancer (P < 0.0001). In case of rs2182644 of Fyn gene, mutant genotype (AA) indicated statistically significant 17-fold increased risk of thyroid cancer (P < 0.0001). Statistically significant threefold increased risk of thyroid cancer was observed in genotype AC (P < 0.0001) of Fyn gene polymorphism rs12910. In SNP rs2668011 of Lyn gene, TT genotype showed statistically significant threefold increased risk of thyroid cancer (P < 0.0001). In case of rs45587541 of Lyn gene, GA genotypes showed statistically significant 11-fold increased risk in thyroid cancer (P < 0.0001). Haplotype analysis revealed that AAATAG*, AGACAG*, AGCCAA*, AGCCAG*, CAATAG*, CGCCAG* and CGCCGA* haplotypes of Fyn and Lyn polymorphisms are associated with increased thyroid cancer risk. These results showed that genotypes and allele distribution of Fyn and Lyn are significantly linked with increased thyroid cancer risk and could be genetic adjuster for said disease.

A Post-GWAS Functional Analysis Confirming Effects of Three BTA13 Genes CACNB2, SLC39A12, and ZEB1 on Dairy Cattle Reproduction

In our previous GWAS of Chinese and Nordic dairy cattle, genes CACNB2, SLC39A12, and ZEB1 locating on BTA 13 were suggested as candidate genes for reproduction. In this study, validation of these associations was performed in an independent population with records of nine reproductive traits. More importantly, functions of these genes in the reproductive process were verified by employing the expression data of ovarian follicles. The potential variants within the three genes were firstly detected in 68 Chinese Holstein bulls, and then screened in 1,588 Chinese Holstein cows using the KASP (Kompetitive allele-specific PCR) method. There were nine variants with polymorphisms in CACNB2, five in SLC39A12, and four in ZEB1, respectively, of which one SNP was in the upstream regulatory region, two in exon region, four in downstream regulatory region, and 11 SNPs in intronic regions. Amongst the 18 variants, g.33267056T/G in CACNB2 explained the largest phenotypic variance for age at first calving (0.011%), interval from first to last insemination (0.004%), and calving ease (0.002%), while g.32751518G/A in SLC39A12 contributed the most to stillbirth in heifers (0.038%). Two haplotype blocks were constructed for CACNB2 while one each for SLC39A12 and ZEB1, which were significantly associated with five reproductive traits, including age at the first service, age at the first calving, calving ease in heifers and cows, and the interval from calving to the first insemination. We then studied the profile of gene expression in granulosa cells isolated from four developmental stages of ovarian follicles from eight dairy cows. All three genes were differentially expressed between ovarian follicles with different sizes (p < 0.05), indicating their potential roles in the reproductive process of dairy cows. This study successfully demonstrated the associations of three BTA 13 genes CACNB2, SLC39A12, and ZEB1 with reproduction and further examined their expression levels in ovarian follicles directly. These findings can be beneficial for the ongoing genomic selection program for reproductive traits which have long been considered as traits that are difficult to achieve genetic improvement due to the lack of efficient genetic markers.

Mitochondrial sirtuins genetic variations and gastric cancer risk: Evidence from retrospective observational study

AIMS

The purpose of the present study was to analyze the role of selected polymorphisms of SIRT3 and SIRT5 in gastric carcinogenesis.

METHODS

For this study, 500 blood samples of GC patients and 500 blood samples of healthy individuals were collected. Six selected polymorphisms of mitochondrial sirtuins were analyzed for analysis using Tetra-Arms PCR followed by DNA sequencing.

RESULTS

Mutant allele frequencies of selected polymorphisms [rs3782116 (p < 0.0001), rs6598072 (p < 0.0001) and rs11246020 (p < 0.0001), rs938222 (p = 0.0136), rs3757261 (p = 0.0005) and rs2841511 (p = 0.0015)] were observed significant higher in GC patients vs controls. Haplotype analysis was performed, and 51 haplotypes were generated using haploview software. Among these haplotypes, eleven haplotypes were found associated with a significantly increased risk of GC. Furthermore, SNP-SNP interaction showed a significant correlation between studied SNPs and GC risk. Kaplan Meier analysis showed that mutant allele frequencies of selected polymorphisms are linked with a significant decrease in survival of GC patients CONCLUSIONS: It can be concluded that selected SNPs may be associated with enhanced risk of GC and hence can be potential prognostic markers for prognosis and predisposition of GC.

Chemerin levels and its genetic variants are associated with Gestational Diabetes Mellitus: a hospital-based study in a Chinese cohort

BACKGROUND

Gestational diabetes mellitus (GDM) is a glucose intolerance condition encounters for the first time in a fraction of pregnant women. The role of different host inflammatory molecules in GDM etiology has been deciphered. Chemerin is a chemoattractant protein primarily associated with the pathogenesis of type 2 diabetes, obesity, and metabolic syndrome. However, the association of chemerin and its genetic variants with the predisposition of GDM is not clear, and our present study is aimed to address the issue.

MATERIALS AND METHODS

A total of 703 Chinese women comprising of GDM (n=303), glucose tolerant pregnant women (n=211), and non-pregnant glucose tolerant controls (n=189) were recruited in the present investigation. GDM was diagnosed according to the World Health Organization recommendation for diagnosis of gestational diabetes during pregnancy. Plasma levels of chemerin were quantified by an Enzyme-linked Immunosorbent Assay (ELISA). Common variants in the chemerin gene (rs4721, rs17173617, rs7806429, and rs17173608) were genotyped by using TaqMan assay.

RESULTS

Plasma chemerin level was found higher in subjects with GDM as compared to glucose tolerant pregnant and non-pregnant women. Further, a positive correlation between plasma chemerin and HOMA-IR index suggesting an essential role of chemerin in mediating insulin resistance. Variants of rs4721 and rs17173608 polymorphisms were associated with lower levels of plasma chemerin and low HOMA-IR index. Furthermore, mutants of rs4721 and rs17173608 polymorphisms were associated with protection against the development of GDM in the Chinese cohort.

CONCLUSIONS

Plasma chemerin is elevated in GDM patients. Genetic variation in chemerin gene associated with lower plasma levels of chemerin, HOMA-IR index and protects against the development of GDM in Chinese.

Genetic Variations in Prostaglandin E2 Pathway Identified as Susceptibility Biomarkers for Gastric Cancer in an Intermediate Risk European Country

The cyclooxygenase-2 (COX-2)/prostaglandin E₂ (PGE₂) pathway exerts deleterious pleiotropic effects in inflammation-induced gastric carcinogenesis. We aimed to assess the association of genetic variants in prostaglandin-endoperoxide synthase 2 (PTGS2), ATP binding cassette subfamily C member 4 (ABCC4), hydroxyprostaglandin dehydrogenase 15-(NAD) (HPGD), and solute carrier organic anion transporter family member 2A1 (SLCO2A1) PGE₂ pathway-related genes with gastric cancer (GC) risk in a European Caucasian population. A hospital-based case-control study gathering 260 GC cases and 476 cancer-free controls was implemented. Using a tagSNP approach, 51 single nucleotide polymorphisms (SNPs) were genotyped through MassARRAY^® iPLEX Gold Technology or allelic discrimination by real-time polymerase chain reaction (PCR). Homozygous carriers of the minor allele for both rs689466 and rs10935090 SNPs were associated with a 2.98 and 4.30-fold increased risk for GC, respectively (95% confidence interval (CI): 1.14–7.74, p = 0.027; 95% CI: 1.22–15.16, p = 0.026), with the latter also being associated with an anticipated diagnosis age. A multifactor dimensionality reduction analysis identified an overall three-factor best interactive model composed of age, rs689466, and rs1678374 that was associated with a 17.6-fold GC increased risk (95% CI: 11.67–26.48, p < 0.0001, (cross-validation) CV consistency of 8/10 and accuracy of 0.807). In this preliminary study, several tagSNPs in PGE₂ pathway-related genes were identified as risk biomarkers for GC development. This approach may help to identify higher-risk individuals and may contribute to the tailoring screening of GC in intermediate-risk European countries.

Linkage disequilibrium and haplotype analysis of Src and Yes1 genes in thyroid cancer

Purpose: This study was planned to examine the effects of Src and Yes1 single nucleotide polymorphism (SNPs) on the risk of thyroid cancer in 499 patients and 500 controls. Materials & methods: Three SNPs of Src gene and three SNPs of Yes1 gene were analyzed using Tetra-primer ARMS-PCR followed by sequencing. Results: rs121913314 of Src gene genotype TT showed 32-fold increased risk of thyroid cancer and rs2305994 of Yes1 genotypes TT and CT showed 2.7-fold and 16-fold increased risk in thyroid cancer (p < 0.0001). Haplotype analysis revealed that CATGCC, CATGCT, CATGTC, CATGTT, TATGCC and TATGTTA haplotypes are associated with thyroid cancer risk. Conclusion: Results showed that genotypes and allele distribution of Src and Yes1 genes are significantly linked with increased risk of thyroid cancer.

Association of intronic polymorphisms (rs1549339, rs13402242) and mRNA expression variations in PSMD1 gene in arsenic-exposed workers

Ubiquitin-proteasome system (UPS) gene, PSMD1, is an important gene for neutralization of damaged and misfolded protein(s). The current study was designed to study the genetic and expression variations of PSMD1 gene as a consequence of arsenic exposure and its potential implications in arsenic induced diseases. In the present study, 250 blood samples of exposed industrial workers along with 250 controls were used. Initially, tetra amplification refractory mutation system-PCR was used to determine the role of PSMD1 gene polymorphisms (rs1549339, rs13402242) in industrial workers and controls. Frequency of homozygous mutant genotype of rs1549339 (OR: 2.23, 95% CI: 1.51-3.32, p = 0.0001) and rs13402242 (OR: 2.96, 95% CI: 1.52-5.75, p = 0.001) was observed significantly higher in exposed individuals vs controls. Secondly, qPCR was performed for expression analysis of PSMD1 gene. Significant down-regulated expression of PSMD1 gene (p < 0.0001) was observed vs controls, and this down-regulation was observed more pronounced in smokers (p < 0.0001) with maximum exposure duration (p < 0.0008). This down-regulated expression was observed significantly more pronounced in welding (p < 0.004) and brick kiln industries (p < 0.04) compared to other selected industries. The obtained results suggest that the exposure to arsenic may have an increased risk of developing disease(s) because of arsenic-induced PSMD1 variations.

SNP Diversity in CD14 Gene Promoter Suggests Adaptation Footprints in Trypanosome Tolerant N'Dama (Bos taurus) but not in Susceptible White Fulani (Bos indicus) Cattle

Immune response to infections has been shown to be mediated by genetic diversity in pattern recognition receptors, leading to disease tolerance or susceptibility. We elucidated naturally occurring variations within the bovine CD14 gene promoter in trypanosome-tolerant (N'Dama) and susceptible (White Fulani) cattle, with genomic and computational approaches. Blood samples were collected from White Fulani and N'Dama cattle, genomic DNA extracted and the entire promoter region of the CD14 gene amplified by PCR. We sequenced this region and performed in silico computation to identify SNP variants, transcription factor binding sites, as well as micro RNAs in the region. CD14 promoter sequences were compared with the reference bovine genome from the Ensembl database to identify various SNPs. Furthermore, we validated three selected N'Dama specific SNPs using custom Taqman SNP genotyping assay for genetic diversity. In all, we identified a total of 54 and 41 SNPs at the CD14 promoter for N'Dama and White Fulani respectively, including 13 unique SNPs present in N'Dama only. The significantly higher SNP density at the CD14 gene promoter region in N'Dama may be responsible for disease tolerance, possibly an evolutionary adaptation. Our genotype analysis of the three loci selected for validation show that mutant alleles (A/A, C/C, and A/A) were adaptation profiles within disease tolerant N'Dama. A similar observation was made for our haplotype analysis revealing that haplotypes H1 (ACA) and H2 (ACG) were significant combinations within the population. The SNP effect prediction revealed 101 and 89 new transcription factor binding sites in N'Dama and White Fulani, respectively. We conclude that disease tolerant N'Dama possessing higher SNP density at the CD14 gene promoter and the preponderance of mutant alleles potentially confirms the significance of this promoter in immune response, which is lacking in susceptible White Fulani. We, therefore, recommend further in vitro and in vivo study of this observation in infected animals, as the next step for understanding genetic diversity relating to varying disease phenotypes in both breeds.

Molecular Imaging of Opioid and Dopamine Systems: Insights Into the Pharmacogenetics of Opioid Use Disorders

Opioid use in the United States has steadily risen since the 1990s, along with staggering increases in addiction and overdose fatalities. With this surge in prescription and illicit opioid abuse, it is paramount to understand the genetic risk factors and neuropsychological effects of opioid use disorder (OUD). Polymorphisms disrupting the opioid and dopamine systems have been associated with increased risk for developing substance use disorders. Molecular imaging studies have revealed how these polymorphisms impact the brain and contribute to cognitive and behavioral differences across individuals. Here, we review the current molecular imaging literature to assess how genetic variations in the opioid and dopamine systems affect function in the brain's reward, cognition, and stress pathways, potentially resulting in vulnerabilities to OUD. Continued research of the functional consequences of genetic variants and corresponding alterations in neural mechanisms will inform prevention and treatment of OUD.

Contribution of Genetic Polymorphisms and Haplotypes in DRD2, BDNF, and Opioid Receptors to Heroin Dependence and Endophenotypes Among the Han Chinese

Heroin and drug dependence are major contributors to global health burden worldwide, but their underlying mechanisms remain elusive and may vary from population to population. Reward- and memory-related candidate genes dopamine D2 receptor (DRD2) and brain-derived neurotrophic factor (BDNF), as well as the opioid receptor genes (OPRM1, OPRD1, and OPRK1), have been implicated in drug dependence, but relatively little is known on their contributions to heroin dependence in populations worldwide. Hence, we evaluated the contributions of the above five candidate genes in heroin dependence and several important related endophenotypes (the onset age of heroin use and subjective response to first heroin use), at single single-nucleotide polymorphism as well as haplotype levels, in a Han Chinese population sample. We genotyped 546 unrelated and heroin-dependent subjects for the candidate genes noted, and 228 sex- and age-matched unrelated controls. The G allele of rs4654327 (OPRD1), DRD2 haplotype block CCGCCGTT (rs6277-rs1076560-rs2283265-rs2734833-rs2075652-rs1079596-rs4436578-rs11214607), and OPRD1 haplotypes TACG (rs6669447-rs2236857-rs508448-rs4654327), CG (rs508448-rs4654327), and TG (rs6669447-rs4654327) were significantly associated with heroin dependence phenotype. Homozygotes AA at rs6265 (BDNF), TT at rs16917234 (BDNF), and CC at rs508448 (OPRD1) also appeared as risk factors for the endophenotype earlier age of onset for heroin use. Two OPRM1 haplotypes, AG (rs1799971-rs1381376) and AT (rs1799971-rs3778151), were observed as potential protective factors. These emerging findings contribute to the literature on genetic biomarkers of drug dependence and related endophenotypes, and call for replication in independent population.

Epigenetic modifications in KDM lysine demethylases associate with survival of early-stage NSCLC

Background

KDM lysine demethylase family members are related to lung cancer clinical outcomes and are potential biomarkers for chemotherapeutics. However, little is known about epigenetic alterations in KDM genes and their roles in lung cancer survival.

Methods

Tumor tissue samples of 1230 early-stage non-small cell lung cancer (NSCLC) patients were collected from the five independent cohorts. The 393 methylation sites in KDM genes were extracted from epigenome-wide datasets and analyzed by weighted random forest (Ranger) in discovery phase and validation dataset, respectively. The variable importance scores (VIS) for the sites in top 5% of both discovery and validation sets were carried forward for Cox regression to further evaluate the association with patient’s overall survival. TCGA transcriptomic data were used to evaluate the correlation with the corresponding DNA methylation.

Results

DNA methylation at sites cg11637544 in KDM2A and cg26662347 in KDM1A were in the top 5% of VIS in both discovery phase and validation for squamous cell carcinomas (SCC), which were also significantly associated with SCC survival (HR_cg11637544 = 1.32, 95%CI, 1.16–1.50, P = 1.1 × 10⁻⁴; HR_cg26662347 = 1.88, 95%CI, 1.37–2.60, P = 3.7 × 10⁻³), and correlated with corresponding gene expression (cg11637544 for KDM2A, P = 1.3 × 10⁻¹⁰; cg26662347 for KDM1A P = 1.5 × 10⁻⁵). In addition, by using flexible criteria for Ranger analysis followed by survival classification tree analysis, we identified four clusters for adenocarcinomas and five clusters for squamous cell carcinomas which showed a considerable difference of clinical outcomes with statistical significance.

Conclusions

These findings highlight the association between somatic DNA methylation in KDM genes and early-stage NSCLC patient survival, which may reveal potential epigenetic therapeutic targets.

Electronic supplementary material

The online version of this article (10.1186/s13148-018-0474-3) contains supplementary material, which is available to authorized users.

Haplotype analysis of XRCC2 gene polymorphisms and association with increased risk of head and neck cancer

We aimed to investigate the effect of hotspot variations of XRCC2 gene on the risk of head and neck cancer (HNC) in 400 patients and 400 controls. Five polymorphisms of XRCC2 gene G4234C (rs3218384), G4088T (rs3218373), G3063A (rs2040639), R188H (rs3218536) and rs7802034 were analyzed using Allele- specific polymerase chain reaction (ARMS-PCR) followed by sequence analysis. For rs3218373, the GG genotype indicated a statistically significant 3-fold increased risk of HNC (P < 0.001) after multivariate adjustment. For rs7802034, the GG genotype suggested statistically significant 2-fold increased risk of HNC (P < 0.001). For SNP of rs3218536, the AA genotype indicated a significant 3-fold increased risk of HNC (P < 0.001). Additionally, haplotype analysis revealed that TACAG, TGGAG, TACGG and TAGGA haplotypes of XRCC2 polymorphisms are associated with HNC risk. Two SNPs in XRCC2 (rs2040639 and rs3218384) were found increased in strong linkage disequilibrium. Furthermore, joint effect model showed 20 fold (OR = 19.89; 95% CI = 2.65–149.36, P = 0.003) increased HNC risk in patients carrying four homozygous risk alleles of selected polymorphisms. These results show that allele distributions and genotypes of XRCC2 SNPs are significantly associated with increased HNC risk and could be a genetic adjuster for the said disease.

Analysis of Genes with Alternatively Spliced Transcripts in the Leaf, Root, Panicle and Seed of Rice Using a Long Oligomer Microarray and RNA-Seq

Pre-mRNA splicing further increases protein diversity acquired through evolution. The underlying driving forces for this phenomenon are unknown, especially in terms of gene expression. A rice alternatively spliced transcript detection microarray (ASDM) and RNA sequencing (RNA-Seq) were applied to differentiate the transcriptome of 4 representative organs of Oryza sativa L. cv. Ilmi: leaves, roots, 1-cm-stage panicles and young seeds at 21 days after pollination. Comparison of data obtained by microarray and RNA-Seq showed a bell-shaped distribution and a co-lineation for highly expressed genes. Transcripts were classified according to the degree of organ enrichment using a coefficient value (CV, the ratio of the standard deviation to the mean values): highly variable (CVI), variable (CVII), and constitutive (CVIII) groups. A higher index of the portion of loci with alternatively spliced transcripts in a group (IAST) value was observed for the constitutive group. Genes of the highly variable group showed the characteristics of the examined organs, and alternatively spliced transcripts tended to exhibit the same organ specificity or less organ preferences, with avoidance of 'organ distinctness'. In addition, within a locus, a tendency of higher expression was found for transcripts with a longer coding sequence (CDS), and a spliced intron was the most commonly found type of alternative splicing for an extended CDS. Thus, pre-mRNA splicing might have evolved to retain maximum functionality in terms of organ preference and multiplicity.

Prioritizing single-nucleotide polymorphisms and variants associated with clinical mastitis

Next-generation sequencing technology has provided resources to easily explore and identify candidate single-nucleotide polymorphisms (SNPs) and variants. However, there remains a challenge in identifying and inferring the causal SNPs from sequence data. A problem with different methods that predict the effect of mutations is that they produce false positives. In this hypothesis, we provide an overview of methods known for identifying causal variants and discuss the challenges, fallacies, and prospects in discerning candidate SNPs. We then propose a three-point classification strategy, which could be an additional annotation method in identifying causalities.

Genomic strategies to understand causes of keratoconus

Keratoconus (KTCN) is a degenerative disorder of the eye characterized by the conical shape and thinning of the cornea. The abnormal structure of KTCN-affected cornea results in loss of visual acuity. While many studies examine how environmental factors influence disease development, finding the genetic triggers has been a major emphasis of KTCN research. This paper focuses on genomic strategies that were implemented for finding candidate genes, including linkage and association studies, and presents different approaches of mutation screening. The advantages and limitations of particular tools are discussed based on literature and personal experience. Since etiology underlying KTCN is complex, numerous findings indicating heterogeneity of genetic factors involved KTCN etiology are presented.

Genomic translational research: Paving the way to individualized cardiac functional analyses and personalized cardiology

For most of Medicine's past, the best that physicians could do to cope with disease prevention and treatment was based on the expected response of an average patient. Currently, however, a more personalized/precise approach to cardiology and medicine in general is becoming possible, as the cost of sequencing a human genome has declined substantially. As a result, we are witnessing an era of precipitous advances in biomedicine and bourgeoning understanding of the genetic basis of cardiovascular and other diseases, reminiscent of the resurgence of innovations in physico-mathematical sciences and biology-anatomy-cardiology in the Renaissance, a parallel time of radical change and reformation of medical knowledge, education and practice. Now on the horizon is an individualized, diverse patient-centered, approach to medical practice that encompasses the development of new, gene-based diagnostics and preventive medicine tactics, and offers the broadest range of personalized therapies based on pharmacogenetics. Over time, translation of genomic and high-tech approaches unquestionably will transform clinical practice in cardiology and medicine as a whole, with the adoption of new personalized medicine approaches and procedures. Clearly, future prospects far outweigh present accomplishments, which are best viewed as a promising start. It is now essential for pluridisciplinary health care providers to examine the drivers and barriers to the clinical adoption of this emerging revolutionary paradigm, in order to expedite the realization of its potential. So, we are not there yet, but we are definitely on our way.

Direct Transcriptional Consequences of Somatic Mutation in Breast Cancer

Disordered transcriptomes of cancer encompass direct effects of somatic mutation on transcription, coordinated secondary pathway alterations, and increased transcriptional noise. To catalog the rules governing how somatic mutation exerts direct transcriptional effects, we developed an exhaustive pipeline for analyzing RNA sequencing data, which we integrated with whole genomes from 23 breast cancers. Using X-inactivation analyses, we found that cancer cells are more transcriptionally active than intermixed stromal cells. This is especially true in estrogen receptor (ER)-negative tumors. Overall, 59% of substitutions were expressed. Nonsense mutations showed lower expression levels than expected, with patterns characteristic of nonsense-mediated decay. 14% of 4,234 rearrangements caused transcriptional abnormalities, including exon skips, exon reusage, fusions, and premature polyadenylation. We found productive, stable transcription from sense-to-antisense gene fusions and gene-to-intergenic rearrangements, suggesting that these mutation classes drive more transcriptional disruption than previously suspected. Systematic integration of transcriptome with genome data reveals the rules by which transcriptional machinery interprets somatic mutation.

•

Greater transcriptional activity in cancer than stromal cells, particularly when ER-ve

•

Intron mutations only infrequently affect splicing, even at essential splice sites

•

Distinctive RNA effects of sense-to-antisense and gene-to-intergenic rearrangements

•

Exhaustive pipeline for identifying aberrant transcripts from RNA-sequencing data

Enhancer-promoter interactions are encoded by complex genomic signatures on looping chromatin

Discriminating the gene target of a distal regulatory element from other nearby transcribed genes is a challenging problem with the potential to illuminate the causal underpinnings of complex diseases. We present TargetFinder, a computational method that reconstructs regulatory landscapes from diverse features along the genome. The resulting models accurately predict individual enhancer-promoter interactions across multiple cell lines with a false discovery rate up to 15 times smaller than that obtained using the closest gene. By evaluating the genomic features driving this accuracy, we uncover interactions between structural proteins, transcription factors, epigenetic modifications, and transcription that together distinguish interacting from non-interacting enhancer-promoter pairs. Most of this signature is not proximal to the enhancers and promoters but instead decorates the looping DNA. We conclude that complex but consistent combinations of marks on the one-dimensional genome encode the three-dimensional structure of fine-scale regulatory interactions.

A Dual Model for Prioritizing Cancer Mutations in the Non-coding Genome Based on Germline and Somatic Events

We address here the issue of prioritizing non-coding mutations in the tumoral genome. To this aim, we created two independent computational models. The first (germline) model estimates purifying selection based on population SNP data. The second (somatic) model estimates tumor mutation density based on whole genome tumor sequencing. We show that each model reflects a different set of constraints acting either on the normal or tumor genome, and we identify the specific genome features that most contribute to these constraints. Importantly, we show that the somatic mutation model carries independent functional information that can be used to narrow down the non-coding regions that may be relevant to cancer progression. On this basis, we identify positions in non-coding RNAs and the non-coding parts of mRNAs that are both under purifying selection in the germline and protected from mutation in tumors, thus introducing a new strategy for future detection of cancer driver elements in the expressed non-coding genome.

Cancer cells undergo a mutation/selection process that resembles that of any living cell. Most mutations in cancer cell DNA occur in the so-called "non-coding" regions that represent 98.5% of the genome length. Pinning down which of these mutations contribute to the fitness of cancer cells would be important for identifying new "cancer drivers", which may in turn lead to future treatments. Unfortunately, predicting the impact of a non-coding DNA alteration remains extremely difficult. In this study, we analyze millions of non-coding cancer mutations and show cancer-specific mutational patterns can be used to predict non-coding regions that are preserved from mutations and may thus be important for cancer cell survival. Combining this information with population data, we propose a new scoring system that should help prioritize important non-coding mutations in future studies.

Molecular Screening of Keratoconus Susceptibility Sequence Variants in VSX1, TGFBI, DOCK9, STK24, and IPO5 Genes in Polish Patients and Novel TGFBI Variant Identification

PURPOSE

Keratoconus (KTCN) is a degenerative disorder of the eye that results in the conical shape and thinning of the cornea and is a leading cause for corneal transplantations. A number of studies suggest that genetic factors play a role in KTCN etiology. Some candidate gene variants have recently been shown to be associated with KTCN. The purpose of our study was to verify the role of VSX1, TGFBI, DOCK9, IPO5, and STK24 sequence variants in Polish KTCN patients.

METHODS

Forty-two Polish patients with sporadic KTCN and 50 control individuals were enrolled into this study. Both affected and unaffected individuals underwent detailed ophthalmic examination. The mutations screening in the candidate genes was performed by the direct sequencing method.

RESULTS

Analysis of VSX1, TGFBI, DOCK9, IPO5, and STK24 genes identified numerous sequence variants. Variants c.-264_-255delGGGGTGGGGT, c.627 + 23G > A, c.809-6_809-5insT, and c.*200G > T in the VSX1 gene, and heterozygous c.1598G > A mutation (Arg533Gln) in exon 12 of TGFBI were detected for the first time in KTCN patients. Two known sequence variants of TGFBI c.1620T > C (Phe540Phe) and c.1678 + 23G > A were observed in KTCN patients and control individuals. The newly reported c.717 + 43A > G substitution in intron 7 of DOCK9 was identified in both KTCN patients and healthy individuals.

CONCLUSIONS

Our investigation showed that KTCN-related sequence variants of analyzed genes were found in a very small proportion of the studied patients indicating that genes other than VSX1, TGFBI, DOCK9, IPO5, and STK24 are involved in the development and progression of KTCN in Polish patients. Our results support the hypothesis about the genetic heterogeneity of KTCN.

Weak negative and positive selection and the drift load at splice sites

Splice sites (SSs) are short sequences that are crucial for proper mRNA splicing in eukaryotic cells, and therefore can be expected to be shaped by strong selection. Nevertheless, in mammals and in other intron-rich organisms, many of the SSs often involve nonconsensus (Nc), rather than consensus (Cn), nucleotides, and beyond the two critical nucleotides, the SSs are not perfectly conserved between species. Here, we compare the SS sequences between primates, and between Drosophila fruit flies, to reveal the pattern of selection acting at SSs. Cn-to-Nc substitutions are less frequent, and Nc-to-Cn substitutions are more frequent, than neutrally expected, indicating, respectively, negative and positive selection. This selection is relatively weak (1 < |4Nes| < 4), and has a similar efficiency in primates and in Drosophila. Within some nucleotide positions, the positive selection in favor of Nc-to-Cn substitutions is weaker than the negative selection maintaining already established Cn nucleotides; this difference is due to site-specific negative selection favoring current Nc nucleotides. In general, however, the strength of negative selection protecting the Cn alleles is similar in magnitude to the strength of positive selection favoring replacement of Nc alleles, as expected under the simple nearly neutral turnover. In summary, although a fraction of the Nc nucleotides within SSs is maintained by selection, the abundance of deleterious nucleotides in this class suggests a substantial genome-wide drift load.

HPA axis genetic variation, cortisol and psychosis in major depression

Genetic variation underlying hypothalamic pituitary adrenal (HPA) axis overactivity in healthy controls (HCs) and patients with severe forms of major depression has not been well explored, but could explain risk for cortisol dysregulation. In total, 95 participants were studied: 40 patients with psychotic major depression (PMD); 26 patients with non-psychotic major depression (NPMD); and 29 HCs. Collection of genetic material was added one third of the way into a larger study on cortisol, cognition and psychosis in major depression. Subjects were assessed using the Brief Psychiatric Rating Scale, the Hamilton Depression Rating Scale and the Structured Clinical Interview for Diagnostic and Statistical Manual of Mental Disorders. Blood was collected hourly for determination of cortisol from 1800 to 0900 h and for the assessment of alleles for six genes involved in HPA axis regulation. Two of the six genes contributed significantly to cortisol levels, psychosis measures or depression severity. After accounting for age, depression and psychosis, and medication status, only allelic variation for the glucocorticoid receptor (GR) gene accounted for a significant variance for mean cortisol levels from 1800 to 0100 h (r(2)=0.288) and from 0100 to 0900 h (r(2)=0.171). In addition, GR and corticotropin-releasing hormone receptor 1 (CRHR1) genotypes contributed significantly to psychosis measures and CRHR1 contributed significantly to depression severity rating.

Association of OPRD1 polymorphisms with heroin dependence in a large case-control series

Genes encoding the opioid receptors (OPRM1, OPRD1 and OPRK1) are obvious candidates for involvement in risk for heroin dependence. Prior association studies commonly had samples of modest size, included limited single nucleotide polymorphism (SNP) coverage of these genes and yielded inconsistent results. Participants for the current investigation included 1459 heroin-dependent cases ascertained from maintenance clinics in New South Wales, Australia, 1495 unrelated individuals selected from an Australian sample of twins and siblings as not meeting DSM-IV criteria for lifetime alcohol or illicit drug dependence (non-dependent controls) and 531 controls ascertained from economically disadvantaged neighborhoods in proximity to the maintenance clinics. A total of 136 OPRM1, OPRD1 and OPRK1 SNPs were genotyped in this sample. After controlling for admixture with principal components analysis, our comparison of cases to non-dependent controls found four OPRD1 SNPs in fairly high linkage disequilibrium for which adjusted P values remained significant (e.g. rs2236857; OR 1.25; P=2.95×10(-4) ) replicating a previously reported association. A post hoc analysis revealed that the two SNP (rs2236857 and rs581111) GA haplotype in OPRD1 is associated with greater risk (OR 1.68; P=1.41×10(-5) ). No OPRM1 or OPRK1 SNPs reached more than nominal significance. Comparisons of cases to neighborhood controls reached only nominal significance. Our results replicate a prior report providing strong evidence implicating OPRD1 SNPs and, in particular, the two SNP (rs2236857 and rs581111) GA haplotype in liability for heroin dependence. Support was not found for similar association involving either OPRM1 or OPRK1 SNPs.

Effects of hydroxy-delta-5-steroid dehydrogenase, 3 beta- and steroid delta-isomerase 1 polymorphisms on fat androstenone level and gene expression in Duroc pigs

A high level of androstenone in porcine adipose tissue is a major factor contributing to boar taint. Porcine hydroxy-delta-5-steroid dehydrogenase, 3 beta- and steroid delta-isomerase 1 (3β-HSD, also known as HSD3B1) plays a key role in the hepatic metabolism that catalyzes androstenone to β-androstenol. Therefore, 3β-HSD is a candidate gene for boar taint. This study aimed to investigate functional 3β-HSD polymorphisms in Duroc pigs. We found eight single nucleotide polymorphisms (SNPs) in the full-length porcine 3β-HSD. Four of the SNPs had restriction enzyme sites, and we genotyped them in 147 uncastrated male Duroc pigs using a polymerase chain reaction-restriction fragment length polymorphism method. Pigs with the GG genotype at the g.165262G>A locus (SNP5) had significantly lower androstenone levels than did those with other genotypes (P = 0.030). SNP5 also was associated with differences in 3β-HSD mRNA levels: pigs with the GG genotype had higher levels than those with other genotypes (P = 0.019). The SNP5 polymorphism could affect the hepatic catabolism of androstenone and consequently impact androstenone accumulation in the adipose tissue. Therefore, SNP5 in the 3β-HSD of Duroc pigs could be a useful selective marker for decreasing boar taint.

From expression pattern to genetic association in asthma and asthma-related phenotypes

BACKGROUND

Asthma is a complex disease characterized by hyperresponsiveness, obstruction and inflammation of the airways. To date, several studies using different approaches as candidate genes approach, genome wide association studies, linkage analysis and genomic expression leaded to the identification of over 300 genes involved in asthma pathophysiology. Combining results from two studies of genomic expression, this study aims to perform an association analysis between genes differently expressed in bronchial biopsies of asthmatics compared to controls and asthma-related phenotypes using the same French-Canadian Caucasian population.

RESULTS

Before correction, 31 of the 85 genes selected were associated with at least one asthma-related phenotype. We found four genes that survived the correction for multiple testing. The rs11630178 in aggrecan gene (AGC1) is associated with atopy (p=0.0003) and atopic asthma (p=0.0001), the rs1247653 in the interferon alpha-inducible protein 6 (IFI6), the rs1119529 in adrenergic, alpha-2A-, receptor (ADRA2A) and the rs13103321 in the alcohol dehydrogenase 1B (class I), beta polypeptide (ADH1B), are associated with asthma (p=0.019; 0.01 and 0.002 respectively).

CONCLUSION

To our knowledge, this is the first time those genes are associated with asthma and related traits. Consequently, our study confirms that genetic and expression studies are complementary to identify new candidate genes and to investigate their role to improve the comprehension of the complexity of asthma pathophysiology.

Evidence of abundant purifying selection in humans for recently acquired regulatory functions

Although only 5% of the human genome is conserved across mammals, a substantially larger portion is biochemically active, raising the question of whether the additional elements evolve neutrally or confer a lineage-specific fitness advantage. To address this question, we integrate human variation information from the 1000 Genomes Project and activity data from the ENCODE Project. A broad range of transcribed and regulatory nonconserved elements show decreased human diversity, suggesting lineage-specific purifying selection. Conversely, conserved elements lacking activity show increased human diversity, suggesting that some recently became nonfunctional. Regulatory elements under human constraint in nonconserved regions were found near color vision and nerve-growth genes, consistent with purifying selection for recently evolved functions. Our results suggest continued turnover in regulatory regions, with at least an additional 4% of the human genome subject to lineage-specific constraint.

Genetic variations and alternative splicing: the Glioma associated oncogene 1, GLI1

Alternative splicing is a post-transcriptional regulatory process that is attaining stronger recognition as a modulator of gene expression. Alternative splicing occurs when the primary RNA transcript is differentially processed into more than one mature RNAs. This is the result of a variable definition/inclusion of the exons, the sequences that are excised from the primary RNA to form the mature RNAs. Consequently, RNA expression can generate a collection of differentially spliced RNAs, which may distinctly influence subsequent biological events, such as protein synthesis or other biomolecular interactions. Still the mechanisms that control exon definition and exon inclusion are not fully clarified. This mini-review highlights advances in this field as well as the impact of single nucleotide polymorphisms in affecting splicing decisions. The Glioma-associated oncogene 1, GLI1, is taken as an example in addressing the role of nucleotide substitutions for splicing regulation.

Mutation in a primate-conserved retrotransposon reveals a noncoding RNA as a mediator of infantile encephalopathy

The human genome is densely populated with transposons and transposon-like repetitive elements. Although the impact of these transposons and elements on human genome evolution is recognized, the significance of subtle variations in their sequence remains mostly unexplored. Here we report homozygosity mapping of an infantile neurodegenerative disease locus in a genetic isolate. Complete DNA sequencing of the 400-kb linkage locus revealed a point mutation in a primate-specific retrotransposon that was transcribed as part of a unique noncoding RNA, which was expressed in the brain. In vitro knockdown of this RNA increased neuronal apoptosis, consistent with the inappropriate dosage of this RNA in vivo and with the phenotype. Moreover, structural analysis of the sequence revealed a small RNA-like hairpin that was consistent with the putative gain of a functional site when mutated. We show here that a mutation in a unique transposable element-containing RNA is associated with lethal encephalopathy, and we suggest that RNAs that harbor evolutionarily recent repetitive elements may play important roles in human brain development.

Novel mutation and three other sequence variants segregating with phenotype at keratoconus 13q32 susceptibility locus

Keratoconus (KTCN), a non-inflammatory corneal disorder characterized by stromal thinning, represents a major cause of corneal transplantations. Genetic and environmental factors have a role in the etiology of this complex disease. Previously reported linkage analysis revealed that chromosomal region 13q32 is likely to contain causative gene(s) for familial KTCN. Consequently, we have chosen eight positional candidate genes in this region: MBNL1, IPO5, FARP1, RNF113B, STK24, DOCK9, ZIC5 and ZIC2, and sequenced all of them in 51 individuals from Ecuadorian KTCN families and 105 matching controls. The mutation screening identified one mutation and three sequence variants showing 100% segregation under a dominant model with KTCN phenotype in one large Ecuadorian family. These substitutions were found in three different genes: c.2262A>C (p.Gln754His) and c.720+43A>G in DOCK9; c.2377-132A>C in IPO5 and c.1053+29G>C in STK24. PolyPhen analyses predicted that c.2262A>C (Gln754His) is possibly damaging for the protein function and structure. Our results suggest that c.2262A>C (p.Gln754His) mutation in DOCK9 may contribute to the KTCN phenotype in the large KTCN-014 family.

SLC6A4 methylation modifies the effect of the number of traumatic events on risk for posttraumatic stress disorder

BACKGROUND

Posttraumatic stress disorder (PTSD) is a common and debilitating mental disorder that occurs following exposure to a traumatic event. However, most individuals do not develop PTSD following even a severe trauma, leading to a search for new variables, such as genetic and other molecular variation, associated with vulnerability and resilience in the face of trauma exposure.

METHOD

We examined whether serotonin transporter (SLC6A4) promoter genotype and methylation status modified the association between number of traumatic events experienced and PTSD in a subset of 100 individuals from the Detroit Neighborhood Health Study.

RESULTS

Number of traumatic events was strongly associated with risk of PTSD. Neither SLC6A4 genotype nor methylation status was associated with PTSD in main effects models. However, SLC6A4 methylation levels modified the effect of the number of traumatic events on PTSD after controlling for SLC6A4 genotype. Persons with more traumatic events were at increased risk for PTSD, but only at lower methylation levels. At higher methylation levels, individuals with more traumatic events were protected from this disorder. This interaction was observed whether the outcome was PTSD diagnosis, symptom severity, or number of symptoms.

CONCLUSIONS

Gene-specific methylation patterns may offer potential molecular signatures of increased risk for and resilience to PTSD.

The impact of cis-acting polymorphisms on the human phenotype

Cis-acting polymorphisms that affect gene expression are now known to be frequent, although the extent and mechanisms by which such variation affects the human phenotype are, as yet, only poorly understood. Key signatures of cis-acting variation are differences in gene expression that are tightly associated with regulatory SNPs or expression Quantitative Trait Loci (eQTL) and an imbalance of allelic expression (AEI) in heterozygous samples. Such cis-acting sequence differences appear often to have been under selection within and between populations and are also thought to be important in speciation. Here we describe the example of lactase persistence. In medical research, variants that affect regulation in cis have been implicated in both monogenic and polygenic disorders, and in the metabolism of drugs. In this review we suggest that by further understanding common regulatory variations and how they interact with other genetic and environmental variables it will be possible to gain insight into important mechanisms behind complex disease, with the potential to lead to new methods of diagnosis and treatments.

Global assessment of genetic variation influencing response to retinoid chemoprevention in head and neck cancer patients

Head and neck squamous cell carcinoma (HNSCC) patients are at an increased risk of developing a second primary tumor (SPT) or recurrence following curative treatment. 13-cis-retinoic acid (13-cRA) has been tested in chemoprevention clinical trials, but the results have been inconclusive. We genotyped 9,465 single nucleotide polymorphisms (SNP) in 450 patients from the Retinoid Head and Neck Second Primary Trial. SNPs were analyzed for associations with SPT/recurrence in patients receiving placebo to identify prognosis markers and further analyzed for effects of 13-cRA in patients with these prognostic loci. Thirteen loci identified a majority subgroup of patients at a high risk of SPT/recurrence and in whom 13-cRA was protective. Patients carrying the common genotype of rs3118570 in the retinoid X receptor (RXRA) were at a 3.33-fold increased risk (95% CI, 1.67-6.67) and represented more than 70% of the study population. This locus also identified individuals who received benefit from chemoprevention with a 38% reduced risk (95% CI, 0.43-0.90). Analyses of cumulative effect and potential gene-gene interactions also implicated CDC25C:rs6596428 and JAK2:rs1887427 as 2 other genetic loci with major roles in prognosis and 13-cRA response. Patients with all 3 common genotypes had a 76% reduction in SPT/recurrence (95% CI, 0.093-0.64) following 13-cRA chemoprevention. Carriers of these common genotypes constituted a substantial percentage of the study population, indicating that a pharmacogenetic approach could help select patients for 13-cRA chemoprevention. The lack of any alternatives for reducing risk in these patients highlights the need for future clinical trials to prospectively validate our findings.

Polymorphism of 3'UTR region of TNFR2 coding gene and its role in clinical tuberculosis in Han Chinese pediatric population

Genetic factors of human susceptibility to tuberculosis (TB) are multiple and their effect may be ethnic- and age-dependent. TNFR2 encoded by the TNFRSF1B gene is one of the important TNF-α receptors; its polymorphisms were previously suggested as potential markers of host susceptibility to TB. Here, genotyping of three SNPs in TNFRSF1B 3'UTR (rs1061624, rs5030792, rs3397) was performed in Han Chinese pediatric population (229 TB patients and 233 control subjects). rs5030792 was found homozygous (TT genotype) in all studied individuals. The rs3397-T allele was almost equally represented in both gender groups in this study; in particular, it was detected in 33.9% and 35.2% in female cases and controls, respectively (P=0.8). This latter result differs strikingly from an African study where rs3397-T was found in only 12.8 and 16.2% of Ghanaian female cases and controls, respectively (P=0.007 [Möller et al., 2010. Am. J. Respir. Crit. Care. Med. 181, 388-393]). In contrast, rs1061624-A allele, acting recessively, was found to be a possible risk factor for clinical TB in females (P=0.03). The rs1061624 heterozygotes were overdominant in controls versus patients (P=0.015) that warrants further study of their hypothetical advantage in TB. Neither of the common haplotypes was associated with susceptibility to TB. Compared to the published contrasting data on African (7-15%) and European (57%) populations, GTT haplotype was found in an intermediate frequency (26%). Further studies on both adult and pediatric populations in ethnically diverse settings are needed to elucidate the functionality of these 3'UTR SNPs of the TNFR2 gene.

Sequence variants in COL4A1 and COL4A2 genes in Ecuadorian families with keratoconus

PURPOSE

Keratoconus (KTCN) is a non-inflammatory, usually bilateral disorder of the eye which results in the conical shape and the progressive thinning of the cornea. Several studies have suggested that genetic factors play a role in the etiology of the disease. Several loci were previously described as possible candidate regions for familial KTCN; however, no causative mutations in any genes have been identified for any of these loci. The purpose of this study was to evaluate role of the collagen genes collagen type IV, alpha-1 (COL4A1) and collagen type IV, alpha-2 (COL4A2) in KTCN in Ecuadorian families.

METHODS

COL4A1 and COL4A2 in 15 Ecuadorian KTCN families were examined with polymerase chain reaction amplification, and direct sequencing of all exons, promoter and intron-exon junctions was performed.

RESULTS

Screening of COL4A1 and COL4A2 revealed numerous alterations in coding and non-coding regions of both genes. We detected three missense substitutions in COL4A1: c.19G>C (Val7Leu), c.1663A>C (Thr555Pro), and c.4002A>C (Gln1334His). Five non-synonymous variants were identified in COL4A2: c.574G>T (Val192Phe), c.1550G>A (Arg517Lys), c.2048G>C (Gly683Ala), c.2102A>G (Lys701Arg), and c.2152C>T (Pro718Ser). None of the identified sequence variants completely segregated with the affected phenotype. The Gln1334His variant was possibly damaging to protein function and structure.

CONCLUSIONS

This is the first mutation screening of COL4A1 and COL4A2 genes in families with KTCN and linkage to a locus close to these genes. Analysis of COL4A1 and COL4A2 revealed no mutations indicating that other genes are involved in KTCN causation in Ecuadorian families.

Predicting functionally important SNP classes based on negative selection

Background

With the advent of cost-effective genotyping technologies, genome-wide association studies allow researchers to examine hundreds of thousands of single nucleotide polymorphisms (SNPs) for association with human disease. Recently, many researchers applying this strategy have detected strong associations to disease with SNP markers that are either not in linkage disequilibrium with any nonsynonymous SNP or large distances from any annotated gene. In such cases, no well-established standard practice for effective SNP selection for follow-up studies exists. We aim to identify and prioritize groups of SNPs that are more likely to affect phenotypes in order to facilitate efficient SNP selection for follow-up studies.

Results

Based on the annotations available in the Ensembl database, we categorized SNPs in the human genome into classes related to regulatory attributes, such as epigenetic modifications and transcription factor binding sites, in addition to classes related to gene structure and cross-species conservation. Using the distribution of derived allele frequencies (DAF) within each class, we assessed the strength of natural selection for each class relative to the genome as a whole. We applied this DAF analysis to Perlegen resequenced SNPs genome-wide. Regulatory elements annotated by Ensembl such as specific histone methylation sites as well as classes defined by cross-species conservation showed negative selection in comparison to the genome as a whole.

Conclusions

These results highlight which annotated classes are under purifying selection, have putative functional importance, and contain SNPs that are strong candidates for follow-up studies after genome-wide association. Such SNP annotation may also be useful in interpreting results of whole-genome sequencing studies.

Mining Functional Elements in Messenger RNAs: Overview, Challenges, and Perspectives

Eukaryotic messenger RNA (mRNA) contains not only protein-coding regions but also a plethora of functional cis-elements that influence or coordinate a number of regulatory aspects of gene expression, such as mRNA stability, splicing forms, and translation rates. Understanding the rules that apply to each of these element types (e.g., whether the element is defined by primary or higher-order structure) allows for the discovery of novel mechanisms of gene expression as well as the design of transcripts with controlled expression. Bioinformatics plays a major role in creating databases and finding non-evident patterns governing each type of eukaryotic functional element. Much of what we currently know about mRNA regulatory elements in eukaryotes is derived from microorganism and animal systems, with the particularities of plant systems lagging behind. In this review, we provide a general introduction to the most well-known eukaryotic mRNA regulatory motifs (splicing regulatory elements, internal ribosome entry sites, iron-responsive elements, AU-rich elements, zipcodes, and polyadenylation signals) and describe available bioinformatics resources (databases and analysis tools) to analyze eukaryotic transcripts in search of functional elements, focusing on recent trends in bioinformatics methods and tool development. We also discuss future directions in the development of better computational tools based upon current knowledge of these functional elements. Improved computational tools would advance our understanding of the processes underlying gene regulations. We encourage plant bioinformaticians to turn their attention to this subject to help identify novel mechanisms of gene expression regulation using RNA motifs that have potentially evolved or diverged in plant species.

Context-dependent robustness to 5' splice site polymorphisms in human populations

There has been growing evidence for extensive diversity of alternative splicing in human populations. Genetic variants within the 5' splice site can cause splicing differences among human individuals and constitute an important class of human disease mutations. In this study, we explored whether natural variations of splicing could reveal important signals of 5' splice site recognition. In seven lymphoblastoid cell lines of Asian, European and African ancestry, we identified 1174 single nucleotide polymorphisms (SNPs) within the consensus 5' splice site. We selected 129 SNPs predicted to significantly alter the splice site activity, and quantitatively examined their splicing impact in the seven individuals. Surprisingly, outside of the essential GT dinucleotide position, only ∼14% of the tested SNPs altered splicing. Bioinformatic and minigene analyses identified signals that could modify the impact of 5' splice site polymorphisms, most notably a strong 3' splice site and the presence of intronic motifs downstream of the 5' splice site. Strikingly, we found that the poly-G run, a known intronic splicing enhancer, was the most significantly enriched motif downstream of exons unaffected by 5' splice site SNPs. In TRIM62, the upstream 3' splice site and downstream intronic poly-G runs functioned redundantly to protect an exon from its 5' splice site polymorphism. Collectively, our study reveals widespread context-dependent robustness to 5' splice site polymorphisms in human transcriptomes. Consequently, certain exons are more susceptible to 5' splice site mutations. Additionally, our work demonstrates that genetic diversity of alternative splicing can provide significant insights into the splicing code of mammalian cells.

RNA editing of nuclear transcripts in Arabidopsis thaliana

BACKGROUND

RNA editing is a transcript-based layer of gene regulation. To date, no systemic study on RNA editing of plant nuclear genes has been reported. Here, a transcriptome-wide search for editing sites in nuclear transcripts of Arabidopsis (Arabidopsis thaliana) was performed.

RESULTS

MPSS (massively parallel signature sequencing) and PARE (parallel analysis of RNA ends) data retrieved from public databases were utilized, focusing on one-base-conversion editing. Besides cytidine (C)-to-uridine (U) editing in mitochondrial transcripts, many nuclear transcripts were found to be diversely edited. Interestingly, a sizable portion of these nuclear genes are involved in chloroplast- or mitochondrion-related functions, and many editing events are tissue-specific. Some editing sites, such as adenosine (A)-to-U editing loci, were found to be surrounded by peculiar elements. The editing events of some nuclear transcripts are highly enriched surrounding the borders between coding sequences (CDSs) and 3' untranslated regions (UTRs), suggesting site-specific editing. Furthermore, RNA editing is potentially implicated in new start or stop codon generation, and may affect alternative splicing of certain protein-coding transcripts. RNA editing in the precursor microRNAs (pre-miRNAs) of ath-miR854 family, resulting in secondary structure transformation, implies its potential role in microRNA (miRNA) maturation.

CONCLUSIONS

To our knowledge, the results provide the first global view of RNA editing in plant nuclear transcripts.

Extending rare-variant testing strategies: analysis of noncoding sequence and imputed genotypes

Next Generation Sequencing Technology has revolutionized our ability to study the contribution of rare genetic variation to heritable traits. However, existing single-marker association tests are underpowered for detecting rare risk variants. A more powerful approach involves pooling methods that combine multiple rare variants from the same gene into a single test statistic. Proposed pooling methods can be limited because they generally assume high-quality genotypes derived from deep-coverage sequencing, which may not be available. In this paper, we consider an intuitive and computationally efficient pooling statistic, the cumulative minor-allele test (CMAT). We assess the performance of the CMAT and other pooling methods on datasets simulated with population genetic models to contain realistic levels of neutral variation. We consider study designs ranging from exon-only to whole-gene analyses that contain noncoding variants. For all study designs, the CMAT achieves power comparable to that of previously proposed methods. We then extend the CMAT to probabilistic genotypes and describe application to low-coverage sequencing and imputation data. We show that augmenting sequence data with imputed samples is a practical method for increasing the power of rare-variant studies. We also provide a method of controlling for confounding variables such as population stratification. Finally, we demonstrate that our method makes it possible to use external imputation templates to analyze rare variants imputed into existing GWAS datasets. As proof of principle, we performed a CMAT analysis of more than 8 million SNPs that we imputed into the GAIN psoriasis dataset by using haplotypes from the 1000 Genomes Project.

Late onset sporadic dilated cardiomyopathy caused by a cardiac troponin T mutation

Mutations in TNNT2, encoding cardiac troponin T, commonly shows early onset, aggressive dilated cardiomyopathy (DCM). This observation may influence the decision of whether to undertake clinical genetic testing for TNNT2 in later onset DCM. Further, the trigger for late onset DCM remains enigmatic. A 70-year-old woman, previously healthy with a left ventricular ejection fraction of 50%-55% at age 69, presented with DCM of unknown cause and a 4-month history progressive heart failure requiring cardiac transplantation. Clinical genetic testing revealed a novel TNNT2 R139H mutation but no relevant variants in 18 other DCM genes. Her explanted heart showed partial fatty replacement in the right ventricle. Sequencing for five arrhythmogenic right ventricular dysplasia genes was negative. Functional studies in porcine cardiac skinned fibers reconstituted with the mutant R139H troponin T protein showed decreased Ca(2+) sensitivity at pH 7, characteristic of DCM. Because fatty infiltration may acidify the myocellular environment, maximal force development examined at pH 6.5 was diminished, suggesting a possible environmental trigger. We conclude that the TNNT2 R139H mutation was likely to be disease causing. Further, later age of onset may not be relevant to exclude genetic testing for TNNT2 mutations.

Genes, mutations, and human inherited disease at the dawn of the age of personalized genomics

The number of reported germline mutations in human nuclear genes, either underlying or associated with inherited disease, has now exceeded 100,000 in more than 3,700 different genes. The availability of these data has both revolutionized the study of the morbid anatomy of the human genome and facilitated "personalized genomics." With approximately 300 new "inherited disease genes" (and approximately 10,000 new mutations) being identified annually, it is pertinent to ask how many "inherited disease genes" there are in the human genome, how many mutations reside within them, and where such lesions are likely to be located? To address these questions, it is necessary not only to reconsider how we define human genes but also to explore notions of gene "essentiality" and "dispensability."Answers to these questions are now emerging from recent novel insights into genome structure and function and through complete genome sequence information derived from multiple individual human genomes. However, a change in focus toward screening functional genomic elements as opposed to genes sensu stricto will be required if we are to capitalize fully on recent technical and conceptual advances and identify new types of disease-associated mutation within noncoding regions remote from the genes whose function they disrupt.

The two stem cell microRNA gene clusters C19MC and miR-371-3 are activated by specific chromosomal rearrangements in a subgroup of thyroid adenomas

Thyroid adenomas are common benign human tumors with a high prevalence of about 5% of the adult population even in iodine sufficient areas. Rearrangements of chromosomal band 19q13.4 represent a frequent clonal cytogenetic deviation in these tumors making them the most frequent non-random chromosomal translocations in human epithelial tumors at all. Two microRNA (miRNA) gene clusters i.e. C19MC and miR-371-3 are located in close proximity to the breakpoint region of these chromosomal rearrangements and have been checked for a possible up-regulation due to the genomic alteration. In 4/5 cell lines established from thyroid adenomas with 19q13.4 rearrangements and 5/5 primary adenomas with that type of rearrangement both the C19MC and miR-371-3 cluster were found to be significantly overexpressed compared to controls lacking that particular chromosome abnormality. In the remaining cell line qRT-PCR revealed overexpression of members of the miR-371-3 cluster only which might be due to a deletion accompanying the chromosomal rearrangement in that case. In depth molecular characterization of the breakpoint in a cell line from one adenoma of this type reveals the existence of large Pol-II mRNA fragments as the most likely source of up-regulation of the C19MC cluster. The up-regulation of the clusters is likely to be causally associated with the pathogenesis of the corresponding tumors. Of note, the expression of miRNAs miR-520c and miR-373 is known to characterize stem cells and in terms of molecular oncology has been implicated in invasive growth of epithelial cells in vitro and in vivo thus allowing to delineate a distinct molecular subtype of thyroid adenomas. Besides thyroid adenomas rearrangements of 19q13.4 are frequently found in other human neoplasias as well, suggesting that activation of both clusters might be a more general phenomenon in human neoplasias.