Neonatal screening for amino acidaemias in Karnataka, south India

Consanguineous marriages are strongly favoured among the peoples of South India. Because of the potential genetic risks resulting from inbreeding, a neonatal screening project was established in 1980 in the state of Karnataka for the identification of amino acidaemias. To date, blood samples obtained by toe-stab from 98,256 neonates have been tested by thin layer chromatography, with 46 single and 70 general amino acidaemias detected. The coefficients of inbreeding (F) for the two groups of neonates were 0.0336 and 0.0350, by comparison with a previously determined F value for the general, new-born population of 0.0298. The most common single abnormality detected was tyrosinaemia, with spontaneous resolution in the majority of cases.

Genetic and epigenetic mechanisms of tumor predisposition in hereditary non-polyposis colorectal carcinoma and sporadic cancers

Individuals with inherited deficiency in DNA mismatch repair(MMR) (Lynch syndrome) LS are predisposed to different cancers in a non-random fashion. Endometrial cancer (EC) is the most common extracolonic malignancy in LS. LS represents the best characterized form of hereditary nonpolyposis colorectal carcinoma (HNPCC). Other forms of familial non-polyposis colon cancer exist, including familial colorectal cancer type X (FCCX). This syndrome resembles LS, but MMR gene defects are excluded and the predisposition genes are unknown so far. To address why different organs are differently susceptible to cancer development, we examined molecular similarities and differences in selected cancers whose frequency varies in LS individuals. Tumors that are common (colorectal, endometrial, gastric) and less common (brain, urological) in LS were characterized for MMR protein expression, microsatellite instability (MSI), and by altered DNA methylation. We also studied samples of histologically normal endometrium, endometrial hyperplasia,and cancer for molecular alterations to identify potential markers that could predict malignant transformation in LS and sporadic cases. Our results suggest that brain and kidney tumors follow a different pathway for cancer development than the most common LS related cancers.Our results suggest also that MMR defects are detectable in endometrial tissues from a proportion of LS mutation carriers prior to endometrial cancer development. Traditionally (complex) atypical hyperplasia has been considered critical for progression to malignancy. Our results suggest that complex hyperplasia without atypia is equally important as a precursor lesion of malignancy. Tumor profiles from Egypt were compared with colorectal tumors from Finland to evaluate if there are differences specific to the ethnic origin (East vs.West). Results showed for the first time a distinct genetic and epigenetic signature in the Egyptian CRC marked by high methylation of microsatellite stable tumors associated with advanced stage, and low frequency of Wnt signaling activation, suggesting a novel pathway. DNA samples from FCCX families were studied with genome wide linkage analysis using microsatellite markers. Selected genes from the linked areas were tested for possible mutations that could explain predisposition to a large number of colon adenomas and carcinomas seen in these families. Based on the results from the linkage analysis, a number of areas with tentative linkage were identified in family 20. We narrowed down these areas by additional microsatellite markers to found a mutation in the BMPR1A gene. Sequencing of an additional 17 FCCX families resulted in a BMPR1A mutation frequency of 2/18 families (11%). Clarification of the mechanisms of the differential tumor susceptibility in LS increases the understanding of gene and organ specific targets of MMR deficiency. While it is generally accepted that widespread MMR deficiency and consequent microsatellite instability (MSI) drives tumorigenesis in LS, the timing of molecular alterations is controversial. In particular, it is important to know that alterations may occur several years before cancer formation, at stages that are still histologically regarded as normal. Identification of molecular markers that could predict the risk of malignant transformation may be used to improve surveillance and cancer prevention in genetically predisposed individuals. Significant fractions of families with colorectal and/or endometrial cancer presently lack molecular definition altogether. Our findings expand the phenotypic spectrum of BMPR1A mutations and, for the first time, link FCCX families to the germline mutation of a specific gene. In particular, our observations encourage screening of additional families with FCCX for BMPR1A mutation, which is necessary in obtaining a reliable estimate of the share of BMPR1A-associated cases among all FCCX families worldwide. Clinically, the identification of predisposing mutations enables targeted cancer prevention in proven mutation carriers and thereby reduces cancer morbidity and mortality in the respective families.

Dissection of the genetic background of childhood onset progressive myoclonic epilepsies

The progressive myoclonic epilepsies (PMEs) are a clinically and etiologically heterogeneous group of symptomatic epilepsies characterized by myoclonus, tonic-clonic seizures, psychomotor regression and ataxia. Different disorders have been classified as PMEs. Of these, the group of neuronal ceroid lipofuscinoses (NCLs) comprise an entity that has onset in childhood, being the most common cause of neurodegeneration in children. The primary aim of this thesis was to dissect the molecular genetic background of patients with childhood onset PME by studying candidate genes and attempting to identify novel PME-associated genes. Another specific aim was to study the primary protein properties of the most recently identified member of the NCL-causing proteins, MFSD8. To dissect the genetic background of a cohort of Turkish patients with childhood onset PME, a screen of the NCL-associated genes PPT1, TPP1, CLN3, CLN5, CLN6, MFSD8, CLN8 and CTSD was performed. Altogether 49 novel mutations were identified, which together with 56 mutations found by collaborators raised the total number of known NCL mutations to 364. Fourteen of the novel mutations affect the recently identified MFSD8 gene, which had originally been identified in a subset of mainly Turkish patients as the underlying cause of CLN7 disease. To investigate the distribution of MFSD8 defects, a total of 211 patients of different ethnic origins were evaluated for mutations in the gene. Altogether 45 patients from nine different countries were provided with a CLN7 molecular diagnosis, denoting the wide geographical occurrence of MFSD8 defects. The mutations are private with only one having been established by a founder-effect in the Roma population from the former Czechoslovakia. All mutations identified except one are associated with the typical clinical picture of variant late-infantile NCL. To address the trafficking properties of MFSD8, lysosomal targeting of the protein was confirmed in both neuronal and non-neuronal cells. The major determinant for this lysosomal sorting was identified to be an N-terminal dileucine based signal (9-EQEPLL-14), recognized by heterotetrameric AP-1 adaptor proteins, suggesting that MFSD8 takes the direct trafficking pathway en route to the lysosomes. Expression studies revealed the neurons as the primary cell-type and the hippocampus and cerebellar granular cell layer as the predominant regions in which MFSD8 is expressed. To identify novel genes associated with childhood onset PME, a single nucleotide polymorphism (SNP) genomewide scan was performed in three small families and 18 sporadic patients followed by homozygosity mapping to determine the candidate loci. One of the families and a sporadic patient were positive for mutations in PLA2G6, a gene that had previously been shown to cause infantile neuroaxonal dystrophy. Application of next-generation sequencing of candidate regions in the remaining two families led to identification of a homozygous missense mutation in USP19 for the first and TXNDC6 for the second family. Analysis of the 18 sporadic cases mapped the best candidate interval in a 1.5 Mb region on chromosome 7q21. Screening of the positional candidate KCTD7 revealed six mutations in seven unrelated families. All patients with mutations in KCTD7 were reported to have early onset PME, rapid disease progression leading to dementia and no pathologic hallmarks. The identification of KCTD7 mutations in nine patients and the clinical delineation of their phenotype establish KCTD7 as a gene for early onset PME. The findings presented in this thesis denote MFSD8 and KCTD7 as genes commonly associated with childhood onset symptomatic epilepsy. The disease-associated role of TXNDC6 awaits verification through identification of additional mutations in patients with similar phenotypes. Completion of the genetic spectrum underlying childhood onset PMEs and understanding of the gene products functions will comprise important steps towards understanding the underlying pathogenetic mechanisms, and will possibly shed light on the general processes of neurodegeneration and nervous system regulation, facilitating the diagnosis, classification and possibly treatment of the affected cases.

Colored polydimethylsiloxane micropillar arrays for high throughput measurements of forces applied by genetic model organisms

Measuring forces applied by multi-cellular organisms is valuable in investigating biomechanics of their locomotion. Several technologies have been developed to measure such forces, for example, strain gauges, micro-machined sensors, and calibrated cantilevers. We introduce an innovative combination of techniques as a high throughput screening tool to assess forces applied by multiple genetic model organisms. First, we fabricated colored Polydimethylsiloxane (PDMS) micropillars where the color enhances contrast making it easier to detect and track pillar displacement driven by the organism. Second, we developed a semiautomated graphical user interface to analyze the images for pillar displacement, thus reducing the analysis time for each animal to minutes. The addition of color reduced the Young's modulus of PDMS. Therefore, the dye-PDMS composite was characterized using Yeoh's hyperelastic model and the pillars were calibrated using a silicon based force sensor. We used our device to measure forces exerted by wild type and mutant Caenorhabditis elegans moving on an agarose surface. Wild type C. elegans exert an average force of similar to 1 mu N on an individual pillar and a total average force of similar to 7.68 mu N. We show that the middle of C. elegans exerts more force than its extremities. We find that C. elegans mutants with defective body wall muscles apply significantly lower force on individual pillars, while mutants defective in sensing externally applied mechanical forces still apply the same average force per pillar compared to wild type animals. Average forces applied per pillar are independent of the length, diameter, or cuticle stiffness of the animal. We also used the device to measure, for the first time, forces applied by Drosophila melanogaster larvae. Peristaltic waves occurred at 0.4Hz applying an average force of similar to 1.58 mu N on a single pillar. Our colored microfluidic device along with its displacement tracking software allows us to measure forces applied by multiple model organisms that crawl or slither to travel through their environment. (C) 2015 AIP Publishing LLC.

The role of androgens in fetal growth: observational study in two genetic models of disordered androgen signalling.

OBJECTIVE: To examine the role of androgens on birth weight in genetic models of altered androgen signalling. SETTING: Cambridge Disorders of Sex Development (DSD) database and the Swedish national screening programme for congenital adrenal hyperplasia (CAH). PATIENTS: (1) 29 girls with XY karyotype and mutation positive complete androgen insensitivity syndrome (CAIS); (2) 43 girls and 30 boys with genotype confirmed CAH. MAIN OUTCOME MEASURES: Birth weight, birth weight-for-gestational-age (birth weight standard deviation score (SDS)) calculated by comparison with national references. RESULTS: Mean birth weight SDS in CAIS XY infants was higher than the reference for girls (mean, 95% CI: 0.4, 0.1 to 0.7; p=0.02) and was similar to the national reference for boys (0.1, -0.2 to 0.4). Birth weight SDS in CAH girls was similar to the national reference for girls (0.0, -0.2 to 0.2) and did not vary by severity of gene mutation. Birth weight SDS in CAH boys was also similar to the national reference for boys (0.2, -0.2 to 0.6). CONCLUSION: CAIS XY infants have a birth weight distribution similar to normal male infants and birth weight is not increased in infants with CAH. Alterations in androgen signalling have little impact on birth weight. Sex dimorphism in birth size is unrelated to prenatal androgen exposure.

The MLH1 c.1852_1853delinsGC (p.K618A) Variant in Colorectal Cancer: Genetic Association Study in 18,723 Individuals

Colorectal cancer is one of the most frequent neoplasms and an important cause of mortality in the developed world. Mendelian syndromes account for about 5% of the total burden of CRC, being Lynch syndrome and familial adenomatous polyposis the most common forms. Lynch syndrome tumors develop mainly as a consequence of defective DNA mismatch repair associated with germline mutations in MLH1, MSH2, MSH6 and PMS2. A significant proportion of variants identified by screening these genes correspond to missense or noncoding changes without a clear pathogenic consequence, and they are designated as "variants of uncertain significance'', being the c.1852_1853delinsGC (p.K618A) variant in the MLH1 gene a clear example. The implication of this variant as a low-penetrance risk variant for CRC was assessed in the present study by performing a case-control study within a large cohort from the COGENT consortium-COST Action BM1206 including 18,723 individuals (8,055 colorectal cancer cases and 10,668 controls) and a case-only genotype-phenotype correlation with several clinical and pathological characteristics restricted to the Epicolon cohort. Our results showed no involvement of this variant as a low-penetrance variant for colorectal cancer genetic susceptibility and no association with any clinical and pathological characteristics including family history for this neoplasm or Lynch syndrome.

Genetic analysis of the gametophytes of Undaria pinnatifida (Phaeophyceae) with ISSR method

Inter-simple sequence repeat (ISSR) analysis was used to assess eleven pairs of Undaria pinnatifida (Harv.) Suringar male and female gametophytes. After screening fifty primers, 18 ISSR primers were selected for final analysis. A total of 104 loci were obtained, of which 77 were polymorphic, among the gametophytes studied. Genetic relationships were analyzed with simple matching (S), Jaccard's (J) and Dice's (D) distance coefficients. Little genetic variations were found among the selected Undaria gametophytes, for instance, the genetic distances ranging from 0.010 to 0.125 with Dice coefficients. UPGMA dendrograms showed that 11 pairs of Undaria gametophytes were distributed into five groups. Most Undaria strains cultivated in China exhibited closely genetic relationships with the strains from Japan. However, gametophytes from Qingdao appeared as distinct clades from other Undaria strains with all three distance coefficients used. Mantel test showed that the three distance measurements generated congruent clustering patterns on the same data. Our results demonstrated the feasibility of applying ISSR markers for genetic analysis of Undaria gametophytes. (c) 2006 Elsevier B.V. All rights reserved.

The genetic analysis and germplasm identification of the gametophytes of Undaria pinnatifida (Phaeophyceae) with RAPD method

Eleven pairs of Undaria pinnatifida (Harv.) Suringar gametophytes were identified with random amplified polymorphic DNA (RAPD) technique. After screening 100 primers, 20 ten-base primers were determined for the RAPD analysis. A total of 312 polymorphic loci were obtained, of which 97.7% were polymorphic. The primer S198 was found to distinguish all the selected Undaria pinnatifida gametophytes. The genetic distances between each two of the twenty-two U. pinnatifida gametophytes ranged from 0.080 to 0.428, while the distances to the Laminaria was 0.497 on average. After reexamination, two sequences characterized amplification region (SCAR) markers were successfully converted, which could be applied to U. pinnatifida germplasm identification. All these results demonstrated the feasibility of applying RAPD markers to germplasm characterization and identification of U. pinnatifida gametophytes, and to provide a molecular basis for Undaria breeding.

Screening of genes related to ovary development in Chinese shrimp Fenneropenaeus chinensis by suppression subtractive hybridization

The ovary of triploid shrimp Fenneropenaeus chinensis was apparently impaired compared to that of the diploid shrimp at the same age. Therefore triploid shrimp ovary is possible to be taken as a model to understand the mechanism of ovary development of shrimp compared to that of the ovary of diploid shrimp at the same age. In the present study, a suppression subtractive hybridization (SSH) technique was applied to identify differentially expressed genes in the ovary between diploid and triploid shrimp. For the forward library (RNA from the ovary of triploid shrimp as the tester), 54 genes were identified. For the reverse library (RNA from the ovary of diploid shrimp as the tester), 16 genes were identified. The identified genes encoded proteins with multiple functions, including extracellular matrix components, cytoskeleton, cell growth and death, metabolism, genetic information processing, signal transduction/transport or immunity related proteins. Eleven differentially expressed genes were selected to be confirmed in the ovaries of triploid and diploid shrimp by semi-quantitative RT-PCR. Genes encoding spermatogonial stem-cell renewal factor, cytochrome c oxidase subunits I and II, clottable protein, antimicrobial peptide and transposase showed up-regulated expressions in the ovary of triploid shrimp. Genes encoding tubulin, cellular apoptosis susceptibility protein, farnesoic acid O-methyltransferase, thrombospondin and heat shock protein 90 genes showed higher expressions in the ovary of diploid shrimp. The differential expressions of the above genes are suggested to be related to the ovary development of shrimp. It will provide a new clue to uncover the molecular mechanisms underlying the ovarian development in penaeid shrimp. (C) 2010 Elsevier Inc. All rights reserved.

Sensitive and Rapid Detection of Genetic Modified Soybean (Roundup Ready) by Loop-Mediated Isothermal Amplification

Using the LAMP method, a highly specific and sensitive detection system for genetically modified soybean (Roundup Ready) was designed. In this detection system, a set of four primers was designed by targeting the exogenous 35S epsps gene. Target DNA was amplified and visualized on agarose gel within 45 min under isothermal conditions at 65 degrees C. Without gel electrophoresis, the LAMP amplicon was visualized directly in the reaction tube by the addition of SYBR Green I for naked-eye inspection. The detection sensitivity of LAMP was 10-fold higher than the nested PCR established in our laboratory. Moreover, the LAMP method was much quicker, taking only 70 min, as compared with 300 min for nested PCR to complete the analysis of the GM soybean. Compared with traditional PCR approaches, the LAMP procedure is faster and more sensitive, and there is no need for a special PCR machine or electrophoresis equipment. Hence, this method can be a very useful tool for GMO detection and is particularly convenient for fast screening.

Molecular genetic analysis of DNA alterations in Irish colorectal cancer

Colorectal cancer is the most common cause of death due to malignancy in nonsmokers in the western world. In 1995 there were 1,757 cases of colon cancer in Ireland. Most colon cancer is sporadic, however ten percent of cases occur where there is a previous family history of the disease. In an attempt to understand the tumorigenic pathway in Irish colon cancer patients, a number of genes associated with colorectal cancer development were analysed in Irish sporadic and HNPCC colon cancer patients. The hereditary forms of colon cancer include Familial adenomatous polyposis coli (FAP) and Hereditary Non-Polyposis Colon Cancer (HNPCC). Genetic analysis of the gene responsible for FAP, (the APC gene) has been previously performed on Irish families, however the genetic analysis of HNPCC families is limited. In an attempt to determine the mutation spectrum in Irish HNPCC pedigrees, the hMSH2 and hMLHl mismatch repair genes were screened in 18 Irish HNPCC families. Using SSCP analysis followed by DNA sequencing, five mutations were identified, four novel and a previously reported mutation. In families where a mutation was detected, younger asyptomatic members were screened for the presence of the predisposing mutation (where possible). Detection of mutations is particularly important for the identification of at risk individuals as the early diagnosis of cancer can vastly improve the prognosis. The sensitive and efficient detection of multiple different mutations and polymorphisms in DNA is of prime importance for genetic diagnosis and the identification of disease genes. A novel mutation detection technique has recently been developed in our laboratory. In order to assess the efficacy and application of the methodology in the analysis of cancer associated genes, a protocol for the analysis of the K-ras gene was developed and optimised. Matched normal and tumour DNA from twenty sporadic colon cancer patients was analysed for K-ras mutations using the Glycosylase Mediated Polymorphism Detection technique. Five mutations of the K-ras gene were detected using this technology. Sequencing analysis verified the presence of the mutations and SSCP analysis of the same samples did not identify any additional mutations. The GMPD technology proved to be highly sensitive, accurate and efficient in the identification of K-ras gene mutations. In order to investigate the role of the replication error phenomenon in Irish colon cancer, 3 polyA tract repeat loci were analysed. The repeat loci included a 10 bp intragenic repeat of the TGF-β-RII gene. TGF-β-RII is involved in the TGF-β epithelial cell growth pathway and mutation of the gene is thought to play a role in cell proliferation and tumorigenesis. Due to the presence of a repeat sequence within the gene, TGFB-RII defects are associated with tumours that display the replication error phenomenon. Analysis of the TGF-β-RII 10 bp repeat failed to identify mutations in any colon cancer patients. Analysis of the Bat26 and Bat 40 polyA repeat sequences in the sporadic and HNPCC families revealed that instability is associated with HNPCC tumours harbouring mismatch repair defects and with 20 % of sporadic colon cancer tumours. No correlation between K-ras gene mutations and the RER+ phenotype was detected in sporadic colon cancer tumours.

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.

The biodiversity, control and genetic characterisation of the lactococcal 936 group phages

Phages belonging to the 936 group represent one of the most prevalent and frequently isolated phages in dairy fermentation processes using Lactococcus lactis as the primary starter culture. In recent years extensive research has been carried out to characterise this phage group at a genomic level in an effort to understand how the 936 group phages dominate this particular niche and cause regular problems during large scale milk fermentations. This thesis describes a large scale screening of industrial whey samples, leading to the isolation of forty three genetically different lactococcal phages. Using multiplex PCR, all phages were identified as members of the 936 group. The complete genome of thirty eight of these phages was determined using next generation sequencing technologies which identified several regions of divergence. These included the structural region surrounding the major tail protein, the replication region as well as the genes involved in phage DNA packing. For a number of phages the latter genomic region was found to harbour genes encoding putative orphan methyltransferases. Using small molecule real time (SMRT) sequencing and heterologous gene expression, the target motifs for several of these MTases were determined and subsequently shown to actively protect phage DNA from restriction endonuclease activity. Comparative analysis of the thirty eight phages with fifty two previously sequenced members of this group showed that the core genome consists of 28 genes, while the non-core genome was found to fluctuate irrespective of geographical location or time of isolation. This study highlights the continued need to perform large scale characterisation of the bacteriophage populations infecting industrial fermentation facilities in effort to further our understanding dairy phages and ways to control their proliferation.

Analysis of educational materials and destruction/opt-out initiatives for storage and use of residual newborn screening samples.

In recent years, the storage and use of residual newborn screening (NBS) samples has gained attention. To inform ongoing policy discussions, this article provides an update of previous work on new policies, educational materials, and parental options regarding the storage and use of residual NBS samples. A review of state NBS Web sites was conducted for information related to the storage and use of residual NBS samples in January 2010. In addition, a review of current statutes and bills introduced between 2005 and 2009 regarding storage and/or use of residual NBS samples was conducted. Fourteen states currently provide information about the storage and/or use of residual NBS samples. Nine states provide parents the option to request destruction of the residual NBS sample after the required storage period or the option to exclude the sample for research uses. In the coming years, it is anticipated that more states will consider policies to address parental concerns about the storage and use of residual NBS samples. Development of new policies regarding storage and use of residual NBS samples will require careful consideration of impact on NBS programs, parent and provider educational materials, and respect for parents among other issues.

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.