Non-coding Regions Of Genome Research Articles

Investigation of the tick-borne encephalitis virus in Norway

Expansion of the range of Ixodes ricinus ticks further north and to higher altitudes affects the spread of tick-borne encephalitis virus (TBEV) in new territories and increases the risk of human infections. Over the past decade, the number of tick-borne encephalitis (TBE) cases has increased in Norway and parts of the southern coast have been established as endemic. The aims of this study were to confirm the existence of TBEV endemic foci, determine the spread of TBE in different localities of Norway, review and compare relevant published records on the prevalence of TBEV in the different developmental stages and sex of the ticks. Ticks were collected from nine locations along the southern coast of Norway during June-July 2009. For the detection of TBEV-specific RNA, a real-time RT-PCR targeting a part of the 3’ non-coding region of the TBEV genome was used. A total of 1542 ticks grouped in 250 pools were examined and TBEV RNA was detected in 5 pools with overall prevalence of 0.32%.

Abstract 1221: Chronic lymphocytic leukemia risk variants are associated with alteration of epigenetic states in regulatory regions

Abstract Chronic lymphocytic leukemia risk variants are associated with alteration of epigenetic state in regulatory regions The vast majority of disease-associated single-nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWASs) reside in noncoding genomic regions. Functional impact of noncoding casual variants on target gene expression often involves regulatory mechanisms, typically through altering local chromatin accessibility, histone modifications, transcription factor (TF) binding affinity and chromatin interactions. In chronic lymphocytic leukemia (CLL), we and others have identified 41 GWAS risk variants, whose mechanism of action in CLL etiology remains largely unknown. In this study, we aim to decipher the regulatory potential of these variants and to better understand the mechanisms contributing to dysregulation of genes implicated in CLL pathogenesis. We performed ChIP-seq for H3K4me3, H3K4me1, H3K27ac and H3K27me3 in CLL. In addition, we also analyzed over 250 public ChIP-seq, ATAC-seq, DNase-seq, FAIRE-seq, and chromatin interaction data in CLL, B cells and/or lymphoblastoid cell lines. Comparing the 41 index SNPs and the nearby SNPs they associated with in linkage disequilibrium (R2 at least 0.5 in the EUR ethnic group in 1000 Genomes Project) to the above epigenetic data sets revealed overlaps with regulatory regions in 38 of the risk loci, with 28 of the loci overlapping super-enhancers. We used two approaches to elucidate the mechanisms of potential SNP causality. First, we used GATK UnifiedGenotyper to infer genotype for the SNPs that overlapped open chromatin regions and ChIP-seq peaks, by which evidence of imbalance between reference and alternative allele was identified at heterozygous sites. Second, for SNPs located within open chromatin regions or ChIP-seq peaks, we used FIMO program to scan the 100bp sequences centered on the SNPs for significant matches (p=1e-05) to the position weight matrices in five TF binding motif databases; SNPs that had TF motifs altered by the alternative alleles were retained. Together, these two approaches revealed evidence for allele imbalance or TF binding motif alterations at 36 of the 38 CLL risk loci, suggesting the generality of epigenetic regulation in CLL pathogenesis. Further analysis of chromatin interaction, gene expression and expression quantitative trait loci data will likely link the casual SNPs to their target genes. Citation Format: Huihuang Yan, Shulan Tian, Geffen Kleinstern, Susan Slager. Chronic lymphocytic leukemia risk variants are associated with alteration of epigenetic states in regulatory regions [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 1221.

Abstract 2281: Flexible discovery of recurrent coding and non-coding mutations in tumor whole genomes

Abstract Analysis of somatic protein-coding mutations, particularly via tumor exome sequencing, has driven discovery of novel cancer genes in recent years. However, ~97% of the human genome is non-coding and not addressable by exome-targeting methods. Tumor whole genome sequencing, by contrast, allows for genome-wide discovery of somatic events and may identify novel non-coding driver mutations. In conjunction, flexible tools employing principled statistical methods are needed for discovery of recurrent somatic driver mutations from tumor whole genomes. Here, we developed a computational toolset that addresses these needs and describe its application to a novel cohort of lung adenocarcinoma tumors generated by the APOLLO (Applied Proteogenomics OrganizationaL Learning and Outcomes) consortium. Our toolset flexibly computes feature-wise and hotspot somatic mutation enrichment statistics in both protein-coding and non-coding (e.g. promoters, enhancers) genomic regions. Additionally, we incorporated an affinity propagation-based clustering procedure that groups test regions by user-specified genomic covariates that can influence local mutation properties. We applied our toolset to several TCGA datasets of tumor exome sequencing. We found strong agreement between our results and those of other tools: among the top 15 genes called at FDR < 0.01 by MutSigCV, we observed a median overlap of 75.7% with genes ranked in the top 15 by our methods. We next applied our toolset to tumor whole genome sequencing data from our APOLLO lung adenocarcinoma cohort. During analysis, we simultaneously controlled for several genomic covariates, including GC/CpG content, replication timing, and proximal mRNA expression. Gene-centric analysis identified several significantly mutated genes, including KRAS, STK11, and TP53 (FDR < 0.005). We examined matched mRNA-seq data from this cohort and found significant over-expression of KRAS mRNA in samples possessing a non-silent coding mutation in this gene (p < 0.02). We additionally assessed somatic mutation enrichment in promoters, 3' UTRs, non-coding genes, and lung-specific enhancers. From these results, we found significant enrichment of somatic mutations in the body of the long non-coding RNA NEAT1. Among the tumors with somatic NEAT1 mutations, a single sample was hypermutated at this locus and corresponded with the lowest observed RNA expression for this lncRNA in our cohort (TPM of 12.2 vs. overall mean TPM of 84.2). In conclusion, we developed a flexible toolset for interrogating both coding and non-coding landscapes from tumor whole genomes. We applied our methods to a novel cohort of lung adenocarcinoma tumors and identified recurrently mutated genes and non-coding regions, including the NEAT1 lncRNA. Citation Format: Anthony R. Soltis, Coralie Viollet, Harvey B. Pollard, Christopher A. Moskaluk, Robert F. Browning, Clifton L. Dalgard, Craig D. Shriver, Matthew D. Wilkerson. Flexible discovery of recurrent coding and non-coding mutations in tumor whole genomes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2281.

The use of chloroplast genome sequences to solve phylogenetic incongruences in Polystachya Hook (Orchidaceae Juss).

BackgroundCurrent evidence suggests that for more robust estimates of species tree and divergence times, several unlinked genes are required. However, most phylogenetic trees for non-model organisms are based on single sequences or just a few regions, using traditional sequencing methods. Techniques for massive parallel sequencing or next generation sequencing (NGS) are an alternative to traditional methods that allow access to hundreds of DNA regions. Here we use this approach to resolve the phylogenetic incongruence found in Polystachya Hook. (Orchidaceae), a genus that stands out due to several interesting aspects, including cytological (polyploid and diploid species), evolutionary (reticulate evolution) and biogeographical (species widely distributed in the tropics and high endemism in Brazil). The genus has a notoriously complicated taxonomy, with several sections that are widely used but probably not monophyletic.MethodsWe generated the complete plastid genome of 40 individuals from one clade within the genus. The method consisted in construction of genomic libraries, hybridization to RNA probes designed from available sequences of a related species, and subsequent sequencing of the product. We also tested how well a smaller sample of the plastid genome would perform in phylogenetic inference in two ways: by duplicating a fast region and analyzing multiple copies of this dataset, and by sampling without replacement from all non-coding regions in our alignment. We further examined the phylogenetic implications of non-coding sequences that appear to have undergone hairpin inversions (reverse complemented sequences associated with small loops).ResultsWe retrieved 131,214 bp, including coding and non-coding regions of the plastid genome. The phylogeny was able to fully resolve the relationships among all species in the targeted clade with high support values. The first divergent species are represented by African accessions and the most recent ones are among Neotropical species.DiscussionOur results indicate that using the entire plastid genome is a better option than screening highly variable markers, especially when the expected tree is likely to contain many short branches. The phylogeny inferred is consistent with the proposed origin of the genus, showing a probable origin in Africa, with later dispersal into the Neotropics, as evidenced by a clade containing all Neotropical individuals. The multiple positions of Polystachya concreta (Jacq.) Garay & Sweet in the phylogeny are explained by allotetraploidy. Polystachya estrellensis Rchb.f. can be considered a genetically distinct species from P. concreta and P. foliosa (Lindl.) Rchb.f., but the delimitation of P. concreta remains uncertain. Our study shows that NGS provides a powerful tool for inferring relationships at low taxonomic levels, even in taxonomically challenging groups with short branches and intricate morphology.

ALUminating the Path of Atherosclerosis Progression: Chaos Theory Suggests a Role for Alu Repeats in the Development of Atherosclerotic Vascular Disease.

Atherosclerosis (ATH) and coronary artery disease (CAD) are chronic inflammatory diseases with an important genetic background; they derive from the cumulative effect of multiple common risk alleles, most of which are located in genomic noncoding regions. These complex diseases behave as nonlinear dynamical systems that show a high dependence on their initial conditions; thus, long-term predictions of disease progression are unreliable. One likely possibility is that the nonlinear nature of ATH could be dependent on nonlinear correlations in the structure of the human genome. In this review, we show how chaos theory analysis has highlighted genomic regions that have shared specific structural constraints, which could have a role in ATH progression. These regions were shown to be enriched with repetitive sequences of the Alu family, genomic parasites that have colonized the human genome, which show a particular secondary structure and are involved in the regulation of gene expression. Here, we show the impact of Alu elements on the mechanisms that regulate gene expression, especially highlighting the molecular mechanisms via which the Alu elements alter the inflammatory response. We devote special attention to their relationship with the long noncoding RNA (lncRNA); antisense noncoding RNA in the INK4 locus (ANRIL), a risk factor for ATH; their role as microRNA (miRNA) sponges; and their ability to interfere with the regulatory circuitry of the (nuclear factor kappa B) NF-κB response. We aim to characterize ATH as a nonlinear dynamic system, in which small initial alterations in the expression of a number of repetitive elements are somehow amplified to reach phenotypic significance.

Integration of hepatitis B virus DNA in chronically infected patients assessed by Alu-PCR.

Hepatitis B virus (HBV) infection is the main risk factor for hepatocellular carcinoma (HCC) worldwide. Integration of HBV DNA into the human genome has been found in >80% of HBV-related HCC cases. Some studies have, however, found similar integration patterns in tumorous and nontumorous tissues. Thus, the role of integrations for the development of HCC as well as the rate of integration in different stages of infection remain unclear. The aim of this study was to investigate integrations in patients without HCC, representing different stages of chronic HBV (CHB) infection. Extracted DNA in liver biopsies from 74 patients (one with 2 available biopsies) with CHB infection was analyzed by Alu-PCR. Amplicons were further analyzed by Sanger sequencing. Integration was detected in 39 biopsies (52%) as an amplicon containing both human and HBV sequences by Alu-PCR with one primer targeting a region in the HBV genome. Integrations were found in patients representing the different stages of CHB infection. A majority of the HBV sequences were located upstream or downstream of nucleotide position 1820, which previously has been identified as a common breakpoint in the HBV genome in integrated sequences. Approximately 60% of the HBV integrations were found in noncoding regions of the human genome. Integrations of HBV DNA into the human genome is an event frequentlyfound in mild phases of chronic hepatitis.

Regulatory variants: from detection to predicting impact.

Variants within non-coding genomic regions can greatly affect disease. In recent years, increasing focus has been given to these variants, and how they can alter regulatory elements, such as enhancers, transcription factor binding sites and DNA methylation regions. Such variants can be considered regulatory variants. Concurrently, much effort has been put into establishing international consortia to undertake large projects aimed at discovering regulatory elements in different tissues, cell lines and organisms, and probing the effects of genetic variants on regulation by measuring gene expression. Here, we describe methods and techniques for discovering disease-associated non-coding variants using sequencing technologies. We then explain the computational procedures that can be used for annotating these variants using the information from the aforementioned projects, and prediction of their putative effects, including potential pathogenicity, based on rule-based and machine learning approaches. We provide the details of techniques to validate these predictions, by mapping chromatin–chromatin and chromatin–protein interactions, and introduce Clustered Regularly Interspaced Short Palindromic Repeats-Associated Protein 9 (CRISPR-Cas9) technology, which has already been used in this field and is likely to have a big impact on its future evolution.We also give examples of regulatory variants associated with multiple complex diseases. This review is aimed at bioinformaticians interested in the characterization of regulatory variants, molecular biologists and geneticists interested in understanding more about the nature and potential role of such variants from a functional point of views, and clinicians who may wish to learn about variants in non-coding genomic regions associated with a given disease and find out what to do next to uncover how they impact on the underlying mechanisms.

Synthetic RNA derived from the foot-and-mouth disease virus genome elicits antiviral responses in bovine and porcine cells through IRF3 activation

Foot-and-mouth disease virus (FMDV) is the causative agent of a highly transmissible disease affecting wild and domestic animals including pigs, cattle and sheep. The ability of synthetic RNA transcripts mimicking distinct domains in the non-coding regions of the FMDV genome (ncRNAs) to induce a potent innate immune response in swine cultured cells and mice has been previously described, as well as their enhancing effect on conventional inactivated FMD vaccines. Here, we provide evidence of the activation of interferon regulatory factor 3 (IRF3), a key transcriptional regulator of type I interferon (IFN)-dependent immune responses after transfection of swine and bovine cells with transcripts corresponding to the FMDV 3´ non-coding region (3´NCR). Induction of IFN-β and Mx1expression, concomitantly with antiviral activity and IRF3 activation was observed in bovine MDBK cells transfected with the 3´NCR. Our results link the stimulation of the innate immune response observed in 3´NCR-transfected cells to the intracellular type I IFN signaling pathway and suggest the potential use of these molecules for antiviral strategies in cattle.

Integrative Analysis Identifies Genetic Variants Associated With Autoimmune Diseases Affecting Putative MicroRNA Binding Sites.

Genome-wide and fine mapping studies have shown that more than 90% of genetic variants associated with autoimmune diseases (AID) are located in non-coding regions of the human genome and especially in regulatory sequences, including microRNAs (miRNA) target sites. MiRNAs are small endogenous noncoding RNAs that modulate gene expression at the post-transcriptional level. Single nucleotide polymorphisms (SNPs) located within the 3′ untranslated region of their target mRNAs (miRSNP) can alter miRNA binding sites. Yet, little is known about their effect on regulation by miRNA and the consequences for AID. Conversely, it is well known that two or more AID may share part of their genetic background. Here, we hypothesized that miRSNPs could be associated with more than one AID. To identify miRSNPs associated with AID, we integrated results from three different prediction tools (Polymirts, miRSNP, and miRSNPscore) using a naïve Bayes classifier approach to identify miRSNPs predicted to affect binding sites of miRNAs. Further, to detect miRSNPs associated with two or more AID, we integrated predictions with summary statistics from 12 AID studies. In addition, to prioritize miRSNPs, miRNAs and AID-associated target genes, we used public expression quantitative trait locus (eQTL) data and mRNA-seq and small RNA-seq data. We identified 34 miRNSPs associated with at least two AID. Furthermore, we found 86 miRNAs predicted to target 18 of the associated gene's mRNAs. Our integrative approach revealed new insights into miRNAs and AID associated target genes. Thus, it helped to prioritize AID noncoding risk SNPs that might be involved in the causal mechanisms, providing valuable information for further functional studies.

Mutation hotspots at CTCF binding sites coupled to chromosomal instability in gastrointestinal cancers

Tissue-specific driver mutations in non-coding genomic regions remain undefined for most cancer types. Here, we unbiasedly analyze 212 gastric cancer (GC) whole genomes to identify recurrently mutated non-coding regions in GC. Applying comprehensive statistical approaches to accurately model background mutational processes, we observe significant enrichment of non-coding indels (insertions/deletions) in three gastric lineage-specific genes. We further identify 34 mutation hotspots, of which 11 overlap CTCF binding sites (CBSs). These CBS hotspots remain significant even after controlling for a genome-wide elevated mutation rate at CBSs. In 3 out of 4 tested CBS hotspots, mutations are nominally associated with expression change of neighboring genes. CBS hotspot mutations are enriched in tumors showing chromosomal instability, co-occur with neighboring chromosomal aberrations, and are common in gastric (25%) and colorectal (19%) tumors but rare in other cancer types. Mutational disruption of specific CBSs may thus represent a tissue-specific mechanism of tumorigenesis conserved across gastrointestinal cancers.

Emerging applications of genome-editing technology to examine functionality of GWAS-associated variants for complex traits.

Over the last decade, genome-wide association studies (GWAS) have propelled the discovery of thousands of loci associated with complex diseases. The focus is now turning toward the function of these association signals, determining the causal variant(s) among those in strong linkage disequilibrium, and identifying their underlying mechanisms, such as long-range gene regulation. Genome-editing techniques utilizing zinc-finger nucleases (ZFN), transcription activator-like effector nucleases (TALENs), and clustered regularly-interspaced short palindromic repeats with Cas9 nuclease (CRISPR-Cas9) are becoming the tools of choice to establish functionality for these variants, due to the ability to assess effects of single variants in vivo. This review will discuss examples of how these technologies have begun to aid functional analysis of GWAS loci for complex traits such as cardiovascular disease, Type 2 diabetes, cancer, obesity, and autoimmune disease. We focus on analysis of variants occurring within noncoding genomic regions, as these comprise the majority of GWAS variants, providing the greatest challenges to determining functionality, and compare editing strategies that provide different levels of evidence for variant functionality. The review describes molecular insights into some of these potentially causal variants and how these may relate to the pathology of the trait and look toward future directions for these technologies in post-GWAS analysis, such as base-editing.

Genomic non-coding regions reveal hidden patterns of mumps virus circulation in Spain, 2005 to 2015.

BackgroundSince mumps vaccination was introduced in 1981 in Spain, the incidence of the disease has dropped significantly. However, cyclic epidemic waves and outbreaks still occur, despite high vaccination coverage. The World Health Organization (WHO) recommends genotyping to trace the pattern of mumps virus (MuV) circulation. Genotype H was predominant in Spain, but was replaced in 2005 by genotype G which has subsequently remained dominant. Of the small hydrophobic protein gene sequences, 78% are identical and belong to the MuVi/ Sheffield.GBR.1.05/[G]-variant. Aim: Our study aimed to investigate whether the circulation of MuV strains in Spain was continuous after the emergence of genotype G in 2005. Method: We obtained 46 samples from Spanish patients infected with MuVi/Sheffield.GBR.1.05/[G] during two epidemic waves and analysed them using new molecular markers based on genomic non-coding regions (NCRs) that discriminate subvariants of this virus strain. Results: Phylogenetic analyses of the nucleoprotein–phosphoprotein and matrix protein–fusion protein NCR indicated strain replacement after a drop in incidence in 2009, which had not been detectable by SH sequencing. Clustering of sequences from patients epidemiologically linked in the same outbreak suggests a potential use for these NCRs in outbreak characterisation. Conclusion: We suggest to consider their use in conjunction with the SH gene in the future WHO recommendations for MuV epidemiological surveillance.

High levels of genetic diversity in marginal populations of the marsh orchid Dactylorhiza majalis subsp. majalis

Because of harsh conditions, suboptimal habitat quality and poor connectivity to other populations, plant populations at the margin of a distribution are expected to be less genetically diverse, but to be more divergent from each other than populations in the centre of a distribution. In northern Europe, northern marginal populations may also be younger than populations further to the south, and may have had less time to accumulate genetic diversity by mutation and gene flow. However, orchids have very small seeds, which are easily dispersed long distances by wind, and orchids are therefore expected to show less differentiation between marginal and central populations than other groups of seed plants. Here, we analysed whether Scandinavian populations of the tetraploid marsh orchid Dactylorhiza majalis subsp. majalis differ from central European populations in genetic diversity patterns. A total of 220 plants from eight central European and ten Scandinavian populations was examined for variation at five nuclear microsatellite loci, nuclear ITS and 13 polymorphic sites in noncoding regions of the plastid genome. The total genetic diversity was slightly lower in Scandinavia than in central Europe, both in plastid and nuclear markers, but the differences were small. Also, the Scandinavian populations were less diverse and somewhat more strongly differentiated from each other than the central European ones. Dactylorhiza majalis subsp. majalis has apparently colonized Scandinavia on multiple independent occasions and from different source areas in the south. Seed flow between Scandinavian populations has still not fully erased the patterns imprinted by early colonization. Our results suggest that marginal populations of orchids may be as important as central ones in preserving genetic diversity through Pleistocene glacial cycles. We also predict that orchids with their light seeds are better adapted than many other plants to respond to future climate changes by dispersing into new suitable areas.

17th Chromosome-Centric Human Proteome Project Symposium in Tehran.

This report describes the 17th Chromosome-Centric Human Proteome Project which was held in Tehran, Iran, April 27 and 28, 2017. A brief summary of the symposium's talks including new technical and computational approaches for the identification of novel proteins from non-coding genomic regions, physicochemical and biological causes of missing proteins, and the close interactions between Chromosome- and Biology/Disease-driven Human Proteome Project are presented. A synopsis of decisions made on the prospective programs to maintain collaborative works, share resources and information, and establishment of a newly organized working group, the task force for missing protein analysis are discussed.

Influence of Electron-Holes on DNA Sequence-Specific Mutation Rates.

Biases in mutation rate can influence molecular evolution, yielding rates of evolution that vary widely in different parts of the genome and even among neighboring nucleotides. Here, we explore one possible mechanism of influence on sequence-specific mutation rates, the electron–hole, which can localize and potentially trigger a replication mismatch. A hole is a mobile site of positive charge created during one-electron oxidation by, for example, radiation, contact with a mutagenic agent, or oxidative stress. Its quantum wavelike properties cause it to localize at various sites with probabilities that vary widely, by orders of magnitude, and depend strongly on the local sequence. We find significant correlations between hole probabilities and mutation rates within base triplets, observed in published mutation accumulation experiments on four species of bacteria. We have also computed hole probability spectra for hypervariable segment I of the human mtDNA control region, which contains several mutational hotspots, and for heptanucleotides in noncoding regions of the human genome, whose polymorphism levels have recently been reported. We observe significant correlations between hole probabilities, and context-specific mutation and substitution rates. The correlation with hole probability cannot be explained entirely by CpG methylation in the heptanucleotide data. Peaks in hole probability tend to coincide with mutational hotspots, even in mtDNA where CpG methylation is rare. Our results suggest that hole-enhanced mutational mechanisms, such as oxidation-stabilized tautomerization and base deamination, contribute to molecular evolution.

Novel GREM1 Variations in Sub-Saharan African Patients With Cleft Lip and/or Cleft Palate.

Cleft lip and/or cleft palate (CL/P) are congenital anomalies of the face and have multifactorial etiology, with both environmental and genetic risk factors playing crucial roles. Though at least 40 loci have attained genomewide significant association with nonsyndromic CL/P, these loci largely reside in noncoding regions of the human genome, and subsequent resequencing studies of neighboring candidate genes have revealed only a limited number of etiologic coding variants. The present study was conducted to identify etiologic coding variants in GREM1, a locus that has been shown to be largely associated with cleft of both lip and soft palate. We resequenced DNA from 397 sub-Saharan Africans with CL/P and 192 controls using Sanger sequencing. Following analyses of the sequence data, we observed 2 novel coding variants in GREM1. These variants were not found in the 192 African controls and have never been previously reported in any public genetic variant database that includes more than 5000 combined African and African American controls or from the CL/P literature. The novel variants include p.Pro164Ser in an individual with soft palate cleft only and p.Gly61Asp in an individual with bilateral cleft lip and palate. The proband with the p.Gly61Asp GREM1 variant is a van der Woude (VWS) case who also has an etiologic variant in IRF6 gene. Our study demonstrated that there is low number of etiologic coding variants in GREM1, confirming earlier suggestions that variants in regulatory elements may largely account for the association between this locus and CL/P.

Selective constraint and adaptive potential of West Nile virus within and among naturally infected avian hosts and mosquito vectors.

Arthropod-borne viruses are among the most genetically constrained RNA viruses, yet they have a remarkable propensity to adapt and emerge. We studied wild birds and mosquitoes naturally infected with West Nile virus (WNV) in a ‘hot spot’ of virus transmission in Chicago, IL, USA. We generated full coding WNV genome sequences from spatiotemporally matched bird and mosquito samples using high-throughput sequencing, allowing a molecular evolutionary assessment with deep coverage. Mean FST among samples was 0.66 (±0.02 SE) and was bimodal, with mean nucleotide diversity being higher between samples (interhost πN = 0.001; πS = 0.024) than within them (intrahost πN < 0.0001; πS < 0.001). Eight genomic sites with FST > 1.01 (in the PrM, NS2a, NS3, NS4b, and 5'-noncoding genomic regions) showed bird versus mosquito variant frequency differences of >30 per cent and/or polymorphisms fixed in ≥5 host or vector individuals, suggesting host tropism for these variants. However, phylogenetic analyses demonstrated a lack of grouping by bird or mosquito, most inter-sample differences were synonymous (mean interhost πN/πS = 0.04), and there was no significant difference between hosts and vectors in either their nucleotide diversities or levels of purifying selection (mean intrahost πN/πS = 0.28 in birds and πN/πS = 0.21 in mosquitoes). This finding contrasts with the ‘trade-off’ and ‘selective sieve’ hypotheses that have been proposed and tested in the laboratory, which predict strong host versus vector effects on WNV genetic variation, with heightened selective constraint in birds alternating with heightened viral diversity in mosquitoes. Overall, our data show WNV to be highly selectively constrained within and between both hosts and vectors but still able to vary at a limited number of sites across the genome. Such site-specific plasticity in the face of overall selective constraint may offer a mechanism whereby highly constrained viruses such as WNV and its relatives can still adapt and emerge.

Mining and analysis of simple sequence repeats in the chloroplast genomes of genus Vigna

The genus Vigna is distributed throughout the world and it is phylogenetically closely related to agriculturally important crops. Simple sequence repeats (SSRs) are defined as sequence repeat units that occur between 1 and 6bp and found in both coding and non-coding regions of genome. The availability of organelle genome sequence helps in understanding the organization of SSRs in them. In this study, chloroplast genome sequences of the genus Vigna were screened for the identification of chloroplast simple sequence repeats (cpSSRs). A total of 21, 24 and 25 cpSSRs in Vigna angularis, Vigna radiata and Vigna unguiculata, respectively, were mined computationally. The density of ∼1SSR/6.1–7.2kb was observed. Depending on the different repeat units, the length of SSRs ranged from 12 to 24 bp. Mononucleotides were the most abundant repeats followed by dinucleotide repeats. Hexanucleotide repeats were completely absent in the chloroplast genomes of genus Vigna. Seven cpSSRs were selected on the basis of repeats to study transferability. The transferable markers can be employed in genome mapping and comparative mapping among Vigna species in future studies.

Invention and Early History of Exon Skipping and Splice Modulation.

Since its discovery in 1977, much has been known about RNA splicing and how it plays a central role in human development, function, and, notably, disease. Defects in RNA splicing account for at least 10% of all genetic disorders, with the number expected to increase as more information is uncovered on the contribution of noncoding genomic regions to disease. Splice modulation through the use of antisense oligonucleotides (AOs) has emerged as a promising avenue for the treatment of these disorders. In fact, two splice-switching AOs have recently obtained approval from the US Food and Drug Administration: eteplirsen (Exondys 51) for Duchenne muscular dystrophy, and nusinersen (Spinraza) for spinal muscular atrophy. These work by exon skipping and exon inclusion, respectively. In this chapter, we discuss the early development of AO-based splice modulation therapy-its invention, first applications, and its evolution into the approach we are now familiar with. We give a more extensive history of exon skipping in particular, as it is the splice modulation approach given the most focus in this book.

An integrative approach to predicting the functional effects of small indels in non-coding regions of the human genome

BackgroundSmall insertions and deletions (indels) have a significant influence in human disease and, in terms of frequency, they are second only to single nucleotide variants as pathogenic mutations. As the majority of mutations associated with complex traits are located outside the exome, it is crucial to investigate the potential pathogenic impact of indels in non-coding regions of the human genome.ResultsWe present FATHMM-indel, an integrative approach to predict the functional effect, pathogenic or neutral, of indels in non-coding regions of the human genome. Our method exploits various genomic annotations in addition to sequence data. When validated on benchmark data, FATHMM-indel significantly outperforms CADD and GAVIN, state of the art models in assessing the pathogenic impact of non-coding variants. FATHMM-indel is available via a web server at indels.biocompute.org.uk.ConclusionsFATHMM-indel can accurately predict the functional impact and prioritise small indels throughout the whole non-coding genome.