An ordered comparative map of the cattle and human genomes

Triamcinolone acetonide has minimal effect on short- and long-term metabolic activities of cartilage

Intra-articular corticosteroid injections, such as triamcinolone acetonide (TA), are commonly used by clinicians to manage joint synovial inflammation. However, due to conflicting evidence in literature, there is a fear among clinicians that the injections may be harmful to otherwise healthy cartilage in young patients. The purpose of this study was to evaluate the effects of TA on young, healthy chondrocytes. Articular cartilage samples were harvested from bovine knee joints (1-2 months old). In both healthy and inflammatory (interleukin-1β) challenged cartilage, samples were treated with TA at doses ranging from 1 nM to 200 μM. Following a short- (2 days) or long-term (10-14 days) treatment, chondrocyte viability, proliferation, and extracellular matrix (ECM) synthesis and degradation were evaluated with a click chemistry-based technique. Chondrocyte viability, proliferation, and anabolic activity were all minimally affected by short-term and long-term TA treatment. After both acute and sustained inflammatory challenges, TA reduced the catabolic activities in cartilage, reducing nascent glycosaminoglycan loss and maintaining cartilage mechanical properties. Overall, at physiologically relevant doses, TA had minimal negative impact on chondrocytes when maintained within their native ECM. Clinical significance: The findings provide new insight for current clinical practices concerning the use of TA in intra-articular injections, especially in young patients, and established a foundation for future investigations into the impact of corticosteroids on joint homeostasis.

The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomic resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomic resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, anti-predator strategies, and resilience and adaptive responses. They also serve as essential models for studying broad genomic traits, such as evolutionary genome expansions and contractions, as they exhibit the widest range of genome sizes among all animal taxa and possess multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The emergence of long-read sequencing technologies, combined with advanced molecular and computational techniques that improve scaffolding and reduce computational workloads, is now making it possible to address some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC, https://mvs.unimelb.edu.au/amphibian-genomics-consortium ) in early 2023. This burgeoning community already has more than 282 members from 41 countries. The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and call on the research and conservation communities to unite as part of the AGC to enable amphibian genomics research to "leap" to the next level.

Association of STAT1 gene with milk fat and protein yield in Holstein Friesian crossbred cattle maintained in the sub-tropical climate of India

Signal transducers and activators of transcription (STAT) genes are involved in signal mediation of various hormones and cytokines. STAT1 located on chromosome number 2 is involved in mammary gland development and is associated with milk composition traits in bovines. This study aimed to find any relationship and impact of STAT1/BspHI gene with milk fat and protein yields in a herd of Holstein Friesian (HF) crossbred cattle of sub-tropical climate of Northern India. Milk composition data of 535 adult HF crossbred cows for a period of 12 years was collected from the records maintained at Livestock Farm, Guru Angad Dev Veterinary and Animal Sciences University. First lactation data of 222 animals was chosen for further analysis. After data correction for non-genetic factors (season of calving, period of calving, interaction effect of season and period of calving and age at first calving) these animals were categorised into two groups based on corrected high and low milk fat and protein yields. Forty animals were then selected for blood collection and further laboratory analysis. Amplified using PCR-RFLP technique, the 314 bp STAT1 gene was digested using BspHI restriction enzyme. C-T polymorphism at nucleotide position 201 and 260 of the STAT1 amplicon was observed. At 201, for genotype AA and Aa, the genotypic frequencies were 0.80 and 0.20%. At 260, for genotype BB and Bb, the genotypic frequencies were 0.25 and 0.75%. Least square analysis showed a significant association of all genotypes with milk fat and protein yields. Hence, STAT1 can be used as a potential candidate gene to aid in better animal selection in breeding programmes.

Prevalence of Autosomal Monosomy and Trisomy Estimated Using Single Nucleotide Polymorphism Genotype Intensity Chip Information in a Large Population of Juvenile Dairy and Beef Cattle

Aneuploidy, a genetic condition characterised by the deletion (monosomy) or duplication (trisomy) of a chromosome, has been extensively studied in humans, particularly in the context of trisomy on chromosome 21, also known as Down syndrome. Research on autosomal aneuploidy in live-born cattle has been limited to case reports, resulting in a lack of prevalence estimates of aneuploidy in cattle. Furthermore, the viability or lethality of aneuploidy on specific autosomes in cattle has not been well documented. The objective of this study was to estimate the prevalence of autosomal aneuploidy in a large population of new-born and juvenile beef and dairy cattle using single nucleotide polymorphism (SNP) chip genotype intensity data. Of the population of 779,138 cattle genotyped when younger than 15 months of age, 139 cattle (i.e., 0.017%) were diagnosed with one case of autosomal trisomy. Trisomy in only 10 different autosomes were detected (BTA 4, 6, 12, 15, 20, 24, 26, 27, 28 and 29) albeit the one case of trisomy detected on Bos taurus autosome (BTA) 4 was in an additional population of 341,927 cattle that were genotyped at > 15 months of age and was therefore excluded from prevalence estimates to minimise bias. The prevalence of trisomy per chromosome was generally inversely related to the length of the chromosome. Although the number of affected individuals was few, there was no evidence of differences in prevalence by breed, inbreeding level or parental age. The parental origin of the detected cases of trisomy was maternal for 92% of the cases. No cases of monosomy were detected despite the large dataset, which included calves genotyped at birth, indicating the potential lethal nature of monosomy in cattle. Cytogenetic testing was used to verify three of the animals with detected autosomal trisomy who were still alive. Eighteen of the 139 animals identified with autosomal trisomy were recorded as being stillborn, resulting in a prevalence of autosomal aneuploidy in live-born cattle of 0.015%. Of the 121 live-born cattle with autosomal trisomy, a total of 68 died on farm at, on average (standard deviation), 6.8 (8.7) months of age. All animals with autosomal trisomy on BTA 6, 12, 15, 20 or 24 were either stillborn or died on farm within 15 days of birth. This study is the first report of trisomy on BTA 4, 6, 15, 20 and 27 in live-born cattle, as well as the first to document fertile cows with trisomy on BTA 4, 27 or 28. Given that genotype intensity SNP data from SNP-chips are readily available, identifying animals affected with autosomal aneuploidy as well as quantifying and monitoring the incidence can be easily undertaken.

The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomic resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomic resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, antipredator strategies, and resilience and adaptive responses. They also serve as essential models for studying broad genomic traits, such as evolutionary genome expansions and contractions, as they exhibit the widest range of genome sizes among all animal taxa and possess multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The emergence of long-read sequencing technologies, combined with advanced molecular and computational techniques that improve scaffolding and reduce computational workloads, is now making it possible to address some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC, https://mvs.unimelb.edu.au/amphibian-genomics-consortium) in early 2023. This burgeoning community already has more than 282 members from 41 countries. The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and call on the research and conservation communities to unite as part of the AGC to enable amphibian genomics research to "leap" to the next level.

Comparative genomics reveals the dynamics of chromosome evolution in Lepidoptera

Chromosomes are a central unit of genome organization. One-tenth of all described species on Earth are butterflies and moths, the Lepidoptera, which generally possess 31 chromosomes. However, some species display dramatic variation in chromosome number. Here we analyse 210 chromosomally complete lepidopteran genomes and show that the chromosomes of extant lepidopterans are derived from 32 ancestral linkage groups, which we term Merian elements. Merian elements have remained largely intact through 250 million years of evolution and diversification. Against this stable background, eight lineages have undergone extensive reorganization either through numerous fissions or a combination of fusion and fission events. Outside these lineages, fusions are rare and fissions are rarer still. Fusions often involve small, repeat-rich Merian elements and the sex-linked element. Our results reveal the constraints on genome architecture in Lepidoptera and provide a deeper understanding of chromosomal rearrangements in eukaryotic genome evolution.

Two novel SNPs identified in STAT1 gene adjoining a QTL for milk production in Holstein Friesian crossbreds of Kerala

We analyzed the effect of a single nucleotide polymorphism, g. C3141T in the 3' UTR of Signal transducer and activator of transcription-1 gene (STAT1) on milk production traits in the Holstein Friesian crossbred cattle of Kerala (n = 144) by association analysis and expression study. The population was genotyped by restriction fragment length polymorphism using Pag1. Association study using the General Linear Model-Analysis of Variance revealed that none of the yield or composition traits analyzed were significantly differed. The expression profile of STAT1 gene in leucocytes of animals bearing homozygous genotypes was compared by quantitative real time PCR using SYBR green chemistry with and relative expression was not found to be significantly differed. The second stage of the study, the STAT1 mRNA spanning 3213 bp was amplified from leucocytes and sequenced (GenBank: MT459802.1). Two novel SNPs were identified; one synonymous mutation in the coding region (g.A1212G) and the other in the 3'UTR (g.T3042C). The novel SNPs might contribute to STAT1 gene regulation mediated by alternate spicing or binding sites for regulatory molecules. The results reiterate the importance of extensive studies of STAT1 gene variants to substantiate the presence of a quantitative trait loci for dairy traits in the vicinity of STAT1 gene.

Runs of homocigosity and its association with productive traits in Mexican Holstein cattle

The objective of this study was to describe the runs of homozygosity (ROH) detected in the Mexican Holstein population and to associate them with milk, fat and protein yields, and conformation final score. After imputation and genomic quality control, 4,227 genotyped animals with 100,806 SNPs markers each were used. ROH with a minimum length of 1 Mb and a minimum of 10 SNPs were included in the analysis. One heterozygous SNP marker and five missing genotypes per ROH were allowed. A total of 425,098 ROH were found in the studied population (71.83 ± 10.73 ROH per animal), with an average length and coverage of 4.80 ± 0.77 Mb, and 276.89 Mb, respectively. The average chromosome length covered by ROH was 10.40 ± 3.70 Mb. ROH between 1 and 2 Mb were the most frequent in the population (51.33%) while those between 14 and 16 Mb were the least frequent (1.20%). Long chromosomes showed a larger number of ROH. Chromosomes 10 and 20, had a greater percentage of their length covered by ROH because they presented a largest number of long ROH (>8 Mb). From the total ROH, 17 were detected in 1,847 animals and distributed among different chromosomes, and were associated with milk, fat and protein yield and percentage, and conformation final score. Of the ROH with effects on production traits, the majority were found with a length between 1 and 4 Mb. These results show evidence of genomic regions preserved by genetic selection and associated with the improvement of the productivity and functionality of dairy cattle.

A telomere-to-telomere assembly of Oscheius tipulae and the evolution of rhabditid nematode chromosomes

Eukaryotic chromosomes have phylogenetic persistence. In many taxa, each chromosome has a single functional centromere with essential roles in spindle attachment and segregation. Fusion and fission can generate chromosomes with no or multiple centromeres, leading to genome instability. Groups with holocentric chromosomes (where centromeric function is distributed along each chromosome) might be expected to show karyotypic instability. This is generally not the case, and in Caenorhabditis elegans, it has been proposed that the role of maintenance of a stable karyotype has been transferred to the meiotic pairing centers, which are found at one end of each chromosome. Here, we explore the phylogenetic stability of nematode chromosomes using a new telomere-to-telomere assembly of the rhabditine nematode Oscheius tipulae generated from nanopore long reads. The 60-Mb O. tipulae genome is resolved into six chromosomal molecules. We find the evidence of specific chromatin diminution at all telomeres. Comparing this chromosomal O. tipulae assembly with chromosomal assemblies of diverse rhabditid nematodes, we identify seven ancestral chromosomal elements (Nigon elements) and present a model for the evolution of nematode chromosomes through rearrangement and fusion of these elements. We identify frequent fusion events involving NigonX, the element associated with the rhabditid X chromosome, and thus sex chromosome-associated gene sets differ markedly between species. Despite the karyotypic stability, gene order within chromosomes defined by Nigon elements is not conserved. Our model for nematode chromosome evolution provides a platform for investigation of the tensions between local genome rearrangement and karyotypic evolution in generating extant genome architectures.

An Evaluation of the Fasciola hepatica miRnome Predicts a Targeted Regulation of Mammalian Innate Immune Responses

Understanding mechanisms by which parasitic worms (helminths) control their hosts’ immune responses is critical to the development of effective new disease interventions. Fasciola hepatica, a global scourge of humans and their livestock, suppresses host innate immune responses within hours of infection, ensuring that host protective responses are quickly incapacitated. This allows the parasite to freely migrate from the intestine, through the liver to ultimately reside in the bile duct, where the parasite establishes a chronic infection that is largely tolerated by the host. The recent identification of micro(mi)RNA, small RNAs that regulate gene expression, within the extracellular vesicles secreted by helminths suggest that these non-coding RNAs may have a role in the parasite-host interplay. To date, 77 miRNAs have been identified in F. hepatica comprising primarily of ancient conserved species of miRNAs. We hypothesized that many of these miRNAs are utilized by the parasite to regulate host immune signaling pathways. To test this theory, we first compiled all of the known published F. hepatica miRNAs and critically curated their sequences and annotations. Then with a focus on the miRNAs expressed by the juvenile worms, we predicted gene targets within human innate immune cells. This approach revealed the existence of targets within every immune cell, providing evidence for the universal management of host immunology by this parasite. Notably, there was a high degree of redundancy in the potential for the parasite to regulate the activation of dendritic cells, eosinophils and neutrophils, with multiple miRNAs predicted to act on singular gene targets within these cells. This original exploration of the Fasciola miRnome offers the first molecular insight into mechanisms by which F. hepatica can regulate the host protective immune response.

Vertebrate Chromosome Evolution

The study of chromosome evolution is undergoing a resurgence of interest owing to advances in DNA sequencing technology that facilitate the production of chromosome-scale whole-genome assemblies de novo. This review focuses on the history, methods, discoveries, and current challenges facing the field, with an emphasis on vertebrate genomes. A detailed examination of the literature on the biology of chromosome rearrangements is presented, specifically the relationship between chromosome rearrangements and phenotypic evolution, adaptation, and speciation. A critical review of the methods for identifying, characterizing, and visualizing chromosome rearrangements and computationally reconstructing ancestral karyotypes is presented. We conclude by looking to the future, identifying the enormous technical and scientific challenges presented by the accumulation of hundreds and eventually thousands of chromosome-scale assemblies.

Potential candidate gene markers for milk fat in bovines: A review

In dairy animals, the principal goal of selection is to improve quality and quantity of milk. Genetic information inferred from single nucleotide polymorphism (SNP) primarily linked to Quantitative Trait Loci (QTL) can be used to improve selection for milk and milk constituent traits in bovines. Selection for a marker allele known to be associated with a beneficial QTL increases the frequency of that allele and hence, dairy performance can be enhanced. One of the potential benefit of selection based on molecular marker is that the marker genotypes can be determined in a dairy animal just after birth. Thus, marker information can be used to predict an animal's genotype before its actual performance recording for a trait is available, which considerably reduces generation interval and thus improves genetic gain in a herd for milk and its constituent traits. This review article is an attempt to comprehend the idea behind marker based selection for milk fat and genes regulating milk fat with significant effects that can be targeted specifically in selection of superior dairy animals. Once an association is established, itcan be utilized in a marker assisted breeding program for improvement of bovines.

Genetic Basis of Blood-Based Traits and Their Relationship With Performance and Environment in Beef Cattle at Weaning

The objectives of this study were to explore the usefulness of blood-based traits as indicators of health and performance in beef cattle at weaning and identify the genetic basis underlying the different blood parameters obtained from complete blood counts (CBCs). Disease costs represent one of the main factors determining profitability in animal production. Previous research has observed associations between blood cell counts and an animal’s health status in some species. CBC were recorded from approximately 570 Angus based, crossbred beef calves at weaning born between 2015 and 2016 and raised on toxic or novel tall fescue. The calves (N = ∼600) were genotyped at a density of 50k SNPs and the genotypes (N = 1160) were imputed to a density of 270k SNPs. Genetic parameters were estimated for 15 blood and 4 production. Finally, with the objective of identifying the genetic basis underlying the different blood-based traits, genome-wide association studies (GWAS) were performed for all traits. Heritability estimates ranged from 0.11 to 0.60, and generally weak phenotypic correlations and strong genetic correlations were observed among blood-based traits only. Genome-wide association study identified ninety-one 1-Mb windows that accounted for 0.5% or more of the estimated genetic variance for at least 1 trait with 21 windows overlapping across two or more traits (explaining more than 0.5% of estimated genetic variance for two or more traits). Five candidate genes have been identified in the most interesting overlapping regions related to blood-based traits. Overall, this study represents one of the first efforts represented in scientific literature to identify the genetic basis of blood cell traits in beef cattle. The results presented in this study allow us to conclude that: (1) blood-based traits have weak phenotypic correlations but strong genetic correlations among themselves. (2) Blood-based traits have moderate to high heritability. (3) There is evidence of an important overlap of genetic control among similar blood-based traits which will allow for their use in improvement programs in beef cattle.

Identification of small and large genomic candidate variants in bovine pulmonary hypoplasia and anasarca syndrome

The pulmonary hypoplasia and anasarca syndrome (PHA) is a congenital lethal disorder, which until now has been reported in cattle and sheep. PHA is characterized by extensive subcutaneous fetal edema combined with hypoplasia or aplasia of the lungs and dysplasia of the lymphatic system. PHA is assumed to be of genetic etiology. This study presents the occurrence of PHA in two different cattle breeds and their genetic causation. Two PHA cases from one sire were observed in Slovenian Cika cattle. Under the assumption of monogenic inheritance, genome-wide homozygosity mapping scaled down the critical regions to 3% of the bovine genome including a 43.6 Mb-sized segment on chromosome 6. Whole-genome sequencing of one case, variant filtering against controls and genotyping of a larger cohort of Cika cattle led to the detection of a likely pathogenic protein-changing variant perfectly associated with the disease: a missense variant on chromosome 6 in ADAMTS3 (NM_001192797.1: c.1222C>T), which affects an evolutionary conserved residue (NP_001179726.1: p.(His408Tyr)). A single PHA case was found in Danish Holstein cattle and was whole-genome sequenced along with its parents. However, as there was no plausible private protein-changing variant, mining for structural variation revealed a likely pathogenic trisomy of the entire chromosome 20. The identified ADAMTS3 associated missense variant and the trisomy 20 are two different genetic causes, which shows a compelling genetic heterogeneity for bovine PHA.

Assessment of genome stability in various breeds of cattle

Chromosomal instability is a type of genome instability involving changes in genetic information at the chromosomal level. The basic tests used to identify this form of instability are sister chromatid exchange (SCE) tests and identification of fragile sites (FS). SCE is the process by which sister chromatids become fragmented as a result of DNA strand breakage and reassembly, followed by exchange of these fragments. FS can be observed in the form of breaks, gaps or constrictions on chromosomes, which often result from multiple nucleotide repeats in DNA that are difficult to replicate. The research material was the peripheral blood of ten breeds of cattle raised in Poland, including four native breeds covered by a genetic resources conservation programme, i.e. Polish Red, Polish Red-and-White, White-Backed, and Polish Black-and-White, as well as Polish Holstein-Friesian, Simmental, Montbéliarde, Jersey, Limousine and Danish Red. Two tests were performed on chromosomes obtained from in vitro cultures: SCE and FS. The average frequency of SCE was 5.08 ± 1.31, while the incidence of FS was 3.45 ± 0.94. Differences in the incidence of SCE and FS were observed between breeds. The least damage was observed in the Polish Red and White-Backed breeds, and the most in Polish Holstein-Friesians. The most damage was observed in the interstitial part of the chromosomes. Age was shown to significantly affect the incidence of SCE and FS. Younger cows showed less damage than older ones (SCE: 4.84 ± 1.25; 5.34 ± 1.24; FS: 3.10 ± 0.88, 3.80 ± 0.92).

Comparative analyses of sperm DNA methylomes among human, mouse and cattle provide insights into epigenomic evolution and complex traits

Sperm DNA methylation is crucial for fertility and viability of offspring but epigenome evolution in mammals is largely understudied. By comparing sperm DNA methylomes and large-scale genome-wide association study (GWAS) signals between human and cattle, we aimed to examine the DNA methylome evolution and its associations with complex phenotypes in mammals. Our analysis revealed that genes with conserved non-methylated promoters (e.g., ANKS1A and WNT7A) among human and cattle were involved in common system and embryo development, and enriched for GWAS signals of body conformation traits in both species, while genes with conserved hypermethylated promoters (e.g., TCAP and CD80) were engaged in immune responses and highlighted by immune-related traits. On the other hand, genes with human-specific hypomethylated promoters (e.g., FOXP2 and HYDIN) were engaged in neuron system development and enriched for GWAS signals of brain-related traits, while genes with cattle-specific hypomethylated promoters (e.g., LDHB and DGAT2) mainly participated in lipid storage and metabolism. We validated our findings using sperm-retained nucleosome, preimplantation transcriptome, and adult tissue transcriptome data, as well as sequence evolutionary features, including motif binding sites, mutation rates, recombination rates and evolution signatures. In conclusion, our results demonstrate important roles of epigenome evolution in shaping the genetic architecture underlying complex phenotypes, hence enhance signal prioritization in GWAS and provide valuable information for human neurological disorders and livestock genetic improvement.

Mapping Genes Is Good for You

I abandoned my original career choice of high school teaching to pursue dentistry and soon abandoned that path for genetics. The latter decision was due to a challenge by a professor that led to me reading Nobel speeches by pioneer geneticists before I had formal exposure to the subject. Even then, I was 15 years into my career before my interest in rodent genomes gave way to mapping cattle genes. Events behind these twists and turns in my career path comprise the first part of this review. The remainder is a review of the development of the field of bovine genomics from my personal perspective. I have had the pleasure of working with outstanding graduate students, postdocs, and colleagues to contribute my small part to a discipline that has evolved from a few individuals mapping an orphan genome to a discipline underlying a revolution in animal breeding.

Developmental Transformation and Reduction of Connective Cavities within the Subchondral Bone

It is widely accepted that the subchondral bone (SCB) plays a crucial role in the physiopathology of osteoarthritis (OA), although its contribution is still debated. Much of the pre-clinical research on the role of SCB is concentrated on comparative evaluations of healthy vs. early OA or early OA vs. advanced OA cases, while neglecting how pure maturation could change the SCB’s microstructure. To assess the transformations of the healthy SCB from young age to early adulthood, we examined the microstructure and material composition of the medial condyle of the femur in calves (three months) and cattle (18 months) for the calcified cartilage (CC) and the subchondral bone plate (SCBP). The entire subchondral zone (SCZ) was significantly thicker in cattle compared to calves, although the proportion of the CC and SCBP thicknesses were relatively constant. The trabecular number (Tb.N.) and the connectivity density (Conn.D) were significantly higher in the deeper region of the SCZ, while the bone volume fraction (BV/TV), and the degree of anisotropy (DA) were more affected by age rather than the region. The mineralization increased within the first 250 µm of the SCZ irrespective of sample type, and became stable thereafter. Cattle exhibited higher mineralization than calves at all depths, with a mean Ca/P ratio of 1.59 and 1.64 for calves and cattle, respectively. Collectively, these results indicate that the SCZ is highly dynamic at early age, and CC is the most dynamic layer of the SCZ.

Characterization of recombination features and the genetic basis in multiple cattle breeds

BACKGROUND

Crossover generated by meiotic recombination is a fundamental event that facilitates meiosis and sexual reproduction. Comparative studies have shown wide variation in recombination rate among species, but the characterization of recombination features between cattle breeds has not yet been performed. Cattle populations in North America count millions, and the dairy industry has genotyped millions of individuals with pedigree information that provide a unique opportunity to study breed-level variations in recombination.

RESULTS

Based on large pedigrees of Jersey, Ayrshire and Brown Swiss cattle with genotype data, we identified over 3.4 million maternal and paternal crossover events from 161,309 three-generation families. We constructed six breed- and sex-specific genome-wide recombination maps using 58,982 autosomal SNPs for two sexes in the three dairy cattle breeds. A comparative analysis of the six recombination maps revealed similar global recombination patterns between cattle breeds but with significant differences between sexes. We confirmed that male recombination map is 10% longer than the female map in all three cattle breeds, consistent with previously reported results in Holstein cattle. When comparing recombination hotspot regions between cattle breeds, we found that 30% and 10% of the hotspots were shared between breeds in males and females, respectively, with each breed exhibiting some breed-specific hotspots. Finally, our multiple-breed GWAS found that SNPs in eight loci affected recombination rate and that the PRDM9 gene associated with hotspot usage in multiple cattle breeds, indicating a shared genetic basis for recombination across dairy cattle breeds.

CONCLUSIONS

Collectively, our results generated breed- and sex-specific recombination maps for multiple cattle breeds, provided a comprehensive characterization and comparison of recombination patterns between breeds, and expanded our understanding of the breed-level variations in recombination features within an important livestock species.

A High-Density Linkage Map Reveals Sexual Dimorphism in Recombination Landscapes in Red Deer (Cervus elaphus)

High-density linkage maps are an important tool to gain insight into the genetic architecture of traits of evolutionary and economic interest, and provide a resource to characterize variation in recombination landscapes. Here, we used information from the cattle genome and the 50 K Cervine Illumina BeadChip to inform and refine a high-density linkage map in a wild population of red deer (Cervus elaphus). We constructed a predicted linkage map of 38,038 SNPs and a skeleton map of 10,835 SNPs across 34 linkage groups. We identified several chromosomal rearrangements in the deer lineage relative to sheep and cattle, including six chromosome fissions, one fusion, and two large inversions. Otherwise, our findings showed strong concordance with map orders in the cattle genome. The sex-averaged linkage map length was 2739.7 cM and the genome-wide autosomal recombination rate was 1.04 cM/Mb. The female autosomal map length was 1.21 longer than that of males (2767.4 cM vs. 2280.8 cM, respectively). Sex differences in map length were driven by high female recombination rates in peri-centromeric regions, a pattern that is unusual relative to other mammal species. This effect was more pronounced in fission chromosomes that would have had to produce new centromeres. We propose two hypotheses to explain this effect: (1) that this mechanism may have evolved to counteract centromeric drive associated with meiotic asymmetry in oocyte production; and/or (2) that sequence and structural characteristics suppressing recombination in close proximity to the centromere may not have evolved at neo-centromeres. Our study provides insight into how recombination landscapes vary and evolve in mammals, and will provide a valuable resource for studies of evolution, genetic improvement, and population management in red deer and related species.

Animal breeding in the (post-) genomic era

One of the benefits of the genomics revolution for animal production will be knowledge of genes that can be used to select more profitable livestock. Although it is possible to use genetic markers linked to genes of economic importance, tests for the genes themselves will be much more successful. Consequently finding genes of economic importance to livestock will be a major research aim for the future. Most traits of economic importance are quantitative traits affected by many genes. Mutations at many genes (e.g. 500) and at many positions within a gene (e.g. 1000 coding and non-coding bases) can affect a typical quantitative trait. The effect of these mutations on phenotype is usually small (e.g. 0·1 standard deviation) but occasionally large. Many mutations are lost from the population through genetic drift and selection, so that polymorphisms exist at only a subset of the relevant genes (e.g. 100 genes). Finding these genes, that have relatively small effects, is more difficult than finding genes for a classical Mendellian trait but, as the genomic tools become more powerful, it is becoming feasible and some successes have already occurred. The standard approach is to map a quantitative trait loci (QTL) to a chromosome region using linkage and linkage disequilibrium. Then test polymorphisms in positional candidate genes for an effect on the trait. Tools such as genomic sequence, EST collections and comparative maps make this approach feasible. Candidate genes can be selected based on functional data such as gene expression obtained from microarrays. At present the gain in rate of genetic improvement from use of DNA-based tests for QTL is small, because selection without them is already quite accurate, not enough QTL have been identified and genotyping is too expensive. However, in the future, with many QTL identified and inexpensive genotyping combined with decreased generation intervals, large gains are possible.

Use of canonical discriminant analysis to study signatures of selection in cattle

BACKGROUND

Cattle include a large number of breeds that are characterized by marked phenotypic differences and thus constitute a valuable model to study genome evolution in response to processes such as selection and domestication. Detection of "signatures of selection" is a useful approach to study the evolutionary pressures experienced throughout history. In the present study, signatures of selection were investigated in five cattle breeds farmed in Italy using a multivariate approach.

METHODS

A total of 4094 bulls from five breeds with different production aptitudes (two dairy breeds: Italian Holstein and Italian Brown Swiss; two beef breeds: Piemontese and Marchigiana; and one dual purpose breed: Italian Simmental) were genotyped using the Illumina BovineSNP50 v.1 beadchip. Canonical discriminant analysis was carried out on the matrix of single nucleotide polymorphisms (SNP) genotyping data, separately for each chromosome. Scores for each canonical variable were calculated and then plotted in the canonical space to quantify the distance between breeds. SNPs for which the correlation with the canonical variable was in the 99th percentile for a specific chromosome were considered to be significantly associated with that variable. Results were compared with those obtained using an FST-based approach.

RESULTS

Based on the results of the canonical discriminant analysis, a large number of signatures of selection were detected, among which several had strong signals in genomic regions that harbour genes known to have an impact on production and morphological bovine traits, including MSTN, LCT, GHR, SCD, NCAPG, KIT, and ASIP. Moreover, new putative candidate genes were identified, such as GCK, B3GALNT1, MGAT1, GALNTL1, PRNP, and PRND. Similar results were obtained with the FST-based approach.

CONCLUSIONS

The use of canonical discriminant analysis on 50 K SNP genotypes allowed the extraction of new variables that maximize the separation between breeds. This approach is quite straightforward, it can compare more than two groups simultaneously, and relative distances between breeds can be visualized. The genes that were highlighted in the canonical discriminant analysis were in concordance with those obtained using the FST index.

Cattle Sex-Specific Recombination and Genetic Control from a Large Pedigree Analysis

Meiotic recombination is an essential biological process that generates genetic diversity and ensures proper segregation of chromosomes during meiosis. From a large USDA dairy cattle pedigree with over half a million genotyped animals, we extracted 186,927 three-generation families, identified over 8.5 million maternal and paternal recombination events, and constructed sex-specific recombination maps for 59,309 autosomal SNPs. The recombination map spans for 25.5 Morgans in males and 23.2 Morgans in females, for a total studied region of 2,516 Mb (986 kb/cM in males and 1,085 kb/cM in females). The male map is 10% longer than the female map and the sex difference is most pronounced in the subtelomeric regions. We identified 1,792 male and 1,885 female putative recombination hotspots, with 720 hotspots shared between sexes. These hotspots encompass 3% of the genome but account for 25% of the genome-wide recombination events in both sexes. During the past forty years, males showed a decreasing trend in recombination rate that coincided with the artificial selection for milk production. Sex-specific GWAS analyses identified PRDM9 and CPLX1 to have significant effects on genome-wide recombination rate in both sexes. Two novel loci, NEK9 and REC114, were associated with recombination rate in both sexes, whereas three loci, MSH4, SMC3 and CEP55, affected recombination rate in females only. Among the multiple PRDM9 paralogues on the bovine genome, our GWAS of recombination hotspot usage together with linkage analysis identified the PRDM9 paralogue on chromosome 1 to be associated in the U.S. Holstein data. Given the largest sample size ever reported for such studies, our results reveal new insights into the understanding of cattle and mammalian recombination.

A clone-free, single molecule map of the domestic cow (Bos taurus) genome

BACKGROUND

The cattle (Bos taurus) genome was originally selected for sequencing due to its economic importance and unique biology as a model organism for understanding other ruminants, or mammals. Currently, there are two cattle genome sequence assemblies (UMD3.1 and Btau4.6) from groups using dissimilar assembly algorithms, which were complemented by genetic and physical map resources. However, past comparisons between these assemblies revealed substantial differences. Consequently, such discordances have engendered ambiguities when using reference sequence data, impacting genomic studies in cattle and motivating construction of a new optical map resource--BtOM1.0--to guide comparisons and improvements to the current sequence builds. Accordingly, our comprehensive comparisons of BtOM1.0 against the UMD3.1 and Btau4.6 sequence builds tabulate large-to-immediate scale discordances requiring mediation.

RESULTS

The optical map, BtOM1.0, spanning the B. taurus genome (Hereford breed, L1 Dominette 01449) was assembled from an optical map dataset consisting of 2,973,315 (439 X; raw dataset size before assembly) single molecule optical maps (Rmaps; 1 Rmap = 1 restriction mapped DNA molecule) generated by the Optical Mapping System. The BamHI map spans 2,575.30 Mb and comprises 78 optical contigs assembled by a combination of iterative (using the reference sequence: UMD3.1) and de novo assembly techniques. BtOM1.0 is a high-resolution physical map featuring an average restriction fragment size of 8.91 Kb. Comparisons of BtOM1.0 vs. UMD3.1, or Btau4.6, revealed that Btau4.6 presented far more discordances (7,463) vs. UMD3.1 (4,754). Overall, we found that Btau4.6 presented almost double the number of discordances than UMD3.1 across most of the 6 categories of sequence vs. map discrepancies, which are: COMPLEX (misassembly), DELs (extraneous sequences), INSs (missing sequences), ITs (Inverted/Translocated sequences), ECs (extra restriction cuts) and MCs (missing restriction cuts).

CONCLUSION

Alignments of UMD3.1 and Btau4.6 to BtOM1.0 reveal discordances commensurate with previous reports, and affirm the NCBI's current designation of UMD3.1 sequence assembly as the "reference assembly" and the Btau4.6 as the "alternate assembly." The cattle genome optical map, BtOM1.0, when used as a comprehensive and largely independent guide, will greatly assist improvements to existing sequence builds, and later serve as an accurate physical scaffold for studies concerning the comparative genomics of cattle breeds.

Deficit of mitonuclear genes on the human X chromosome predates sex chromosome formation

Two taxa studied to date, the therian mammals and Caenorhabditis elegans, display underrepresentations of mitonuclear genes (mt-N genes, nuclear genes whose products are imported to and act within the mitochondria) on their X chromosomes. This pattern has been interpreted as the result of sexual conflict driving mt-N genes off of the X chromosome. However, studies in several other species have failed to detect a convergent biased distribution of sex-linked mt-N genes, leading to questions over the generality of the role of sexual conflict in shaping the distribution of mt-N genes. Here we tested whether mt-N genes moved off of the therian X chromosome following sex chromosome formation, consistent with the role of sexual conflict, or whether the paucity of mt-N genes on the therian X is a chance result of an underrepresentation on the ancestral regions that formed the X chromosome. We used a synteny-based approach to identify the ancestral regions in the platypus and chicken genomes that later formed the therian X chromosome. We then quantified the movement of mt-N genes on and off of the X chromosome and the distribution of mt-N genes on the human X and ancestral X regions. We failed to find an excess of mt-N gene movement off of the X. The bias of mt-N genes on ancestral therian X chromosomes was also not significantly different from the biases on the human X. Together our results suggest that, rather than conflict driving mt-N genes off of the mammalian X, random biases on chromosomes that formed the X chromosome could explain the paucity of mt-N genes in the therian lineage.

Molecular cloning, sequence identification, and gene expression analysis of bovine ADCY2 gene

Adenylyl cyclase 2 (ADCY2), a class B member of adenylyl cyclases, is important in accelerating phosphor-acidification as well as glycogen synthesis and breakdown. Given its distinct role in flesh tenderization after butchering, we cloned and sequenced the ADCY2 gene from Yanbian cattle and assessed its expression in bovine tissues. A 2947 bp nucleotide sequence representing the full-length cDNA of bovine ADCY2 gene was obtained by 5' and 3' remote analysis computations for gene expression. Analyses of the putative protein sequence showed that ADCY2 had high homology among species, except with the non-mammal Oreochromis niloticus. Gene structural domain analyses in humans and rats indicated that the ADCY2 protein had no flaw; only the transmembrane domain was reduced and the CYCc structure domain was shortened. Assessment of ADCY2 expression in bovine tissues by real-time PCR showed that the highest expression was in the testes, followed by the longissimus dorsi, tensor fasciae latae, and latissimus dorsi. These data will serve as a foundation for further insight into the cattle ADCY2 gene.

Use of locally weighted scatterplot smoothing (LOWESS) regression to study selection signatures in Piedmontese and Italian Brown cattle breeds

Selection is the major force affecting local levels of genetic variation in species. The availability of dense marker maps offers new opportunities for a detailed understanding of genetic diversity distribution across the animal genome. Over the last 50 years, cattle breeds have been subjected to intense artificial selection. Consequently, regions controlling traits of economic importance are expected to exhibit selection signatures. The fixation index (Fst ) is an estimate of population differentiation, based on genetic polymorphism data, and it is calculated using the relationship between inbreeding and heterozygosity. In the present study, locally weighted scatterplot smoothing (LOWESS) regression and a control chart approach were used to investigate selection signatures in two cattle breeds with different production aptitudes (dairy and beef). Fst was calculated for 42 514 SNP marker loci distributed across the genome in 749 Italian Brown and 364 Piedmontese bulls. The statistical significance of Fst values was assessed using a control chart. The LOWESS technique was efficient in removing noise from the raw data and was able to highlight selection signatures in chromosomes known to harbour genes affecting dairy and beef traits. Examples include the peaks detected for BTA2 in the region where the myostatin gene is located and for BTA6 in the region harbouring the ABCG2 locus. Moreover, several loci not previously reported in cattle studies were detected.

Independent specialization of the human and mouse X chromosomes for the male germ line

We compared the human and mouse X chromosomes to systematically test Ohno’s law, which states that the gene content of X chromosomes is conserved across placental mammals¹. First, we improved the accuracy of the human X-chromosome reference sequence through single-haplotype sequencing of ampliconic regions. This closed gaps in the reference sequence, corrected previously misassembled regions, and identified new palindromic amplicons. Our subsequent analysis led us to conclude that the evolution of human and mouse X chromosomes was bimodal. In accord with Ohno’s law, 94–95% of X-linked single-copy genes are shared between human and mouse; most are expressed in both sexes. Strikingly, most X-ampliconic genes are exceptions to Ohno’s law: only 31% of human and 22% of mouse X-ampliconic genes share orthologs. X-ampliconic genes are expressed predominantly in testicular germ cells, and many were independently acquired since the common ancestor of humans and mice, specializing portions of their X chromosomes for sperm production.

First steps: bovine genomics in historical perspective

The vision of Morris Soller was instrumental in launching the field of bovine genomics. This study is a review of the early years of bovine gene mapping leading up to the sequencing and assembly of the bovine genome in 2009. A historical perspective of parasexual, linkage and physical mapping is provided with a focus on the contribution of these maps to the eventual assignment and orientation of genes and sequence to cattle chromosomes.

Quantitative trait loci associated with the immune response to a bovine respiratory syncytial virus vaccine

Infectious disease is an important problem for animal breeders, farmers and governments worldwide. One approach to reducing disease is to breed for resistance. This linkage study used a Charolais-Holstein F2 cattle cross population (n = 501) which was genotyped for 165 microsatellite markers (covering all autosomes) to search for associations with phenotypes for Bovine Respiratory Syncytial Virus (BRSV) specific total-IgG, IgG1 and IgG2 concentrations at several time-points pre- and post-BRSV vaccination. Regions of the bovine genome which influenced the immune response induced by BRSV vaccination were identified, as well as regions associated with the clearance of maternally derived BRSV specific antibodies. Significant positive correlations were detected within traits across time, with negative correlations between the pre- and post-vaccination time points. The whole genome scan identified 27 Quantitative Trait Loci (QTL) on 13 autosomes. Many QTL were associated with the Thymus Helper 1 linked IgG2 response, especially at week 2 following vaccination. However the most significant QTL, which reached 5% genome-wide significance, was on BTA 17 for IgG1, also 2 weeks following vaccination. All animals had declining maternally derived BRSV specific antibodies prior to vaccination and the levels of BRSV specific antibody prior to vaccination were found to be under polygenic control with several QTL detected.

Heifers from the same population (n = 195) were subsequently immunised with a 40-mer Foot-and-Mouth Disease Virus peptide (FMDV) in a previous publication. Several of these QTL associated with the FMDV traits had overlapping peak positions with QTL in the current study, including the QTL on BTA23 which included the bovine Major Histocompatibility Complex (BoLA), and QTL on BTA9 and BTA24, suggesting that the genes underlying these QTL may control responses to multiple antigens. These results lay the groundwork for future investigations to identify the genes underlying the variation in clearance of maternal antibody and response to vaccination.

Marsupial genome sequences: providing insight into evolution and disease

Marsupials (metatherians), with their position in vertebrate phylogeny and their unique biological features, have been studied for many years by a dedicated group of researchers, but it has only been since the sequencing of the first marsupial genome that their value has been more widely recognised. We now have genome sequences for three distantly related marsupial species (the grey short-tailed opossum, the tammar wallaby, and Tasmanian devil), with the promise of many more genomes to be sequenced in the near future, making this a particularly exciting time in marsupial genomics. The emergence of a transmissible cancer, which is obliterating the Tasmanian devil population, has increased the importance of obtaining and analysing marsupial genome sequence for understanding such diseases as well as for conservation efforts. In addition, these genome sequences have facilitated studies aimed at answering questions regarding gene and genome evolution and provided insight into the evolution of epigenetic mechanisms. Here I highlight the major advances in our understanding of evolution and disease, facilitated by marsupial genome projects, and speculate on the future contributions to be made by such sequences.

RNA interference-mediated knockdown of DGAT1 decreases triglyceride content of bovine mammary epithelial cell line

Diacylglyceroltransferase-1 (DGAT1) expresses in nearly all tissues, including the mammary gland. Mice lacking DGAT1 exhibit decreased triglyceride content in mammary tissue, and are resistant to diet-induced obesity and diabetes mellitus. Thus, DGAT1 has received considerable attention. In the present study, the function of DGAT1 was examined by liposome mediated RNA interference (RNAi) to knockdown the expression of endogenous DGAT1 expression in bovine mammary epithelial cells (BMEC) and the changes of the biological functions of cells were analyzed. The mRNA and protein levels, intracellular triglyceride (TG) content, and total protein of BMECs were analyzed by real-time PCR, Western blot, TG kit, and ultraviolet spectrophotometer, respectively, before and after RNAi treatment. The results indicated that knockdown of DGAT1 expression significantly reduced TG content in BMECs. This study further confirmed the importance of DGAT1 in triglyceride synthesis in bovine mammary tissue.

Cattle genomics and its implications for future nutritional strategies for dairy cattle

The recently sequenced cattle (Bos taurus) genome unraveled the unique genomic features of the species and provided the molecular basis for applying a systemic approach to systematically link genomic information to metabolic traits. Comparative analysis has identified a variety of evolutionary adaptive features in the cattle genome, such as an expansion of the gene families related to the rumen function, large number of chromosomal rearrangements affecting regulation of genes for lactation, and chromosomal rearrangements that are associated with segmental duplications and copy number variations. Metabolic reconstruction of the cattle genome has revealed that core metabolic pathways are highly conserved among mammals although five metabolic genes are deleted or highly diverged and seven metabolic genes are present in duplicate in the cattle genome compared to their human counter parts. The evolutionary loss and gain of metabolic genes in the cattle genome may reflect metabolic adaptations of cattle. Metabolic reconstruction also provides a platform for better understanding of metabolic regulation in cattle and ruminants. A substantial body of transcriptomics data from dairy and beef cattle under different nutritional management and across different stages of growth and lactation are already available and will aid in linking the genome with metabolism and nutritional physiology of cattle. Application of cattle genomics has great potential for future development of nutritional strategies to improve efficiency and sustainability of beef and milk production. One of the biggest challenges is to integrate genomic and phenotypic data and interpret them in a biological and practical platform. Systems biology, a holistic and systemic approach, will be very useful in overcoming this challenge.

Genome Scan for Parent-of-Origin QTL Effects on Bovine Growth and Carcass Traits

Parent-of-origin effects (POE) such as genomic imprinting influence growth and body composition in livestock, rodents, and humans. Here, we report the results of a genome scan to detect quantitative trait loci (QTL) with POE on growth and carcass traits in Angus × Brahman cattle crossbreds. We identified 24 POE–QTL on 15 Bos taurus autosomes (BTAs) of which six were significant at 5% genome-wide (GW) level and 18 at the 5% chromosome-wide (CW) significance level. Six QTL were paternally expressed while 15 were maternally expressed. Three QTL influencing post-weaning growth map to the proximal end of BTA2 (linkage region of 0–9 cM; genomic region of 5.0–10.8 Mb), for which only one imprinted ortholog is known so far in the human and mouse genomes, and therefore may potentially represent a novel imprinted region. The detected QTL individually explained 1.4 ∼ 5.1% of each trait’s phenotypic variance. Comparative in silico analysis of bovine genomic locations show that 32 out of 1,442 known mammalian imprinted genes from human and mouse homologs map to the identified QTL regions. Although several of the 32 genes have been associated with quantitative traits in cattle, only two (GNAS and PEG3) have experimental proof of being imprinted in cattle. These results lend additional support to recent reports that POE on quantitative traits in mammals may be more common than previously thought, and strengthen the need to identify and experimentally validate cattle orthologs of imprinted genes so as to investigate their effects on quantitative traits.

Genetic mapping of quantitative trait loci for meat quality and muscle metabolic traits in cattle

A whole-genome scan was carried out in New Zealand and Australia to detect quantitative trait loci (QTL) for live animal and carcass composition traits and meat quality attributes in cattle. Backcross calves (385 heifers and 398 steers) were generated, with Jersey and Limousin backgrounds. The New Zealand cattle were reared and finished on pasture, whilst Australian cattle were reared on grass and finished on grain for at least 180 days. This paper reports on meat quality traits (tenderness measured as shear force at 4-5 ages on two muscles as well as associated traits of meat colour, pH and cooking loss) and a number of metabolic traits. For meat quality traits, 18 significant QTL (P < 0.05), located in nine linkage groups, were detected on a genome-wise basis, in combined-sire (seven QTL) or within-sire analyses (11 QTL). For metabolic traits, 11 significant QTL (P < 0.05), located in eight linkage groups, were detected on a genome-wise basis, in combined-sire (five QTL) or within-sire analyses (six QTL). BTA2 and BTA3 had QTL for both metabolic traits and meat quality traits. Six significant QTL for meat quality and metabolic traits were found at the proximal end of chromosome 2. BTA2 and BTA29 were the most common chromosomes harbouring QTL for meat quality traits; QTL for improved tenderness were associated with Limousin-derived and Jersey-derived alleles on these two chromosomes, respectively.

Status of the cattle genome map

During the last 30 years, the cattle genome map has been expanded from 4 genes linked on chromosome X to over 22,000 genes identified in the cattle genome sequence assembly. This progress has been achieved due to numerous projects on linkage and physical mapping of the cattle genome driven by its agricultural and scientific significance. Indeed, the high-resolution mapping and functional analysis of the genome led to the discovery of major quantitative trait loci (QTL) regions and several quantitative trait nucleotides (QTNs), as well as some disease genes in the cow population. In addition, a comparison of the cattle genome to the genomes of other mammals has revealed its unique features gained during the speciation and adaptation. With the development of non-expensive sequencing techniques, the analysis of the cattle genome will shift towards the identification of differences between breeds or individuals within breeds that account for the unique features of each breed. This approach holds promise for the development of effective tools for the marker assistant selection and disease diagnostics in cattle.

A 4 Mb high resolution BAC contig on bovine chromosome 1q12 and comparative analysis with human chromosome 21q22

The bovine RPCI-42 BAC library was screened to construct a sequence-ready ~4 Mb single contig of 92 BAC clones on BTA 1q12. The contig covers the region between the genes KRTAP8P1 and CLIC6. This genomic segment in cattle is of special interest as it contains the dominant gene responsible for the hornless or polled phenotype in cattle. The construction of the BAC contig was initiated by screening the bovine BAC library with heterologous cDNA probes derived from 12 human genes of the syntenic region on HSA 21q22. Contig building was facilitated by BAC end sequencing and chromosome walking. During the construction of the contig, 165 BAC end sequences and 109 single-copy STS markers were generated. For comparative mapping of 25 HSA 21q22 genes, genomic PCR primers were designed from bovine EST sequences and the gene-associated STSs mapped on the contig. Furthermore, bovine BAC end sequence comparisons against the human genome sequence revealed significant matches to HSA 21q22 and allowed the in silico mapping of two new genes in cattle. In total, 31 orthologues of human genes located on HSA 21q22 were directly mapped within the bovine BAC contig, of which 16 genes have been cloned and mapped for the first time in cattle. In contrast to the existing comparative bovine–human RH maps of this region, these results provide a better alignment and reveal a completely conserved gene order in this 4 Mb segment between cattle, human and mouse. The mapping of known polled linked BTA 1q12 microsatellite markers allowed the integration of the physical contig map with existing linkage maps of this region and also determined the exact order of these markers for the first time. Our physical map and transcript map may be useful for positional cloning of the putative polled gene in cattle. The nucleotide sequence data reported in this paper have been submitted to EMBL and have been assigned Accession Numbers AJ698510–AJ698674.

Cloning and expression of activation induced cytidine deaminase from Bos taurus

Activation induced cytidine deaminase is an enzyme crucial to somatic hypermutation and gene conversion, processes that are essential for the diversification of Ig V genes. The bovine Ig repertoire appears to be diversified by mechanisms that are significantly different to those that operate in humans and mice. This study set out to test the hypothesis that differences in the organization, coding sequence, expression or genomic location of the bovine AICDA gene enables the encoded enzyme to catalyse the unusual Ig diversification mechanism seen in cattle as well as conventional antigen-driven mutation. Characterization of bovine AICDA excluded the first two possibilities. AICDA expression was detected in lymphoid tissues from neonatal and older cattle, but AICDA cDNA could not be detected in muscle tissue. The pattern of gene expression did not therefore differ from that in other vertebrates. The AICDA cDNA was cloned and expressed successfully in Escherichia coli generating a phenotype consistent with the mutating action of this deaminase. Using a whole genome radiation hybrid panel, bovine AICDA was mapped to a region of bovine chromosome 5 syntenic with the location of human AICDA on chromosome 12. We conclude that the unusual nature of Ig diversification in cattle is unlikely to be attributable to the structure, sequence, activity or genomic location of bovine AICDA.

A high density linkage map of the bovine genome

BACKGROUND

Recent technological advances have made it possible to efficiently genotype large numbers of single nucleotide polymorphisms (SNPs) in livestock species, allowing the production of high-density linkage maps. Such maps can be used for quality control of other SNPs and for fine mapping of quantitative trait loci (QTL) via linkage disequilibrium (LD).

RESULTS

A high-density bovine linkage map was constructed using three types of markers. The genotypic information was obtained from 294 microsatellites, three milk protein haplotypes and 6769 SNPs. The map was constructed by combining genetic (linkage) and physical information in an iterative mapping process. Markers were mapped to 3,155 unique positions; the 6,924 autosomal markers were mapped to 3,078 unique positions and the 123 non-pseudoautosomal and 19 pseudoautosomal sex chromosome markers were mapped to 62 and 15 unique positions, respectively. The linkage map had a total length of 3,249 cM. For the autosomes the average genetic distance between adjacent markers was 0.449 cM, the genetic distance between unique map positions was 1.01 cM and the average genetic distance (cM) per Mb was 1.25.

CONCLUSION

There is a high concordance between the order of the SNPs in our linkage map and their physical positions on the most recent bovine genome sequence assembly (Btau 4.0). The linkage maps provide support for fine mapping projects and LD studies in bovine populations. Additionally, the linkage map may help to resolve positions of unassigned portions of the bovine genome.

A first generation whole genome RH map of the river buffalo with comparison to domestic cattle

Background

The recently constructed river buffalo whole-genome radiation hybrid panel (BBURH₅₀₀₀) has already been used to generate preliminary radiation hybrid (RH) maps for several chromosomes, and buffalo-bovine comparative chromosome maps have been constructed. Here, we present the first-generation whole genome RH map (WG-RH) of the river buffalo generated from cattle-derived markers. The RH maps aligned to bovine genome sequence assembly Btau_4.0, providing valuable comparative mapping information for both species.

Results

A total of 3990 markers were typed on the BBURH₅₀₀₀ panel, of which 3072 were cattle derived SNPs. The remaining 918 were classified as cattle sequence tagged site (STS), including coding genes, ESTs, and microsatellites. Average retention frequency per chromosome was 27.3% calculated with 3093 scorable markers distributed in 43 linkage groups covering all autosomes (24) and the X chromosomes at a LOD ≥ 8. The estimated total length of the WG-RH map is 36,933 cR₅₀₀₀. Fewer than 15% of the markers (472) could not be placed within any linkage group at a LOD score ≥ 8. Linkage group order for each chromosome was determined by incorporation of markers previously assigned by FISH and by alignment with the bovine genome sequence assembly (Btau_4.0).

Conclusion

We obtained radiation hybrid chromosome maps for the entire river buffalo genome based on cattle-derived markers. The alignments of our RH maps to the current bovine genome sequence assembly (Btau_4.0) indicate regions of possible rearrangements between the chromosomes of both species. The river buffalo represents an important agricultural species whose genetic improvement has lagged behind other species due to limited prior genomic characterization. We present the first-generation RH map which provides a more extensive resource for positional candidate cloning of genes associated with complex traits and also for large-scale physical mapping of the river buffalo genome.

Evolutionary-new centromeres preferentially emerge within gene deserts

BACKGROUND

Evolutionary-new centromeres (ENCs) result from the seeding of a centromere at an ectopic location along the chromosome during evolution. The novel centromere rapidly acquires the complex structure typical of eukaryote centromeres. This phenomenon has played an important role in shaping primate karyotypes. A recent study on the evolutionary-new centromere of macaque chromosome 4 (human 6) showed that the evolutionary-new centromere domain was deeply restructured, following the seeding, with respect to the corresponding human region assumed as ancestral. It was also demonstrated that the region was devoid of genes. We hypothesized that these two observations were not merely coincidental and that the absence of genes in the seeding area constituted a crucial condition for the evolutionary-new centromere fixation in the population.

RESULTS

To test our hypothesis, we characterized 14 evolutionary-new centromeres selected according to conservative criteria. Using different experimental approaches, we assessed the extent of genomic restructuring. We then determined the gene density in the ancestral domain where each evolutionary-new centromere was seeded.

CONCLUSIONS

Our study suggests that restructuring of the seeding regions is an intrinsic property of novel evolutionary centromeres that could be regarded as potentially detrimental to the normal functioning of genes embedded in the region. The absence of genes, which was found to be of high statistical significance, appeared as a unique favorable scenario permissive of evolutionary-new centromere fixation in the population.

Anti-major histocompatibility complex class I-related chain A antibodies in organ transplantation

Candidate trigger antigens for alloreactive responses have been appearing continuously in the organ transplant scenario. Major histocompatibility complex class I-related chain A (MICA) is a polymorphic gene family, located near the HLA-B locus on chromosome 6, that encode a 62-kd cell surface glycoprotein. Endothelial cells, in addition to many cell lines, express MICA, whereas resting lymphocytes do not, making this polymorphic molecule a target for both cellular and humoral immune responses. Major histocompatibility complex class I-related chain A antigens are able to elicit the synthesis of alloantibodies in transplant recipients. These antibodies have been found in association with irreversible allograft rejection, an increased frequency of acute rejection episodes, and a significantly lower deceased donor graft survival, as well as in the eluates from rejected grafts. This review summarizes currently available information on MICA in the transplant setting. Undoubtedly, the questions that have surfaced surpass in excess the currently available answers.

Molecular characterization of the bovine chromodomain Y-like genes

The human chromodomain protein, Y-like (CDYL) gene family consists of three members, one on the Y chromosome (CDY) and two on autosomes (CDYL and CDYL2). Studies in the human and mouse showed that genes in the CDYL family are abundantly expressed in testis and play an important role in spermatogenesis. In this study, we have characterized the bovine CDYL (bCDYL) and CDYL2 (bCDYL2) genes. We found that bCDYL and bCDYL2 are very similar to the human orthologues at both mRNA (79% and 85%) and protein (89% and 93%) levels. However, the similarity between the bCDYL and bCDYL2 proteins is low (41%). The bCDYL gene is composed of nine exons, and the bCDYL2 has seven exons. The bCDYL and bCDYL2 genes were mapped by radiation hybrid mapping to bovine chromosomes (BTA) 24 and 18 respectively. The bCDYL gene has four transcript variants that produce four protein isoforms. RT-PCR expression analysis in 12 bovine tissues showed that bCDYL variant 2 was expressed in the testis only, bCDYL variants 1, 3 and 4 were expressed predominantly in the testis and at very low or undetectable levels in the remaining tissues and bCDYL2 was expressed ubiquitously. Examination of bovine testis with in situ hybridization revealed that the bCDYL and bCDYL2 transcripts were found mainly in spermatids, though the amounts of transcripts varied among genes/variants. In addition, antisense transcripts were detected in bCDYL variants 2/3 and 4, as well as in the bCDYL2 gene.

Construction of a dense SNP map for bovine chromosome 6 to assist the assembly of the bovine genome sequence

A linkage map was constructed for bovine chromosome 6 (BTA6), using 399 single nucleotide polymorphisms (SNPs) detected primarily from PCR-resequencing. The efficiency of SNP detection was highly dependent on the source of sequence information chosen for primer design (BAC-end sequences, introns or promoters). The SNPs were used to build a linkage map comprising 104 cM on BTA6. The SNP order in the linkage map corresponded very well with radiation hybrid (RH) maps available for BTA6 as well as with expected positions in the human comparative map, but diverged significantly from the current assembly of the bovine genome (Btau_3.1). When performing linkage analysis with the marker order suggested from the Btau_3.1 we observed an expansion of the genetic map from 104 cM to 137 cM, strongly suggesting a reordering of scaffolds in the current version of the bovine genome assembly. The extent of LD on BTA6 was evaluated by calculating the average r(2) for SNP pairs separated by given distances. The decline of LD was rapid with distance, such that r(2) was 0.1 at 100 kb. Our results indicate that linkage mapping will be a valuable source of information for correcting errors in the current bovine assembly. These errors were sufficiently frequent to be of concern for the accuracy of mapping QTL with panels of SNPs whose positions are based on the current assembly.