The genomic substrate for adaptive radiation in African cichlid fish

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

Activation of common signaling pathways during remodeling of the heart and the bladder

The heart and the urinary bladder are hollow muscular organs, which can be afflicted by pressure overload injury due to pathological conditions such as hypertension and bladder outlet obstruction. This increased outflow resistance induces hypertrophy, marked by dramatic changes in the organs' phenotype and function. The end result in both the heart and the bladder can be acute organ failure due to advanced fibrosis and the subsequent loss of contractility. There is emerging evidence that microRNAs (miRNAs) play an important role in the pathogenesis of heart failure and bladder dysfunction. MiRNAs are endogenous non-coding single-stranded RNAs, which regulate gene expression and control adaptive and maladaptive organ remodeling processes. This Review summarizes the current knowledge of molecular alterations in the heart and the bladder and highlights common signaling pathways and regulatory events. The miRNA expression analysis and experimental target validation done in the heart provide a valuable source of information for investigators working on the bladder and other organs undergoing the process of fibrotic remodeling. Aberrantly expressed miRNA are amendable to pharmacological manipulation, offering an opportunity for development of new therapies for cardiac and bladder hypertrophy and failure.

Regulation of two P450 isoforms: CYP2D6 and CYP4F11

Cytochromes P450 catalyze a monooxygenase reaction in which molecular oxygen is split and one oxygen atom is incorporated into the substrate. As a whole, P450 researchers have focused most of their attention on substrate metabolism and relatively little on how these enzymes are regulated. This study will focus on the regulation of two P450 isoforms known as, CYP2D6 and CYP4F11. ^ The human CYP2D gene locus contains two pseudogenes and one functional gene known as CYP2D6. This locus is highly polymorphic and produces several alternatively spliced transcripts from the pseudogene CYP2D7. My objective was to understand the role of SV5-in (splice variant 5), one of several alternative splice variants transcribed from the CYP2D7 pseudogene. My results indicate that SV5-in mRNA causes an increase in CYP2D6 protein levels and suggest that there is a role for SV5-in in regulation of CYP2D6 expression. ^ Second, CYP4F11 is a recently discovered and uncharacterized isoform, derived from the CYP4F subfamily. It metabolizes several clinically relevant drugs (i.e.—erythromycin and benzphetamine) and some endogenous inflammatory mediators (i.e.—LTB4). After evaluation of microarray data, I observed an increase in CYP4F11 mRNA levels from wild-type HCT116 cells compared to p53-null cells. Our objectives were to explore and understand this connection between p53 and CYP4F11. Microarray data were confirmed by Q-PCR, after which this effect was again observed at the protein level via Western blot and again at the promoter level via luciferase assay and chromatin immunoprecipitation. Our results indicate that p53 protein regulates expression of CYP4F11 mRNA and protein through CYP4F11 promoter binding (note that p53 binding to CYP4F11 DNA was not shown to be direct). These results signify a whole new level of regulation of drug metabolizing enzymes by p53. ^ An understanding of CYP4F11 regulation by p53 could help us understand another pathway leading to apoptosis or cell growth arrest. This can aid future drug studies and discover new drug metabolism pathways under the control of a tumor suppressor protein. An understanding of the CYP2D6 regulation pathway could illuminate the role of non-coding RNAs in the P450 field and potentially explain several inter-individual drug response variations observed in clinical medicine that are not yet completely explained by genotyping analysis. ^

Association of folate metabolic pathway genes with human spina bifida meningomyelocele

Neural tube defects including spina bifida meningomyelocele (SBMM) are common malformations of the brain and spinal cord, and include all abnormalities resulting from lack of closure of the developing neural tube during embryological development.^ The specific aims of this study were to determine if single nucleotide polymorphic variants (SNPs) in the folate/homocysteine metabolic pathway genes confer a risk for NTD susceptibility within this SBMM population.^ In completion of the first specific aim, two novel SNPs were identified in the FOLR1 gene in Chromosome 11of patients including one in non-coding exon 1 with a C → T transition at nucleotide position 71578317 and another in non-coding exon 3 with a T → G transversion at nucleotide position 71579123. It will be important to determine if these variants are present in the respective parents of these individuals. If they are in fact de novo variants, then these SNPs may be more likely to contribute to the birth defect.^ The second project aim was to analyze genotypes associated with SBMM risk by transmission disequilibrium tests (TDT) and association was detected on several SNPs across the folate metabolic pathway genes in this population. SNPs with significant RC-TDT values were found within the DHFR gene (rs1650723), the MTRR gene (rs327592), the FOLR2 gene (rs13908), four tightly linked variants in the FOLR3 gene (rs7925545, rs7926875, rs7926987, rs7926360) and a variant in the SLC19A1 gene (rs1888530). The product of each of these genes performs a vital function in the folate metabolic pathway. It is conceivable, therefore, that if the individual SNP or SNPs can be proven to perturb the function in some way that they may be involved in the disruption of folate metabolism and in the resulting birth defect. Validating the results of this study in other independent populations will further strengthen the evidence that dysfunction of folate enzymes and receptors may confer SBMM risk in humans. ^

Arabidopsis thaliana bZIP44: a transcription factor affecting seed germination and expression of the mannanase-encoding gene AtMAN7.

Endo-β-mannanases (MAN; EC. 3.2.1.78) catalyze the cleavage of β1[RIGHTWARDS ARROW]4 bonds in mannan polymers and have been associated with the process of weakening the tissues surrounding the embryo during seed germination. In germinating Arabidopsis thaliana seeds, the most highly expressed MAN gene is AtMAN7 and its transcripts are restricted to the micropylar endosperm and to the radicle tip just before radicle emergence. Mutants with a T-DNA insertion in AtMAN7 have a slower germination than the wild type. To gain insight into the transcriptional regulation of the AtMAN7 gene, a bioinformatic search for conserved non-coding cis-elements (phylogenetic shadowing) within the Brassicaceae MAN7 gene promoters has been done, and these conserved motifs have been used as bait to look for their interacting transcription factors (TFs), using as a prey an arrayed yeast library from A. thaliana. The basic-leucine zipper TF AtbZIP44, but not the closely related AtbZIP11, has thus been identified and its transcriptional activation upon AtMAN7 has been validated at the molecular level. In the knock-out lines of AtbZIP44, not only is the expression of the AtMAN7 gene drastically reduced, but these mutants have a significantly slower germination than the wild type, being affected in the two phases of the germination process, both in the rupture of the seed coat and in the breakage of the micropylar endosperm cell walls. In the over-expression lines the opposite phenotype is observed.

Transcription factors of the bzip factors of the BZIP family affecting arabidopsis thaliana germination sensu stricto

During seed germination, the endosperm cell walls (CWs) suffer an important weakening process mainly driven by hydrolytic enzymes, such are endo-?- mannanases (MAN; EC. 3.2.1.78) that catalyze the cleavage of ?1?4 bonds in the mannan-polymers. In Arabidopsis thaliana seeds, endo-?-mannanase activity increases during seed imbibition, decreasing after radicle emergence1. AtMAN7 is the most highly expressed MAN gene in seeds upon germination and their transcripts are restricted to the micropylar endosperm and to the radicle tip just before radicle emergence. Mutants with a T-DNA insertion in this gene (K.O. MAN7) have a slower germination rate than the wild type (t50=34 h versus t50=25 h). To gain insight into the transcriptional regulation of the AtMAN7 gene, a bioinformatic search for conserved non-coding cis-elements (phylogenetic shadowing) within the Brassicaceae orthologous MAN7 gene promoters has been done and these conserved motives have been used as baits to look for their interacting transcription factors (TFs), using as a prey an arrayed yeast library of circa 1,200 TFs from A. thaliana. The basic leucine zipper AtbZIP44, but not its closely related ortholog AtbZIP11, has been thus identified and its regulatory function upon AtMAN7 during seed germination validated by different molecular and physiological techniques, such are RT-qPCR analyses, mRNA Fluorescence in situ Hybridization (FISH) experiments, and by the establishment of the germination kinetics of both over-expression (oex) lines and TDNA insertion mutants in AtbZIP44. The transcriptional combinatorial network through which AtbZIP44 regulates AtMAN7 gene expression during seed germination has been further explored through protein-protein interactions between AtbZIP44 and other bZIP members. In such a way, AtbZIP9 has been identified by yeast two-hybrid experiments and its physiological implication in the control of AtMAN7 expression similarly established.

Nucleotide sequence of the Kaposi sarcoma-associated herpesvirus (HHV8)

The genome of the Kaposi sarcoma-associated herpesvirus (KSHV or HHV8) was mapped with cosmid and phage genomic libraries from the BC-1 cell line. Its nucleotide sequence was determined except for a 3-kb region at the right end of the genome that was refractory to cloning. The BC-1 KSHV genome consists of a 140.5-kb-long unique coding region flanked by multiple G+C-rich 801-bp terminal repeat sequences. A genomic duplication that apparently arose in the parental tumor is present in this cell culture-derived strain. At least 81 ORFs, including 66 with homology to herpesvirus saimiri ORFs, and 5 internal repeat regions are present in the long unique region. The virus encodes homologs to complement-binding proteins, three cytokines (two macrophage inflammatory proteins and interleukin 6), dihydrofolate reductase, bcl-2, interferon regulatory factors, interleukin 8 receptor, neural cell adhesion molecule-like adhesin, and a D-type cyclin, as well as viral structural and metabolic proteins. Terminal repeat analysis of virus DNA from a KS lesion suggests a monoclonal expansion of KSHV in the KS tumor.

A host-specific function is required for ligation of a wide variety of ribozyme-processed RNAs

Hepatitis δ virus (HDV) replicates its circular RNA genome via a rolling circle mechanism. During this process, cis-acting ribozymes cleave adjacent upstream sequences and thereby resolve replication intermediates to unit-length RNA. The subsequent ligation of these 5′OH and 2′,3′-cyclic phosphate termini to form circular RNA is an essential step in the life cycle of the virus. Here we present evidence for the involvement of a host activity in the ligation of HDV RNA. We used both HDV and hammerhead ribozymes to generate a panel of HDV and non-HDV RNA substrates that bear 5′ hydroxyl and 2′,3′- cyclic phosphate termini. We found that ligation of these substrates occurred in host cells, but not in vitro or in Escherichia coli. The host-specific ligation activity was capable of joining RNA in both bimolecular and intramolecular reactions and functioned in a sequence-independent manner. We conclude that mammalian cells contain a default pathway that efficiently circularizes ribozyme processed RNAs. This pathway could be exploited in the delivery of stable antisense and decoy RNA to the nucleus.

Directional cDNA library construction assisted by the in vitro recombination reaction

We report here a new directional cDNA library construction method using an in vitro site-specific recombination reaction, based on the integrase–excisionase system of bacteriophage λ. Preliminary experiments revealed that in vitro recombinational cloning (RC) provided important advantages over conventional ligation-assisted cloning: it eliminated restriction digestion for directional cloning, generated low levels of chimeric clones, reduced size bias and, in our hands, gave a higher cloning efficiency than conventional ligation reactions. In a cDNA cloning experiment using an in vitro synthesized long poly(A)+ RNA (7.8 kb), the RC gave a higher full-length cDNA clone content and about 10 times more transformants than conventional ligation-assisted cloning. Furthermore, characterization of rat brain cDNA clones yielded by the RC method showed that the frequency of cDNA clones >2 kb having internal NotI sites was ∼6%, whereas these cDNAs could not be cloned at all or could be isolated only in a truncated form by conventional methods. Taken together, these results indicate that the RC method makes it possible to prepare cDNA libraries better representing the entire population of cDNAs, without sacrificing the simplicity of current conventional ligation-assisted methods.

Regulation of c-maf gene expression by Pax6 in cultured cells

c-Maf is a bZip transcription factor expressed in developmental and cellular differentiation processes. Recently, a c-maf knockout mouse model, showing abnormal lens development, has been reported. In order to study the regulation mechanisms of c-maf gene expression during the differentiation process we have cloned and functionally characterized the rat c-maf (maf-2) gene. The rat c-maf gene is an intronless gene, covering a length of 3.5 kb. Transient transfection analysis of the 5′-flanking region of the c-maf gene using luciferase as the reporter gene shows that Pax6, a master transcription factor for lens development, strongly activates the c-maf promoter construct. Endogenous c-maf is also activated by the Pax6 expression vector. Electrophoresis mobility shift assay and DNase I footprinting analysis show that at least three Pax6-binding sites are located in the 5′-flanking and 5′-non-coding regions of the rat c-maf gene. The c-maf gene was also markedly activated by its own product, c-Maf, through the MARE (Maf recognition element), suggesting that a positive autoregulatory mechanism controls this gene. In situ hybridization histochemical detection of Pax6 and c-Maf in the E14 lens showed that both mRNAs are expressed in the lens equator where lens epithelial cells are differentiating to lens fiber cells. These results suggest that a Pax6/c-Maf transcription factor cascade is working in lens development.

Ruolo di Clusterina nella progressione del tumore prostatico

Clusterina (CLU) è una proteina ubiquitaria, presente nella maggior parte dei fluidi corporei e implicata in svariati processi fisiologici. Dalla sua scoperta fino ad oggi, CLU è risultata essere una proteina enigmatica, la cui funzione non è ancora stata compresa appieno. Il gene codifica per 3 varianti trascrizionali identificate nel database NCBI con i codici: NM_001831 (CLU 1 in questo lavoro di tesi), NR_038335 (CLU 2 in questo lavoro di tesi) e NR_045494 (CLU 3 in questo lavoro di tesi). Tutte le varianti sono trascritte come pre-mRNA contenenti 9 esoni e 8 introni e si differenziano per l’esone 1, la cui sequenza è unica e caratteristica di ogni variante. Sebbene in NCBI sia annotato che le varianti CLU 2 e CLU 3 non sono codificanti, tramite analisi bioinformatica è stato predetto che da tutti e tre i trascritti possono generarsi proteine di differente lunghezza e localizzazione cellulare. Tra tutte le forme proteiche ipotizzate, l’unica a essere stata isolata e sequenziata è quella tradotta dall’AUG presente sull’esone 2 che dà origine a una proteina di 449 aminoacidi. Il processo di maturazione prevede la formazione di un precursore citoplasmatico (psCLU) che subisce modificazioni post-traduzionali tra cui formazione di ponti disolfuro, glicosilazioni, taglio in due catene denominate β e α prima di essere secreta come eterodimero βα (sCLU) nell’ambiente extracellulare, dove esercita la sua funzione di chaperone ATP-indipendente. Oltre alla forma extracellulare, è possibile osservare una forma intracellulare con localizzazione citosolica la cui funzione non è stata ancora completamente chiarita. Questo lavoro di tesi si è prefissato lo scopo di incrementare le conoscenze in merito ai trascritti CLU 1 e CLU 2 e alla loro regolazione, oltre ad approfondire il ruolo della forma citosolica della proteina in relazione al signaling di NF-kB che svolge un ruolo importante nel processo di sviluppo e metastatizzazione del tumore. Nella prima parte, uno screening di differenti linee cellulari, quali cellule epiteliali di prostata e di mammella, sia normali sia tumorali, fibroblasti di origine polmonare e linfociti di tumore non-Hodgkin, ha permesso di caratterizzare i trascritti CLU 1 e CLU 2. Dall’analisi è emerso che la sequenza di CLU 1 è più corta al 5’ rispetto a quella depositata in NCBI con l’identificativo NM_001831 e il primo AUG disponibile per l’inizio della traduzione è localizzato sull’esone 2. È stato dimostrato che CLU 2, al contrario di quanto riportato in NCBI, è tradotto in proteina a partire dall’AUG presente sull’esone 2, allo stesso modo in cui viene tradotto CLU 1. Inoltre, è stato osservato che i livelli d’espressione dei trascritti variano notevolmente tra le diverse linee cellulari e nelle cellule epiteliali CLU 2 è espressa sempre a bassi livelli. In queste cellule, l’espressione di CLU 2 è silenziata per via epigenetica e la somministrazione di farmaci capaci di rendere la cromatina più accessibile, quali tricostatina A e 5-aza-2’-deossicitidina, è in grado di incrementarne l’espressione. Nella seconda parte, un’analisi bioinformatica seguita da saggi di attività in vitro in cellule epiteliali prostatiche trattate con farmaci epigenetici, hanno permesso di identificare, per la prima volta in uomo, una seconda regione regolatrice denominata P2, capace di controllare l’espressione di CLU 2. Rispetto a P1, il classico promotore di CLU già ampiamente studiato da altri gruppi di ricerca, P2 è un promotore debole, privo di TATA box, che nelle cellule epiteliali prostatiche è silente in condizioni basali e la cui attività incrementa in seguito alla somministrazione di farmaci epigenetici capaci di alterare le modificazioni post-traduzionali delle code istoniche nell’intorno di P2. Ne consegue un rilassamento della cromatina e un successivo aumento di trascrizione di CLU 2. La presenza di un’isola CpG differentemente metilata nell’intorno di P1 spiegherebbe, almeno in parte, i differenti livelli di espressione di CLU che si osservano tra le diverse linee cellulari. Nella terza parte, l’analisi del pathway di NF-kB in un modello sperimentale di tumore prostatico in cui CLU è stata silenziata o sovraespressa, ha permesso di capire come la forma citosolica di CLU abbia un ruolo inibitorio nei confronti dell’attività del fattore trascrizionale NF-kB. CLU inibisce la fosforilazione e l’attivazione di p65, il membro più rappresentativo della famiglia NF-kB, con conseguente riduzione della trascrizione di alcuni geni da esso regolati e coinvolti nel rimodellamento della matrice extracellulare, quali l’urochinasi attivatrice del plasminogeno, la catepsina B e la metallo proteinasi 9. È stato dimostrato che tale inibizione non è dovuta a un’interazione fisica diretta tra CLU e p65, per cui si suppone che CLU interagisca con uno dei componenti più a monte della via di segnalazione responsabile della fosforilazione ed attivazione di p65.

Análise do papel da via miR156/SQUAMOSA Promoter-Binding Protein-Like (SPL) na organogênese in vitro a partir de raízes de Arabidopsis thaliana

Os microRNAs (miRNAs) são pequenos RNAs endógenos não codantes de 21-24 nucleotídeos (nt) que regulam a expressão gênica de genes-alvos. Eles estão envolvidos em diversos aspectos de desenvolvimento da planta, tanto na parte aérea, quanto no sistema radicular. Entre os miRNAs, o miRNA156 (miR156) regula a família de fatores de transcrição SQUAMOSA Promoter-Binding Protein-Like (SPL) afetando diferentes processos do desenvolvimento vegetal. Estudos recentes mostram que a via gênica miR156/SPL apresenta efeito positivo tanto no aumento da formação de raízes laterais, quanto no aumento de regeneração de brotos in vitro a partir de folhas e hipocótilos em Arabidopsis thaliana. Devido ao fato de que a origem da formação de raiz lateral e a regeneração in vitro de brotos a partir de raiz principal compartilham semelhanças anatômicas e moleculares, avaliou-se no presente estudo se a via miR156/SPL, da mesma forma que a partir de explantes aéreos, também é capaz de influenciar na regeneração de brotos in vitro a partir de explantes radiculares. Para tanto foram comparados taxa de regeneração, padrão de distribuição de auxina e citocinina, análises histológicas e histoquímicas das estruturas regeneradas em plantas com via miR156/SPL alterada, incluindo planta mutante hyl1, na qual a produção desse miRNA é severamente reduzida. Além disso, foi avaliado o padrão de expressão do miR156 e específicos genes SPL durante a regeneração de brotos in vitro a partir da raiz principal de Arabidopsis thaliana. No presente trabalho observou-se que a alteração da via gênica miR156/SPL é capaz de modular a capacidade de regeneração de brotos in vitro a partir de raiz principal de Arabidopsis thaliana e a distribuição de auxina e citocinina presente nas células e tecidos envolvidos no processo de regeneração. Plantas superexpressando o miR156 apresentaram redução no número de brotos regenerados, além de ter o plastochron reduzido quando comparado com plantas controle. Adicionalmente, plantas contento o gene SPL9 resistente à clivagem pelo miR156 (rSPL9) apresentaram severa redução na quantidade de brotos, além de terem o plastochron alongado. Interessantemente, plantas mutantes hyl1-2 e plantas rSPL10 não apresentaram regeneração de brotos ao longo da raiz principal, mas sim intensa formação de raízes laterais e protuberâncias, respectivamente, tendo essa última apresentado indícios de diferenciação celular precoce. Tomados em conjunto os dados sugerem que o miR156 apresenta importante papel no controle do processo de regeneração de brotos in vitro. Entretanto, esse efeito é mais complexo em regeneração in vitro a partir de raízes do que a partir de cotilédones ou hipocótilos.

MicroRNAs in pathology and as therapeutic target in gene therapy?

Dissertação para obtenção do grau de Mestre no Instituto Superior de Ciências da Saúde Egas Moniz

The prion protein gene: Identifying regulatory signals using marsupial sequence

The function of the prion protein gene (PRNP) and its normal product PrPC is elusive. We used comparative genomics as a strategy to understand the normal function of PRNP. As the reliability of comparisons increases with the number of species and increased evolutionary distance, we isolated and sequenced a 66.5 kb BAC containing the PRNP gene from a distantly related mammal, the model Australian marsupial Macropus eugenii (tammar wallaby). Marsupials are separated from eutherians such as human and mouse by roughly 180 million years of independent evolution. We found that tammar PRNP, like human PRNP, has two exons. Prion proteins encoded by the tammar wallaby and a distantly related marsupial, Monodelphis domestica (Brazilian opossum) PRNP contain proximal PrP repeats with a distinct, marsupial-specific composition and a variable number. Comparisons of tammar wallaby PRNP with PRNPs from human, mouse, bovine and ovine allowed us to identify non-coding gene regions conserved across the marsupial-eutherian evolutionary distance, which are candidates for regulatory regions. In the PRNP 3' UTR we found a conserved signal for nuclear-specific polyadenylation and the putative cytoplasmic polyadenylation element (CPE), indicating that post-transcriptional control of PRNP mRNA activity is important. Phylogenetic footprinting revealed conserved potential binding sites for the MZF-1 transcription factor in both upstream promoter and intron/intron 1, and for the MEF2, MyTI, Oct-1 and NFAT transcription factors in the intron(s). The presence of a conserved NFAT-binding site and CPE indicates involvement of PrPC in signal transduction and synaptic plasticity. (c) 2004 Elsevier B.V. All rights reserved.

Genomic organization of human arylamine N-acetyltransferase Type I reveals alternative promoters that generate different 5'-UTR splice variants with altered translational activities

In humans, a polymorphic gene encodes the drug-metabolizing enzyme NATI (arylamine N-acetyltransferase Type 1), which is widely expressed throughout the body. While the protein-coding region of NATI is contained within a single exon, examination of the human EST (expressed sequence tag) database at the NCBI revealed the presence of nine separate exons, eight of which were located in the 5'non-coding region of NATI. Differential splicing produced at least eight unique mRNA isoforms that could be grouped according to the location of the first exon, which suggested that NATI expression occurs from three alternative promoters. Using RT (reverse transcriptase)-PCR, we identified one major transcript in various epithelial cells derived from different tissues. In contrast, multiple transcripts were observed in blood-derived cell lines (CEM, THP-1 and Jurkat), with a novel variant, not identified in the EST database, found in CEM cells only. The major splice variant increased gene expression 9-11-fold in a luciferase reporter assay, while the other isoforrns were similar or slightly greater than the control. We examined the upstream region of the most active splice variant in a promoter-reporter assay, and isolated a 257 bp sequence that produced maximal promoter activity. This sequence lacked a TATA box, but contained a consensus Sp1 site and a CAAT box, as well as several other putative transcription-factor-binding sites. Cell-specific expression of the different NATI transcripts may contribute to the variation in NATI activity in vivo.