Análise do grau de conservação de resíduos em proteínas com estrutura 3D resolvida utilizando o SMS.

HSSP e entropia relativa. Módulos do SMS para análise de conservação. Discussão e trabalhos futuros.

Analysis of the elements of catabolite repression in Clostridium acetobutylicum.

The ptsH gene, encoding the phosphotransferase protein HPr, from Clostridium acetobutylicum ATCC 824 was identified from the genome sequence, cloned and shown to complement a ptsH mutant of Escherichia coli. The deduced protein sequence shares significant homology with HPr proteins from other low-GC gram-positive bacteria, although the highly conserved sequence surrounding the Ser-46 phosphorylation site is not well preserved in the clostridial protein. Nevertheless, the HPr was phosphorylated in an ATP-dependent manner in cell-free extracts of C. acetobutylicum. Furthermore, purified His-tagged HPr from Bacillus subtilis was also a substrate for the clostridial HPr kinase/phosphorylase. This phosphorylation reaction is a key step in the mechanism of carbon catabolite repression proposed to operate in B. subtilis and other low-GC gram-positive bacteria. Putative genes encoding the HPr kinase/phosphorylase and the other element of this model, namely the catabolite control protein CcpA, were identified from the C. acetobutylicum genome sequence, suggesting that a similar mechanism of carbon catabolite repression may operate in this industrially important organism.

Scalable Multi-Relational Association Mining

Clare, A., Williams, H. E. and Lester, N. M. (2004) Scalable Multi-Relational Association Mining. In proceedings of the 4th International Conference on Data Mining ICDM '04.

Confirmation of Data Mining Based Predictions of Protein Function

King, R. D. and Wise, P. H. and Clare, A. (2004) Confirmation of Data Mining Based Predictions of Protein Function. Bioinformatics 20(7), 1110-1118

Bioinformatyczna analiza struktury i oddziaływań z substratami metylotransferaz RNA odpowiedzialnych za oporność bakterii na antybiotyki

Wydział Biologii: Instytut Biologii Molekularnej i Biotechnologii

Visible Volume: a Robust Measure for Protein Structure Characterization

We propose a new characterization of protein structure based on the natural tetrahedral geometry of the β carbon and a new geometric measure of structural similarity, called visible volume. In our model, the side-chains are replaced by an ideal tetrahedron, the orientation of which is fixed with respect to the backbone and corresponds to the preferred rotamer directions. Visible volume is a measure of the non-occluded empty space surrounding each residue position after the side-chains have been removed. It is a robust, parameter-free, locally-computed quantity that accounts for many of the spatial constraints that are of relevance to the corresponding position in the native structure. When computing visible volume, we ignore the nature of both the residue observed at each site and the ones surrounding it. We focus instead on the space that, together, these residues could occupy. By doing so, we are able to quantify a new kind of invariance beyond the apparent variations in protein families, namely, the conservation of the physical space available at structurally equivalent positions for side-chain packing. Corresponding positions in native structures are likely to be of interest in protein structure prediction, protein design, and homology modeling. Visible volume is related to the degree of exposure of a residue position and to the actual rotamers in native proteins. In this article, we discuss the properties of this new measure, namely, its robustness with respect to both crystallographic uncertainties and naturally occurring variations in atomic coordinates, and the remarkable fact that it is essentially independent of the choice of the parameters used in calculating it. We also show how visible volume can be used to align protein structures, to identify structurally equivalent positions that are conserved in a family of proteins, and to single out positions in a protein that are likely to be of biological interest. These properties qualify visible volume as a powerful tool in a variety of applications, from the detailed analysis of protein structure to homology modeling, protein structural alignment, and the definition of better scoring functions for threading purposes.

Space, Time and Learning in the Hippocampus: How Fine Spatial and Temporal Scales Are Expanded into Population Codes for Behavioral Control

The hippocampus participates in multiple functions, including spatial navigation, adaptive timing, and declarative (notably, episodic) memory. How does it carry out these particular functions? The present article proposes that hippocampal spatial and temporal processing are carried out by parallel circuits within entorhinal cortex, dentate gyrus, and CA3 that are variations of the same circuit design. In particular, interactions between these brain regions transform fine spatial and temporal scales into population codes that are capable of representing the much larger spatial and temporal scales that are needed to control adaptive behaviors. Previous models of adaptively timed learning propose how a spectrum of cells tuned to brief but different delays are combined and modulated by learning to create a population code for controlling goal-oriented behaviors that span hundreds of milliseconds or even seconds. Here it is proposed how projections from entorhinal grid cells can undergo a similar learning process to create hippocampal place cells that can cover a space of many meters that are needed to control navigational behaviors. The suggested homology between spatial and temporal processing may clarify how spatial and temporal information may be integrated into an episodic memory.

Investigation of polyhydroxyalkanoate production by an activated sludge microbial consortium treating artificial dairy wastewater

Petrochemical plastics/polymers are a common feature of day to day living as they occur in packaging, furniture, mobile phones, computers, construction equipment etc. However, these materials are produced from non-renewable materials and are resistant to microbial degradation in the environment. Considerable research has therefore been carried out into the production of sustainable, biodegradable polymers, amenable to microbial catabolism to CO2 and H2O. A key group of microbial polyesters, widely considered as optimal replacement polymers, are the Polyhydroxyalkaonates (PHAs). Primary research in this area has focused on using recombinant pure cultures to optimise PHA yields, however, despite considerable success, the high costs of pure culture fermentation have thus far hindered the commercial viability of PHAs thus produced. In more recent years work has begun to focus on mixed cultures for the optimisation of PHA production, with waste incorporations offering optimal production cost reductions. The scale of dairy processing in Ireland, and the high organic load wastewaters generated, represent an excellent potential substrate for bioconversion to PHAs in a mixed culture system. The current study sought to investigate the potential for such bioconversion in a laboratory scale biological system and to establish key operational and microbial characteristics of same. Two sequencing batch reactors were set up and operated along the lines of an enhanced biological phosphate removal (EBPR) system, which has PHA accumulation as a key step within repeated rounds of anaerobic/aerobic cycling. Influents to the reactors varied only in the carbon sources provided. Reactor 1 received artificial wastewater with acetate alone, which is known to be readily converted to PHA in the anaerobic step of EBPR. Reactor 2 wastewater influent contained acetate and skim milk to imitate a dairy processing effluent. Chemical monitoring of nutrient remediation within the reactors as continuously applied and EBPR consistent performances observed. Qualitative analysis of the sludge was carried out using fluorescence microscopy with Nile Blue A lipophillic stain and PHA production was confirmed in both reactors. Quantitative analysis via HPLC detection of crotonic acid derivatives revealed the fluorescence to be short chain length Polyhydroxybutyrate, with biomass dry weight accumulations of 11% and 13% being observed in reactors 1 and 2, respectively. Gas Chromatography-Mass Spectrometry for medium chain length methyl ester derivatives revealed the presence of hydroxyoctanoic, -decanoic and -dodecanoic acids in reactor 1. Similar analyses in reactor 2 revealed monomers of 3-hydroxydodecenoic and 3-hydroxytetradecanoic acids. Investigation of the microbial ecology of both reactors as conducted in an attempt to identify key species potentially contributing to reactor performance. Culture dependent investigations indicated that quite different communities were present in both reactors. Reactor 1 isolates demonstrated the following species distributions Pseudomonas (82%), Delftia acidovorans (3%), Acinetobacter sp. (5%) Aminobacter sp., (3%) Bacillus sp. (3%), Thauera sp., (3%) and Cytophaga sp. (3%). Relative species distributions among reactor 2 profiled isolates were more evenly distributed between Pseudoxanthomonas (32%), Thauera sp (24%), Acinetobacter (24%), Citrobacter sp (8%), Lactococcus lactis (5%), Lysinibacillus (5%) and Elizabethkingia (2%). In both reactors Gammaproteobacteria dominated the cultured isolates. Culture independent 16S rRNA gene analyses revealed differing profiles for both reactors. Reactor 1 clone distribution was as follows; Zooglea resiniphila (83%), Zooglea oryzae (2%), Pedobacter composti (5%), Neissericeae sp. (2%) Rhodobacter sp. (2%), Runella defluvii (3%) and Streptococcus sp. (3%). RFLP based species distribution among the reactor 2 clones was as follows; Runella defluvii (50%), Zoogloea oryzae (20%), Flavobacterium sp. (9%), Simplicispira sp. (6%), Uncultured Sphingobacteria sp. (6%), Arcicella (6%) and Leadbetterella bysophila (3%). Betaproteobacteria dominated the 16S rRNA gene clones identified in both reactors. FISH analysis with Nile Blue dual staining resolved these divergent findings, identifying the Betaproteobacteria as dominant PHA accumulators within the reactor sludges, although species/strain specific allocations could not be made. GC analysis of the sludge had indicated the presence of both medium chain length as well short chain length PHAs accumulating in both reactors. In addition the cultured isolates from the reactors had been identified previously as mcl and scl PHA producers, respectively. Characterisations of the PHA monomer profiles of the individual isolates were therefore performed to screen for potential novel scl-mcl PHAs. Nitrogen limitation driven PHA accumulation in E2 minimal media revealed a greater propensity among isoates for mcl-pHA production. HPLC analysis indicated that PHB production was not a major feature of the reactor isolates and this was supported by the low presence of scl phaC1 genes among PCR screened isolates. A high percentage distribution of phaC2 mcl-PHA synthase genes was recorded, with the majority sharing high percentage homology with class II synthases from Pseudomonas sp. The common presence of a phaC2 homologue was not reflected in the production of a common polymer. Considerable variation was noted in both the monomer composition and ratios following GC analysis. While co-polymer production could not be demonstrated, potentially novel synthase substrate specificities were noted which could be exploited further in the future.

Zinc finger nuclease gene repair as a treatment for cystinosis

Cystinosis is a multi-system autosomal recessive disorder caused by mutations and/or deletions in both alleles of CTNS, a gene encoding for the low pH dependent lysosomal cystine exporter cystinosin. Cystinosis occurs in approximately 1:200,000 newborns worldwide and is characterised by an accumulation of cystine in the lysosomes. The most severe form of the disorder is nephropathic cystinosis presenting Fanconi syndrome and leads without treatment to an end-stage renal failure before the age of ten. The only treatment available so far is cysteamine therapy, which delays disease progression by five years, but does not provide a cure for cystinosis patients. Current gene and cell based therapeutic approaches have not yet provided a suitable alternative. A potentially approach for a long-term treatment could be to generate autologous gene–modified stem cells by repairing the gene. Zinc Finger Nucleases (ZFNs) serve as a tool to increase HDR up to a 200,000-fold by introducing a double-stranded break (DSB). Thus, simple mutations in the CTNS gene could be corrected by introduction of a double-stranded break using ZFNs to boost the process of HDR with a suitable donor DNA sequence. A permanent repair of the most common lesion CTNS, a 57 kb deletion, could be achieved by ZFN-mediated HDR using a minigene CTNS promoter/cDNA construct. The thesis describes the design and testing of seven zinc finger nuclease pairs for their cleavage activity in vitro and in cellulo.. A highly sensitive assay to detect even low levels of ZFN-mediated HDR was also developed. Finally, to further investigate the role of autophagy in tissue injury in cystinotic cells an assay to monitor autophagy levels in the cells was successfully developed. This assay provides the opportunity to demonstrate functional restoration of CTNS after successful ZFN-HDR in cystinotic cells.

Genetic analysis of bacteriophages from clinical and environmental samples

Bacteriophages, viruses infecting bacteria, are uniformly present in any location where there are high numbers of bacteria, both in the external environment and the human body. Knowledge of their diversity is limited by the difficulty to culture the host species and by the lack of the universal marker gene present in all viruses. Metagenomics is a powerful tool that can be used to analyse viral communities in their natural environments. The aim of this study was to investigate diverse populations of uncultured viruses from clinical (a sputum of patient with cystic fibrosis, CF) and environmental samples (a sludge from a dairy food wastewater treatment plant) containing rich bacterial populations using genetic and metagenomic analyses. Metagenomic sequencing of viruses obtained from these samples revealed that the majority of the metagenomic reads (97-99%) were novel when compared to the NCBI protein database using BLAST. A large proportion of assembled contigs were assignable as novel phages or uncharacterised prophages, the next largest assignable group being single-stranded eukaryotic virus genomes. Sputum from a cystic fibrosis patient contained DNA typical of phages of bacteria that are traditionally involved in CF lung infections and other bacteria that are part of the normal oral flora. The only eukaryotic virus detected in the CF sputum was Torque Teno virus (TTV). A substantial number of assigned sequences from dairy wastewater could be affiliated with phages of bacteria that are typically found in the soil and aquatic environments, including wastewater. Eukaryotic viral sequences were dominated by plant pathogens from the Geminiviridae and Nanoviridae families, and animal pathogens from the Circoviridae family. Antibiotic resistance genes were detected in both metagenomes suggesting phages could be a source for transmissible antimicrobial resistance. Overall, diversity of viruses in the CF sputum was low, with 89 distinct viral genotypes predicted, and higher (409 genotypes) in the wastewater. Function-based screening of a metagenomic library constructed from DNA extracted from dairy food wastewater viruses revealed candidate promoter sequences that have ability to drive expression of GFP in a promoter-trap vector in Escherichia coli. The majority of the cloned DNA sequences selected by the assay were related to ssDNA circular eukaryotic viruses and phages which formed a minority of the metagenome assembly, and many lacked any significant homology to known database sequences. Natural diversity of bacteriophages in wastewater samples was also examined by PCR amplification of the major capsid protein sequences, conserved within T4-type bacteriophages from Myoviridae family. Phylogenetic analysis of capsid sequences revealed that dairy wastewater contained mainly diverse and uncharacterized phages, while some showed a high level of similarity with phages from geographically distant environments.

A study of diet and health in the elderly; the gut microbiota as a source of bioactive agents

The global proportion of older persons is increasing rapidly. Diet and the intestinal microbiota independently and jointly contribute to health in the elderly. The habitual dietary patterns and functional microbiota components of elderly subjects were investigated in order to identify specific effector mechanisms. A study of the dietary intake of Irish community-dwelling elderly subjects showed that the consumption of foods high in fat and/or sugar was excessive, while consumption of dairy foods was inadequate. Elderly females typically had a more nutrient- dense diet than males and a considerable proportion of subjects, particularly males, had inadequate intakes of calcium, magnesium, vitamin D, folate, zinc and vitamin C. The association between dietary patterns, glycaemic index and cognitive function was also investigated. Elderly subjects consuming ‘prudent’ dietary patterns had better cognitive function compared to those consuming ‘Western’ dietary patterns. Furthermore, fully-adjusted regression models revealed that a high glycaemic diet was associated with poor cognitive function, demonstrating a new link between nutrition and cognition. An extensive screening study of the elderly faecal-derived microbiota was also undertaken to examine the prevalence of antimicrobial production by intestinal bacteria. A number of previously characterised bacteriocins were isolated (gassericin T, ABP-118, mutacin II, enterocin L-50 and enterocin P) in this study. Interestingly, a Lactobacillus crispatus strain was found to produce a potentially novel antimicrobial compound. Full genome sequencing of this strain revealed the presence of three loci which exhibited varying degrees of homology with the genes responsible for helveticin J production in Lb. helveticus. An additional study comparing the immunomodulatory capacity of ‘viable’ and ‘non-viable’ Bifidobacterium strains found that Bifidobacterium-fermented milks (BFMs) containing ‘non-viable’ cells could stimulate levels of IL-10 and TNF-α in a manner similar to those stimulated by BFMs containing ‘viable’ cells in vitro.

Cystic fibrosis gene repair: correction of ΔF508 using ZFN and CRISPR/Cas9 guide RNA gene editing tools

Cystic Fibrosis (CF) is an autosomal recessive monogenic disorder caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene with the ΔF508 mutation accounting for approximately 70% of all CF cases worldwide. This thesis investigates whether existing zinc finger nucleases designed in this lab and CRISPR/gRNAs designed in this thesis can mediate efficient homology-directed repair (HDR) with appropriate donor repair plasmids to correct CF-causing mutations in a CF cell line. Firstly, the most common mutation, ΔF508, was corrected using a pair of existing ZFNs, which cleave in intron 9, and the donor repair plasmid pITR-donor-XC, which contains the correct CTT sequence and two unique restriction sites. HDR was initially determined to be <1% but further analysis by next generation sequencing (NGS) revealed HDR occurred at a level of 2%. This relatively low level of repair was determined to be a consequence of distance from the cut site to the mutation and so rather than designing a new pair of ZFNs, the position of the existing intron 9 ZFNs was exploited and attempts made to correct >80% of CF-causing mutations. The ZFN cut site was used as the site for HDR of a mini-gene construct comprising exons 10-24 from CFTR cDNA (with appropriate splice acceptor and poly A sites) to allow production of full length corrected CFTR mRNA. Finally, the ability to cleave closer to the mutation and mediate repair of CFTR using the latest gene editing tool CRISPR/Cas9 was explored. Two CRISPR gRNAs were tested; CRISPR ex10 was shown to cleave at an efficiency of 15% and CRISPR in9 cleaved at 3%. Both CRISPR gRNAs mediated HDR with appropriate donor plasmids at a rate of ~1% as determined by NGS. This is the first evidence of CRISPR induced HDR in CF cell lines.

Characterisation of bacteriophage-host interactions in Lactococcus lactis

Lactococcus lactis is used extensively world-wide for the production of fermented dairy products. Bacteriophages (phages) infecting L. lactis can result in slow or incomplete fermentations, or may even cause total fermentation failure. Therefore, bacteriophages disrupting L. lactis fermentation are of economic concern. This thesis employed a multifaceted approach to investigate various molecular aspects of phage-host interaction in L. lactis. The genome sequence of an Irish dairy starter strain, the prophage-cured L. lactis subsp. cremoris UC509.9, was studied. The 2,250,427 bp circular chromosome represents the smallest among its sequenced lactococcal equivalents. The genome displays clear genetic adaptation to the dairy niche in the form of extensive reductive evolution. Gene prediction identified 2066 protein-encoding genes, including 104 which showed significant homology to transposase-specifying genes. Over 9 % of the identified genes appear to be inactivated through stop codons or frame shift mutations. Many pseudogenes were found in genes that are assigned to carbohydrate and amino acid transport and metabolism orthologous groups, reflecting L. lactis UC509.9’s adaptation to the lactose and casein-rich dairy environment. Sequence analysis of the eight plasmids of L. lactis revealed extensive adaptation to the dairy environment. Key industrial phenotypes were mapped and novel lactococcal plasmid-associated genes highlighted. In addition to chromosomally-encoded bacteriophage resistance systems, six functional such systems were identified, including two abortive infection systems, AbiB and AbiD1, explaining the observed phage resistance of L. lactis UC509.9 Molecular analysis suggests that the constitutive expression of AbiB is not lethal to cells, suggesting the protein is expressed in an un/inactivated form. Analysis of 936 species phage sk1-escape mutants of AbiB revealed that all such mutants harbour mutations in orf6, which encodes the major capsid protein. Results suggest that the major capsid protein is required for activation of the AbiB system, although this requires furrther investigations. Temporal transcriptomes of L. lactis UC509.9 undergoing lytic infection with either one of two distinct bacteriophages, Tuc2009 and c2, was determined and compared to the transcriptome of uninfected UC509.9 cells. Whole genome microarrays performed at various time-points post-infection demonstrated a rather modest impact on host transcription. Alterations in the UC509.9 transcriptome during lytic infection appear phage-specific, with a relatively small number of differentially transcribed genes shared between infection with either Tuc2009 or c2. Transcriptional profiles of both bacteriophages during lytic infection was shown to generally correlate with previous studies and allowed the confirmation of previously predicted promoter sequences. Bioinformatic analysis of genomic regions encoding the presumed cell wall polysaccharide (CW PS) biosynthesis gene cluster of several strains of L. lactis was performed. Results demonstrate the presence of three dominant genetic types of this gene cluster, termed type A, B and C. These regions were used for the development of a multiplex PCR to identify CW PS genotype of various lactococcal strains. Analysis of 936 species phage receptor binding protein phylogeny (RBP) and CW PS genotype revealed an apparent correlation between RBP phylogeny and CW PS type, thereby providing a partial explanation for the observed narrow host range of 936 phages. Further analysis of the genetic locus encompassing the presumed CW PS biosynthesis operon of eight strains identified as belonging to the CW PS C (geno)type, revealed the presence of a variable region among the examined strains. The obtained comparative analysis allowed for the identification of five subgroups of the C type, named C1 to C5. We purified an acidic polysaccharide from the cell wall of L. lactis 3107 (C2 subtype) and confirmed that it is structurally different from the CW PS of the C1 subtype L. lactis MG1363. Combinations of genes from the variable region of C2 subtype were amplified from L. lactis 3107 and introduced into a mutant of the C1 subtype L. lactis NZ9000 (a direct derivative of MG1363) deficient in CW PS biosynthesis. The resulting recombinant mutant synthesized a CW PS with a composition characteristic for that of the C2 subtype L. lactis 3107 and not the wildtype C1 L. lactis NZ9000. The recombinant mutant exhibited a changed phage resistance/sensitivity profile consistent with that of L. lactis 3107, which unambiguously demonstrated that L. lactis 3107 CW PS is the host cell surface receptor of two bacteriophages belonging to the P335 species as well as phages that are member of the 936 species. The research presented in this thesis has significantly advanced our understanding of L. lactis bacteriophage-host interactions in several ways. Firstly, the examination of plasmidencoded bacteriophage resistance systems has allowed inferences to be made regarding the mode of action of AbiB, thereby providing a platform for further elucidation of the molecular trigger of this system. Secondly, the phage infection transcriptome data presented, in addition to previous work, has made L. lactis a model organism in terms of transcriptomic studies of bacteriophage-host interactions. And finally, the research described in this thesis has for the first time explicitly revealed the nature of a carbohydrate bacteriophage receptor in L. lactis, while also providing a logical explanation for the observed narrow host ranges exhibited by 936 and P335 phages. Future research in discerning the structures of other L. lactis CW PS, combined with the determination of the molecular interplay between receptor binding proteins of these phages and CW PS will allow an in depth understanding of the mechanism by which the most prevalent lactococcal phages identify and adsorb to their specific host.

Chiropteran types I and II interferon genes inferred from genome sequencing traces by a statistical gene-family assembler.

BACKGROUND: The rate of emergence of human pathogens is steadily increasing; most of these novel agents originate in wildlife. Bats, remarkably, are the natural reservoirs of many of the most pathogenic viruses in humans. There are two bat genome projects currently underway, a circumstance that promises to speed the discovery host factors important in the coevolution of bats with their viruses. These genomes, however, are not yet assembled and one of them will provide only low coverage, making the inference of most genes of immunological interest error-prone. Many more wildlife genome projects are underway and intend to provide only shallow coverage. RESULTS: We have developed a statistical method for the assembly of gene families from partial genomes. The method takes full advantage of the quality scores generated by base-calling software, incorporating them into a complete probabilistic error model, to overcome the limitation inherent in the inference of gene family members from partial sequence information. We validated the method by inferring the human IFNA genes from the genome trace archives, and used it to infer 61 type-I interferon genes, and single type-II interferon genes in the bats Pteropus vampyrus and Myotis lucifugus. We confirmed our inferences by direct cloning and sequencing of IFNA, IFNB, IFND, and IFNK in P. vampyrus, and by demonstrating transcription of some of the inferred genes by known interferon-inducing stimuli. CONCLUSION: The statistical trace assembler described here provides a reliable method for extracting information from the many available and forthcoming partial or shallow genome sequencing projects, thereby facilitating the study of a wider variety of organisms with ecological and biomedical significance to humans than would otherwise be possible.

Screening the human exome: a comparison of whole genome and whole transcriptome sequencing.

BACKGROUND: There is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important. RESULTS: Here we provide a systematic exploration of how well RNA-Seq can identify human coding variants by comparing variants identified through high coverage whole-genome sequencing to those identified by high coverage RNA-Seq in the same individual. This comparison allowed us to directly evaluate the sensitivity and specificity of RNA-Seq in identifying coding variants, and to evaluate how key parameters such as the degree of coverage and the expression levels of genes interact to influence performance. We find that although only 40% of exonic variants identified by whole genome sequencing were captured using RNA-Seq; this number rose to 81% when concentrating on genes known to be well-expressed in the source tissue. We also find that a high false positive rate can be problematic when working with RNA-Seq data, especially at higher levels of coverage. CONCLUSIONS: We conclude that as long as a tissue relevant to the trait under study is available and suitable quality control screens are implemented, RNA-Seq is a fast and inexpensive alternative approach for finding coding variants in genes with sufficiently high expression levels.