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

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

Mitochondrial genotype segregation and effects during mammalian development: Applications to biotechnology

Mitochondria are endosymbiotic organelles responsible for energy production in practically every eukaryotic cell. Their uniparental fashion of inheritance, maternally inherited in mammals, and the homogeneity of mitochondrial DNA (mtDNA) within individuals and matrilineages, are biological phenomena that remain unexplained. This paper reviews some of the recent findings on mitochondrial influences on the manner in which embryos develop and how their genotypes are inherited in mammals, with particular emphasis on the genetic bottleneck effect. Animal models carrying a mix of mtDNAs (heteroplasmic) have been produced by karyoplast and cytoplast transplantation to analyze the segregation patterns at different stages during embryogenesis, in fetuses and offspring. Comparisons performed between murine and bovine reveal interesting changes in segregation and replication of transplanted mtDNAs. We have recently obtained Bos indicus and Bos taurus fetuses and calves from embryos reconstructed using enucleated polymorphic oocytes of Bos taurus origin. These and other findings on mitochondrial biology will have important implications in determining the cytoplasmic genotype of clones and in the preservation of endangered breeds and species. (C) 1999 by Elsevier B.V.

New Loci Of Lychnophora Ericoides And Transferability To Lychnophora Pinaster, Endangered Medicinal Species From Brazil.

Lychnophora ericoides and Lychnophora pinaster are species used in popular medicine as analgesic or anti-inflammatory agents to treat contusions, rheumatism, and insect bites. In this study, 21 simple sequence repeat loci of L. ericoides were developed and transferred to L. pinaster. Three populations of L. ericoides and 2 populations of L. pinaster were evaluated; they were collected in the State of Minas Gerais. Population parameters were estimated, and the mean values of observed and expected heterozygosity were 0.297 and 0.408 (L. ericoides) and 0.228 and 0.310 (L. pinaster), respectively. Greater genetic variability was observed within populations than between populations of L. ericoides (62 and 37%) and L. pinaster (97 and 2.8%). These results provide information for genetic conservation and taxonomic studies of these endangered species.

Induced spawning of the endangered Neotropical species Steindachneridion parahybae (Siluriformes: Pimelodidae)

The "surubim do Paraíba" (Steindachneridion parahybae) is a freshwater catfish endemic to the Paraíba do Sul River basin, Brazil. This species has been seriously threatened by environmental disturbances in the last several decades. Wild Steindachneridion parahybae males and females were collected in 2003 and taken to the hatchery of a power plant of the Companhia Energética de São Paulo (CESP). Steindachneridion parahybae broodstocks were artificially induced to reproduce in December 2003 using a combination of carp pituitary extract (CPE) and human chorionic gonadotropin (hCG). Oocytes and milt were stripped; the fertilized eggs were transferred to 60-liter conical incubators and hatched larvae distributed in nine horizontal trays. Exogenous feed was started just after yolk sac absorption. A high rate of cannibalism and photophobia were observed during the larval period, resulting in a 26% survival rate from larvae to fingerlings.

Karyotypic analyses and morphological comments on the endemic and endangered Brazilian painted tree rat Callistomys pictus (Rodentia, Echimyidae)

The genus Callistomys belongs to the rodent family Echimyidae, subfamily Echimyinae, and its only living representative is Callistomys pictus, a rare and vulnerable endemic species of the state of Bahia, Brazil. Callistomys has been previously classified as Nelomys, Loncheres, Isothrix and Echimys. In this paper we present the karyotype of Callistomys pictus, including CBG and GTG-banding patterns and silver staining of the nucleolus organizer regions (Ag-NORs). Comments on Callistomys pictus morphological traits and a compilation of Echimyinae chromosomal data are also included. Our analyses revealed that Callistomys can be recognized both by its distintinctive morphology and by its karyotype.

An assessment of population structure in eight breeds of cattle using a whole genome SNP panel

Background: Analyses of population structure and breed diversity have provided insight into the origin and evolution of cattle. Previously, these studies have used a low density of microsatellite markers, however, with the large number of single nucleotide polymorphism markers that are now available, it is possible to perform genome wide population genetic analyses in cattle. In this study, we used a high-density panel of SNP markers to examine population structure and diversity among eight cattle breeds sampled from Bos indicus and Bos taurus. Results: Two thousand six hundred and forty one single nucleotide polymorphisms ( SNPs) spanning all of the bovine autosomal genome were genotyped in Angus, Brahman, Charolais, Dutch Black and White Dairy, Holstein, Japanese Black, Limousin and Nelore cattle. Population structure was examined using the linkage model in the program STRUCTURE and Fst estimates were used to construct a neighbor-joining tree to represent the phylogenetic relationship among these breeds. Conclusion: The whole-genome SNP panel identified several levels of population substructure in the set of examined cattle breeds. The greatest level of genetic differentiation was detected between the Bos taurus and Bos indicus breeds. When the Bos indicus breeds were excluded from the analysis, genetic differences among beef versus dairy and European versus Asian breeds were detected among the Bos taurus breeds. Exploration of the number of SNP loci required to differentiate between breeds showed that for 100 SNP loci, individuals could only be correctly clustered into breeds 50% of the time, thus a large number of SNP markers are required to replace the 30 microsatellite markers that are currently commonly used in genetic diversity studies.

Influence of Chinese breeds on pork quality of commercial pig lines

We compared carcass and meat quality of pigs from the same sire line and two different dam lines, one that included Chinese breeds and one that did not. Line A consisted of 1/4 Landrace, 1/2 Large White, 1/8 Chinese breeds (Meishan, Fengjing, Jiaxing), and 1/8 Large White, Duroc and Pietrain, and line B consisted of 1/2 Large White and 1/2 Pietrain. The animals (N = 144) were slaughtered at a live weight of 108 kg. Backfat thickness, percentage of lean meat, pH 24 h after slaughter, meat color, percentage of drip loss, and percentage of intramuscular fat were measured and compared using analysis of variance in a completely randomized design; the BioEstat 5.0 test was applied for the comparison of means at a significance level of 5% for all analyses. Backfat was significantly lower for line A (12.78 mm) than for line B (15.90 mm). The pH measured 24 h after slaughter was significantly lower in line A (5.68) compared to line B (5.84). Percent lean meat was significantly higher for line A (61.21%) compared to line B (59.72%). Percentage drip loss was significantly higher in line A (2.73%) than in line B (2.23%). Percentage intramuscular fat and meat color were not significantly different between the lines. The inclusion of Chinese breeds produced a higher percentage of lean meat and reduced fat thickness, along with increased heterosis, which are important characteristics for breeding programs.

The taxonomic status of the endangered thin-spined porcupine, Chaetomys subspinosus (Olfers, 1818), based on molecular and karyologic data

Background: The thin-spined porcupine, also known as the bristle-spined rat, Chaetomys subspinosus (Olfers, 1818), the only member of its genus, figures among Brazilian endangered species. In addition to being threatened, it is poorly known, and even its taxonomic status at the family level has long been controversial. The genus Chaetomys was originally regarded as a porcupine in the family Erethizontidae, but some authors classified it as a spiny-rat in the family Echimyidae. Although the dispute seems to be settled in favor of the erethizontid advocates, further discussion of its affinities should be based on a phylogenetic framework. In the present study, we used nucleotide-sequence data from the complete mitochondrial cytochrome b gene and karyotypic information to address this issue. Our molecular analyses included one individual of Chaetomys subspinosus from the state of Bahia in northeastern Brazil, and other hystricognaths. Results: All topologies recovered in our molecular phylogenetic analyses strongly supported Chaetomys subspinosus as a sister clade of the erethizontids. Cytogenetically, Chaetomys subspinosus showed 2n = 52 and FN = 76. Although the sexual pair could not be identified, we assumed that the X chromosome is biarmed. The karyotype included 13 large to medium metacentric and submetacentric chromosome pairs, one small subtelocentric pair, and 12 small acrocentric pairs. The subtelocentric pair 14 had a terminal secondary constriction in the short arm, corresponding to the nucleolar organizer region (Ag-NOR), similar to the erethizontid Sphiggurus villosus, 2n = 42 and FN = 76, and different from the echimyids, in which the secondary constriction is interstitial. Conclusion: Both molecular phylogenies and karyotypical evidence indicated that Chaetomys is closely related to the Erethizontidae rather than to the Echimyidae, although in a basal position relative to the rest of the Erethizontidae. The high levels of molecular and morphological divergence suggest that Chaetomys belongs to an early radiation of the Erethizontidae that may have occurred in the Early Miocene, and should be assigned to its own subfamily, the Chaetomyinae.

Characterization of unisexual flower development in the endangered mahogany tree Swietenia macrophylla King. (Meliaceae)

The selection of candidate plus trees of desirable phenotypes from tropical forest trees and the rapid devastation of the natural environments in which these trees are found have created the need for a more detailed knowledge of the floral and reproductive biology of tropical tree species. In this article, the organogenic processes related to unisexual flower development in tropical mahogany, Swietenia macrophylla, are described. Mahogany inflorescences at different developmental stages were evaluated using scanning electron microscopy or optical microscopy of histological sections. The unisexual flowers of S. macrophylla are usually formed in a thyrse, in which the positions of the female and male flowers are not random. Differences between male and female flowers arise late during development. Both female and male flowers can only be structurally distinguished after stage 9, where ovule primordia development is arrested in male flowers and microspore development is aborted in female flower anthers. After this stage, male and female flowers can be distinguished by the naked eye as a result of differences in the dimensions of the gynoecium. The floral characteristics of S. macrophylla (distribution of male and female flowers within the inflorescence, and the relative number of male to female flowers) have practical implications for conservation strategies of this endangered species. (c) 2008 The Linnean Society of London, Botanical Journal of the Linnean Society, 2008, 156, 529-535.

Fruit flies (Diptera, Tephritidae) and their associations with native host plants in a remnant area of the highly endangered Atlantic Rain Forest in the State of Espirito Santo, Brazil

The results presented in this paper refer to a host survey, lasting approximately three and a half years (February 2003-july 2006), undertaken in the Vale do Rio Doce Natural Reserve, a remnant area of the highly endangered Atlantic Rain Forest located in Linhares County, State of Espirito Santo, Brazil. A total of 330 fruit samples were collected from native plants, representing 248 species and 51 plant families. Myrtaceae was the most diverse family with 54 sampled species. Twenty-eight plant species, from ten families, are hosts of ten Anastrepha species and of Ceratitis capitata (Wiedemann). Among 33 associations between host plants and fruit flies, 20 constitute new records, including the records of host plants for A. fumipennis Lima and A. nascimentoi Zucchi. The findings were discussed in the light of their implications for rain forest conservation efforts and the study of evolutionary relationships between fruit flies and their hosts.

Chemical composition of whole body and carcass of Bos indicus and tropically adapted Bos taurus breeds

Relationships between the chemical composition of the 9th- to 11th-rib section and the chemical composition of the carcass and empty body were evaluated for Bos indicus (108 Nellore and 36 Guzerah; GuS) and tropically adapted Bos taurus (56 Caracu; CaS) bulls, averaging 20 to 24 mo of age at slaughter. Nellore cattle were represented by 56 animals from the selected herd (NeS) and 52 animals from the control herd (NeC). The CaS and GuS bulls were from selected herds. Selected herds were based on 20 yr of selection for postweaning BW. Carcass composition was obtained after grinding, homogenizing, sampling, and analyzing soft tissue and bones. Similarly, empty body composition was obtained after grinding, homogenizing, sampling, analyzing, and combining blood, hide, head + feet, viscera, and carcass. Bulls were separated into 2 groups. Group 1 was composed of 36 NeS, 36 NeC, 36 CaS, and 36 GuS bulls and had water, ether extract (EE), protein, and ash chemically determined in the 9th- to 11th-rib section and in the carcass. Group 2 was composed of 20 NeS, 16 NeC, and 20 CaS bulls and water, EE, protein, and ash were determined in the 9th-to 11th-rib section, carcass, and empty body. Linear regressions were developed between the carcass and the 9th-to 11th-rib section compositions for group 1 and between carcass and empty body compositions for group 2. The 9th-to 11th-rib section percentages of water (RWt) and EE (RF) predicted the percentages of carcass water (CWt) and carcass fat (CF) with high precision: CWt, % = 29.0806 + 0.4873 x RWt, % (r(2) = 0.813, SE = 1.06) and CF, % = 10.4037 + 0.5179 x RF, % (r(2) = 0.863, SE = 1.26), respectively. Linear regressions between percentage of CWt and CF and empty body water (EBWt) and empty body fat (EBF) were also predicted with high precision: EBWt, % = -9.6821 + 1.1626 x CWt, % (r(2) = 0.878, SE = 1.43) and EBF, % = 0.3739 + 1.0386 x CF, % (r(2) = 0.982, SE = 0.65), respectively. Chemical composition of the 9th-to 11th-rib section precisely estimated carcass percentages of water and EE. These regressions can accurately predict carcass and empty body compositions for Nellore, Guzerah, and Caracu breeds.

Evaluation of carcass characteristics of Bos indicus and tropically adapted Bos taurus breeds selected for postweaning weight

Data from 9 studies were compiled to evaluate the effects of 20 yr of selection for postweaning weight (PWW) on carcass characteristics and meat quality in experimental herds of control Nellore (NeC) and selected Nellore (NeS), Caracu (CaS), Guzerah (GuS), and Gir (GiS) breeds. These studies were conducted with animals from a genetic selection program at the Experimental Station of Sertaozinho, Sao Paulo State, Brazil. After the performance test (168 d postweaning), bulls (n = 490) from the calf crops born between 1992 and 2000 were finished and slaughtered to evaluate carcass traits and meat quality. Treatments were different across studies. A meta-analysis was conducted with a random coefficients model in which herd was considered a fixed effect and treatments within year and year were considered as random effects. Either calculated maturity degree or initial BW was used interchangeably as the covariate, and least squares means were used in the multiple-comparison analysis. The CaS and NeS had heavier (P = 0.002) carcasses than the NeC and GiS; GuS were intermediate. The CaS had the longest carcass (P < 0.001) and heaviest spare ribs (P < 0.001), striploin (P < 0.001), and beef plate (P = 0.013). Although the body, carcass, and quarter weights of NeS were similar to those of CaS, NeS had more edible meat in the leg region than did CaS bulls. Selection for PWW increased rib-eye area in Nellore bulls. Selected Caracu had the lowest (most favorable) shear force values compared with the NeS (P = 0.003), NeC (P = 0.005), GuS (P = 0.003), and GiS (P = 0.008). Selection for PWW increased body, carcass, and meat retail weights in the Nellore without altering dressing percentage and body fat percentage.

Evolutionary distinctiveness and status of the endangered Lake Eacham rainbowfish (Melanotaenia eachamensis)

The Lake Eacham rainbowfish (Melanotaenia eachamensis) was declared extinct in the wild in the late 1980s after it disappeared from its only known locality, an isolated crater lake in northeast Queensland. Doubts have been raised about whether this taxon is distinct from surrounding populations of the eastern rainbowfish (Melanotaenia splendida splendida). We examined the evolutionary distinctiveness of M. eachamensis, obtained from captive stocks, relative to M. s. splendida through analysis of variation in mtDNA sequences, nuclear microsatellites, and morphometric characters Captive M. eachamensis had mtDNAs that were highly divergent from those in most populations of M. s. splendida. A broader geographic survey using RFLPs revealed some populations initially identified as M. s. splendida, that carried eachamensis mtDNA, whereas some others had mixtures of eachamensis and splendida mtDNA. The presence of eachamensis-like mtDNA in these populations could in principle be due to (1) sorting of ancestral polymorphisms, (2) introgression of M. eachamensis mtDNA into M. s. splendida, or (3) incorrect species boundaries, such that some populations currently assigned to M. s. splendida are M. eachamensis or are mixtures of the two species. These alternatives hypotheses were evaluated through comparisons of four nuclear microsatellite loci and morphometrics and meristics. In analyses of both data sets, populations of M. s. splendida with eachamensis mtDNA were more similar to captive M. eachamensis than to M. s. splendida with splendida mtDNA, supporting hypothesis 3. These results are significant for the management of M. eachamensis in several respects. First the combined molecular and morphological evidence indicates that M. eachamensis is a distinct species and a discrete evolutionarily significant unit worthy of conservation effort. Second it appears that the species boundary between M. eachamensis and M. s. splendida has been misdiagnosed such that there are extant populations on the Atherton Tableland as well as areas where both forms coexist. Accordingly we suggest that M. eachamensis be listed as vulnerable, rather than critical (or extinct in the wild). Third, the discovery of extant but genetically divergent populations of M. eachamensis on the Atherton Tableland broadens the options for future reintroductions to Lake Eacham.

Microsatellite diversity and effective population size in a germplasm bank of Hymenaea courbaril var. stilbocarpa (Leguminosae), an endangered tropical tree: recommendations for conservation

Deforestation in southeast Brazil has led to the extinction of Hymenaea courbaril var. stilbocarpa and ex situ conservation has been established. In this study, the levels of genetic diversity and the effective population size of H. courbaril in a germplasm bank were investigated using six nuclear microsatellite loci. A total of 79 and 91 alleles were found in 65 seed-trees and their 176 offspring, respectively. Offspring have a higher average number of alleles per locus (A = 15.2) than seed-trees (A = 13.2), but lower observed heterozygosity (offspring: H (o) = 0.566; seed-trees: H (o) = 0.607). The estimate of outcrossing rate shows that the study population is perfectly outcrossed (t (m) = 0.978, P > 0.05). Significant deviations from random mating were detected through mating among relatives and correlated matings. The average variance in effective population size for each family was 2.63, with a total effective population size retained in the bank of 170.1. These results confirm that the preserved population of H. courbaril retains substantial genetic variability.

Development and characterization of 15 microsatellite loci for Cariniana estrellensis and transferability to Cariniana legalis, two endangered tropical tree species

From a genomic library enriched for GA/CA repeats, 15 highly polymorphic microsatellite markers were developed for Cariniana estrellensis, a tropical forest tree. The microsatellite loci were screened in 49 mature trees found between Pardo river and Mogi-Gua double dagger u river basins, in the state of So Paulo, Brazil. A total of 140 alleles were detected with an average of 9.33 alleles per locus. The expected heterozygosity ranged from 0.37 to 0.88. These loci showed a high probability of paternity exclusion. Additionally, 12 loci were effectively transferred to Cariniana legalis. High levels of polymorphism make the present SSR markers useful for population genetic studies.