Comparative analysis of the cranial osteology and evolutionary traits of micrastur gilvicollis, micrastur ruficollis and micrastur semitorquatus (aves: falconiformes)

The cranial osteology of Micrastur gilvicollis (Vieillot, 1817), Micrastur ruficollis (Vieillot, 1817) and Micrastur semitorquatus (Vieillot, 1817) is comparatively and meticulously described to characterize each of the species and to determine which traits the species have in common and which are distinct. These traits will be used a posteriori for phylogenetic analysis. Our results indicate that M. gilvicollis and M. ruficollis are closely related, as they share a large number of traits, including a lacrimal bone with a distal portion that is approximately half as long as the proximal portion and a parasphenoid rostrum that covers 50% of the distance between the occipital condyle and pterygoid. Similarly, M. gilvicollis and M. semitorquatus both have a partially fused craniofacial flexion zone. In both M. ruficollis and M. semitorquatus, the symphyseal region of the mandible is 1/5 the total length of the mandible. The diagnostic traits for each of these species are as follows: a) in M. gilvicollis, the interorbital distance is 1/3 the length of the parietal, and the zygomatic process stretches 1/5 of the distance from the orbital arch to the jugal arch; b) in M. ruficollis, the interorbital distance is 2/5 of the length of the parietal and the zygomatic process extends 1/4 of the distance from the orbital arch to the jugal arch; and c) in M. semitorquatus, the interorbital distance is 3/7 the length of the parietal and the distal portion of the lacrimal is 1/3 the length of the proximal portion. Among the three species, M. gilvicolis and M. ruficollis share the most traits, which leads us to infer that these species are more closely related to one another than they are to M. semitorquatus. Phylogenetic analysis performed a posteriori may confirm the relationship between these three species.

Pedigree analysis and inbreeding depression on growth traits in Brazilian Marchigiana and Bonsmara breeds

The study of population structure by pedigree analysis is useful to identify important circumstances that affect the genetic history of populations. The intensive use of a small number of superior individuals may reduce the genetic diversity of populations. This situation is very common for the beef cattle breeds. Therefore, the objectives of the present study were to analyze the pedigree and possible inbreeding depression on traits of economic interest in the Marchigiana and Bonsmara breeds and to test the inclusion of the individual inbreeding coefficient (F-i) or individual increases in inbreeding coefficient (Delta F-i) in the genetic evaluation model for the quantification of inbreeding depression. The complete pedigree file of the Marchigiana breed included 29,411 animals born between 1950 and 2003. For the Bonsmara breed, the pedigree file included 18,695 animals born between 1988 and 2006. Only animals with at least 2 equivalent generations of known pedigree were kept in the analyses of inbreeding effect on birth weight, weaning weight measured at about 205 d, and BW at 14 mo in the Marchigiana breed, and on birth weight, weaning weight, and scro-tal circumference measured at 12 mo in the Bonsmara breed. The degree of pedigree knowledge was greater for Marchigiana than for Bonsmara animals. The average generation interval was 7.02 and 3.19 for the Marchigiana and Bonsmara breed, respectively. The average inbreeding coefficient was 1.33% for Marchigiana and 0.26% for Bonsmara. The number of ancestors explaining 50% of the gene pool and effective population size computed via individual increase in coancestry were 13 and 97.79 for Marchigiana and 41 and 54.57 for Bonsmara, respectively. These estimates indicate reduction in genetic variability in both breeds. Inbreeding depression was observed for most of the growth traits. The model including Delta F-i can be considered more adequate to quantify inbreeding depression. The inclusion of F-i or Delta F-i in the genetic evaluation model may not result in better fit to the data. A genetic evaluation with simultaneous estimation of inbreeding depression can be performed in Marchigiana and Bonsmara breeds, providing additional information to producers and breeders.

Bayesian model accounting for within-class biological variability in Serial Analysis of Gene Expression (SAGE)

Abstract Background An important challenge for transcript counting methods such as Serial Analysis of Gene Expression (SAGE), "Digital Northern" or Massively Parallel Signature Sequencing (MPSS), is to carry out statistical analyses that account for the within-class variability, i.e., variability due to the intrinsic biological differences among sampled individuals of the same class, and not only variability due to technical sampling error. Results We introduce a Bayesian model that accounts for the within-class variability by means of mixture distribution. We show that the previously available approaches of aggregation in pools ("pseudo-libraries") and the Beta-Binomial model, are particular cases of the mixture model. We illustrate our method with a brain tumor vs. normal comparison using SAGE data from public databases. We show examples of tags regarded as differentially expressed with high significance if the within-class variability is ignored, but clearly not so significant if one accounts for it. Conclusion Using available information about biological replicates, one can transform a list of candidate transcripts showing differential expression to a more reliable one. Our method is freely available, under GPL/GNU copyleft, through a user friendly web-based on-line tool or as R language scripts at supplemental web-site.

Heritability of physical activity traits in Brazilian families: the Baependi Heart Study

Background It is commonly recognized that physical activity has familial aggregation; however, the genetic influences on physical activity phenotypes are not well characterized. This study aimed to (1) estimate the heritability of physical activity traits in Brazilian families; and (2) investigate whether genetic and environmental variance components contribute differently to the expression of these phenotypes in males and females. Methods The sample that constitutes the Baependi Heart Study is comprised of 1,693 individuals in 95 Brazilian families. The phenotypes were self-reported in a questionnaire based on the WHO-MONICA instrument. Variance component approaches, implemented in the SOLAR (Sequential Oligogenic Linkage Analysis Routines) computer package, were applied to estimate the heritability and to evaluate the heterogeneity of variance components by gender on the studied phenotypes. Results The heritability estimates were intermediate (35%) for weekly physical activity among non-sedentary subjects (weekly PA_NS), and low (9-14%) for sedentarism, weekly physical activity (weekly PA), and level of daily physical activity (daily PA). Significant evidence for heterogeneity in variance components by gender was observed for the sedentarism and weekly PA phenotypes. No significant gender differences in genetic or environmental variance components were observed for the weekly PA_NS trait. The daily PA phenotype was predominantly influenced by environmental factors, with larger effects in males than in females. Conclusions Heritability estimates for physical activity phenotypes in this sample of the Brazilian population were significant in both males and females, and varied from low to intermediate magnitude. Significant evidence for heterogeneity in variance components by gender was observed. These data add to the knowledge of the physical activity traits in the Brazilian study population, and are concordant with the notion of significant biological determination in active behavior.

Different approaches for dealing with rare variants in family-based genetic studies: application of a Genetic Analysis Workshop 17 problem

Rare variants are becoming the new candidates in the search for genetic variants that predispose individuals to a phenotype of interest. Their low prevalence in a population requires the development of dedicated detection and analytical methods. A family-based approach could greatly enhance their detection and interpretation because rare variants are nearly family specific. In this report, we test several distinct approaches for analyzing the information provided by rare and common variants and how they can be effectively used to pinpoint putative candidate genes for follow-up studies. The analyses were performed on the mini-exome data set provided by Genetic Analysis Workshop 17. Eight approaches were tested, four using the trait’s heritability estimates and four using QTDT models. These methods had their sensitivity, specificity, and positive and negative predictive values compared in light of the simulation parameters. Our results highlight important limitations of current methods to deal with rare and common variants, all methods presented a reduced specificity and, consequently, prone to false positive associations. Methods analyzing common variants information showed an enhanced sensibility when compared to rare variants methods. Furthermore, our limited knowledge of the use of biological databases for gene annotations, possibly for use as covariates in regression models, imposes a barrier to further research.

Innovative analytical methods for Central Nervous System Drug analysis in biological fluids

During recent years a consistent number of central nervous system (CNS) drugs have been approved and introduced on the market for the treatment of many psychiatric and neurological disorders, including psychosis, depression, Parkinson disease and epilepsy. Despite the great advancements obtained in the treatment of CNS diseases/disorders, partial response to therapy or treatment failure are frequent, at least in part due to poor compliance, but also genetic variability in the metabolism of psychotropic agents or polypharmacy, which may lead to sub-therapeutic or toxic plasma levels of the drugs, and finally inefficacy of the treatment or adverse/toxic effects. With the aim of improving the treatment, reducing toxic/side effects and patient hospitalisation, Therapeutic Drug Monitoring (TDM) is certainly useful, allowing for a personalisation of the therapy. Reliable analytical methods are required to determine the plasma levels of psychotropic drugs, which are often present at low concentrations (tens or hundreds of nanograms per millilitre). The present PhD Thesis has focused on the development of analytical methods for the determination of CNS drugs in biological fluids, including antidepressants (sertraline and duloxetine), antipsychotics (aripiprazole), antiepileptics (vigabatrin and topiramate) and antiparkinsons (pramipexole). Innovative methods based on liquid chromatography or capillary electrophoresis coupled to diode-array or laser-induced fluorescence detectors have been developed, together with the suitable sample pre-treatment for interference removal and fluorescent labelling in case of LIF detection. All methods have been validated according to official guidelines and applied to the analysis of real samples obtained from patients, resulting suitable for the TDM of psychotropic drugs.

A candidate gene approach to identify DNA markers associated with meat quality and production traits: investigation of several cathepsin genes in Italian heavy pigs.

The cathepsin enzymes represent an important family of lysosomal proteinases with a broad spectrum of functions in many, if not in all, tissues and cell types. In addition to their primary role during the normal protein turnover, they possess highly specific proteolytic activities, including antigen processing in the immune response and a direct role in the development of obesity and tumours. In pigs, the involvement of cathepsin enzymes in proteolytic processes have important effects during the conversion of muscle to meat, due to their influence on meat texture and sensory characteristics, mainly in seasoned products. Their contribution is fundamental in flavour development of dry-curing hams. However, several authors have demonstrated that high cathepsin activity, in particular of cathepsin B, is correlated to defects of these products, such as an excessive meat softness together with abnormal free tyrosine content, astringent or metallic aftertastes and formation of a white film on the cut surface. Thus, investigation of their genetic variability could be useful to identify DNA markers associated with these dry cured hams parameters, but also with meat quality, production and carcass traits in Italian heavy pigs. Unfortunately, no association has been found between cathepsin markers and meat quality traits so far, in particular with cathepsin B activity, suggesting that other genes, besides these, affect meat quality parameters. Nevertheless, significant associations were observed with several carcass and production traits in pigs. A recent study has demonstrated that different single nucleotide polymorphisms (SNPs) localized in cathepsin D (CTSD), F (CTSF), H and Z genes were highly associated with growth, fat deposition and production traits in an Italian Large White pig population. The aim of this thesis was to confirm some of these results in other pig populations and identify new cathepsin markers in order to evaluate their effects on cathepsin activity and other production traits. Furthermore, starting from the data obtained in previous studies on CTSD gene, we also analyzed the known polymorphism located in the insulin-like growth factor 2 gene (IGF2 intron3-g.3072G>A). This marker is considered the causative mutation for the quantitative trait loci (QTL) affecting muscle mass and fat deposition in pigs. Since IGF2 maps very close to CTSD on porcine chromosome (SSC) 2, we wanted to clarify if the effects of the CTSD marker were due to linkage disequilibrium with the IGF2 intron3-g.3072G>A mutation or not. In the first chapter, we reported the results from these two SSC2 gene markers. First of all, we evaluated the effects of the IGF2 intron3-g.3072G>A polymorphism in the Italian Large White breed, for which no previous studies have analysed this marker. Highly significant associations were identified with all estimated breeding values for production and carcass traits (P<0.00001), while no effects were observed for meat quality traits. Instead, the IGF2 intron3-g.3072G>A mutation did not show any associations with the analyzed traits in the Italian Duroc pigs, probably due to the low level of variability at this polymorphic site for this breed. In the same Duroc pig population, significant associations were obtained for the CTSD marker for all production and carcass traits (P < 0.001), after excluding possible confounding effects of the IGF2 mutation. The effects of the CTSD g.70G>A polymorphism were also confirmed in a group of Italian Large White pigs homozygous for the IGF2 intron3-g.3072G allele G (IGF2 intron3-g.3072GG) and by haplotype analysis between the markers of the two considered genes. Taken together, all these data indicated that the IGF2 intron3-g.3072G>A mutation is not the only polymorphism affecting fatness and muscle deposition in pigs. In the second chapter, we reported the analysis of two new SNPs identified in cathepsin L (CTSL) and cathepsin S (CTSS) genes and the association results with meat quality parameters (including cathepsin B activity) and several production traits in an Italian Large White pig population. Allele frequencies of these two markers were evaluated in 7 different pig breeds. Furthermore, we mapped using a radiation hybrid panel the CTSS gene on SSC4. Association studies with several production traits, carried out in 268 Italian Large White pigs, indicated positive effects of the CTSL polymorphism on average daily gain, weight of lean cuts and backfat thickness (P<0.05). The results for these latter traits were also confirmed using a selective genotype approach in other Italian Large White pigs (P<0.01). In the 268 pig group, the CTSS polymorphism was associated with feed:gain ratio and average daily gain (P<0.05). Instead, no association was observed between the analysed markers and meat quality parameters. Finally, we wanted to verify if the positive results obtained for the cathepsin L and S markers and for other previous identified SNPs (cathepsin F, cathepsin Z and their inhibitor cystatin B) were confirmed in the Italian Duroc pig breed (third chapter). We analysed them in two groups of Duroc pigs: the first group was made of 218 performance-tested pigs not selected by any phenotypic criteria, the second group was made of 100 Italian Duroc pigs extreme and divergent for visible intermuscular fat trait. In the first group, the CTSL polymorphism was associated with weight of lean cuts (P<0.05), while suggestive associations were obtained for average daily gain and backfat thickness (P<0.10). Allele frequencies of the CTSL gene marker also differed positively among the visible intermuscular extreme tails. Instead, no positive effects were observed for the other DNA markers on the analysed traits. In conclusion, in agreement with the present data and for the biological role of these enzymes, the porcine CTSD and CTSL markers: a) may have a direct effect in the biological mechanisms involved in determining fat and lean meat content in pigs, or b) these markers could be very close to the putative functional mutation(s) present in other genes. These findings have important practical applications, in particular the CTSD and CTSL mutations could be applied in a marker assisted selection (MAS) both in the Italian Large White and Italian Duroc breeds. Marker assisted selection could also increase in efficiency by adding information from the cathepsin S genotype, but only in the Italian Large White breed.

Different approaches to identify markers for meat quality and production traits in pigs: analysis of the transcriptome in post mortem skeletal muscle and investigation of candidate genes

Due to the growing attention of consumers towards their food, improvement of quality of animal products has become one of the main focus of research. To this aim, the application of modern molecular genetics approaches has been proved extremely useful and effective. This innovative drive includes all livestock species productions, including pork. The Italian pig breeding industry is unique because needs heavy pigs slaughtered at about 160 kg for the production of high quality processed products. For this reason, it requires precise meat quality and carcass characteristics. Two aspects have been considered in this thesis: the application of the transcriptome analysis in post mortem pig muscles as a possible method to evaluate meat quality parameters related to the pre mortem status of the animals, including health, nutrition, welfare, and with potential applications for product traceability (chapters 3 and 4); the study of candidate genes for obesity related traits in order to identify markers associated with fatness in pigs that could be applied to improve carcass quality (chapters 5, 6, and 7). Chapter three addresses the first issue from a methodological point of view. When we considered this issue, it was not obvious that post mortem skeletal muscle could be useful for transcriptomic analysis. Therefore we demonstrated that the quality of RNA extracted from skeletal muscle of pigs sampled at different post mortem intervals (20 minutes, 2 hours, 6 hours, and 24 hours) is good for downstream applications. Degradation occurred starting from 48 h post mortem even if at this time it is still possible to use some RNA products. In the fourth chapter, in order to demonstrate the potential use of RNA obtained up to 24 hours post mortem, we present the results of RNA analysis with the Affymetrix microarray platform that made it possible to assess the level of expression of more of 24000 mRNAs. We did not identify any significant differences between the different post mortem times suggesting that this technique could be applied to retrieve information coming from the transcriptome of skeletal muscle samples not collected just after slaughtering. This study represents the first contribution of this kind applied to pork. In the fifth chapter, we investigated as candidate for fat deposition the TBC1D1 [TBC1 (tre-2/USP6, BUB2, cdc16) gene. This gene is involved in mechanisms regulating energy homeostasis in skeletal muscle and is associated with predisposition to obesity in humans. By resequencing a fragment of the TBC1D1 gene we identified three synonymous mutations localized in exon 2 (g.40A>G, g.151C>T, and g.172T>C) and 2 polymorphisms localized in intron 2 (g.219G>A and g.252G>A). One of these polymorphisms (g.219G>A) was genotyped by high resolution melting (HRM) analysis and PCR-RFLP. Moreover, this gene sequence was mapped by radiation hybrid analysis on porcine chromosome 8. The association study was conducted in 756 performance tested pigs of Italian Large White and Italian Duroc breeds. Significant results were obtained for lean meat content, back fat thickness, visible intermuscular fat and ham weight. In chapter six, a second candidate gene (tribbles homolog 3, TRIB3) is analyzed in a study of association with carcass and meat quality traits. The TRIB3 gene is involved in energy metabolism of skeletal muscle and plays a role as suppressor of adipocyte differentiation. We identified two polymorphisms in the first coding exon of the porcine TRIB3 gene, one is a synonymous SNP (c.132T> C), a second is a missense mutation (c.146C> T, p.P49L). The two polymorphisms appear to be in complete linkage disequilibrium between and within breeds. The in silico analysis of the p.P49L substitution suggests that it might have a functional effect. The association study in about 650 pigs indicates that this marker is associated with back fat thickness in Italian Large White and Italian Duroc breeds in two different experimental designs. This polymorphisms is also associated with lactate content of muscle semimembranosus in Italian Large White pigs. Expression analysis indicated that this gene is transcribed in skeletal muscle and adipose tissue as well as in other tissues. In the seventh chapter, we reported the genotyping results for of 677 SNPs in extreme divergent groups of pigs chosen according to the extreme estimated breeding values for back fat thickness. SNPs were identified by resequencing, literature mining and in silico database mining. analysis, data reported in the literature of 60 candidates genes for obesity. Genotyping was carried out using the GoldenGate (Illumina) platform. Of the analyzed SNPs more that 300 were polymorphic in the genotyped population and had minor allele frequency (MAF) >0.05. Of these SNPs, 65 were associated (P<0.10) with back fat thickness. One of the most significant gene marker was the same TBC1D1 SNPs reported in chapter 5, confirming the role of this gene in fat deposition in pig. These results could be important to better define the pig as a model for human obesity other than for marker assisted selection to improve carcass characteristics.

Molecular and biological analysis on the horizontal transfer of insect transposons in Baculoviruses

Summary During the infection of Lepidoptera larvae with baculoviruses the horizontal escape of Tc1-like transposons, termed TCl4.7 and TCp3.2, from the genome of the host Cryptophlebia leucotreta and Cydia pomonella into the genome of Cydia pomonella granulovirus was observed. In this study we addressed the question whether the transposon harboring viruses had a replication advantage over the wild-type and became dominant in the virus population or whether the activity of the host transposable elements is stimulated by virus infection. Biological characterization studies demonstrated that the transposon containing viruses killed C. pomonella larvae slower than CpGV-M. In co-infection experiments of C. pomonella larvae using a mixture of CpGV-M and mutant viruses as inoculum, it was shown that the transposon carrying mutants had a significant selection disadvantage compared to CpGV-M. Transcription levels of the transposase gene of TCp3.2 were investigated in virus infected and uninfected larvae. These experiments demonstrated that a higher level of transposase transcription was detectable in CpGV-M infected than in mock infected control larvae. This observation gave strong evidence that CpGV-M infection might trigger the activity of transposon TCp3.2 within the genome of Cydia pomonella. Our results suggest that the horizontal transfer of insect host transposons into baculovirus genomes might be induced by virus infection.

Candidate genes affecting fat deposition, carcass composition and meat quality traits in pigs

Pig meat quality is determined by several parameters, such as lipid content, tenderness, water-holding capacity, pH, color and flavor, that affect consumers’ acceptance and technological properties of meat. Carcass quality parameters are important for the production of fresh and dry-cure high-quality products, in particular the fat deposition and the lean cut yield. The identification of genes and markers associated with meat and carcass quality traits is of prime interest, for the possibility of improving the traits by marker-assisted selection (MAS) schemes. Therefore, the aim of this thesis was to investigate seven candidate genes for meat and carcass quality traits in pigs. In particular, we focused on genes belonging to the family of the lipid droplet coat proteins perilipins (PLIN1 and PLIN2) and to the calpain/calpastatin system (CAST, CAPN1, CAPN3, CAPNS1) and on the gene encoding for PPARg-coactivator 1A (PPARGC1A). In general, the candidate genes investigation included the protein localization, the detection of polymorphisms, the association analysis with meat and carcass traits and the analysis of the expression level, in order to assess the involvement of the gene in pork quality. Some of the analyzed genes showed effects on various pork traits that are subject to selection in genetic improvement programs, suggesting a possible involvement of the genes in controlling the traits variability. In particular, significant association results have been obtained for PLIN2, CAST and PPARGC1A genes, that are worthwhile of further validation. The obtained results contribute to a better understanding of biological mechanisms important for pig production as well as for a possible use of pig as animal model for studies regarding obesity in humans.

Molecular strategies for genetic diversity analysis and development of markers linked to resistance traits in apple

The goal of many plant scientists’ research is to explain natural phenotypic variation in term of simple changes in DNA sequence. DNA-based molecular markers are extensively used for the construction of genome-wide molecular maps and to perform genetic analysis for simple and complex traits. The PhD thesis was divided into two main research lines according to the different approaches adopted. The first research line is to analyze the genetic diversity in an Italian apple germplasm collection for the identification of markers tightly linked to targeted genes by an association genetic method. This made it possible to identify synomym and homonym accessions and triploids. The fruit red skin color trait has been used to test the reliability of the genetic approaches in this species. The second line is related to the development of molecular markers closely linked to the Rvi13 and Rvi5 scab resistance genes, previously mapped on apple’s chromosome 10 and 17 respectively by using the traditional linkage mapping method. Both region have been fine-mapped with various type of markers that could be used for marker-assisted selection in future breeding programs and to isolate the two resistance genes.

Analysis of the pig genome for the identification of genomic regions affecting production traits

The aim of this work was to identify markers associated with production traits in the pig genome using different approaches. We focused the attention on Italian Large White pig breed using Genome Wide Association Studies (GWAS) and applying a selective genotyping approach to increase the power of the analyses. Furthermore, we searched the pig genome using Next Generation Sequencing (NSG) Ion Torrent Technology to combine selective genotyping approach and deep sequencing for SNP discovery. Other two studies were carried on with a different approach. Allele frequency changes for SNPs affecting candidate genes and at Genome Wide level were analysed to identify selection signatures driven by selection program during the last 20 years. This approach confirmed that a great number of markers may affect production traits and that they are captured by the classical selection programs. GWAS revealed 123 significant or suggestively significant SNP associated with Back Fat Thickenss and 229 associated with Average Daily Gain. 16 Copy Number Variant Regions resulted more frequent in lean or fat pigs and showed that different copies of those region could have a limited impact on fat. These often appear to be involved in food intake and behavior, beside affecting genes involved in metabolic pathways and their expression. By combining NGS sequencing with selective genotyping approach, new variants where discovered and at least 54 are worth to be analysed in association studies. The study of groups of pigs undergone to stringent selection showed that allele frequency of some loci can drastically change if they are close to traits that are interesting for selection schemes. These approaches could be, in future, integrated in genomic selection plans.

Multiplex network analysis with application to biological high-throughput data

In questa tesi vengono studiate alcune caratteristiche dei network a multiplex; in particolare l'analisi verte sulla quantificazione delle differenze fra i layer del multiplex. Le dissimilarita sono valutate sia osservando le connessioni di singoli nodi in layer diversi, sia stimando le diverse partizioni dei layer. Sono quindi introdotte alcune importanti misure per la caratterizzazione dei multiplex, che vengono poi usate per la costruzione di metodi di community detection . La quantificazione delle differenze tra le partizioni di due layer viene stimata utilizzando una misura di mutua informazione. Viene inoltre approfondito l'uso del test dell'ipergeometrica per la determinazione di nodi sovra-rappresentati in un layer, mostrando l'efficacia del test in funzione della similarita dei layer. Questi metodi per la caratterizzazione delle proprieta dei network a multiplex vengono applicati a dati biologici reali. I dati utilizzati sono stati raccolti dallo studio DILGOM con l'obiettivo di determinare le implicazioni genetiche, trascrittomiche e metaboliche dell'obesita e della sindrome metabolica. Questi dati sono utilizzati dal progetto Mimomics per la determinazione di relazioni fra diverse omiche. Nella tesi sono analizzati i dati metabolici utilizzando un approccio a multiplex network per verificare la presenza di differenze fra le relazioni di composti sanguigni di persone obese e normopeso.

A unifying approach for haplotype analysis of quantitative traits in family-based association studies: Testing and estimating gene-environment interactions with complex exposure variables

We propose robust and e±cient tests and estimators for gene-environment/gene-drug interactions in family-based association studies. The methodology is designed for studies in which haplotypes, quantitative pheno- types and complex exposure/treatment variables are analyzed. Using causal inference methodology, we derive family-based association tests and estimators for the genetic main effects and the interactions. The tests and estimators are robust against population admixture and strati¯cation without requiring adjustment for confounding variables. We illustrate the practical relevance of our approach by an application to a COPD study. The data analysis suggests a gene-environment interaction between a SNP in the Serpine gene and smok- ing status/pack years of smoking that reduces the FEV1 volume by about 0.02 liter per pack year of smoking. Simulation studies show that the pro- posed methodology is su±ciently powered for realistic sample sizes and that it provides valid tests and effect size estimators in the presence of admixture and stratification.