Natural genetic variation in Arabidopsis: tools, traits and prospects for evolutionary ecology.

BACKGROUND: The model plant Arabidopsis thaliana (Arabidopsis) shows a wide range of genetic and trait variation among wild accessions. Because of its unparalleled biological and genomic resources, the potential of Arabidopsis for molecular genetic analysis of this natural variation has increased dramatically in recent years. SCOPE: Advanced genomics has accelerated molecular phylogenetic analysis and gene identification by quantitative trait loci (QTL) mapping and/or association mapping in Arabidopsis. In particular, QTL mapping utilizing natural accessions is now becoming a major strategy of gene isolation, offering an alternative to artificial mutant lines. Furthermore, the genomic information is used by researchers to uncover the signature of natural selection acting on the genes that contribute to phenotypic variation. The evolutionary significance of such genes has been evaluated in traits such as disease resistance and flowering time. However, although molecular hallmarks of selection have been found for the genes in question, a corresponding ecological scenario of adaptive evolution has been difficult to prove. Ecological strategies, including reciprocal transplant experiments and competition experiments, and utilizing near-isogenic lines of alleles of interest will be a powerful tool to measure the relative fitness of phenotypic and/or allelic variants. CONCLUSIONS: As the plant model organism, Arabidopsis provides a wealth of molecular background information for evolutionary genetics. Because genetic diversity between and within Arabidopsis populations is much higher than anticipated, combining this background information with ecological approaches might well establish Arabidopsis as a model organism for plant evolutionary ecology.

Temporal variation in glucocorticoid levels during the resting phase is associated in opposite way with maternal and paternal melanic coloration.

Sex-dependent selection can help maintain sexual dimorphism. When the magnitude of selection exerted on a heritable sex trait differs between the sexes, it may prevent each sex to reach its phenotypic optimum. As a consequence, the benefit of expressing a sex trait to a given value may differ between males and females favouring sex-specific adaptations associated with different values of a sex trait. The level of metabolites regulated by genes that are under sex-dependent selection may therefore covary with the degree of ornamentation differently in the two sexes. We investigated this prediction in the barn owl, a species in which females display on average larger black spots on the plumage than males, a heritable ornament. This melanin-based colour trait is strongly selected in females and weakly counter-selected in males indicating sex-dependent selection. In nestling barn owls, we found that daily variation in baseline corticosterone levels, a key hormone that mediates life history trade-offs, covaries with spot diameter displayed by their biological parents. When their mother displayed larger spots, nestlings had lower corticosterone levels in the morning and higher levels in the evening, whereas the opposite pattern was found with the size of paternal spots. Our study suggests a link between daily regulation of glucocorticoids and sex-dependent selection exerted on sexually dimorphic melanin-based ornaments.

Genetic variation in GIPR influences the glucose and insulin responses to an oral glucose challenge.

Glucose levels 2 h after an oral glucose challenge are a clinical measure of glucose tolerance used in the diagnosis of type 2 diabetes. We report a meta-analysis of nine genome-wide association studies (n = 15,234 nondiabetic individuals) and a follow-up of 29 independent loci (n = 6,958-30,620). We identify variants at the GIPR locus associated with 2-h glucose level (rs10423928, beta (s.e.m.) = 0.09 (0.01) mmol/l per A allele, P = 2.0 x 10(-15)). The GIPR A-allele carriers also showed decreased insulin secretion (n = 22,492; insulinogenic index, P = 1.0 x 10(-17); ratio of insulin to glucose area under the curve, P = 1.3 x 10(-16)) and diminished incretin effect (n = 804; P = 4.3 x 10(-4)). We also identified variants at ADCY5 (rs2877716, P = 4.2 x 10(-16)), VPS13C (rs17271305, P = 4.1 x 10(-8)), GCKR (rs1260326, P = 7.1 x 10(-11)) and TCF7L2 (rs7903146, P = 4.2 x 10(-10)) associated with 2-h glucose. Of the three newly implicated loci (GIPR, ADCY5 and VPS13C), only ADCY5 was found to be associated with type 2 diabetes in collaborating studies (n = 35,869 cases, 89,798 controls, OR = 1.12, 95% CI 1.09-1.15, P = 4.8 x 10(-18)).

The branch-site test of positive selection is surprisingly robust but lacks power under synonymous substitution saturation and variation in GC

Positive selection is widely estimated from protein coding sequence alignments by the nonsynonymous-to-synonymous ratio omega. Increasingly elaborate codon models are used in a likelihood framework for this estimation. Although there is widespread concern about the robustness of the estimation of the omega ratio, more efforts are needed to estimate this robustness, especially in the context of complex models. Here, we focused on the branch-site codon model. We investigated its robustness on a large set of simulated data. First, we investigated the impact of sequence divergence. We found evidence of underestimation of the synonymous substitution rate for values as small as 0.5, with a slight increase in false positives for the branch-site test. When dS increases further, underestimation of dS is worse, but false positives decrease. Interestingly, the detection of true positives follows a similar distribution, with a maximum for intermediary values of dS. Thus, high dS is more of a concern for a loss of power (false negatives) than for false positives of the test. Second, we investigated the impact of GC content. We showed that there is no significant difference of false positives between high GC (up to similar to 80%) and low GC (similar to 30%) genes. Moreover, neither shifts of GC content on a specific branch nor major shifts in GC along the gene sequence generate many false positives. Our results confirm that the branch-site is a very conservative test.

Local adaptation maintains clinal variation in melanin-based coloration of European barn owls (Tyto alba).

Ecological parameters vary in space, and the resulting heterogeneity of selective forces can drive adaptive population divergence. Clinal variation represents a classical model to study the interplay of gene flow and selection in the dynamics of this local adaptation process. Although geographic variation in phenotypic traits in discrete populations could be remainders of past adaptation, maintenance of adaptive clinal variation requires recurrent selection. Clinal variation in genetically determined traits is generally attributed to adaptation of different genotypes to local conditions along an environmental gradient, although it can as well arise from neutral processes. Here, we investigated whether selection accounts for the strong clinal variation observed in a highly heritable pheomelanin-based color trait in the European barn owl by comparing spatial differentiation of color and of neutral genes among populations. Barn owl's coloration varies continuously from white in southwestern Europe to reddish-brown in northeastern Europe. A very low differentiation at neutral genetic markers suggests that substantial gene flow occurs among populations. The persistence of pronounced color differentiation despite this strong gene flow is consistent with the hypothesis that selection is the primary force maintaining color variation among European populations. Therefore, the color cline is most likely the result of local adaptation.

Lipoprotein distribution and biological variation of 24S- and 27-hydroxycholesterol in healthy volunteers.

24S- and 27-hydroxycholesterol are obligatory intermediates of cholesterol catabolism and play an important role in the maintenance of whole-body cholesterol homeostasis. Using an HPLC-MS method for oxysterol quantification, the distribution of esterified and unesterified oxysterols in lipoprotein subfractions as well as the influence of daytime, food intake and menstrual cycle on oxysterol concentrations were investigated in healthy volunteers. Moreover, reference intervals for 24S- and 27-hydroxycholesterol in plasma as well as the corresponding levels for 27-hydroxycholesterol in the HDL subfraction were established in 100 healthy volunteers. Both circulating oxysterols are mainly transported in association with HDL and LDL--primarily in the esterified form. No significant diurnal changes and no variations during menstrual cycle of either absolute or cholesterol-related plasma levels were detected. In contrast to 24S-hydroxycholesterol in plasma and 27-hydroxycholesterol in the HDL subfraction, the 95% reference intervals of 27-hydroxycholesterol both in plasma and the non-HDL subfraction were higher in males than in females. The concentrations of 27-hydroxycholesterol in plasma and the non-HDL subfraction showed strong positive correlations with the concentrations of cholesterol, non-HDL cholesterol and triglycerides. Our data on the lipoprotein distribution of oxysterols as well as on their intra- and inter-individual variation set the stage for future clinical studies.

Pharmacogenetic Study on Antidementia Drugs

Differences in efficacy and safety of drugs among patients are a recognized problem in pharmacotherapy. The reasons are multifactorial and, therefore, the choice of a drug and its dosage for a particular patient based on different clinical and genetic factors is suggested to improve the clinical outcome. Four drugs are currently used for the treatment of Alzheimer's disease: three acetylcholinesterase inhibitors (donepezil, galantamine, rivastigmine) and the N-methyl-D-aspartate-antagonist memantine. For these drugs, a high interindividual variability in plasma levels was observed, which might influence the response to treatment. The main objective of this thesis was to provide a better understanding of clinical and genetic factors affecting the plasma levels of antidementia drugs. Furthermore, the relationship between plasma levels, genetic variations and side effects was assessed. For this purpose, a pharmacogenetic study was conducted including 300 patients from a naturalistic clinical setting. Analytical methods for the simultaneous measurement of antidementia drugs in plasma have been developed and validated using liquid chromatography methods coupled with mass spectrometry detection. Presently, these methods are used in the therapeutic drug monitoring service of our laboratory. The routine use of therapeutic drug monitoring for antidementia drugs cannot yet be recommended with the available data, but it may be beneficial for some patients in special clinical cases such as insufficient treatment response, side effects or drug interactions. Donepezil and galantamine are extensively metabolized by the liver enzymes cytochromes P450 (CYP) 2D6 and 3A and are substrates of the drug transporter P-glycoprotein. The relationship of variations in genes affecting the activity of these metabolic enzymes and drug transporter (CYP2D6, CYP3A, POR, NR1I2, ABCB1) with donepezil and galantamine plasma levels was investigated. The CYP2D6 genotype appeared to be the major genetic factor involved in the pharmacokinetics of these two drugs. Thus, CYP2D6 poor metabolizers demonstrated significantly higher drug plasma levels than extensive metabolizers. Additionally, in the donepezil study population, the frequency of side effects was significantly increased in poor metabolizers. Lower donepezil plasma levels were observed in ultra rapid metabolizers, which might expose those patients to the risk of non-response. Memantine is mainly eliminated unchanged by the kidney, with implication of tubular secretion by renal transporters. A population pharmacokinetic model was developed to quantify the effects of clinical factors and genetic variations in renal cation transporters (SLC22A1/2/5, SLC47A1, ABCB1), and nuclear receptors (NR1I2, NR1I3, PPARG) involved in transporter expression, on memantine plasma levels. In addition to the renal function and gender, a genetic variation in the nuclear receptor Pregnane-X-Receptor (NR1I2) significantly affected memantine elimination. These findings suggest that an individualized therapy approach for antidementia drugs, taking into account clinical characteristics and genetic background of a patient, might increase efficacy and safety of the treatment. - Les différences interindividuelles dans l'efficacité et la tolérance des médicaments sont un problème connu en pharmacothérapie. Les raisons sont multiples, et le choix du médicament et de la dose, basé sur des facteurs cliniques et génétiques spécifiques au patient, peut contribuer à améliorer la réponse clinique. Quatre médicaments sont couramment utilisés dans le traitement de la maladie d'Alzheimer : trois inhibiteurs de l'acétylcholinestérase (donépézil, galantamine, rivastigmine) et un antagoniste du récepteur N-méthyl-D-aspartate, la mémantine. Une forte variabilité interindividuelle dans les taux plasmatiques de ces quatre composés a été observée, ce qui pourrait influencer la réponse au traitement. L'objectif principal de ce travail de thèse est de mieux comprendre les facteurs cliniques et génétiques influençant les taux des médicaments pro-cognitifs. En outre, des associations entre les taux, la variabilité génétique et les effets secondaires ont été recherchées. Dans ce but, 300 patients sous traitement avec un médicament pro-cognitif ont été recrutés pour une étude pharmacogénétique. Des méthodes de dosage simultané de médicaments pro-cognitifs par chromatographie liquide couplée à la spectrométrie de masse ont été développées et validées. Ces méthodes sont actuellement utilisées dans le service de suivi thérapeutique de notre unité. Malgré le fait qu'un suivi des taux sanguins des pro-cognitifs ne puisse pas encore être recommandé en routine, un dosage peut être utile dans des cas cliniques spécifiques, comme une réponse insuffisante, une intolérance ou une interaction médicamenteuse. Le donépézil et la galantamine sont fortement métabolisés par les cytochromes P450 (CYP) 2D6 et 3A, et sont également substrats du transporteur P-glycoprotéine. Les associations entre les polymorphismes génétiques de ces enzymes, cofacteur, récepteur nucléaire et transporteur (CYP2D6, CYP3A, POR, NR1I2, ABCB1) et les taux de donépézil et de galantamine ont été étudiées. Le génotype du CYP2D6 a été montré comme le facteur génétique majeur impliqué dans la pharmacocinétique de ces deux médicaments. Ainsi, les métaboliseurs déficients du CYP2D6 ont démontré des taux plasmatiques significativement plus élevés comparé aux bons métaboliseurs. De plus, dans la population traitée avec le donépézil, la fréquence des effets secondaires était plus élevée chez les métaboliseurs déficients. Des taux plasmatiques bas ont été mesurés chez les métaboliseurs ultra-rapides traités avec le donépézil, ce qui pourrait être un facteur de risque à une non-réponse au traitement. La mémantine est principalement éliminée sous forme inchangée par les reins, et partiellement par sécrétion tubulaire grâce à des transporteurs rénaux. Un modèle de cinétique de population a été développé pour quantifier les effets des différents facteurs cliniques et de la variabilité génétique des transporteurs rénaux (SLC22A1/2/5, SLC47A1, ABCB1) et des récepteurs nucléaires (NR1I2, NR1I3, PPARG, impliqués dans l'expression des transporteurs) sur les taux plasmatiques de mémantine. En plus de la fonction rénale et du genre, une variation génétique dans le récepteur nucléaire Pregnane-X-Receptor (NR1I2) a montré une influence significative sur l'élimination de la mémantine. Ces résultats suggèrent qu'une approche thérapeutique individualisée, prenant en compte des facteurs cliniques et génétiques du patient, pourrait améliorer l'efficacité et la sécurité du traitement pro-cognitif.

Long sequence duplications, repeats, and palindromes in HIV-1 gp120: length variation in V4 as the product of misalignment mechanism.

We have shown that indels in gp120 V4 are associated to the presence of duplicated and palindromic sequences, suggesting that they may be produced by strand-slippage misalignment mechanism. Indels in V4 involved region-specific duplications 9 to 15 bp long, and repeats of various lengths, associated to trinucleotides AAT. No duplications were found in V3 and C3. The frequency of palindromic sequences in individual genes was found to be significantly higher in gp120 (p < or = 3.00E-7), and significantly lower in Tat (p < or = 9.00E-7) than the average frequency calculated over the full genome. The finding of elements of misalignment in association with indels in V4 suggests that these mutations may occur in proviral DNA after integration of HIV into the host genome. It also implies that occurrence of large indels in gp120 is not random but is directed by the presence and distribution of elements of misalignment in the HIV genome.

Pheomelanin-based coloration and the ability to cope with variation in food supply and parasitism.

Although gene by environment interactions may play a key role in the maintenance of genetic polymorphisms, little is known about the ecological factors involved in these interactions. We investigated whether food supply and parasites can mediate covariation between the degree of adult pheomelanin-based coloration, a heritable trait, and offspring body mass in the tawny owl (Strix aluco). We swapped clutches between nests to allocate genotypes randomly among environments. Three weeks after hatching, we challenged the immune system of 80 unrelated nestlings with either a phytohemagglutinin (PHA) or a lipopolysaccharide, surrogates of alternative parasites, and then fed them ad lib. or food-restricted them during the following 6 days in the laboratory. Whatever the immune challenge, nestlings fed ad lib. converted food more efficiently into body mass when their biological mother was dark pheomelanic. In contrast, food-restricted nestlings challenged with PHA lost less body mass when their biological mother was pale pheomelanic. Nestling tawny owls born from differently melanic mothers thus show differing reaction norms relative to food availability and parasitism. This suggests that dark and pale pheomelanic owls reflect alternative adaptations to food availability and parasites, factors known to vary in space and time.

Novel method to estimate the phenotypic variation explained by genome-wide association studies reveals large fraction of the missing heritability.

Genome-wide association studies (GWAS) are conducted with the promise to discover novel genetic variants associated with diverse traits. For most traits, associated markers individually explain just a modest fraction of the phenotypic variation, but their number can well be in the hundreds. We developed a maximum likelihood method that allows us to infer the distribution of associated variants even when many of them were missed by chance. Compared to previous approaches, the novelty of our method is that it (a) does not require having an independent (unbiased) estimate of the effect sizes; (b) makes use of the complete distribution of P-values while allowing for the false discovery rate; (c) takes into account allelic heterogeneity and the SNP pruning strategy. We applied our method to the latest GWAS meta-analysis results of the GIANT consortium. It revealed that while the explained variance of genome-wide (GW) significant SNPs is around 1% for waist-hip ratio (WHR), the observed P-values provide evidence for the existence of variants explaining 10% (CI=[8.5-11.5%]) of the phenotypic variance in total. Similarly, the total explained variance likely to exist for height is estimated to be 29% (CI=[28-30%]), three times higher than what the observed GW significant SNPs give rise to. This methodology also enables us to predict the benefit of future GWA studies that aim to reveal more associated genetic markers via increased sample size.

Natural gene-expression variation in Down syndrome modulates the outcome of gene-dosage imbalance.

Down syndrome (DS) is characterized by extensive phenotypic variability, with most traits occurring in only a fraction of affected individuals. Substantial gene-expression variation is present among unaffected individuals, and this variation has a strong genetic component. Since DS is caused by genomic-dosage imbalance, we hypothesize that gene-expression variation of human chromosome 21 (HSA21) genes in individuals with DS has an impact on the phenotypic variability among affected individuals. We studied gene-expression variation in 14 lymphoblastoid and 17 fibroblast cell lines from individuals with DS and an equal number of controls. Gene expression was assayed using quantitative real-time polymerase chain reaction on 100 and 106 HSA21 genes and 23 and 26 non-HSA21 genes in lymphoblastoid and fibroblast cell lines, respectively. Surprisingly, only 39% and 62% of HSA21 genes in lymphoblastoid and fibroblast cells, respectively, showed a statistically significant difference between DS and normal samples, although the average up-regulation of HSA21 genes was close to the expected 1.5-fold in both cell types. Gene-expression variation in DS and normal samples was evaluated using the Kolmogorov-Smirnov test. According to the degree of overlap in expression levels, we classified all genes into 3 groups: (A) nonoverlapping, (B) partially overlapping, and (C) extensively overlapping expression distributions between normal and DS samples. We hypothesize that, in each cell type, group A genes are the most dosage sensitive and are most likely involved in the constant DS traits, group B genes might be involved in variable DS traits, and group C genes are not dosage sensitive and are least likely to participate in DS pathological phenotypes. This study provides the first extensive data set on HSA21 gene-expression variation in DS and underscores its role in modulating the outcome of gene-dosage imbalance.

Length variation in HIV-1 gp120 as the product of DNA misalignment mechanism

Background: To determine whether misalignment structures such as duplications, repeats, and palindromes are associated to insertions/deletions (indels) in gp120, indicating that indels are indeed frameshift mutations generated by DNA misalignment mechanism. Methods: Cloning and sequencing of a fragment of HIV-1 gp120 spanning C2-C4 derived from plasma RNA in 12 patients with early chronic disease and naïve to antiretroviral therapy. Results: Indels in V4 involved always insertion and deletion of duplicated nucleotide segments, and AAT repeats, and were associated to the presence of palindromic sequences. No duplications were detected in V3 and C3. Palindromic sequences occurred with similar frequencies in V3, C3 and V4; the frequency of palindromes in individual genes was found to be significantly higher in structural (gp120, p ≤ 3.00E-7) and significantly lower in regulatory (Tat, p ≤ 9.00E-7) genes, as compared to the average frequency calculated over the full genome. Discussion: Indels in V4 are associated to misalignment structures (i.e. duplications repeat and palindromes) indicating DNA misalignment as the mechanism underlying length variation in V4. The finding that indels in V4 are caused by DNA misalignment has some very important implications: 1) indels in V4 are likely to occur in proviral DNA (and not in RNA), after integration of HIV into the host genome; 2) they are likely to occur as progressive modifications of the early founder virus during chronic infection, as more and more cells get infected; 3) frameshift mutations involving any number of base pairs are likely to occur evenly across gp120; however, only those mutants carrying a functional gp120 (indels as multiples of three base pairs) will be able to perpetuate the virus cycle and to keep spreading through the population.

Pharmacogenetics of Antiviral Treatment

RESUME - FRANÇAISRésuméDans ce travail de thèse, l'importance de la pharmacogénétique des traitements antiviraux a été évaluée en déterminant, au moyen de trois différentes approches, l'impact de variations génétiques sur la pharmacocinétique de deux traitements antirétroviraux (à savoir l'efavirenz et le lopinavir) ainsi que sur la capacité de pouvoir éliminer le virus de l'hépatite C de façon naturelle ou suite à un traitement médicamenteux.L'influence des variations génétiques sur les taux plasmatiques de l'efavirenz et de ses métabolites primaires a été évaluée par l'analyse d'un seul gène candidat : le cytochrome P450 (CYP) 2A6, impliqué dans une voie métabolique accessoire de l'efavirenz. Cette étude a permis de démontrer que le génotype du CYP2A6 devient cliniquement déterminant en l'absence de fonction du CYP2B6, impliqué dans la voie métabolique principale, et que la perte simultanée des voies métaboliques principales et accessoires entraine une augmen¬tation du risque d'interruption du traitement, soulignant la valeur prédictive du génotypage.L'influence de la génétique sur la clairance du lopinavir a été évaluée par l'analyse à grande échelle de gènes candidats, à savoir les gènes potentiellement impliqués dans l'absorption, le métabolisme, la distribution et l'élimination d'un médicament. Cette étude a permis l'identification de 4 polymorphismes, dans des transporteurs et des enzymes métaboliques, associés à la clairance du lopinavir et expliquant 5% de la variabilité inter¬individuelle de ce phénotype.L'influence de la génétique sur la capacité d'éliminer le virus de l'hépatite C, de façon naturelle ou à la suite d'un traitement, a été évaluée par l'analyse du génome entier. Cette étude a permis l'identification d'un polymorphisme situé à proximité de l'interféron-X3. Quatre variations génétiques potentiellement causales ont ensuite pu être identifiées par reséquencage. Finalement, la contribution nette de ce gène sur l'élimination du virus a pu être évaluée dans une cohorte infectée par une seule et même source, permettant ainsi de contrôler l'effet de la diversité virale, du genre et de la présence de co-infections.Cette thèse a permis de mettre en évidence les diverses méthodes disponibles pour la recherche en pharmacogénétique, ainsi que l'importance du reséquencage pour l'identification de variations génétiques causales.SUMMARY - ENGLISHSummaryIn this thesis work the relevance of pharmacogenetics of antiviral treatment has been assessed by investigating, through three different approaches, the impact of host genetic variation on antiretroviral drug disposition (namely efavirenz and lopinavir) and on natural or treatment-induced clearance of hepatitis C virus.The influence of host genetic variation on efavirenz and its primary metabolite plasma levels was assessed by single candidate gene approach, through comprehensive analysis of cytochrome P450 (CYP) 2A6 - involved in efavirenz accessory metabolic pathway. The study could demonstrate that CYP2A6 genotype became increasingly relevant in the setting of limited CYP2B6 function - involved in efavirenz main metabolic pathway - and that individuals with both main and accessory metabolic pathways impaired were at higher risk for treatment discontinuation, overall emphasizing the predictive power of genotyping.The influence of host genetic variation on lopinavir clearance was assessed by large scale candidate gene approach, through analysis of genes involved in the absorption, distribution, metabolism and elimination. The study identified four genetic variants in drug transporters and metabolizing enzymes that explained 5% of the interindividual variability in lopinavir clearance.The influence of host genetic variation on hepatitis C virus (HCV) natural or treatment- induced clearance was assessed through genome-wide association study approach. This study identified an intergenic polymorphism, part of a linkage disequilibrium block encompassing the interferon-3 gene, as highly associated with treatment-induced and spontaneous HCV clearance. Resequencing and recombinant mapping lead to the identification of four potentially causal genetic variants. Finally, we could assess the net contribution of genetic variants in interferon-3 to clearance by controlling for viral diversity, gender and co-infection status in a single source infected cohort.This thesis highlights the various genetic tools available to pharmacogenetic discovery (candidate gene, pathway or and genome-wide approaches), and the importance of resequencing for mapping of causal variants.

Copy number variation modifies expression time courses.

A preliminary understanding into the phenotypic effect of DNA segment copy number variation (CNV) is emerging. These rearrangements were demonstrated to influence, in a somewhat dose-dependent manner, the expression of genes that map within them. They were also shown to modify the expression of genes located on their flanks and sometimes those at a great distance from their boundary. Here we demonstrate, by monitoring these effects at multiple life stages, that these controls over expression are effective throughout mouse development. Similarly, we observe that the more specific spatial expression patterns of CNV genes are maintained through life. However, we find that some brain-expressed genes mapping within CNVs appear to be under compensatory loops only at specific time points, indicating that the effect of CNVs on these genes is modulated during development. Notably, we also observe that CNV genes are significantly enriched within transcripts that show variable time courses of expression between strains. Thus, modifying the copy number of a gene may potentially alter not only its expression level, but also the timing of its expression.