The brown coat colour of Coppernecked goats is associated with a non-synonymous variant at the TYRP1 locus on chromosome 8

The recent development of a goat SNP genotyping microarray enables genome-wide association studies in this important livestock species. We investigated the genetic basis of the black and brown coat colour in Valais Blacknecked and Coppernecked goats. A genome-wide association analysis using goat SNP50 BeadChip genotypes of 22 cases and 23 controls allowed us to map the locus for the brown coat colour to goat chromosome 8. The TYRP1 gene is located within the associated chromosomal region, and TYRP1 variants cause similar coat colour phenotypes in different species. We thus considered TYRP1 as a strong positional and functional candidate. We resequenced the caprine TYRP1 gene by Sanger and Illumina sequencing and identified two non-synonymous variants, p.Ile478Thr and p.Gly496Asp, that might have a functional impact on the TYRP1 protein. However, based on the obtained pedigree and genotype data, the brown coat colour in these goats is not due to a single recessive loss-of-function allele. Surprisingly, the genotype distribution and the pedigree data suggest that the (496) Asp allele might possibly act in a dominant manner. The (496) Asp allele was present in 77 of 81 investigated Coppernecked goats and did not occur in black goats. This strongly suggests heterogeneity underlying the brown coat colour in Coppernecked goats. Functional experiments or targeted matings will be required to verify the unexpected preliminary findings.

A comparative radiation hybrid map of sheep chromosome 10

Comparative radiation hybrid (RH) maps of individual ovine chromosomes are essential to identify genes governing traits of economic importance in sheep, a livestock species for which whole genome sequence data are not yet available. The USUoRH5000 radiation hybrid panel was used to generate a RH map of sheep chromosome 10 (OAR10) with 59 markers that span 1,422 cR over an estimated 92 Mb of the chromosome, thus providing markers every 2 Mb (equivalent to every 24 cR). The markers were derived from 46 BAC end sequences (BESs), a single EST, and 12 microsatellites. Comparative analysis showed that OAR10 shares remarkable conservation of gene order along the entire length of cattle chromosome 12 and that OAR10 contains four major homologous synteny blocks, each related to segments of the homologous human chromosome 13. Extending the comparison to the horse, dog, mouse, and chicken genome showed that these blocks share conserved synteny across species.

Molecular characterization and chromosomal assignment of the bovine glycinamide ribonucleotide formyltransferase (GART) gene on cattle chromosome 1q12.1-q12.2

The mammalian glycinamide ribonucleotide formyltransferase (GART) genes encode a trifunctional polypeptide involved in the de novo purine biosynthesis. We isolated a bacterial artificial chromosome (BAC) clone containing the bovine GART gene and determined the complete DNA sequence of the BAC clone. Cloning and characterization of the bovine GART gene revealed that the bovine gene consists of 23 exons spanning approximately 27 kb. RT-PCR amplification of bovine GART in different organs showed the expression of two GART transcripts in cattle similar to human and mouse. The GART transcripts encode two proteins of 1010 and 433 amino acids, respectively. Eleven single nucleotide polymorphisms (SNPs) were detected in a mutation scan of 24 unrelated animals of three different cattle breeds, including one SNP that affects the amino acid sequence of GART. The chromosomal localization of the gene was determined by fluorescence in situ hybridization. Comparative genome analysis between cattle, human and mouse indicates that the chromosomal location of the bovine GART gene is in agreement with a previously published mapping report.

A Multi-Breed Genome-Wide Association Analysis for Canine Hypothyroidism Identifies a Shared Major Risk Locus on CFA12.

Hypothyroidism is a complex clinical condition found in both humans and dogs, thought to be caused by a combination of genetic and environmental factors. In this study we present a multi-breed analysis of predisposing genetic risk factors for hypothyroidism in dogs using three high-risk breeds-the Gordon Setter, Hovawart and the Rhodesian Ridgeback. Using a genome-wide association approach and meta-analysis, we identified a major hypothyroidism risk locus shared by these breeds on chromosome 12 (p = 2.1x10-11). Further characterisation of the candidate region revealed a shared ~167 kb risk haplotype (4,915,018-5,081,823 bp), tagged by two SNPs in almost complete linkage disequilibrium. This breed-shared risk haplotype includes three genes (LHFPL5, SRPK1 and SLC26A8) and does not extend to the dog leukocyte antigen (DLA) class II gene cluster located in the vicinity. These three genes have not been identified as candidate genes for hypothyroid disease previously, but have functions that could potentially contribute to the development of the disease. Our results implicate the potential involvement of novel genes and pathways for the development of canine hypothyroidism, raising new possibilities for screening, breeding programmes and treatments in dogs. This study may also contribute to our understanding of the genetic etiology of human hypothyroid disease, which is one of the most common endocrine disorders in humans.

The gene for histone RNA hairpin binding protein is located on human chromosome 4 and encodes a novel type of RNA binding protein.

The hairpin structure at the 3' end of animal histone mRNAs controls histone RNA 3' processing, nucleocytoplasmic transport, translation and stability of histone mRNA. Functionally overlapping, if not identical, proteins binding to the histone RNA hairpin have been identified in nuclear and polysomal extracts. Our own results indicated that these hairpin binding proteins (HBPs) bind their target RNA as monomers and that the resulting ribonucleoprotein complexes are extremely stable. These features prompted us to select for HBP-encoding human cDNAs by RNA-mediated three-hybrid selection in Saccharomyces cerevesiae. Whole cell extract from one selected clone contained a Gal4 fusion protein that interacted with histone hairpin RNA in a sequence- and structure-specific manner similar to a fraction enriched for bovine HBP, indicating that the cDNA encoded HBP. DNA sequence analysis revealed that the coding sequence did not contain any known RNA binding motifs. The HBP gene is composed of eight exons covering 19.5 kb on the short arm of chromosome 4. Translation of the HBP open reading frame in vitro produced a 43 kDa protein with RNA binding specificity identical to murine or bovine HBP. In addition, recombinant HBP expressed in S. cerevisiae was functional in histone pre-mRNA processing, confirming that we have indeed identified the human HBP gene.

Meta-analysis identifies seven susceptibility loci involved in the atopic march.

Eczema often precedes the development of asthma in a disease course called the 'atopic march'. To unravel the genes underlying this characteristic pattern of allergic disease, we conduct a multi-stage genome-wide association study on infantile eczema followed by childhood asthma in 12 populations including 2,428 cases and 17,034 controls. Here we report two novel loci specific for the combined eczema plus asthma phenotype, which are associated with allergic disease for the first time; rs9357733 located in EFHC1 on chromosome 6p12.3 (OR 1.27; P=2.1 × 10(-8)) and rs993226 between TMTC2 and SLC6A15 on chromosome 12q21.3 (OR 1.58; P=5.3 × 10(-9)). Additional susceptibility loci identified at genome-wide significance are FLG (1q21.3), IL4/KIF3A (5q31.1), AP5B1/OVOL1 (11q13.1), C11orf30/LRRC32 (11q13.5) and IKZF3 (17q21). We show that predominantly eczema loci increase the risk for the atopic march. Our findings suggest that eczema may play an important role in the development of asthma after eczema.

CYTOGENETIC AND DNA CONTENT ANALYSIS IN MULTIPLE MYELOMA

In this dissertation, the cytogenetic characteristics of bone marrow cells from 41 multiple myeloma patients were investigated. These cytogenetic data were correlated with the total DNA content as measured by flow cytometry. Both the cytogenetic information and DNA content were then correlated with clinical data to determine if diagnosis and prognosis of multiple myeloma could be improved.^ One hundred percent of the patients demonstrated abnormal chromosome numbers per metaphase. The average chromosome number per metaphase ranged from 42 to 49.9, with a mean of 44.99. The percent hypodiploidy ranged from 0-100% and the percent hyperdiploidy from 0-53%. Detailed cytogenetic analyses were very difficult to perform because of the paucity of mitotic figures and the poor chromosome morphology. Thus, detailed chromosome banding analysis on these patients was impossible.^ Thirty seven percent of the patients had normal total DNA content, whereas 63% had abnormal amounts of DNA (one patient with less than normal amounts and 25 patients with greater than normal amounts of DNA).^ Several clinical parameters were used in the statistical analyses: tumor burden, patient status at biopsy, patient response status, past therapy, type of treatment and percent plasma cells. Only among these clinical parameters were any statistically significant correlations found: pretreatment tumor burden versus patient response, patient biopsy status versus patient response and past therapy versus patient response.^ No correlations were found between percent hypodiploid, diploid, hyperdiploid or DNA content, and the patient response status, nor were any found between those patients with: (a) normal plasma cells, low pretreatment tumor mass burden and more than 50% of the analyzed metaphases with 46 chromosomes; (b) normal amounts of DNA, low pretreatment tumor mass burden and more than 50% of the metaphases with 46 chromosomes; (c) normal amounts of DNA and normal quantities of plasma cells; (d) abnormal amounts of DNA, abnormal amounts of plasma cells, high pretreatment tumor mass burden and less than 50% of the metaphases with 46 chromosomes.^ Technical drawbacks of both cytogenetic and DNA content analysis in these multiple myeloma patients are discussed along with the lack of correlations between DNA content and chromosome number. Refined chromosome banding analysis awaits technical improvements before we can understand which chromosome material (if any) makes up the "extra" amounts of DNA in these patients. None of the correlations tested can be used as diagnostic or prognostic aids for multiple myeloma. ^

STUDIES OF THE MOLECULAR MECHANISMS OF CHROMOSOME ABERRATION FORMATION (DNA REPAIR, CROSSLINKS, MITOMYCIN C)

The objective of this research has been to study the molecular basis for chromosome aberration formation. Predicated on a variety of data, Mitomycin C (MMC)-induced DNA damage has been postulated to cause the formation of chromatid breaks (and gaps) by preventing the replication of regions of the genome prior to mitosis. The basic protocol for these experiments involved treating synchronized Hela cells in G(,1)-phase with a 1 (mu)g/ml dose of MMC for one hour. After removing the drug, cells were then allowed to progress to mitosis and were harvested for analysis by selective detachment. Utilizing the alkaline elution assay for DNA damage, evidence was obtained to support the conclusion that Hela cells can progress through S-phase into mitosis with intact DNA-DNA interstrand crosslinks. A higher level of crosslinking was observed in those cells remaining in interphase compared to those able to reach mitosis at the time of analysis. Dual radioisotope labeling experiments revealed that, at this dose, these crosslinks were associated to the same extent with both parental and newly replicated DNA. This finding was shown not to be the result of a two-step crosslink formation mechanism in which crosslink levels increase with time after drug treatment. It was also shown not to be an artefact of the double-labeling protocol. Using neutral CsCl density gradient ultracentrifugation of mitotic cells containing BrdU-labeled newly replicated DNA, control cells exhibited one major peak at a heavy/light density. However, MMC-treated cells had this same major peak at the heavy/light density, in addition to another minor peak at a density characteristic for light/light DNA. This was interpreted as indicating either: (1) that some parental DNA had not been replicated in the MMC treated sample or; (2) that a recombination repair mechanism was operational. To distinguish between these two possibilities, flow cytometric DNA fluorescence (i.e., DNA content) measurements of MMC-treated and control cells were made. These studies revealed that the mitotic cells that had been treated with MMC while in G(,1)-phase displayed a 10-20% lower DNA content than untreated control cells when measured under conditions that neutralize chromosome condensation effects (i.e., hypotonic treatment). These measurements were made under conditions in which the binding of the drug, MMC, was shown not to interfere with the stoichiometry of the ethidium bromide-mithramycin stain. At the chromosome level, differential staining techniques were used in an attempt to visualize unreplicated regions of the genome, but staining indicative of large unreplicated regions was not observed. These results are best explained by a recombinogenic mechanism. A model consistent with these results has been proposed.^

Physical and functional mapping of a novel tumor suppressor locus for prostate cancer within chromosome 10p12.31-q11

Prostate cancer remains the second leading cause of male cancer deaths in the United States, yet the molecular mechanisms underlying this disease remain largely unknown. Cytogenetic and molecular analyses of prostate tumors suggest a consistent association with the loss of chromosome 10. Previously, we have defined a novel tumor suppressor locus PAC-1 within chromosome 10pter-q11. Introduction of the short arm of chromosome 10 into a prostatic adenocarcinoma cell line PC-3H resulted in dramatic tumor suppression and restoration of a programmed cell death pathway. Using a combined approach of comparative genomic hybridization and microsatellite analysis of PC-3H, I have identified a region of hemizygosity within 10p12-p15. This region has been shown to be involved in frequent loss of heterozygosity in gliomas and melanoma. To functionally dissect the region within chromosome 10p containing PAC-1, we developed a strategy of serial microcell fusion, a technique that allows the transfer of defined fragments of chromosome 10p into PC-3H. Serial microcell fusion was used to transfer defined 10p fragments into a mouse A9 fibrosarcoma cell line. Once characterized by FISH and microsatellite analyses, the 10p fragments were subsequently transferred into PC-3H to generate a panel of microcell hybrid clones containing overlapping deletions of chromosome 10p. In vivo and microsatellite analyses of these PC hybrids identified a small chromosome 10p fragment (an estimated 31 Mb in size inclusive of the centromere) that when transferred into the PC-3H background, resulted in significant tumor suppression and limited a region of functional tumor suppressor activity to chromosome 10p12.31-q11. This region coincides with a region of LOH demonstrated in prostate cancer. These studies demonstrate the utility of this approach as a powerful tool to limit regions of functional tumor suppressor activity. Furthermore, these data used in conjunction with data generated by the Human Genome Project lent a focused approach to identify candidate tumor suppressor genes involved in prostate cancer. ^

Identification and analysis of novel mitotic inhibitors of the Caenorhabditis elegans Aurora B kinase

Proper execution of mitosis requires the accurate segregation of replicated DNA into each daughter cell. The highly conserved mitotic kinase AIR-2/Aurora B is a dynamic protein that interacts with subsets of cofactors and substrates to coordinate chromosome segregation and cytokinesis in Caenorhabdiris elegans. To identify components of the AIR-2 regulatory pathway, a genome-wide RNAi-based screen for suppressors of air-2 temperature-sensitive mutant lethality was conducted. Here, I present evidence that two classes of suppressors identified in this screen are bona fide regulators of the AIR-2 kinase. The strongest suppressor cdc-48.3, encodes an Afg2/Spaf-related Cdc48-like AAA+ ATPase that regulates AIR-2 kinase activity and stability during C. elegans embryogenesis. Loss of CDC-48.3 suppresses the lethality of air-2 mutant embryos, marked by the restoration of the dynamic behavior of AIR-2 and rescue of chromosome segregation and cytokinesis defects. Loss of CDC-48.3 leads to mitotic delays and abnormal accumulation of AIR-2 during late telophase/mitotic exit. In addition, AIR-2 kinase activity is significantly upregulated from metaphase through mitotic exit in CDC-48.3 depleted embryos. Inhibition of the AIR-2 kinase is dependent on (1) a direct physical interaction between CDC-48.3 and AIR-2, and (2) CDC-48.3 ATPase activity. Importantly, the increase in AIR-2 kinase activity does not correlate with the stabilization of AIR-2 in late mitosis. Hence, CDC-48.3 is a bi-functional inhibitor of AIR-2 that is likely to act via distinct mechanisms. The second class of suppressors consists of psy-2/smk-1 and pph-4.1, which encode two components of the conserved PP4 phosphatase complex that is essential for spindle assembly, chromosome segregation, and overall mitotic progression. AIR-2 and its substrates are likely to be targets of this complex since mitotic AIR-2 kinase activity is significantly increased during mitosis when either PSY-2/SMK-1 or PPH-4.l is depleted. Altogether, this study demonstrates that during the C. elegans embryonic cell cycle, regulators including the CDC-48.3 ATPase and PP4 phosphatase complex interact with and control the kinase activity, targeting behavior and protein stability of the Aurora B kinase to ensure accurate and timely progression of mitosis. ^

IDENTIFICATION AND ANALYSIS OF A NOVEL ROLE FOR THE TOUSLED-LIKE KINASE IN REGULATING MITOTIC SPINDLE DYNAMICS

Deregulation of kinase activity is one example of how cells become cancerous by evading evolutionary constraints. The Tousled kinase (Tsl) was initially identified in Arabidopsis thaliana as a developmentally important kinase. There are two mammalian orthologues of Tsl and one orthologue in C. elegans, TLK-1, which is essential for embryonic viability and germ cell development. Depletion of TLK-1 leads to embryonic arrest large, distended nuclei, and ultimately embryonic lethality. Prior to terminal arrest, TLK-1-depleted embryos undergo aberrant mitoses characterized by poor metaphase chromosome alignment, delayed mitotic progression, lagging chromosomes, and supernumerary centrosomes. I discovered an unanticipated requirement for TLK-1 in mitotic spindle assembly and positioning. Normally, in the newly-fertilized zygote (P0) the maternal pronucleus migrates toward the paternal pronucleus at the posterior end of the embryo. After pronuclear meeting, the pronuclear-centrosome complex rotates 90° during centration to align on the anteroposterior axis followed by nuclear envelope breakdown (NEBD). However, in TLK-1-depleted P0 embryos, the centrosome-pronuclear complex rotation is significantly delayed with respect to NEBD and chromosome congression, Additionally, centrosome positions over time in tlk-1(RNAi) early embryos revealed a defect in posterior centrosome positioning during spindle-pronuclear centration, and 4D analysis of centrosome positions and movement in newly fertilized embryos showed aberrant centrosome dynamics in TLK-1-depleted embryos. Several mechanisms contribute to spindle rotation, one of which is the anchoring of astral microtubules to the cell cortex. Attachment of these microtubules to the cortices is thought to confer the necessary stability and forces in order to rotate the centrosome-pronuclear complex in a timely fashion. Analysis of a microtubule end-binding protein revealed that TLK-1-depleted embryos exhibit a more stochastic distribution of microtubule growth toward the cell cortices, and the types of microtubule attachments appear to differ from wild-type embryos. Additionally, fewer astral microtubules are in the vicinity of the cell cortex, thus suggesting that the delayed spindle rotation could be in part due to a lack of appropriate microtubule attachments to the cell cortex. Together with recently published biochemical data revealing the Tousled-like kinases associate with components of the dynein microtubule motor complex in humans, these data suggest that Tousled-like kinases play an important role in mitotic spindle assembly and positioning.

Seed image analysis and taxonomy of Diplotaxis DC. (Brassicaceae, Brassiceae)

The genus Diplotaxis, comprising 32 or 34 species, plus several additional infraspecific taxa, displays a considerable degree of heterogeneity in the morphology, molecular markers, chromosome numbers and geographical amplitude of the species. The taxonomic relationships within the genus Diplotaxis were investigated by phenetic characterisation of germplasm belonging to 27 taxa of the genus, because there is an increasing interest in Diplotaxis, since some of its species (D. tenuifolia, D. muralis) are gathered or cultivated for human consumption, whereas others are frequent arable weeds (D. erucoides) in many European vineyards. Using a computer-aided vision system, 33 morpho-colorimetric features of seeds were electronically measured. The data were used to implement a statistical classifier, which is able to discriminate the taxa within the genus Diplotaxis, in order to compare the resulting species grouping with the current infrageneric systematics of this genus. Despite the high heterogeneity of the samples, due to the great intra-population variability, the stepwise Linear Discriminant Analysis method, applied to distinguish the groups, was able to reach over 80% correct identification. The results obtained allowed us to confirm the current taxonomic position of most taxa and suggested the taxonomic position of others for reconsideration.

A cluster of oppositely imprinted transcripts at the Gnas locus in the distal imprinting region of mouse chromosome 2

Imprinted genes tend to occur in clusters. We have identified a cluster in distal mouse chromosome (Chr) 2, known from early genetic studies to contain both maternally and paternally imprinted, but unspecified, genes. Subsequently, one was identified as Gnas, which encodes a G protein α subunit, and there is clinical and biochemical evidence that the human homologue GNAS1, mutated in patients with Albright hereditary osteodystrophy, is also imprinted. We have used representational difference analysis, based on parent-of-origin methylation differences, to isolate candidate imprinted genes in distal Chr 2 and found two oppositely imprinted genes, Gnasxl and Nesp. Gnasxl determines a variant G protein α subunit associated with the trans-Golgi network and Nesp encodes a secreted protein of neuroendocrine tissues. Gnasxl is maternally methylated in genomic DNA and encodes a paternal-specific transcript, whereas Nesp is paternally methylated with maternal-specific expression. Their reciprocal imprinting may offer insight into the distal Chr 2 imprinting phenotypes. Remarkably, Gnasxl, Nesp, and Gnas are all part of the same transcription unit; transcripts for Gnasxl and Nesp are alternatively spliced onto exon 2 of Gnas. This demonstrates an imprinting mechanism in which two oppositely imprinted genes share the same downstream exons.

Genetic analysis of the mouse X inactivation center defines an 80-kb multifunction domain

Dosage compensation in mammals occurs by X inactivation, a silencing mechanism regulated in cis by the X inactivation center (Xic). In response to developmental cues, the Xic orchestrates events of X inactivation, including chromosome counting and choice, initiation, spread, and establishment of silencing. It remains unclear what elements make up the Xic. We previously showed that the Xic is contained within a 450-kb sequence that includes Xist, an RNA-encoding gene required for X inactivation. To characterize the Xic further, we performed deletional analysis across the 450-kb region by yeast-artificial-chromosome fragmentation and phage P1 cloning. We tested Xic deletions for cis inactivation potential by using a transgene (Tg)-based approach and found that an 80-kb subregion also enacted somatic X inactivation on autosomes. Xist RNA coated the autosome but skipped the Xic Tg, raising the possibility that X chromosome domains escape inactivation by excluding Xist RNA binding. The autosomes became late-replicating and hypoacetylated on histone H4. A deletion of the Xist 5′ sequence resulted in the loss of somatic X inactivation without abolishing Xist expression in undifferentiated cells. Thus, Xist expression in undifferentiated cells can be separated genetically from somatic silencing. Analysis of multiple Xic constructs and insertion sites indicated that long-range Xic effects can be generalized to different autosomes, thereby supporting the feasibility of a Tg-based approach for studying X inactivation.

The FEZ1 gene at chromosome 8p22 encodes a leucine-zipper protein, and its expression is altered in multiple human tumors

Alterations of human chromosome 8p occur frequently in many tumors. We identified a 1.5-Mb common region of allelic loss on 8p22 by allelotype analysis. cDNA selection allowed isolation of several genes, including FEZ1. The predicted Fez1 protein contained a leucine-zipper region with similarity to the DNA-binding domain of the cAMP-responsive activating-transcription factor 5. RNA blot analysis revealed that FEZ1 gene expression was undetectable in more than 60% of epithelial tumors. Mutations were found in primary esophageal cancers and in a prostate cancer cell line. Transcript analysis from several FEZ1-expressing tumors revealed truncated mRNAs, including a frameshift. Alteration and inactivation of the FEZ1 gene may play a role in various human tumors.