Association between polymorphisms in the progesterone receptor gene and endometriosis

The progesterone receptor (PR) is a candidate gene for the development of endometriosis, a complex disease with strong hormonal features, common in women of reproductive age. We typed the 306 base pair Alu insertion (AluIns) polymorphism in intron G of PR in 101 individuals, estimated linkage disequilibrium (LD) between five single-nucleotide polymorphisms (SNPs) across the PR locus in 980 Australian triads (endometriosis case and two parents) and used transmission disequilibrium testing (TDT) for association with endometriosis. The five SNPs showed strong pairwise LD, and the AluIns was highly correlated with proximal SNPs rs1042839 ({Delta}2 = 0.877, D9 = 1.00, P < 0.0001) and rs500760 ({Delta}2 = 0.438, D9 = 0.942, P < 0.0001). TDT showed weak evidence of allelic association between endometriosis and rs500760 (P = 0.027) but not in the expected direction. We identified a common susceptibility haplotype GGGCA across the five SNPs (P = 0.0167) in the whole sample, but likelihood ratio testing of haplotype transmission and non-transmission of the AluIns and flanking SNPs showed no significant pattern. Further, analysis of our results pooled with those from two previous studies suggested that neither the T2 allele of the AluIns nor the T1/T2 genotype was associated with endometriosis.

Ultraconserved elements in insect genomes: A highly conserved intronic sequence implicated in the control of homothorax mRNA splicing

Recently, we identified a large number of ultraconserved (uc) sequences in noncoding regions of human, mouse, and rat genomes that appear to be essential for vertebrate and amniote ontogeny. Here, we used similar methods to identify ultraconserved genomic regions between the insect species Drosophila melanogaster and Drosophila pseudoobscura, as well as the more distantly related Anopheles gambiae. As with vertebrates, ultraconserved sequences in insects appear to Occur primarily in intergenic and intronic sequences, and at intron-exon junctions. The sequences are significantly associated with genes encoding developmental regulators and transcription factors, but are less frequent and are smaller in size than in vertebrates. The longest identical, nongapped orthologous match between the three genomes was found within the homothorax (hth) gene. This sequence spans an internal exon-intron junction, with the majority located within the intron, and is predicted to form a highly stable stem-loop RNA structure. Real-time quantitative PCR analysis of different hth splice isoforms and Northern blotting showed that the conserved element is associated with a high incidence of intron retention in hth pre-mRNA, suggesting that the conserved intronic element is critically important in the post-transcriptional regulation of hth expression in Diptera.

Characterization of two polymorphisms in the leukotriene C4 synthase gene in an Australian population of subjects with mild, moderate, and severe asthma

Background: The cysteinyl-leukotrienes (cys-LTs) are proinflammatory mediators that are important in the pathophysiology of asthma. LTC4 synthase is a key enzyme in the cys-LT biosynthetic pathway, and studies in small populations have suggested that a promoter polymorphism (A(-444)C) in the gene might be associated with asthma severity and aspirin intolerance. Objective: We sought to screen the LTC4 synthase gene for polymorphisms and to determine whether there is an association between these polymorphisms and asthma severity or aspirin sensitivity in a large, well-phenotyped population and to determine whether this polymorphism is functionally relevant. Methods: The coding regions of the LTC4 synthase gene were screened for polymorphisms and the A(-444)C polymorphism was analyzed in a large Australian white adult population of mild (n = 282), moderate (n = 236), and severe asthmatic subjects (n = 86) and nonasthmatic subjects (n = 458), as well as in aspirin-intolerant asthmatic subjects (n = 67). The functional activity of the promoter polymorphism was investigated by transient transfection of HL-60 cells with a promoter construct. Results: A new polymorphism was identified in intron 1 of the gene (IVS1-10c>a) but was not associated with asthma. Association studies showed that the A(-444)C polymorphism was weakly associated with asthma per se, but there was no association between the C-444 allele and chronic asthma severity or aspirin intolerance. A meta-analysis of all the genetic studies conducted to date found significant between-study heterogeneity in C-444 allele frequencies within different clinical subgroups. In vitro functional studies showed no significant differences in transcription efficiency between constructs containing the A(-444) allele or the C-444 allele. Conclusions: Our data confirm that, independent of transcriptional activity, the C-444 allele in the LTC4 synthase gene is weakly associated with the asthma phenotype, but it is not related to disease severity or aspirin intolerance.

Cytokine gene polymorphisms in idiopathic pulmonary fibrosis

Pro- and anti-fibrotic cytokine gene polymorphisms may affect expression of idiopathic pulmonary fibrosis (IPF). The aims of the present case-control study were to examine polymorphisms in the IL-6, transforming growth factor (TGF)-beta1, tumour necrosis factor (TNF)-alpha and interleukin-1 (IL-1)Ra genes in patients with IPF (n=22)-compared to healthy controls (n=140). Genotyping was performed on DNA extracted from peripheral blood lymphocytes, using polymerase chain reaction-restriction fragment length polymorphism with gene polymorphisms determined according to-published techniques. The following sites were examined: (i) IL-1Ra*1-5 (86 bp variable tandem repeat intron 2), (ii) IL-6 (-174G>C), (iii) TNF-alpha (-308G>A) and (iv) TGF-beta1 (Arg25Pro). The TNF-alpha (-308 A) allele was over-represented in the IPF (p(corr)=0.004) group compared to controls. Risk of IPF was significant for heterozygotes for: (i) the TNF-alpha (-308 A) allele (A/G) (odds ratio (OR) 2.9; 95% confidence interval (CI) 1.2-7.2; P=0.02), (ii) homozygotes (A/A) (OR 13.9; 95%CI 1.2-160; P=0.04) and (iii) carriage of the allele (A/A+A/G) (OR 4; 95%CI 1.6-10.2; P=0.003). The distribution of alleles and genotypes for IL-6, TGF-beta1 and IL-1Ra between the two groups was not significantly different. This is the third study to independently confirm that there is a significant association of the TNF-alpha (-308 A) allele with IPF. Further research is needed to assess the utility of cytokine gene polymorphisms as markers of disease-susceptibility.

Comparative phylogeographic summary statistics for testing simultaneous vicariance

Testing for simultaneous vicariance across comparative phylogeographic data sets is a notoriously difficult problem hindered by mutational variance, the coalescent variance, and variability across pairs of sister taxa in parameters that affect genetic divergence. We simulate vicariance to characterize the behaviour of several commonly used summary statistics across a range of divergence times, and to characterize this behaviour in comparative phylogeographic datasets having multiple taxon-pairs. We found Tajima's D to be relatively uncorrelated with other summary statistics across divergence times, and using simple hypothesis testing of simultaneous vicariance given variable population sizes, we counter-intuitively found that the variance across taxon pairs in Nei and Li's net nucleotide divergence (pi(net)), a common measure of population divergence, is often inferior to using the variance in Tajima's D across taxon pairs as a test statistic to distinguish ancient simultaneous vicariance from variable vicariance histories. The opposite and more intuitive pattern is found for testing more recent simultaneous vicariance, and overall we found that depending on the timing of vicariance, one of these two test statistics can achieve high statistical power for rejecting simultaneous vicariance, given a reasonable number of intron loci (> 5 loci, 400 bp) and a range of conditions. These results suggest that components of these two composite summary statistics should be used in future simulation-based methods which can simultaneously use a pool of summary statistics to test comparative the phylogeographic hypotheses we consider here.

Further evidence of linkage at 7p22 with familial hyperaldosteronism type II

Primary aldosteronism (PAL) is caused by the autonomous over-production of aldosterone. Once thought rare, it is now reported to be responsible for 5–10% of hypertension. Familial hyperaldosteronism type II (FH-II), unlike familial hyperaldosteronism type I, is not glucocorticoid-remediable and not associated with the hybrid CYP11B1/CYP11B2 gene mutation. At least five times more common than FH-I, FH-II is clinically, biochemically and morphologically indistinguishable from apparently sporadic PAL, suggesting that its incidence maybe even higher. Studies performed in collaboration with C Stratakis (NIH, Bethesda) on our largest Australian FH-II family (eight affected members) demonstrated linkage at chromosome 7p22. Similar linkage at this region was also found in a South American FH-II family (DNA provided by MI New, Presbyterian Hospital, New York). Mutations in the exons and intron/exon boundaries of the PRKARIB gene (which resides at 7p22 and is closely related to PRKARIA gene mutated in Carney complex) have been excluded in our largest Australian FH-II family. Using more finely spaced markers, we have confirmed linkage at 7p22 in these 2 families, and identified a second Australian family with evidence of linkage at this locus. The combined multipoint LOD score for these 3 families is 4.87 (θ=0) with markers D7S462 and D7S2424, which exceeds the critical threshold for genome-wide significance suggested by Lander and Kruglyak (1995), providing strong support for this locus harbouring mutations responsible for FH-II. A newly identified recombination event in our largest Australian family has narrowed the region of linkage by 1.8 Mb, permitting exclusion of approximately half the genes residing in the original reported 5Mb linked locus. In addition, we have strongly excluded linkage to these key markers in two Australian families (maximum multipoint LOD scores −3.51 and −2.77), supporting the notion that FH-II may be genetically heterogeneous. In order to identify candidate genes at 7p22, more closely spaced markers will be used to refine the locus, as well as single nucleotide polymorphism analysis.

Detection of self-complementary inverted repeats by single forward primer driven PCR

Inverted repeat gene structures designed for silencing functional genes have been widely used both in academic and applied research. The correct orientations of such structures are usually validated with restriction analysis and/or sequencing. We speculated that the inverted repeat nature of such constructs can be shown by a simple PCR reaction with a single forward primer. To test this hypothesis five different constructs were established from grapevine sequences in a hairpin-intron style silencing system. We were able to amplify the appropriate products in each case. Thus a forward-primed PCR alone may be sufficient to prove the inverted repeat nature of the desired constructs.

The isolation and characterization of tescalcin, a novel gene differentially expressed in the mouse testis during the early stages of gonadal differentiation

Gonadal development is an ideal model to study organogenesis because a variety of developmental processes can be studied during the differentiation of the bipotential primordium into testis or ovary. To better understand this process, Representational Difference Analysis of cDNA was used to identify genes that are differentially expressed in mouse gonads at 13.5 days post-coitus. The analysis led to the identification of three testis specific genes and a sequence that was only expressed in the ovary. The male genes identified: renin, Col9a3, and a novel gene termed tescalcin had patterns of expression that suggested a role in testis determination. ^ Studies of the tescalcin gene revealed that it is organized into eight exons and seven introns. The gene was located at 64 cM in mouse chromosome 5, where it spans approximately 35 Kb. Three mRNA variants resulting from alternative splicing of intron 5 were identified in mouse tissues. Gel mobility shift assays demonstrated that Sp1 and Sp3 from Y-1, msc-1, and MIN-6 cells nuclear extracts bind the GC-boxes within the tescalcin proximal promoter. Bisulfite sequencing analysis of tescalcin CpG island revealed that it is differentially methylated in male and female mouse embryonic gonads, and that hypermethylation of this region represses expression of tescalcin in the β-TC3 cell line. ^ The major tescalcin mRNA encodes a protein with 214 amino acids that contains a consensus EF-hand Ca2+-binding domain and an N-myristoylation motif. The amino acid sequence of tescalcin is highly conserved among various species, and it showed the highest homology with calcineurin B homologous proteins 1 and 2, and calcineurin B. Western blot analysis using antibodies generated against the tescalcin protein confirmed its presence in specific mouse tissues and cell lines. Immunohistochemical analysis of mouse embryos confirmed the pattern of expression of tescalcin mRNA in fetal testis. Using pull-down assays, glyceraidehydes-3-phosphate dehydrogenase was identified as an interacting and potential functional partner of tescalcin. ^ The identification and characterization of tescalcin as a novel embryonic testicular marker will contribute to the elucidation of the genetic pathways involved in testis development and likely to the understanding of pathological conditions such as sex reversal and infertility. ^

Dioon Lindl. (Zamiaceae): Perspectives from phylogeny and a population genetic study of Dioon edule

Dioon Lindl. (Zamiaceae) is a small genus restricted to Mexico (12 species) and Honduras (one species). Previous systematic studies have been unable to fully resolve species relationships within the genus. Phylogenetic analyses were conducted with data from several sources, including Restriction Fragment Length Polymorphisms from the chloroplast genome, morphology, two introns of the low copy nuclear gene S-adenosyl-L-homocysteine hydrolase (SAHH) and the 5.8S/ITS2 regions of the nuclear ribosomal DNA. The goals of the study were to construct a total evidence species level phylogeny and to explore current biogeographical hypotheses. None of the analyses performed produced a fully resolved topology. Dioon is comprised of two main lineages (the Edule and Spinulosum Clades), which represents an ancient divergence within the genus. The two introns of the nuclear gene SAHH offer additional evidence for the split into two lineages. Intron 2 contains a 18 bp deletion in the Spinulosum Clade, providing a synapomorphy for that group. The 5.8S/ITS2 regions were highly polymorphic and subsequently omitted from the combined analyses. In order to visualize congruence between morphology and molecular data, morphological characters were mapped onto the combined molecular tree. Current biogeographical hypotheses of a general northward pattern of migration and speciation are supported here. However, sister relationships within the Edule Clade are not fully resolved. Seven DNA microsatellite markers were developed to investigate patterns of genetic variation of seven populations of D. edule, a species restricted to Eastern Mexico. We found that most of the genetic variation lies within populations (Ho = 0.2166–0.3657) and that levels of population differentiation are low (Fst = 0.088); this finding is congruent with the breeding system of this species, dioicy. Four of the populations deviate from Hardy Weinberg Equilibrium and have a high number of identical genotypes, we suggest that this unexpected pattern is due to the life-history strategy of the species coupled with the few number of polymorphic loci detected in these populations. Our results are not congruent with earlier evidence from morphology and allozyme markers that suggest that the two northernmost populations represent a distinct entity that is recognized by some taxonomists as D. angustifolium.

The role of splicing factors and small nuclear RNAs in spliceosomal formation

Protein coding genes are comprised of protein-coding exons and non-protein-coding introns. The process of splicing involves removal of the introns and joining of the exons to form a mature messenger RNA, which subsequently undergoes translation into polypeptide. The spliceosome is a large, RNA/protein assembly of five small nuclear RNAs as well as over 300 proteins, which catalyzes intron removal and exon ligation. The selection of specific exons for inclusion in the mature messenger RNA is spatiotemporally regulated and results in production of an enormous diversity of polypeptides from a single gene locus. This phenomenon, known as alternative splicing, is regulated, in part, by protein splicing factors, which target the spliceosome to exon/intron boundaries. The first part of my dissertation (Chapters II and III) focuses on the discovery and characterization of the 45 kilodalton FK506 binding protein (FKBP45), which I discovered in the silk moth, Bombyx mori, as a U1 small nuclear RNA binding protein. This protein family binds the immunosuppressants FK506 and rapamycin and contains peptidyl-prolyl cis-trans isomerase activity, which converts polypeptides from cis to trans about a proline residue. This is the first time that an FKBP has been identified in the spliceosome. The second section of my dissertation (Chapters IV, V, VI and VII) is an investigation of the potential role of small nuclear RNA sequence variants in the control of splicing. I identified 46 copies of small nuclear RNAs in the 6X whole genome shotgun of the Bombyx mori p50T strain. These variants may play a role in differential binding of specific proteins that mediate alternative splicing. Along these lines, further investigation of U2 snRNA sequence variants in Bombyx mori demonstrated that some U2 snRNAs preferentially assemble into high molecular weight spliceosomal complexes over others. Expression of snRNA variants may represent another mechanism by which the cell is able to fine tune the splicing process.

High natural gene expression variation in the reef-building coral Acropora millepora: potential for acclimative and adaptive plasticity

Background: Ecosystems worldwide are suffering the consequences of anthropogenic impact. The diverse ecosystem of coral reefs, for example, are globally threatened by increases in sea surface temperatures due to global warming. Studies to date have focused on determining genetic diversity, the sequence variability of genes in a species, as a proxy to estimate and predict the potential adaptive response of coral populations to environmental changes linked to climate changes. However, the examination of natural gene expression variation has received less attention. This variation has been implicated as an important factor in evolutionary processes, upon which natural selection can act. Results: We acclimatized coral nubbins from six colonies of the reef-building coral Acropora millepora to a common garden in Heron Island (Great Barrier Reef, GBR) for a period of four weeks to remove any site-specific environmental effects on the physiology of the coral nubbins. By using a cDNA microarray platform, we detected a high level of gene expression variation, with 17% (488) of the unigenes differentially expressed across coral nubbins of the six colonies (jsFDR-corrected, p < 0.01). Among the main categories of biological processes found differentially expressed were transport, translation, response to stimulus, oxidation-reduction processes, and apoptosis. We found that the transcriptional profiles did not correspond to the genotype of the colony characterized using either an intron of the carbonic anhydrase gene or microsatellite loci markers. Conclusion: Our results provide evidence of the high inter-colony variation in A. millepora at the transcriptomic level grown under a common garden and without a correspondence with genotypic identity. This finding brings to our attention the importance of taking into account natural variation between reef corals when assessing experimental gene expression differences. The high transcriptional variation detected in this study is interpreted and discussed within the context of adaptive potential and phenotypic plasticity of reef corals. Whether this variation will allow coral reefs to survive to current challenges remains unknown.

Isolation and characterization of the U1 and U2 snRNAs in the silk gland of the 5th instar larval stage silk moth Bombyx mori

In isolation and characterization studies, expression level U1 and U2 snRNA isoforms were obtained from the 5th instar larval stage silk gland (SG). The DNA content of the SG cells is approximately 200,000-fold higher compared to the usual (2N) somatic cells of B. mori due to endoreduplication. In this study, the existence of U1 and U2 snRNA isoforms in the SG of the organism is investigated. Bombyx mori U1 and U2-specific RT-PCR libraries from the silk gland were generated. Five U1 and eight U2 isoforms were isolated and characterized. Nucleotide differences, structural alterations, as well as protein and RNA interaction sites were analyzed in these variants. For the U1 snRNA variants, they were compared to the previously reported BmN isoforms. In all these U-snRNA variants, polymorphic sites do not predominate at the core of known functional sequences, which were interspecifically conserved. Variant sites and inter-species differences are located in moderately conserved regions. Free energy (ΔG) values for the entire U1 and U2 snRNA secondary structures and for the individual stem/loops domains of the isoforms were generated and compared to determine their structural stability. This will be the first time that U1 and U2 variants are shown specific for a development stage (larval) other than embryonic or adult. ^ Using phylogenetic analysis, evolutionary trees were generated for the U1 and U2 snRNAs using animal, plant, protista and fungal species. The resulting trees were boostrapped for robustness and rooted with the self-splicing RNA group II intron sequence from the cyanobacterium Calothrix. Using phylogenetic analyses, possible structural and functional evolutionary interdependence between the U1 and U2 snRNAs was investigated. ^

Transgenic Expression of Endothelin Receptor-B Rescues the Aganglionosis Phenotype of Piebald Lethal Mice

Neural crest cells (NCC) are a unique population of cells in vertebrates that arise between the presumptive epidermis and the dorsal most region of the neural tube. During neurulation, NCC migrate to many regions of the body to give rise to a wide variety of cell types. NCC that originate from the neural tube at the levels of somite 1-7 colonize the gut and give rise to the enteric ganglia. The endothelin signaling pathway has been shown to be crucial for proper development of some neural crest derivatives. Mice and humans with mutations in the Endothelin receptor b (Ednrb) gene exhibit similar phenotypes characterized by hypopigmentation, hearing loss, and megacolon. Thesephenotypes are due to lack of melanocytes in the skin, inner ear and enteric ganglia in the distal portion of the colon, respectively. It is well established that Ednrb is required early during the embryonic development for normal innervation of the gut. However, it is not clear if Ednrb acts on enteric neuron precursor cells or in pre-committed NC precursors. Additionally, it is controversial whether the action of Ednrb is cell autonomous or non- autonomous. We generated transgenic mice that express Ednrb under the control of the Nestin second intron enhancer (Nes) which drives expression to pre-migrating NCC. These mice were crosses to the spontaneous mouse mutant piebald lethal, which carriers a null mutation in Ednrb and exhibits enteric aganglionosis. The Nes-Ednrb was capable of rescuing the aganglianosis phenotype of piebald lethal mutants demonstrating that expression of Ednrb in pre-committed precursors is sufficient for normal enteric ganglia development. This study provides insight in early embryonic development of NCC and could eventually have potential use in cellular therapies for Hirschsprung's disease.

The Role of Splicing Factors and Small Nuclear RNAS in Spliceosomal Formation

Protein coding genes are comprised of protein-coding exons and non-protein-coding introns. The process of splicing involves removal of the introns and joining of the exons to form a mature messenger RNA, which subsequently undergoes translation into polypeptide. The spliceosome is a large, RNA/protein assembly of five small nuclear RNAs as well as over 300 proteins, which catalyzes intron removal and exon ligation. The selection of specific exons for inclusion in the mature messenger RNA is spatio-temporally regulated and results in production of an enormous diversity of polypeptides from a single gene locus. This phenomenon, known as alternative splicing, is regulated, in part, by protein splicing factors, which target the spliceosome to exon/intron boundaries. The first part of my dissertation (Chapters II and III) focuses on the discovery and characterization of the 45 kilodalton FK506 binding protein (FKBP45), which I discovered in the silk moth, Bombyx mori, as a U1 small nuclear RNA binding protein. This protein family binds the immunosuppressants FK506 and rapamycin and contains peptidyl-prolyl cis-trans isomerase activity, which converts polypeptides from cis to trans about a proline residue. This is the first time that an FKBP has been identified in the spliceosome. The second section of my dissertation (Chapters IV, V, VI and VII) is an investigation of the potential role of small nuclear RNA sequence variants in the control of splicing. I identified 46 copies of small nuclear RNAs in the 6X whole genome shotgun of the Bombyx mori p50T strain. These variants may play a role in differential binding of specific proteins that mediate alternative splicing. Along these lines, further investigation of U2 snRNA sequence variants in Bombyx mori demonstrated that some U2 snRNAs preferentially assemble into high molecular weight spliceosomal complexes over others. Expression of snRNA variants may represent another mechanism by which the cell is able to fine tune the splicing process.