Genomic analysis for managing small and endangered populations: a case study in Tyrol Grey cattle

Analysis of genomic data is increasingly becoming part of the livestock industry. Therefore, the routine collection of genomic information would be an invaluable resource for effective management of breeding programs in small, endangered populations. The objective of the paper was to demonstrate how genomic data could be used to analyse (1) linkage disequlibrium (LD), LD decay and the effective population size (NeLD); (2) Inbreeding level and effective population size (NeROH) based on runs of homozygosity (ROH); (3) Prediction of genomic breeding values (GEBV) using small within-breed and genomic information from other breeds. The Tyrol Grey population was used as an example, with the goal to highlight the potential of genomic analyses for small breeds. In addition to our own results we discuss additional use of genomics to assess relatedness, admixture proportions, and inheritance of harmful variants. The example data set consisted of 218 Tyrol Grey bull genotypes, which were all available AI bulls in the population. After standard quality control restrictions 34,581 SNPs remained for the analysis. A separate quality control was applied to determine ROH levels based on Illumina GenCall and Illumina GenTrain scores, resulting into 211 bulls and 33,604 SNPs. LD was computed as the squared correlation coefficient between SNPs within a 10 mega base pair (Mb) region. ROHs were derived based on regions covering at least 4, 8, and 16 Mb, suggesting that animals had common ancestors approximately 12, 6, and 3 generations ago, respectively. The corresponding mean inbreeding coefficients (F ROH) were 4.0% for 4 Mb, 2.9% for 8 Mb and 1.6% for 16 Mb runs. With an average generation interval of 5.66 years, estimated NeROH was 125 (NeROH>16 Mb), 186 (NeROH>8 Mb) and 370 (NeROH>4 Mb) indicating strict avoidance of close inbreeding in the population. The LD was used as an alternative method to infer the population history and the Ne. The results show a continuous decrease in NeLD, to 780, 120, and 80 for 100, 10, and 5 generations ago, respectively. Genomic selection was developed for and is working well in large breeds. The same methodology was applied in Tyrol Grey cattle, using different reference populations. Contrary to the expectations, the accuracy of GEBVs with very small within breed reference populations were very high, between 0.13-0.91 and 0.12-0.63, when estimated breeding values and deregressed breeding values were used as pseudo-phenotypes, respectively. Subsequent analyses confirmed the high accuracies being a consequence of low reliabilities of pseudo-phenotypes in the validation set, thus being heavily influenced by parent averages. Multi-breed and across breed reference sets gave inconsistent and lower accuracies. Genomic information may have a crucial role in management of small breeds, even if its primary usage differs from that of large breeds. It allows to assess relatedness between individuals, trends in inbreeding and to take decisions accordingly. These decisions would be based on the real genome architecture, rather than conventional pedigree information, which can be missing or incomplete. We strongly suggest the routine genotyping of all individuals that belong to a small breed in order to facilitate the effective management of endangered livestock populations.

Ferramenta computacional para identificação de micro-organismos com base em assinaturas genômicas

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

DNA fragmentation by gamma radiation and electron beams using atomic force microscopy

Double-stranded pBS plasmid DNA was irradiated with gamma rays at doses ranging from 1 to 12 kGy and electron beams from 1 to 10 kGy. Fragment-size distributions were determined by direct visualization, using atomic force microscopy with nanometer-resolution operating in non-tapping mode, combined with an improved methodology. The fragment distributions from irradiation with gamma rays revealed discrete-like patterns at all doses, suggesting that these patterns are modulated by the base pair composition of the plasmid. Irradiation with electron beams, at very high dose rates, generated continuous distributions of highly shattered DNA fragments, similar to results at much lower dose rates found in the literature. Altogether, these results indicate that AFM could supplement traditional methods for high-resolution measurements of radiation damage to DNA, while providing new and relevant information.

Synthesis, base pairing properties and trans-lesion synthesis by reverse transcriptases of oligoribonucleotides containing the oxidatively damaged base 5-hydroxycytidine

The synthesis of a caged RNA phosphoramidite building block containing the oxidatively damaged base 5-hydroxycytidine (5-HOrC) has been accomplished. To determine the effect of this highly mutagenic lesion on complementary base recognition and coding properties, this building block was incorporated into a 12-mer oligoribonucleotide for Tm and CD measurements and a 31-mer template strand for primer extension experiments with HIV-, AMV- and MMLV-reverse transcriptase (RT). In UV-melting experiments, we find an unusual biphasic transition with two distinct Tm's when 5-HOrC is paired against a DNA or RNA complement with the base guanine in opposing position. The higher Tm closely matches that of a C-G base pair while the lower is close to that of a C-A mismatch. In single nucleotide extension reactions, we find substantial misincorporation of dAMP and to a lesser extent dTMP, with dAMP almost equaling that of the parent dGMP in the case of HIV-RT. A working hypothesis for the biphasic melting transition does not invoke tautomeric variability of 5-HOrC but rather local structural perturbations of the base pair at low temperature induced by interactions of the 5-HO group with the phosphate backbone. The properties of this RNA damage is discussed in the context of its putative biological function.

Aromatase deficiency in a female who is compound heterozygote for two new point mutations in the P450arom gene: impact of estrogens on hypergonadotropic hypogonadism, multicystic ovaries, and bone densitometry in childhood

We report on a female who is compound heterozygote for two new point mutations in the CYP19 gene. The allele inherited from her mother presented a base pair deletion (C) occurring at P408 (CCC, exon 9), causing a frameshift that results in a nonsense codon 111 bp (37 aa) further down in the CYP19 gene. The allele inherited from her father showed a point mutation from G-->A at the splicing point (canonical GT to mutational AT) between exon and intron 3. This mutation ignores the splice site and a stop codon 3 bp downstream occurs. Aromatase deficiency was already suspected because of the marked virilization occurring prepartum in the mother, and the diagnosis was confirmed shortly after birth. Extremely low levels of serum estrogens were found in contrast to high levels of androgens. Ultrasonographic follow-up studies revealed persistently enlarged ovaries (19.5-22 mL) during early childhood (2 to 4 yr) which contained numerous large cysts up to 4.8 x 3.7 cm and normal-appearing large tertiary follicles already at the age of 2 yr. In addition, both basal and GnRH-induced FSH levels remained consistently strikingly elevated. Low-dose estradiol (E2) (0.4 mg/day) given for 50 days at the age of 3 6/12 yr resulted in normalization of serum gonadotropin levels, regression of ovarian size, and increase of whole body and lumbar spine (L1-L4) bone mineral density. The FSH concentration and ovarian size returned to pretreatment levels shortly (150 days) after cessation of E2 therapy. Therefore, we recommend that affected females be treated with low-dose E2 in amounts sufficient to result in physiological prepubertal E2 concentrations using an ultrasensitive estrogen assay. However, E2 replacement needs to be adjusted throughout childhood and puberty to ensure normal skeletal maturation and adequate adolescent growth spurt, normal accretion of bone mineral density, and, at the appropriate age, female secondary sex maturation.

Analysis of mutations in $\beta$-tubulin genes affecting tubulin stability and assembly

Cmd4 is a colcemid-sensitive CHO cell line that is temperature sensitive for growth and expresses an altered $\beta$-tubulin, $\beta\sb1$. One revertant of this cell line, D2, exhibits a further alteration in $\beta\sb1$ resulting in an acidic shift in its isoelectric point and a decrease in its molecular weight to 40 kD, as measured by two dimensional gel electrophoresis. This $\beta$-tubulin variant has been shown to be assembly-defective and unstable. Characterization of the mutant $\beta\sb1$ in D2 by high pressure liquid chromatography (HPLC) revealed the loss of methionine containing tryptic peptides 7,8,9, and 10. Southern analysis of the genomic DNA digested with several different restriction enzymes resulted in the appearance of new restriction fragments 250 base pairs shorter than the corresponding fragments from the wild-type $\beta\sb1$-tubulin gene. Northern analysis on mRNA from D2 revealed two new message products that also differed by 250 bases from the corresponding wild type $\beta$-tubulin transcripts. To precisely define the region of the alteration, cloning and sequencing of the mutant and wild type genomic $\beta$-tubulin genes were conducted. A size-selected EcoRI genomic library was prepared using the Stratagene lambda Zap II phage cloning system. Using subclones of CHO $\beta$-tubulin cDNA as probes, a 2.5 kb wild type clone and a 2.3 kb mutant clone were identified from this library. Each of these was shown to contain a portion of the gene extending from intron 3 through the end of the coding sequence in exon 4 and into the 3$\sp\prime$ untranslated region on the basis of alignment with the published human $\beta$-tubulin sequence. Sequencing of the mutant 2.3 kb clone revealed that the mutation is due to a 246 base pair internal deletion in exon 4 (base pair 756-1001) that encodes amino acids 253-334. This deletion results in the loss of a putative binding site for GTP which could potentially explain the phenotype of this mutant $\beta$-tubulin. Also sequence comparison of the 3$\sp\prime$ untranslated region between different species revealed the conservation of 200 base pairs with 78% homology. It is proposed that this region could play an important role in the regulation of $\beta$-tubulin gene expression. ^

Transition metal ligands as novel DNA-base substitutes

The base modified nucleoside dBP, carrying a non-hydrogen-bonding non-shape complementary base was incorporated into oligonucleotides (Brotschi, C.; Haberli, A.; Leumann C.J. Angew. Chem. Int. Ed. 2001, 40, 3012-3014). This base was designed to coordinate transition metal ions into well defined positions within a DNA double helix. Melting experiments revealed that the stability of a dBP:dBP base couple in a DNA duplex is similar to a dG:dC base pair even in the absence of transition metal ions. In the presence of transition metal ions, melting experiments revealed a decrease in duplex stability which is on a similar order for all metal ions (Mn2+, Cu2+, Zn2+, Ni2+) tested

Watson-Crick base-pairing properties of tricyclo-DNA

Tricyclo-DNA belongs to the family of conformationally restricted oligodeoxynucleotide analogues. It differs structurally from DNA by an additional ethylene bridge between the centers C(3') and C(5') of the nucleosides, to which a cyclopropane unit is fused for further enhancement of structural rigidity. The synthesis of the hitherto unknown tricyclodeoxynucleosides containing the bases cytosine and guanine and of the corresponding phosphoramidite building blocks is described, as well as a structural description of a representative of an alpha- and a beta-tricyclodeoxynucleoside by X-ray analysis. Tricyclodeoxynucleoside building blocks of all four bases were used for the synthesis of fully modified mixed-base oligonucleotides. Their Watson-Crick pairing properties with complementary DNA, RNA, and with itself were investigated by UV melting curves, CD spectroscopy, and molecular modeling. Tricyclo-DNA was found to be a very stable Watson-Crick base-pairing system. A UV melting curve analysis of the decamers tcd(pcgtgacagtt) and tcd(paactgtcacg) showed increased thermal stabilities of up to DeltaT(m)/mod. = +1.2 degrees C with complementary DNA and +2.4 degrees C with complementary RNA. With itself, tricyclo-DNA showed an increase in stability of +3.1 degrees C/base pair relative to DNA. Investigations into the thermodynamic properties of these decamers revealed an entropic stabilization and an enthalpic destabilization for the tricyclo-DNA/DNA duplexes. CD spectroscopic structural investigations indicated that tricyclo-DNA containing duplexes preferrably exist in an A-conformation, a fact which is in agreement with results from molecular modeling

X-ray structure of a lectin-bound DNA duplex containing an unnatural phenanthrenyl pair

DNA duplexes containing unnatural base-pair surrogates are attractive biomolecular nanomaterials with potentially beneficial photophysical or electronic properties. Herein we report the first X-ray structure of a duplex containing a phen-pair in the center of the double helix in a zipper like stacking arrangement.

GLUCAGON GENE EXPRESSION: ANALYSIS OF PREPROGLUCAGON

Glucagon is a 29 amino acid polypeptide hormone produced in the (alpha) cells of the pancreatic islets. The purpose of this research was to understand better the role of glucagon in the regulation of metabolic processes. As with other polypeptide hormones, the synthesis of glucagon is thought to involve a larger precursor, which is then enzymatically cleaved to the functional form. The specific research objectives were to obtain cloned copies of the messenger RNA (mRNA) for pancreatic glucagon, to determine their primary sequences, and from this coding information to deduce the amino acid sequence of the initial glucagon precursor. From this suggested preproglucagon sequence and prior information on possible proglucagon intermediate processing products, the overall objective of this research is to propose a possible pathway for the biosynthesis of pancreatic glucagon.^ Synthetic oligodeoxynucleotide probes of 14-nucleotides (14-mer) and 17-nucleotides (a 17-mer) complementary to codons specifying a unique sequence of mature glucagon were synthesized. The ('32)P-labeled-14-mer was hybridized with size-fractionated fetal bovine pancreatic poly(A('+))RNA bound to nitrocellulose. RNA fractions of (TURN)14S were found to hybridize specifically, resulting in an (TURN)10-fold enrichment for these sequences. These poly(A('+))RNAs were translated in a cell-free system and the products analyzed by gel electrophoresis. The translation products were found to be enriched for a protein of the putative size of mammalian preproglucagon ((TURN)21 kd). These enriched RNA fractions were used to construct a complementary DNA (cDNA) library is plasmid pBR322.^ Screening of duplicate colony filters with the ('32)P-labeled-17-mer and a ('32)P-labeled-17-mer-primed cDNA probe indicated 25 possible glucagon clones from 3100 colonies screened. Restriction mapping of 6 of these clones suggested that they represented a single mRNA species. Primary sequence analysis of one clone containing a 1200 base pair DNA insert revealed that it contained essentially a full-length copy of glucagon cDNA.^ Analaysis of the cDNA suggested that it encoded an initial translation product of 180 amino acids with an M(,r) = 21 kd. The first initiation codon (ATG, methionine) followed by the longest open reading frame of 540 nucleotides was preceded by a 5'-untranslated region of 90 nucleotides, and was followed by a longer 3'-untranslated region of 471 nucleotides, resulting in a total of 1101 nucleotides. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author.) UMI ^

The crystal structure of the Rev binding element of HIV-1 reveals novel base pairing and conformational variability

The crystal and molecular structure of an RNA duplex corresponding to the high affinity Rev protein binding element (RBE) has been determined at 2.1-Å resolution. Four unique duplexes are present in the crystal, comprising two structural variants. In each duplex, the RNA double helix consists of an annealed 12-mer and 14-mer that form an asymmetric internal loop consisting of G-G and G-A noncanonical base pairs and a flipped-out uridine. The 12-mer strand has an A-form conformation, whereas the 14-mer strand is distorted to accommodate the bulges and noncanonical base pairing. In contrast to the NMR model of the unbound RBE, an asymmetric G-G pair with N2-N7 and N1-O6 hydrogen bonding, is formed in each helix. The G-A base pairing agrees with the NMR structure in one structural variant, but forms a novel water-mediated pair in the other. A backbone flip and reorientation of the G-G base pair is required to assume the RBE conformation present in the NMR model of the complex between the RBE and the Rev peptide.

Recognition of nonhybridizing base pairs during nucleotide excision repair of DNA

Nondistorting C4′ backbone adducts serve as molecular tools to analyze the strategy by which a limited number of human nucleotide excision repair (NER) factors recognize an infinite variety of DNA lesions. We have constructed composite DNA substrates containing a noncomplementary site adjacent to a nondistorting C4′ adduct to show that the loss of hydrogen bonding contacts between partner strands is an essential signal for the recruitment of NER enzymes. This specific conformational requirement for excision is mediated by the affinity of xeroderma pigmentosum group A (XPA) protein for nonhybridizing sites in duplex DNA. XPA recognizes defective Watson–Crick base pair conformations even in the absence of DNA adducts or other covalent modifications, apparently through detection of hydrophobic base components that are abnormally exposed to the double helical surface. This recognition function of XPA is enhanced by replication protein A (RPA) such that, in combination, XPA and RPA constitute a potent molecular sensor of denatured base pairs. Our results indicate that the XPA–RPA complex may promote damage recognition by monitoring Watson–Crick base pair integrity, thereby recruiting the human NER system preferentially to sites where hybridization between complementary strands is weakened or entirely disrupted.

Unnatural base pairs for specific transcription

An unnatural base pair of 2-amino-6-(N,N-dimethylamino)purine (designated as x) and pyridin-2-one (designated as y) has been developed for specific transcription. The ribonucleoside triphosphates of y and a modified y, 5-methylpyridin-2-one, are selectively incorporated into RNA opposite x in the templates by T7 RNA polymerase. In addition, the sequences of the DNA templates containing x can be confirmed by a dideoxynucleotide chain-terminator method supplemented with the deoxynucleoside triphosphate of y. The bulky dimethylamino group of x in the templates effectively eliminates noncognate pairing with the natural bases. These results enable RNA biosynthesis for the specific incorporation of unnatural nucleotides at the desired positions.

Two new blood flukes (Digenea: Sanguinicolidae) from epinephelinae (Perciformes: Serranidae) of the Pacific Ocean

Pearsonellum pygmaeus n. sp. is described from Cromileptes altivelis (Serranidae), the Barramundi Cod, from Heron Island (southern Great Barrier Reef) and Lizard Island (northern Great Bat-Her Reef). This new species differs from Pearsonellum eorventum (type and only species) in the combination of smaller overall body size, the relative distance of the brain from the anterior end, the relative lengths of both the oesophagus and the testis, the degree to which the testis extends outside the intercaecal field, the shape of the testis, the shape and size of the ovary and the extent to which the uterzus loops around the ovary. There are in addition, 20 base pair differences between the ITS2 rDNA sequence of P. pygmaeus n. sp. and that of P corventum. Three new host records for P. corventum are reported. Adelomyllos teenae n. g., n. sp. is described from Epinephelus coioides (Serranidae), the Estuary Cod, from Moreton Bay, southeast Queensland. The new genus differs from the 22 other sanguinicolid genera in the combined possession of two testes, a cirrus-sac, separate genital pores, a post-ovarian uterus and an H-shaped intestine. A. teenae n. sp. is the third sanguinicolid described from the Epinephelinae. Sanguinicolids have now been reported from 11 species of Serranidae. (C) 2004 Elsevier Ireland Ltd. All rights reserved.

STRUCTURAL ANALYSIS OF RIBOSOME BINDING ELEMENTS OF THE 3´ UTR IN TWO POSITIVE-STRAND PLANT RNA VIRUSES

Turnip crinkle virus (TCV) and Pea enation mosaic virus (PEMV) are two positive (+)-strand RNA viruses that are used to investigate the regulation of translation and replication due to their small size and simple genomes. Both viruses contain cap-independent translation elements (CITEs) within their 3´ untranslated regions (UTRs) that fold into tRNA-shaped structures (TSS) according to nuclear magnetic resonance and small angle x-ray scattering analysis (TCV) and computational prediction (PEMV). Specifically, the TCV TSS can directly associate with ribosomes and participates in RNA-dependent RNA polymerase (RdRp) binding. The PEMV kissing-loop TSS (kl-TSS) can simultaneously bind to ribosomes and associate with the 5´ UTR of the viral genome. Mutational analysis and chemical structure probing methods provide great insight into the function and secondary structure of the two 3´ CITEs. However, lack of 3-D structural information has limited our understanding of their functional dynamics. Here, I report the folding dynamics for the TCV TSS using optical tweezers (OT), a single molecule technique. My study of the unfolding/folding pathways for the TCV TSS has provided an unexpected unfolding pathway, confirmed the presence of Ψ3 and hairpin elements, and suggested an interconnection between the hairpins and pseudoknots. In addition, this study has demonstrated the importance of the adjacent upstream adenylate-rich sequence for the formation of H4a/Ψ3 along with the contribution of magnesium to the stability of the TCV TSS. In my second project, I report on the structural analysis of the PEMV kl-TSS using NMR and SAXS. This study has re-confirmed the base-pair pattern for the PEMV kl-TSS and the proposed interaction of the PEMV kl-TSS with its interacting partner, hairpin 5H2. The molecular envelope of the kl-TSS built from SAXS analysis suggests the kl-TSS has two functional conformations, one of which has a different shape from the previously predicted tRNA-shaped form. Along with applying biophysical methods to study the structural folding dynamics of RNAs, I have also developed a technique that improves the production of large quantities of recombinant RNAs in vivo for NMR study. In this project, I report using the wild-type and mutant E.coli strains to produce cost-effective, site-specific labeled, recombinant RNAs. This technique was validated with four representative RNAs of different sizes and complexity to produce milligram amounts of RNAs. The benefit of using site-specific labeled RNAs made from E.coli was demonstrated with several NMR techniques.