A short critical history of the application of genomics to animal breeding

Two scientific schools have been in coexistence from the beginning of genetics, one of them searching for factors of inheritance and the other one applying biometrical models to study the relationships between relatives. With the development of molecular genetics, the possibilities of detecting genes having a noticeable effect in traits augmented. Some genes with large or medium effects were localized in animals, although the most common result was to detect markers linked to these genes, allowing the possibility of assisting selection programs with markers. When a large amount of simple and inexpensive markers were available, the SNPs, new possibilities were opened since they did not need the presence of genes of large or medium effect controlling a trait, because the whole genome was scanned. Using a large amount of SNPs permits having a prediction of the breeding value at birth accurate enough to be used in some cases, like dairy cattle, to halve its generation interval. In other animal breeding programs, the implementation of genomic selection is less clear and the way in which it can be useful should be carefully studied. The need for large populations for associating phenotypic data and markers, plus the need for repeating the process continuously, complicates its application in some cases. The implementation of the information provided by the SNPs in current genetic programs has led to the development of complex statistical tools, joining the efforts of the two schools, factorial and biometrical, that nowadays work closely related.

Mouse mammary tumor virus/v-Ha-ras transgene-induced mammary tumors exhibit strain-specific allelic loss on mouse chromosome 4

Hybrid mice carrying oncogenic transgenes afford powerful systems for investigating loss of heterozygosity (LOH) in tumors. Here, we apply this approach to a neoplasm of key importance in human medicine: mammary carcinoma. We performed a whole genome search for LOH using the mouse mammary tumor virus/v-Ha-ras mammary carcinoma model in female (FVB/N × Mus musculus castaneus)F1 mice. Mammary tumors developed as expected, as well as a few tumors of a second type (uterine leiomyosarcoma) not previously associated with this transgene. Genotyping of 94 anatomically independent tumors revealed high-frequency LOH (≈38%) for markers on chromosome 4. A marked allelic bias was observed, with M. musculus castaneus alleles almost exclusively being lost. No evidence of genomic imprinting effects was noted. These data point to the presence of a tumor suppressor gene(s) on mouse chromosome 4 involved in mammary carcinogenesis induced by mutant H-ras expression, and for which a significant functional difference may exist between the M. musculus castaneus and FVB/N alleles. Provisional subchromosomal localization of this gene, designated Loh-3, can be made to a distal segment having syntenic correspondence to human chromosome 1p; LOH in this latter region is observed in several human malignancies, including breast cancers. Evidence was also obtained for a possible second locus associated with LOH with less marked allele bias on proximal chromosome 4.

Gene content phylogeny of herpesviruses

Clusters of orthologous groups [COGs; Tatusov, R. L., Koonin, E. V. & Lipman, D. J. (1997) Science 278, 631–637] were identified for a set of 13 completely sequenced herpesviruses. Each COG represented a family of gene products conserved across several herpes genomes. These families were defined without using an arbitrary threshold criterion based on sequence similarity. The COG technique was modified so that variable stringency in COG construction was possible. High stringencies identify a core set of highly conserved genes. Varying COG stringency reveals differences in the degree of conservation between functional classes of genes. The COG data were used to construct whole-genome phylogenetic trees based on gene content. These trees agree well with trees based on other methods and are robust when tested by bootstrap analysis. The COG data also were used to construct a reciprocal tree that clustered genes with similar phylogenetic profiles. This clustering may give clues to genes with related functions or with related histories of acquisition and loss during herpesvirus evolution.

Y chromosomal fertility factors kl-2 and kl-3 of Drosophila melanogaster encode dynein heavy chain polypeptides

The molecular identity and function of the Drosophila melanogaster Y-linked fertility factors have long eluded researchers. Although the D. melanogaster genome sequence was recently completed, the fertility factors still were not identified, in part because of low cloning efficiency of heterochromatic Y sequences. Here we report a method for iterative blast searching to assemble heterochromatic genes from shotgun assemblies, and we successfully identify kl-2 and kl-3 as 1β- and γ-dynein heavy chains, respectively. Our conclusions are supported by formal genetics with X-Y translocation lines. Reverse transcription–PCR was successful in linking together unmapped sequence fragments from the whole-genome shotgun assembly, although some sequences were missing altogether from the shotgun effort and had to be generated de novo. We also found a previously undescribed Y gene, polycystine-related (PRY). The closest paralogs of kl-2, kl-3, and PRY (and also of kl-5) are autosomal and not X-linked, suggesting that the evolution of the Drosophila Y chromosome has been driven by an accumulation of male-related genes arising de novo from the autosomes.

The Comprehensive Microbial Resource

One challenge presented by large-scale genome sequencing efforts is effective display of uniform information to the scientific community. The Comprehensive Microbial Resource (CMR) contains robust annotation of all complete microbial genomes and allows for a wide variety of data retrievals. The bacterial information has been placed on the Web at http://www.tigr.org/CMR for retrieval using standard web browsing technology. Retrievals can be based on protein properties such as molecular weight or hydrophobicity, GC-content, functional role assignments and taxonomy. The CMR also has special web-based tools to allow data mining using pre-run homology searches, whole genome dot-plots, batch downloading and traversal across genomes using a variety of datatypes.

Complete genomic sequence of Pasteurella multocida,Pm70

We present here the complete genome sequence of a common avian clone of Pasteurella multocida, Pm70. The genome of Pm70 is a single circular chromosome 2,257,487 base pairs in length and contains 2,014 predicted coding regions, 6 ribosomal RNA operons, and 57 tRNAs. Genome-scale evolutionary analyses based on pairwise comparisons of 1,197 orthologous sequences between P. multocida, Haemophilus influenzae, and Escherichia coli suggest that P. multocida and H. influenzae diverged ≈270 million years ago and the γ subdivision of the proteobacteria radiated about 680 million years ago. Two previously undescribed open reading frames, accounting for ≈1% of the genome, encode large proteins with homology to the virulence-associated filamentous hemagglutinin of Bordetella pertussis. Consistent with the critical role of iron in the survival of many microbial pathogens, in silico and whole-genome microarray analyses identified more than 50 Pm70 genes with a potential role in iron acquisition and metabolism. Overall, the complete genomic sequence and preliminary functional analyses provide a foundation for future research into the mechanisms of pathogenesis and host specificity of this important multispecies pathogen.

Grass genomes

For the most part, studies of grass genome structure have been limited to the generation of whole-genome genetic maps or the fine structure and sequence analysis of single genes or gene clusters. We have investigated large contiguous segments of the genomes of maize, sorghum, and rice, primarily focusing on intergenic spaces. Our data indicate that much (>50%) of the maize genome is composed of interspersed repetitive DNAs, primarily nested retrotransposons that insert between genes. These retroelements are less abundant in smaller genome plants, including rice and sorghum. Although 5- to 200-kb blocks of methylated, presumably heterochromatic, retrotransposons flank most maize genes, rice and sorghum genes are often adjacent. Similar genes are commonly found in the same relative chromosomal locations and orientations in each of these three species, although there are numerous exceptions to this collinearity (i.e., rearrangements) that can be detected at the levels of both the recombinational map and cloned DNA. Evolutionarily conserved sequences are largely confined to genes and their regulatory elements. Our results indicate that a knowledge of grass genome structure will be a useful tool for gene discovery and isolation, but the general rules and biological significance of grass genome organization remain to be determined. Moreover, the nature and frequency of exceptions to the general patterns of grass genome structure and collinearity are still largely unknown and will require extensive further investigation.

Comparative organization of cattle chromosome 5 revealed by comparative mapping by annotation and sequence similarity and radiation hybrid mapping

A whole genome cattle-hamster radiation hybrid cell panel was used to construct a map of 54 markers located on bovine chromosome 5 (BTA5). Of the 54 markers, 34 are microsatellites selected from the cattle linkage map and 20 are genes. Among the 20 mapped genes, 10 are new assignments that were made by using the comparative mapping by annotation and sequence similarity strategy. A LOD-3 radiation hybrid framework map consisting of 21 markers was constructed. The relatively low retention frequency of markers on this chromosome (19%) prevented unambiguous ordering of the other 33 markers. The length of the map is 398.7 cR, corresponding to a ratio of ≈2.8 cR5,000/cM. Type I genes were binned for comparison of gene order among cattle, humans, and mice. Multiple internal rearrangements within conserved syntenic groups were apparent upon comparison of gene order on BTA5 and HSA12 and HSA22. A similarly high number of rearrangements were observed between BTA5 and MMU6, MMU10, and MMU15. The detailed comparative map of BTA5 should facilitate identification of genes affecting economically important traits that have been mapped to this chromosome and should contribute to our understanding of mammalian chromosome evolution.

Regulation of the Mycobacterium tuberculosis hypoxic response gene encoding α-crystallin

Unlike many pathogens that are overtly toxic to their hosts, the primary virulence determinant of Mycobacterium tuberculosis appears to be its ability to persist for years or decades within humans in a clinically latent state. Since early in the 20th century latency has been linked to hypoxic conditions within the host, but the response of M. tuberculosis to a hypoxic signal remains poorly characterized. The M. tuberculosis α-crystallin (acr) gene is powerfully and rapidly induced at reduced oxygen tensions, providing us with a means to identify regulators of the hypoxic response. Using a whole genome microarray, we identified >100 genes whose expression is rapidly altered by defined hypoxic conditions. Numerous genes involved in biosynthesis and aerobic metabolism are repressed, whereas a high proportion of the induced genes have no known function. Among the induced genes is an apparent operon that includes the putative two-component response regulator pair Rv3133c/Rv3132c. When we interrupted expression of this operon by targeted disruption of the upstream gene Rv3134c, the hypoxic regulation of acr was eliminated. These results suggest a possible role for Rv3132c/3133c/3134c in mycobacterial latency.

Use of cDNA subtraction and RNA interference screens in combination reveals genes required for germ-line development in Caenorhabditis elegans

Caenorhabditis elegans is an ideal organism for the study of the molecular basis of fundamental biological processes such as germ-line development, especially because of availability of the whole genome sequence and applicability of the RNA interference (RNAi) technique. To identify genes involved in germ-line development, we produced subtracted cDNA pools either enriched for or deprived of the cDNAs from germ-line tissues. We then performed differential hybridization on the high-density cDNA grid, on which about 7,600 nonoverlapping expressed sequence tag (EST) clones were spotted, to identify a set of genes specifically expressed in the germ line. One hundred and sixty-eight clones were then tested with the RNAi technique. Of these, 15 clones showed sterility with a variety of defects in germ-line development. Seven of them led to the production of unfertilized eggs, because of defects in spermatogenesis (4 clones), or defects in the oocytes (3 clones). The other 8 clones led to failure of oogenesis. These failures were caused by germ-line proliferation defect (Glp phenotype), meiotic arrest, and defects in sperm–oocyte switch (Mog phenotype) among others. These results demonstrate the efficacy of the screening strategy using the EST library combined with the RNAi technique in C. elegans.

DNA repeats identify novel virulence genes in Haemophilus influenzae.

The whole genome sequence (1.83 Mbp) of Haemophilus influenzae strain Rd was searched to identify tandem oligonucleotide repeat sequences. Loss or gain of one or more nucleotide repeats through a recombination-independent slippage mechanism is known to mediate phase variation of surface molecules of pathogenic bacteria, including H. influenzae. This facilitates evasion of host defenses and adaptation to the varying microenvironments of the host. We reasoned that iterative nucleotides could identify novel genes relevant to microbe-host interactions. Our search of the Rd genome sequence identified 9 novel loci with multiple (range 6-36, mean 22) tandem tetranucleotide repeats. All were found to be located within putative open reading frames and included homologues of hemoglobin-binding proteins of Neisseria, a glycosyltransferase (IgtC gene product) of Neisseria, and an adhesin of Yersinia. These tetranucleotide repeat sequences were also shown to be present in two other epidemiologically different H. influenzae type b strains, although the number and distribution of repeats was different. Further characterization of the IgtC gene showed that it was involved in phenotypic switching of a lipopolysaccharide epitope and that this variable expression was associated with changes in the number of tetranucleotide repeats. Mutation of IgtC resulted in attenuated virulence of H. influenzae in an infant rat model of invasive infection. These data indicate the rapidity, economy, and completeness with which whole genome sequences can be used to investigate the biology of pathogenic bacteria.

Repair of x-ray-induced DNA double-strand breaks in specific Not I restriction fragments in human fibroblasts: joining of correct and incorrect ends.

An assay that allows measurement of absolute induction frequencies for DNA double-strand breaks (dsbs) in defined regions of the genome and that quantitates rejoining of correct DNA ends has been used to study repair of dsbs in normal human fibroblasts after x-irradiation. The approach involves hybridization of single-copy DNA probes to Not I restriction fragments separated according to size by pulsed-field gel electrophoresis. Induction of dsbs is quantitated from the decrease in the intensity of the hybridizing restriction fragment and an accumulation of a smear below the band. Rejoining of dsbs results in reconstitution of the intact restriction fragment only if correct DNA ends are joined. By comparing results from this technique with results from a conventional electrophoresis assay that detects all rejoining events, it is possible to quantitate the misrejoining frequency. Three Not I fragments on the long arm of chromosome 21 were investigated with regard to dsb induction, yielding an identical induction rate of 5.8 X 10(-3) break per megabase pair per Gy. Correct dsb rejoining was measured for two of these Not I fragments after initial doses of 80 and 160 Gy. The misrejoining frequency was about 25% for both fragments and was independent of dose. This result appears to be representative for the whole genome as shown by analysis of the entire Not I fragment distribution. The correct rejoining events primarily occurred within the first 2 h, while the misrejoining kinetics included a much slower component, with about half of the events occurring between 2 and 24 h. These misrejoining kinetics are similar to those previously reported for production of exchange aberrations in interphase chromosomes.

Assessment of fight outcome is needed to activate socially driven transcriptional changes in the zebrafish brain

Group living animals must be able to express different behavior profiles depending on their social status. Therefore, the same genotype may translate into different behavioral phenotypes through socially driven differential gene expression. However, how social information is translated into a neurogenomic response and what are the specific cues in a social interaction that signal a change in social status are questions that have remained unanswered. Here, we show for the first time, to our knowledge, that the switch between status-specific neurogenomic states relies on the assessment of fight outcome rather than just on self- or opponent-only assessment of fighting ability. For this purpose, we manipulated the perception of fight outcome in male zebrafish and measured its impact on the brain transcriptome using a zebrafish whole genome gene chip. Males fought either a real opponent, and a winner and a loser were identified, or their own image on a mirror, in which case, despite expressing aggressive behavior, males did not experience either a victory or a defeat. Massive changes in the brain transcriptome were observed in real opponent fighters, with losers displaying both a higher number of differentially expressed genes and of coexpressed gene modules than winners. In contrast, mirror fighters expressed a neurogenomic state similar to that of noninteracting fish. The genes that responded to fight outcome included immediate early genes and genes involved in neuroplasticity and epigenetic modifications. These results indicate that, even in cognitively simple organisms such as zebrafish, neurogenomic responses underlying changes in social status rely on mutual assessment of fighting ability.

Initial characterization of stiff skin-like syndrome in West Highland white terriers

BACKGROUND Stiff skin syndrome and systemic or localized scleroderma are cutaneous disorders characterized by dermal fibrosis and present clinically with induration of the skin, with or without joint, internal organ or vascular involvement. OBJECTIVES To provide clinical, histological and preliminary genetic analysis of two West Highland white terrier siblings presenting with indurated skin resembling stiff skin syndrome in humans. ANIMALS Two client owned full sibling West Highland white terriers from two different litters. METHODS Clinical examination, histopathological examination and whole genome sequencing analysis of affected and unaffected West Highland white terriers. RESULTS Affected dogs exhibited markedly indurated skin that was attached firmly to the underlying tissue and incomplete closure of the mouth and eyes. No abnormalities were found by neurological or orthopaedic examination, radiographs of the head or whole body computed tomography. Histologically, the dermis and pannicular septa were thickened by a marked increase in coarse collagen fibres and a mild to moderate increase in collagen fibre diameter. The syndrome most likely follows an autosomal recessive mode of inheritance. The sequence analysis did not reveal any obvious causative variant in the investigated candidate genes ADAMTSL2 and FBN1. CONCLUSION AND CLINICAL IMPORTANCE The clinical phenotype and histopathological features of two West Highland white terrier siblings resembled stiff skin syndrome in humans. Unlike in humans, or previously described beagles with stiff skin, there was no restriction of joint mobility. Genetic analysis did not detect a candidate causative variant and warrants further research.

Insight into the evolution of the Solanaceae from the parental genomes of Petunia hybrida

Petunia hybrida is a popular bedding plant that has a long history as a genetic model system. We report the whole-genome sequencing and assembly of inbred derivatives of its two wild parents, P. axillaris N and P. inflata S6. The assemblies include 91.3% and 90.2% coverage of their diploid genomes (1.4 Gb; 2n = 14) containing 32,928 and 36,697 protein-coding genes, respectively. The genomes reveal that the Petunia lineage has experienced at least two rounds of hexaploidization: the older gamma event, which is shared with most Eudicots, and a more recent Solanaceae event that is shared with tomato and other solanaceous species. Transcription factors involved in the shift from bee to moth pollination reside in particularly dynamic regions of the genome, which may have been key to the remarkable diversity of floral colour patterns and pollination systems. The high-quality genome sequences will enhance the value of Petunia as a model system for research on unique biological phenomena such as small RNAs, symbiosis, self-incompatibility and circadian rhythms.