Mapping quantitative trait loci (QTL) in sheep. I. A new male framework linkage map and QTL for growth rate and body weight

A male sheep linkage map comprising 191 microsatellites was generated from a single family of 510 Awassi-Merino backcross progeny. Except for ovine chromosomes 1, 2, 10 and 17, all other chromosomes yielded a LOD score difference greater than 3.0 between the best and second-best map order. The map is on average 11% longer than the Sheep Linkage Map v4.7 male-specific map. This map was employed in quantitative trait loci (QTL) analyses on body-weight and growth-rate traits between birth and 98 weeks of age. A custom maximum likelihood program was developed to map QTL in half-sib families for non-inbred strains (QTL-MLE) and is freely available on request. The new analysis package offers the advantage of enabling QTL × fixed effect interactions to be included in the model. Fifty-four putative QTL were identified on nine chromosomes. Significant QTL with sex-specific effects (i.e. QTL × sex interaction) in the range of 0.4 to 0.7 SD were found on ovine chromosomes 1, 3, 6, 11, 21, 23, 24 and 26.

Extent of genome-wide linkage disequilibrium in Australian Holstein-Friesian cattle based on a high-density SNP panel

BACKGROUND

The extent of linkage disequilibrium (LD) within a population determines the number of markers that will be required for successful association mapping and marker-assisted selection. Most studies on LD in cattle reported to date are based on microsatellite markers or small numbers of single nucleotide polymorphisms (SNPs) covering one or only a few chromosomes. This is the first comprehensive study on the extent of LD in cattle by analyzing data on 1,546 Holstein-Friesian bulls genotyped for 15,036 SNP markers covering all regions of all autosomes. Furthermore, most studies in cattle have used relatively small sample sizes and, consequently, may have had biased estimates of measures commonly used to describe LD. We examine minimum sample sizes required to estimate LD without bias and loss in accuracy. Finally, relatively little information is available on comparative LD structures including other mammalian species such as human and mouse, and we compare LD structure in cattle with public-domain data from both human and mouse.

RESULTS

We computed three LD estimates, D', Dvol and r2, for 1,566,890 syntenic SNP pairs and a sample of 365,400 non-syntenic pairs. Mean D' is 0.189 among syntenic SNPs, and 0.105 among non-syntenic SNPs; mean r2 is 0.024 among syntenic SNPs and 0.0032 among non-syntenic SNPs. All three measures of LD for syntenic pairs decline with distance; the decline is much steeper for r2 than for D' and Dvol. The value of D' and Dvol are quite similar. Significant LD in cattle extends to 40 kb (when estimated as r2) and 8.2 Mb (when estimated as D'). The mean values for LD at large physical distances are close to those for non-syntenic SNPs. Minor allelic frequency threshold affects the distribution and extent of LD. For unbiased and accurate estimates of LD across marker intervals spanning < 1 kb to > 50 Mb, minimum sample sizes of 400 (for D') and 75 (for r2) are required. The bias due to small samples sizes increases with inter-marker interval. LD in cattle is much less extensive than in a mouse population created from crossing inbred lines, and more extensive than in humans.

CONCLUSION

For association mapping in Holstein-Friesian cattle, for a given design, at least one SNP is required for each 40 kb, giving a total requirement of at least 75,000 SNPs for a low power whole-genome scan (median r2 > 0.19) and up to 300,000 markers at 10 kb intervals for a high power genome scan (median r2 > 0.62). For estimation of LD by D' and Dvol with sufficient precision, a sample size of at least 400 is required, whereas for r2 a minimum sample of 75 is adequate.

A primary assembly of a bovine haplotype block map based on a 15,036-single-nucleotide polymorphism panel genotyped in holstein-friesian cattle

Analysis of data on 1000 Holstein-Friesian bulls genotyped for 15,036 single-nucleotide polymorphisms (SNPs) has enabled genomewide identification of haplotype blocks and tag SNPs. A final subset of 9195 SNPs in Hardy-Weinberg equilibrium and mapped on autosomes on the bovine sequence assembly (release Btau 3.1) was used in this study. The average intermarker spacing was 251.8 kb. The average minor allele frequency (MAF) was 0.29 (0.05-0.5). Following recent precedents in human HapMap studies, a haplotype block was defined where 95% of combinations of SNPs within a region are in very high linkage disequilibrium. A total of 727 haplotype blocks consisting of > or =3 SNPs were identified. The average block length was 69.7 +/- 7.7 kb, which is approximately 5-10 times larger than in humans. These blocks comprised a total of 2964 SNPs and covered 50,638 kb of the sequence map, which constitutes 2.18% of the length of all autosomes. A set of tag SNPs, which will be useful for further fine-mapping studies, has been identified. Overall, the results suggest that as many as 75,000-100,000 tag SNPs would be needed to track all important haplotype blocks in the bovine genome. This would require approximately 250,000 SNPs in the discovery phase.

Genome-wide genetic diversity of Holstein Friesian cattle reveals new insights into Australian and global population variability, including impact of selection

Past breeding strategies for dairy cattle have been very effective in producing rapid genetic gain to achieve industry targets and raise profitability. Such gains have been largely facilitated by intense selection of sires combined with the use of artificial insemination. However, this practice can potentially limit the level of genetic diversity through inbreeding and selection plateaus. The rate of inbreeding in Australia is increasing, primarily as a result of semen importation from a small number of prominent bulls from the USA. The effect of this genetic influx in the Australian dairy cattle population is poorly understood both in terms of diversity and local adaptation/divergence. This study uses 845 genome-wide SNP genetic markers and 431 bulls to characterize the level of genetic diversity and genetic divergence within the Australian and international Holstein Friesian dairy population. No significant differences in genetic diversity (as measured by heterozygosity [H(o)] and allelic richness [A]) were observed over the 25-year time period (1975-1999) for bulls used in Australia. The importation of foreign semen into Australia has increased the effective population size until it was in effect a sub-sample of the global population. Our data indicate that most individuals are equally closely related to one another, regardless of country of origin and year of birth. In effect, the global population can be considered as one single population unit. These results indicate that inbreeding, genetic drift and selection has had little effect at reducing genetic diversity and differentiating the Australian Holstein Friesian population at a genome-wide level.

Widespread Utility of Highly Informative AFLP Molecular Markers across Divergent Shark Species

Population numbers of many shark species are declining rapidly around the world. Despite the commercial and conservation significance, little is known on even the most fundamental aspects of their population biology. Data collection that relies on direct observation can be logistically challenging with sharks. Consequently, molecular methods are becoming increasingly important to obtain knowledge that is critical for conservation and management. Here we describe an amplified fragment length polymorphism method that can be applied universally to sharks to identify highly informative genome-wide polymorphisms from 12 primer pairs. We demonstrate the value of our method on 15 divergent shark species within the superorder Galeomorphii, including endangered species which are notorious for low levels of genetic diversity. Both the endangered sand tiger shark (Carcharodon taurus, N = 18) and the great white shark (Carcharodon carcharias, N = 7) displayed relatively high levels of allelic diversity. A total of 59 polymorphic loci (H(e) = 0.373) and 78 polymorphic loci (H(e) = 0.316) were resolved in C. taurus and C. carcharias, respectively. Results from other sharks (e.g., Orectolobus ornatus, Orectolobus sp., and Galeocerdo cuvier) produced remarkably high numbers of polymorphic loci (106, 94, and 86, respectively) from a limited sample size of only 2. A major constraint to obtaining much needed genetic data from sharks is the time-consuming process of developing molecular markers. Here we demonstrate the general utility of a technique that provides large numbers of informative loci in sharks.

Technical note: Whole-genome amplification of DNA extracted from cattle semen samples

The bovine genome sequence project and the discovery of many thousands of bovine single nucleotide polymorphisms has opened the door for large-scale genotyping studies to identify genes that contribute to economically important traits with relevance to the beef and dairy industries. Large amounts of DNA will be required for these research projects. This study reports the use of the whole-genome amplification (WGA) method to create an unlimited supply of DNA for use in genotyping studies and long-term storage for future gene discovery projects. Two commercial WGA kits (GenomiPhi, Amersham Biosciences, Sydney, Australia, and REPLI-g, Qiagen, Doncaster, Australia) were used to amplify DNA from straws of bull semen, resulting in an average of 7.2 and 67 microg of DNA per reaction, respectively. The comparison of 3.5 kb of sequences from the amplified and unamplified DNA indicated no detectable DNA differences. Similarly, gene marker analysis conducted on genomic DNA and DNA after WGA indicated no difference in marker amplification or clarity and accuracy of scoring for approximately 10,000 single nucleotide polymorphism markers when compared with WGA samples genotyped in duplicate. These results illustrate that WGA is a suitable method for the amplification and recovery of DNA from bull semen samples for routine genomic investigations.

A first-generation metric linkage disequilibrium map of bovine chromosome 6

We constructed a metric linkage disequilibrium (LD) map of bovine chromosome 6 (BTA6) on the basis of data from 220 SNPs genotyped on 433 Australian dairy bulls. This metric LD map has distances in LD units (LDUs) that are analogous to centimorgans in linkage maps. The LD map of BTA6 has a total length of 8.9 LDUs. Within the LD map, regions of high LD (represented as blocks) and regions of low LD (steps) are observed, when plotted against the integrated map in kilobases. At the most stringent block definition, namely a set of loci with zero LDU increase over the span of these markers, BTA6 comprises 40 blocks, accounting for 41% of the chromosome. At a slightly lower stringency of block definition (a set of loci covering a maximum of 0.2 LDUs on the LD map), up to 81% of BTA6 is spanned by 46 blocks and with 13 steps that are likely to reflect recombination hot spots. The mean swept radius (the distance over which LD is likely to be useful for mapping) is 13.3 Mb, confirming extensive LD in Holstein-Friesian dairy cattle, which makes such populations ideal for whole-genome association studies.

Quantitative trait loci for steady-state platelet count in mice

Platelet count in humans is a strongly genetically regulated trait, with approximately 85% of the interindividual variance in platelet numbers attributable to genetic factors. Inbred mouse strains also have strain-specific platelet count ranges. As part of a project to identify novel factors that regulate platelet count, we identified two inbred mouse strains, CBA/CaH and QSi5, with substantial differences in platelet count (mean values of 581 vs. 1062 x 10(9)/L). An F(2) intercross resource of 1126 animals was bred from these two parental strains for a genomewide scan for quantitative trait loci (QTL) for platelet count. QTL were identified on MMU1 (LOD 6.8, p < 0.0005) and MMU11 (LOD 11.2, p < 0.0005) by selectively genotyping animals from the extremes of the F(2) platelet count distribution. Three other QTL of suggestive statistical significance were also detected on MMU7, 13, and 17. It is noteworthy that no QTL were detected in the vicinity of the genes encoding thrombopoietin ( Thpo), and its receptor ( c-Mpl), both known to influence platelet production. Comparison of gene expression levels between the parental mouse strains by microarrays also showed little difference in the mRNA levels of these known candidate genes. These results represent the first published use of a genetic linkage-based approach in a mouse model toward the identification of genetic factors that regulate platelet count.

Genetic analysis of a population crash in brush-tailed rock-wallabies (Petrogale penicillata) from Jenolan Caves, south-eastern Australia

Although the theoretical effects of a severe reduction in effective population size (i.e. a bottleneck) are well known, relatively few empirical studies of bottlenecks have been based on extensive temporally spaced samples of a population both before and after a bottleneck. Here we describe the results of one such study, utilising the Jenolan Caves (JC) population of the brush-tailed rock-wallaby (Petrogale penicillata). When first sampled in 1985 (n = 20) the JC population comprised ~90 individuals. Subsequently the population crashed, and by 1992 only seven individuals remained. In 1996 the entire population (n = 10) was again sampled. Genetic diversity in the pre- and post-crash JC population was compared using 11 polymorphic microsatellite loci and PCR–SSCP analysis of the mitochondrial DNA control region. Only a single unique control region haplotype was detected in the pre- and post-crash JC population, although variant haplotypes were present in other P. penicillata populations. Of the 35 microsatellite alleles present in the pre-crash population, nine (26%) were lost during the bottleneck. The average number of rare alleles declined by 72%, allelic diversity was reduced by 30% and average heterozygosity declined by 10%. These observations are consistent with theoretical predictions. Additional analyses revealed that a P.�penicillata female at Wombeyan Caves was the only survivor of a 1990/91 reintroduction attempt using animals from JC. Of the microsatellite alleles detected in this female, 21% (4/19) were no longer present in the post-crash JC population. Furthermore, the genetic profiles of animals from the recently discovered Taralga population indicate that they are not derived from JC stock, but represent a threatened remnant of a hitherto undetected natural P. penicillata population.

Intraspecific variation, sex-biased dispersal and phylogeography of the eastern grey kangaroo (Macropus giganteus)

Genetic information has played an important role in the development of management units by focusing attention on the evolutionary properties and genetics of populations. Wildlife authorities cannot hope to manage species effectively without knowledge of geographical boundaries and demic structure. The present investigation provides an analysis of mitochondrial DNA and microsatellite data, which is used to infer both historical and contemporary patterns of population structuring and dispersal in the eastern grey kangaroo (Macropus giganteus) in Australia. The average level of genetic variation across sample locations was one of the highest observed for marsupials (h=0.95, HE=0.82). Contrary to ecological studies, both genic and genotypic analyses reveal weak genetic structure of populations, where high levels of dispersal may be inferred up to 230 km. The movement of individuals was predominantly male-biased (average Nem=22.61, average Nfm=2.73). However, neither sex showed significant isolation by distance. On a continental scale, there was strong genetic differentiation and phylogeographic distinction between southern (TAS, VIC and NSW) and northern (QLD) populations, indicating a current and/or historical restriction of gene flow. In addition, it is evident that northern populations are historically more recent, and were derived from a small number of southern founders. Phylogenetic comparisons between M. g. giganteus and M. g. tasmaniensis indicated that the current taxonomic status of these subspecies should be revised as there was a lack of genetic differentiation between the populations sampled.

A rapid population expansion retains genetic diversity within European rabbits in Australia

The well documented historical translocations of the European rabbit (Oryctolagus cuniculus) offer an excellent framework to test the genetic effects of reductions in effective population size. It has been proposed that rabbits went through an initial bottleneck at the time of their establishment in Australia, as well as multiple founder events during the rabbit's colonization process. To test these hypotheses, genetic variation at seven microsatellite loci was measured in 252 wild rabbits from five populations across Australia. These populations were compared to each other and to data from Europe. No evidence of a genetic bottleneck was observed with the movement of 13 rabbits from Europe to Australia when compared to French data. Within Australia the distribution of genetic diversity did not reflect the suggested pattern of sequential founder effects. In fact, the current pattern of genetic variation in Australia is most likely a result of multiple factors including mutation, genetic drift and geographical differentiation. The absence of reduced genetic diversity is almost certainly a result of the rabbit's rapid population expansion at the time of establishment in Australia. These results highlight the importance of population growth following a demographic bottleneck, which largely determines the severity of genetic loss.

Echidna IgA supports mammalian unity and traditional Therian relationship

IgA is found only in birds and mammals where it is the principal immunoglobulin class found in secretions, providing protection at mucosal surfaces. The structure of IgA in birds is different from that of marsupials and eutherians. The avian heavy-chain constant region of IgA (Ca) consists of four domains, while marsupial and eutherian Ca consists of three domains plus a hinge. Here we describe the cloning and characterization of the heavy chain of IgA from the short-beaked echidna, Tachyglossus aculeatus, and report that monotreme Ca is composed of three domains plus a hinge, making it similar to its therian counterparts. The amino acid sequence identity of echidna Ca is approximately 47% with the therians and 30% with birds. Phylogenetic analysis of the Ca sequences provides strong support for the Theria hypothesis, which proposes that monotremes diverged prior to the separation of marsupial and eutherians, and directly contradicts the results of the mitochondrial data, which support a "Marsupionta" relationship which has marsupials and monotremes closer to each other. The characterization of the heavy chain of IgA from monotremes, in conjunction with the recent description of monotreme IgG and IgE nucleotide sequence, confirms that the "second big bang" of immunoglobulin evolution predated the divergence of extant mammals.

The First Comprehensive Genetic Linkage Map of a Marsupial: The Tammar Wallaby (Macropus eugenii)

The production of a marsupial genetic linkage map is perhaps one of the most important objectives in marsupial research. This study used a total of 353 informative meioses and 64 genetic markers to construct a framework genetic linkage map for the tammar wallaby (Macropus eugenii). Nearly all markers (93.8%) formed a significant linkage (LOD &gt; 3.0) with at least one other marker, indicating that the majority of the genome had been mapped. In fact, when compared with chiasmata data, &gt;70% (828 cM) of the genome has been covered. Nine linkage groups were identified, with all but one (LG7; X-linked) allocated to the autosomes. These groups ranged in size from 15.7 to 176.5 cM and have an average distance of 16.2 cM between adjacent markers. Of the autosomal linkage groups (LGs), LG2 and LG3 were assigned to chromosome 1 and LG4 localized to chromosome 3 on the basis of physical localization of genes. Significant sex-specific distortions toward reduced female recombination rates were revealed in 22% of comparisons. When comparing the X chromosome data to closely related species it is apparent that they are conserved in both synteny and gene order.

Ontogeny of immunoglobulin expression in the brushtail possum (Trichosurus vulpecula)

Marsupials, unlike eutherians, are born immunologically immature, without circulating lymphocytes or organised lymphoid tissue. Their immune response develops while they are in the pouch not in the uterus. In this study, the onset time of immunoglobulin expression in Trichosurus vulpecula pouch young was estimated by reverse transcription polymerase chain reaction. As in eutherian species, IgM heavy chain transcripts were detected first, at day 10 post partum. The first switched transcript, detected at day 18, was Calpha. Cgamma and Cvarepsilon transcripts were not present at day 72, but were seen at day 103, approximately corresponding to the time of release of the teat and exposure to new antigens, as well as the time of the loss of capacity to absorb maternal Igs through the gut.