Cattle-Derived Microsatellite Markers to Study the Water Buffalo Genome

Microsatellites are well-known DNA markers used in a variety of studies such as genome mapping, genetic diversity analysis, genetic conservation and phylogenetic studies. Although microsatellites are important markers, their development and characterization demands extensive time and high cost. Thus, before new markers are developed for a particular species, it is worthwhile to test the available markers from related species. In the present study, we evaluate cattle-derived microsatellite markers for genetic studies of water buffalo. Eighty-five percents of a total of 120 microsatellite markers were optimized using buffalo DNA (Bubalus bubalis). The results showed in this paper were also deposited in the National Center for Biological Information database (NCBI) (ProbeDB and UniSTS) for use in population genetic studies of buffalo by the scientific community. The use of heterologous primers significantly reduces the cost of developing specific markers for buffalo, providing a useful short cut for the genetic population analysis and gene mapping studies.

Construction and preliminary characterization of a river buffalo bacterial artificial chromosome library

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Precision of distances and ordering of microsatellite markers in consensus linkage maps of chromosomes 1, 3 and 4 from two reciprocal chicken populations using bootstrap sampling

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genome wide scan for quantitative trait loci affecting tick resistance in cattle (Bos taurus × Bos indicus)

Abstract Background In tropical countries, losses caused by bovine tick Rhipicephalus (Boophilus) microplus infestation have a tremendous economic impact on cattle production systems. Genetic variation between Bos taurus and Bos indicus to tick resistance and molecular biology tools might allow for the identification of molecular markers linked to resistance traits that could be used as an auxiliary tool in selection programs. The objective of this work was to identify QTL associated with tick resistance/susceptibility in a bovine F2 population derived from the Gyr (Bos indicus) × Holstein (Bos taurus) cross. Results Through a whole genome scan with microsatellite markers, we were able to map six genomic regions associated with bovine tick resistance. For most QTL, we have found that depending on the tick evaluation season (dry and rainy) different sets of genes could be involved in the resistance mechanism. We identified dry season specific QTL on BTA 2 and 10, rainy season specific QTL on BTA 5, 11 and 27. We also found a highly significant genome wide QTL for both dry and rainy seasons in the central region of BTA 23. Conclusions The experimental F2 population derived from Gyr × Holstein cross successfully allowed the identification of six highly significant QTL associated with tick resistance in cattle. QTL located on BTA 23 might be related with the bovine histocompatibility complex. Further investigation of these QTL will help to isolate candidate genes involved with tick resistance in cattle.

Mitochondrial Genomics: structure, inheritance and phylogenetic utility of the Mitochondrial genome in Bivalvia and Insecta

This PhD Thesis is the result of my research activity in the last three years. My main research interest was centered on the evolution of mitochondrial genome (mtDNA), and on its usefulness as a phylogeographic and phylogenetic marker at different taxonomic levels in different taxa of Metazoa. From a methodological standpoint, my main effort was dedicated to the sequencing of complete mitochondrial genomes, and the approach to whole-genome sequencing was based on the application of Long-PCR and shotgun sequences. Moreover, this research project is a part of a bigger sequencing project of mtDNAs in many different Metazoans’ taxa, and I mostly dedicated myself to sequence and analyze mtDNAs in selected taxa of bivalves and hexapods (Insecta). Sequences of bivalve mtDNAs are particularly limited, and my study contributed to extend the sampling. Moreover, I used the bivalve Musculista senhousia as model taxon to investigate the molecular mechanisms and the evolutionary significance of their aberrant mode of mitochondrial inheritance (Doubly Uniparental Inheritance, see below). In Insects, I focused my attention on the Genus Bacillus (Insecta Phasmida). A detailed phylogenetic analysis was performed in order to assess phylogenetic relationships within the genus, and to investigate the placement of Phasmida in the phylogenetic tree of Insecta. The main goal of this part of my study was to add to the taxonomic coverage of sequenced mtDNAs in basal insects, which were only partially analyzed.

Transcriptional functions of DNA Topoisomerases at a genome-wide scale and a single gene levels

The DNA topology is an important modifier of DNA functions. Torsional stress is generated when right handed DNA is either over- or underwound, producing structural deformations which drive or are driven by processes such as replication, transcription, recombination and repair. DNA topoisomerases are molecular machines that regulate the topological state of the DNA in the cell. These enzymes accomplish this task by either passing one strand of the DNA through a break in the opposing strand or by passing a region of the duplex from the same or a different molecule through a double-stranded cut generated in the DNA. Because of their ability to cut one or two strands of DNA they are also target for some of the most successful anticancer drugs used in standard combination therapies of human cancers. An effective anticancer drug is Camptothecin (CPT) that specifically targets DNA topoisomerase 1 (TOP 1). The research project of the present thesis has been focused on the role of human TOP 1 during transcription and on the transcriptional consequences associated with TOP 1 inhibition by CPT in human cell lines. Previous findings demonstrate that TOP 1 inhibition by CPT perturbs RNA polymerase (RNAP II) density at promoters and along transcribed genes suggesting an involvement of TOP 1 in RNAP II promoter proximal pausing site. Within the transcription cycle, promoter pausing is a fundamental step the importance of which has been well established as a means of coupling elongation to RNA maturation. By measuring nascent RNA transcripts bound to chromatin, we demonstrated that TOP 1 inhibition by CPT can enhance RNAP II escape from promoter proximal pausing site of the human Hypoxia Inducible Factor 1 (HIF-1) and c-MYC genes in a dose dependent manner. This effect is dependent from Cdk7/Cdk9 activities since it can be reversed by the kinases inhibitor DRB. Since CPT affects RNAP II by promoting the hyperphosphorylation of its Rpb1 subunit the findings suggest that TOP 1inhibition by CPT may increase the activity of Cdks which in turn phosphorylate the Rpb1 subunit of RNAP II enhancing its escape from pausing. Interestingly, the transcriptional consequences of CPT induced topological stress are wider than expected. CPT increased co-transcriptional splicing of exon1 and 2 and markedly affected alternative splicing at exon 11. Surprisingly despite its well-established transcription inhibitory activity, CPT can trigger the production of a novel long RNA (5’aHIF-1) antisense to the human HIF-1 mRNA and a known antisense RNA at the 3’ end of the gene, while decreasing mRNA levels. The effects require TOP 1 and are independent from CPT induced DNA damage. Thus, when the supercoiling imbalance promoted by CPT occurs at promoter, it may trigger deregulation of the RNAP II pausing, increased chromatin accessibility and activation/derepression of antisense transcripts in a Cdks dependent manner. A changed balance of antisense transcripts and mRNAs may regulate the activity of HIF-1 and contribute to the control of tumor progression After focusing our TOP 1 investigations at a single gene level, we have extended the study to the whole genome by developing the “Topo-Seq” approach which generates a map of genome-wide distribution of sites of TOP 1 activity sites in human cells. The preliminary data revealed that TOP 1 preferentially localizes at intragenic regions and in particular at 5’ and 3’ ends of genes. Surprisingly upon TOP 1 downregulation, which impairs protein expression by 80%, TOP 1 molecules are mostly localized around 3’ ends of genes, thus suggesting that its activity is essential at these regions and can be compensate at 5’ ends. The developed procedure is a pioneer tool for the detection of TOP 1 cleavage sites across the genome and can open the way to further investigations of the enzyme roles in different nuclear processes.

IgE, IgGa, IgGb and IgG(T) serum antibody levels in offspring of two sires affected with equine recurrent airway obstruction

Equine recurrent airway obstruction (RAO) is a chronic lower airway disease of the horse caused by hypersensitivity reactions to inhaled stable dust, including mould spores such as Aspergillus fumigatus. The goals of this study were to investigate whether total serum IgE levels and allergen-specific IgE and IgG subclasses are influenced by genetic factors and/or RAO and whether quantitative trait loci (QTL) could be identified for these parameters. The offspring of two RAO-affected sires (S1: n=56 and S2: n=65) were grouped by stallion and disease status, and total serum IgE levels and specific IgE, IgGa, IgGb and IgG(T) levels against recombinant Aspergillus fumigatus 7 (rAspf7) were measured by ELISA. A panel of 315 microsatellite markers covering the 31 equine autosomes were used to genotype the stallions and their offspring. A whole-genome scan using half-sib regression interval mapping was performed for each of the IgG and IgE subclasses. There was no significant effect of disease status or sire on total IgE levels, but there was a significant effect of gender and age. rAspf7-specific IgGa levels were significantly higher in RAO-affected than in healthy horses. The offspring of S1 had significantly higher rAspf7-specific IgGa and IgE levels than those of S2. Five QTLs were significant chromosome-wide (P<0.01). QTLs for rAspf7-specific IgGa and IgE were identified on ECA 1, for rAspf7-specific IgGa and IgGb on ECA 24 and for rAspf7 IgGa on ECA 26. These results provide evidence for effects of disease status and genetics on allergen-specific IgGa and IgE.

Genome-wide analysis of rare copy number variations reveals PARK2 as a candidate gene for attention-deficit/hyperactivity disorder

Attention-deficit/hyperactivity disorder (ADHD) is a common, highly heritable neurodevelopmental disorder. Genetic loci have not yet been identified by genome-wide association studies. Rare copy number variations (CNVs), such as chromosomal deletions or duplications, have been implicated in ADHD and other neurodevelopmental disorders. To identify rare (frequency 1%) CNVs that increase the risk of ADHD, we performed a whole-genome CNV analysis based on 489 young ADHD patients and 1285 adult population-based controls and identified one significantly associated CNV region. In tests for a global burden of large (>500 kb) rare CNVs, we observed a nonsignificant (P=0.271) 1.126-fold enriched rate of subjects carrying at least one such CNV in the group of ADHD cases. Locus-specific tests of association were used to assess if there were more rare CNVs in cases compared with controls. Detected CNVs, which were significantly enriched in the ADHD group, were validated by quantitative (q)PCR. Findings were replicated in an independent sample of 386 young patients with ADHD and 781 young population-based healthy controls. We identified rare CNVs within the parkinson protein 2 gene (PARK2) with a significantly higher prevalence in ADHD patients than in controls (P=2.8 × 10(-4) after empirical correction for genome-wide testing). In total, the PARK2 locus (chr 6: 162 659 756-162 767 019) harboured three deletions and nine duplications in the ADHD patients and two deletions and two duplications in the controls. By qPCR analysis, we validated 11 of the 12 CNVs in ADHD patients (P=1.2 × 10(-3) after empirical correction for genome-wide testing). In the replication sample, CNVs at the PARK2 locus were found in four additional ADHD patients and one additional control (P=4.3 × 10(-2)). Our results suggest that copy number variants at the PARK2 locus contribute to the genetic susceptibility of ADHD. Mutations and CNVs in PARK2 are known to be associated with Parkinson disease.Molecular Psychiatry advance online publication, 20 November 2012; doi:10.1038/mp.2012.161.

Genome-wide search for markers associated with osteochondrosis in Hanoverian warmblood horses

A genome-wide scan was performed to detect quantitative trait loci (QTLs) for osteochondrosis (OC) and osteochondrosis dissecans (OCD) in horses. The marker set comprised 260 microsatellites. We collected data from 211 Hanoverian warmblood horses consisting of 14 paternal half-sib families. Traits used were OC (fetlock and/or hock joints affected), OCD (fetlock and/or hock joints affected), fetlock OC, fetlock OCD, hock OC, and hock OCD. The first genome scan included 172 microsatellite markers. In a second step 88 additional markers were chosen to refine putative QTLs found in the first scan. Genome-wide significant QTLs were located on equine chromosomes 2, 4, 5, and 16. QTLs for fetlock OC and hock OC partly overlapped on the same chromosomes, indicating that these traits may be genetically related. QTLs reached the chromosome-wide significance level on eight different equine chromosomes: 2, 3, 4, 5, 15, 16, 19, and 21. This whole-genome scan was a first step toward the identification of candidate genome regions harboring genes responsible for equine OC. Further investigations are necessary to refine the map positions of the QTLs already identified for OC.

A radiation hybrid map of sheep chromosome 23 based on ovine BAC-end sequences

More than 375,000 BAC-end sequences (BES) of the CHORI-243 ovine BAC library have been deposited in public databases. blastn searches with these BES against HSA18 revealed 1806 unique and significant hits. We used blastn-anchored BES for an in silico prediction of gene content and chromosome assignment of comparatively mapped ovine BAC clones. Ovine BES were selected at approximately 1.3-Mb intervals of HSA18 and incorporated into a human-sheep comparative map. An ovine 5000-rad whole-genome radiation hybrid panel (USUoRH5000) was typed with 70 markers, all of which mapped to OAR23. The resulting OAR23 RH map included 43 markers derived from BES with high and unique BLAST hits to the sequence of the orthologous HSA18, nine EST-derived markers, 16 microsatellite markers taken from the ovine linkage map and two bovine microsatellite markers. Six new microsatellite markers derived from the 43 mapped BES and the two bovine microsatellite markers were linkage-mapped using the International Mapping Flock (IMF). Thirteen additional microsatellite markers were derived from other ovine BES with high and unique BLAST hits to the sequence of the orthologous HSA18 and also positioned on the ovine linkage map but not incorporated into the OAR23 RH map. This resulted in 24 markers in common and in the same order between the RH and linkage maps. Eight of the BES-derived markers were mapped using fluorescent in situ hybridization (FISH), to thereby align the RH and cytogenetic maps. Comparison of the ovine chromosome 23 RH map with the HSA18 map identified and localized three major breakpoints between HSA18 and OAR23. The positions of these breakpoints were equivalent to those previously shown for syntenic BTA24 and HSA18. This study presents evidence for the usefulness of ovine BES when constructing a high-resolution comprehensive map for a single sheep chromosome. The comparative analysis confirms and refines knowledge about chromosomal conservation and rearrangements between sheep, cattle and human. The constructed RH map demonstrates the resolution and utility of the newly constructed ovine RH panel.

The locus for bovine dilated cardiomyopathy maps to chromosome 18

Bovine dilated cardiomyopathy (BDCMP) is a severe and terminal disease of the heart muscle observed in Holstein-Friesian cattle over the last 30 years. There is strong evidence for an autosomal recessive mode of inheritance for BDCMP. The objective of this study was to genetically map BDCMP, with the ultimate goal of identifying the causative mutation. A whole-genome scan using 199 microsatellite markers and one SNP revealed an assignment of BDCMP to BTA18. Fine-mapping on BTA18 refined the candidate region to the MSBDCMP06-BMS2785 interval. The interval containing the BDCMP locus was confirmed by multipoint linkage analysis using the software loki. The interval is about 6.7 Mb on the bovine genome sequence (Btau 3.1). The corresponding region of HSA19 is very gene-rich and contains roughly 200 genes. Although telomeric of the marker interval, TNNI3 is a possible positional and a functional candidate for BDCMP given its involvement in a human form of dilated cardiomyopathy. Sequence analysis of TNNI3 in cattle revealed no mutation in the coding sequence, but there was a G-to-A transition in intron 6 (AJ842179:c.378+315G>A). The analysis of this SNP using the study's BDCMP pedigree did not conclusively exclude TNNI3 as a candidate gene for BDCMP. Considering the high density of genes on the homologous region of HSA19, further refinement of the interval on BTA18 containing the BDCMP locus is needed.

Functional Diversification of Dicer-like Proteins and Small RNAs Required for Genome Sculpting.

In eukaryotes, small RNAs (sRNAs) have key roles in development, gene expression regulation, and genome integrity maintenance. In ciliates, such as Paramecium, sRNAs form the heart of an epigenetic system that has evolved from core eukaryotic gene silencing components to selectively target DNA for deletion. In Paramecium, somatic genome development from the germline genome accurately eliminates the bulk of typically gene-interrupting, noncoding DNA. We have discovered an sRNA class (internal eliminated sequence [IES] sRNAs [iesRNAs]), arising later during Paramecium development, which originates from and precisely delineates germline DNA (IESs) and complements the initial sRNAs ("scan" RNAs [scnRNAs]) in targeting DNA for elimination. We show that whole-genome duplications have facilitated successive differentiations of Paramecium Dicer-like proteins, leading to cooperation between Dcl2 and Dcl3 to produce scnRNAs and to the production of iesRNAs by Dcl5. These innovations highlight the ability of sRNA systems to acquire capabilities, including those in genome development and integrity.

Deletion in the EVC2 gene causes chondrodysplastic dwarfism in Tyrolean Grey cattle.

During the summer of 2013 seven Italian Tyrolean Grey calves were born with abnormally short limbs. Detailed clinical and pathological examination revealed similarities to chondrodysplastic dwarfism. Pedigree analysis showed a common founder, assuming autosomal monogenic recessive transmission of the defective allele. A positional cloning approach combining genome wide association and homozygosity mapping identified a single 1.6 Mb genomic region on BTA 6 that was associated with the disease. Whole genome re-sequencing of an affected calf revealed a single candidate causal mutation in the Ellis van Creveld syndrome 2 (EVC2) gene. This gene is known to be associated with chondrodysplastic dwarfism in Japanese Brown cattle, and dwarfism, abnormal nails and teeth, and dysostosis in humans with Ellis-van Creveld syndrome. Sanger sequencing confirmed the presence of a 2 bp deletion in exon 19 (c.2993_2994ACdel) that led to a premature stop codon in the coding sequence of bovine EVC2, and was concordant with the recessive pattern of inheritance in affected and carrier animals. This loss of function mutation confirms the important role of EVC2 in bone development. Genetic testing can now be used to eliminate this form of chondrodysplastic dwarfism from Tyrolean Grey cattle.

Looking the cow in the eye: deletion in the NID1 gene is associated with recessive inherited cataract in Romagnola cattle

Cataract is a known condition leading to opacification of the eye lens causing partial or total blindness. Mutations are known to cause autosomal dominant or recessive inherited forms of cataracts in humans, mice, rats, guinea pigs and dogs. The use of large-sized animal models instead of those using mice for the study of this condition has been discussed due to the small size of rodent lenses. Four juvenile-onset cases of bilateral incomplete immature nuclear cataract were recently observed in Romagnola cattle. Pedigree analysis suggested a monogenic autosomal recessive inheritance. In addition to the cataract, one of the cases displayed abnormal head movements. Genome-wide association and homozygosity mapping and subsequent whole genome sequencing of a single case identified two perfectly associated sequence variants in a critical interval of 7.2 Mb on cattle chromosome 28: a missense point mutation located in an uncharacterized locus and an 855 bp deletion across the exon 19/intron 19 border of the bovine nidogen 1 (NID1) gene (c.3579_3604+829del). RT-PCR showed that NID1 is expressed in bovine lenses while the transcript of the second locus was absent. The NID1 deletion leads to the skipping of exon 19 during transcription and is therefore predicted to cause a frameshift and premature stop codon (p.1164fs27X). The truncated protein lacks a C-terminal domain essential for binding with matrix assembly complexes. Nidogen 1 deficient mice show neurological abnormalities and highly irregular crystal lens alterations. This study adds NID1 to the list of candidate genes for inherited cataract in humans and is the first report of a naturally occurring mutation leading to non-syndromic catarct in cattle provides a potential large animal model for human cataract.

A 4,103 marker integrated physical and comparative map of the horse genome

A comprehensive second-generation whole genome radiation hybrid (RH II), cytogenetic and comparative map of the horse genome (2n = 64) has been developed using the 5000rad horse x hamster radiation hybrid panel and fluorescence in situ hybridization (FISH). The map contains 4,103 markers (3,816 RH; 1,144 FISH) assigned to all 31 pairs of autosomes and the X chromosome. The RH maps of individual chromosomes are anchored and oriented using 857 cytogenetic markers. The overall resolution of the map is one marker per 775 kilobase pairs (kb), which represents a more than five-fold improvement over the first-generation map. The RH II incorporates 920 markers shared jointly with the two recently reported meiotic maps. Consequently the two maps were aligned with the RH II maps of individual autosomes and the X chromosome. Additionally, a comparative map of the horse genome was generated by connecting 1,904 loci on the horse map with genome sequences available for eight diverse vertebrates to highlight regions of evolutionarily conserved syntenies, linkages, and chromosomal breakpoints. The integrated map thus obtained presents the most comprehensive information on the physical and comparative organization of the equine genome and will assist future assemblies of whole genome BAC fingerprint maps and the genome sequence. It will also serve as a tool to identify genes governing health, disease and performance traits in horses and assist us in understanding the evolution of the equine genome in relation to other species.