Alexandre R Caetano

Resolving the evolution of extant and extinct ruminants with high-throughput phylogenomics

The Pecorans (higher ruminants) are believed to have rapidly speciated in the Mid-Eocene, resulting in five distinct extant families: Antilocapridae, Giraffidae, Moschidae, Cervidae, and Bovidae. Due to the rapid radiation, the Pecoran phylogeny has proven difficult to resolve, and 11 of the 15 possible rooted phylogenies describing ancestral relationships among the Antilocapridae, Giraffidae, Cervidae, and Bovidae have each been argued as representations of the true phylogeny. Here we demonstrate that a genome-wide single nucleotide polymorphism (SNP) genotyping platform designed for one species can be used to genotype ancient DNA from an extinct species and DNA from species diverged up to 29 million years ago and that the produced genotypes can be used to resolve the phylogeny for this rapidly radiated infraorder. We used a high-throughput assay with 54,693 SNP loci developed for Bos taurus taurus to rapidly genotype 678 individuals representing 61 Pecoran species. We produced a highly resolved phylogeny for this diverse group based upon 40,843 genome-wide SNP, which is five times as many informative characters as have previously been analyzed. We also establish a method to amplify and screen genomic information from extinct species, and place Bison priscus within the Bovidae. The quality of genotype calls and the placement of samples within a well-supported phylogeny may provide an important test for validating the fidelity and integrity of ancient samples. Finally, we constructed a phylogenomic network to accurately describe the relationships between 48 cattle breeds and facilitate inferences concerning the history of domestication and breed formation.

A physical map of the bovine genome

BackgroundCattle are important agriculturally and relevant as a model organism. Previously described genetic and radiation hybrid (RH) maps of the bovine genome have been used to identify genomic regions and genes affecting specific traits. Application of these maps to identify influential genetic polymorphisms will be enhanced by integration with each other and with bacterial artificial chromosome (BAC) libraries. The BAC libraries and clone maps are essential for the hybrid clone-by-clone/whole-genome shotgun sequencing approach taken by the bovine genome sequencing project.ResultsA bovine BAC map was constructed with HindIII restriction digest fragments of 290,797 BAC clones from animals of three different breeds. Comparative mapping of 422,522 BAC end sequences assisted with BAC map ordering and assembly. Genotypes and pedigree from two genetic maps and marker scores from three whole-genome RH panels were consolidated on a 17,254-marker composite map. Sequence similarity allowed integrating the BAC and composite maps with the bovine draft assembly (Btau3.1), establishing a comprehensive resource describing the bovine genome. Agreement between the marker and BAC maps and the draft assembly is high, although discrepancies exist. The composite and BAC maps are more similar than either is to the draft assembly.ConclusionFurther refinement of the maps and greater integration into the genome assembly process may contribute to a high quality assembly. The maps provide resources to associate phenotypic variation with underlying genomic variation, and are crucial resources for understanding the biology underpinning this important ruminant species so closely associated with humans.

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

BackgroundThe recently constructed river buffalo whole-genome radiation hybrid panel (BBURH5000) 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.ResultsA total of 3990 markers were typed on the BBURH5000 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 cR5000. 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).ConclusionWe 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.

Genetic diversity analysis of two commercial breeds of pigs using genomic and pedigree data.

BackgroundGenetic improvement in livestock populations can be achieved without significantly affecting genetic diversity if mating systems and selection decisions take genetic relationships among individuals into consideration. The objective of this study was to examine the genetic diversity of two commercial breeds of pigs. Genotypes from 1168 Landrace (LA) and 1094 Large White (LW) animals from a commercial breeding program in Brazil were obtained using the Illumina PorcineSNP60 Beadchip. Inbreeding estimates based on pedigree (Fx) and genomic information using runs of homozygosity (FROH) and the single nucleotide polymorphisms (SNP) by SNP inbreeding coefficient (FSNP) were obtained. Linkage disequilibrium (LD), correlation of linkage phase (r) and effective population size (Ne) were also estimated.ResultsEstimates of inbreeding obtained with pedigree information were lower than those obtained with genomic data in both breeds. We observed that the extent of LD was slightly larger at shorter distances between SNPs in the LW population than in the LA population, which indicates that the LW population was derived from a smaller Ne. Estimates of Ne based on genomic data were equal to 53 and 40 for the current populations of LA and LW, respectively. The correlation of linkage phase between the two breeds was equal to 0.77 at distances up to 50xa0kb, which suggests that genome-wide association and selection should be performed within breed. Although selection intensities have been stronger in the LA breed than in the LW breed, levels of genomic and pedigree inbreeding were lower for the LA than for the LW breed.ConclusionsThe use of genomic data to evaluate population diversity in livestock animals can provide new and more precise insights about the effects of intense selection for production traits. Resulting information and knowledge can be used to effectively increase response to selection by appropriately managing the rate of inbreeding, minimizing negative effects of inbreeding depression and therefore maintaining desirable levels of genetic diversity.

Linkage disequilibrium, persistence of phase and effective population size estimates in Hereford and Braford cattle

BackgroundThe existence of moderate to high levels of linkage disequilibrium (LD) between genetic markers and quantitative trait loci (QTL) affecting traits of interest is fundamental for the success of genome-wide association (GWAS) and genomic selection (GS) studies. Knowledge about the extent and the pattern of LD in livestock populations is essential to determine the density of single nucleotide polymorphisms (SNP) required for accurate GWAS and GS. Moreover, observed LD is related to historical effective population sizes (Ne), and can provide insights into the genetic diversity history of populations. Estimates of the consistency of linkage phase across breeds (RH,B) can be used to determine if there is sufficient relationship to use pooled reference populations in multi-breed GS programs. The objective of this study was to estimate LD levels, persistence of phase and effective population size in Hereford and Braford cattle populations sampled in Brazil.ResultsMean LD estimates, measured using the squared correlation of alleles at two loci (r2), obtained between adjacent SNP across all chromosomes were 0.21u2009±u20090.27 for Herefords (391 samples with 41,241 SNP) and 0.16u2009±u20090.22 for Brafords (2044 samples and 41,207 SNP). Estimated r2 wasu2009>u20090.2 and 0.3, respectively, for 34 and 25xa0% of adjacent markers in Herefords, and 26 and 17xa0% in Brafords. Estimated Ne for Brafords and Herefords at the current generation was 220 and 153 individuals, respectively. The two breeds demonstrated moderate to strong persistence of phase at all distances (RH,Bu2009=u20090.53 to 0.97). The largest phase correlations were found in the 0 to 50xa0Kb bins (RH,Bu2009=u20090.92 to 0.97). Estimated LD decreased rapidly with increasing distance between SNP, however, useful linkage for GWAS and GS (r2u2009>u20090.2) was found spanning to ~50xa0Kb.ConclusionsPanels containing about 50,000 and 150,000 SNP markers are necessary to detect minimal levels of LD between adjacent markers that would be useful for GWAS and GS studies to Hereford and Braford breeds, respectively. Markers are expected to be linked to the same QTL alleles in distancesu2009<u200950xa0Kb in both populations due to observed high persistence of phase levels.

Molecular Phylogeny Inferred from the Concatenated Genes of Two Neotropical Catfish Species and Implications for Conservation

The Neotropics host the most diverse ichthyofauna in the world, with catfish species forming one of the most diverse groups in the region. Nuclear (RAG1) and mitochondrial (ATPase and Cytb) markers were analyzed to identify genetic variability in populations of Pseudoplatystoma reticulatum and Pseudoplatystoma corruscans from the La Plata and Sao Francisco Basins. Bayesian topology identified the division of P. corruscansinto two main clades. One of these clades was formed of samples from the Sao Francisco Basin and the other was formed of samples from the Parana+Paraguay Basins. P. reticulatum was grouped together without any clear geographic or taxonomic patterns in Bayesian topology. While only a few common nuclear haplotypes were widely identified in both species, there was great variability in the mitochondrial sequences. The genetic and geographical distance correlations were tested using the Mantel permutation, which detected no significant relationships.The results of the present study suggest a panmitic population for both species (excluding P. corruscansin the Sao Francisco Basin, which is suggested as a new species), with the greatest diversity concentrated in the region covered by the Pantanal biome, and the lowest diversity in Mogi Guacu, in the Parana Basin. These findings support the establishment of public conservation policies and provide information regarding genetic diversity and population differentiation patterns for these ecological and economically important species.

Complete mitochondrial genome from South American catfish Pseudoplatystoma reticulatum (Eigenmann & Eigenmann) and its impact in Siluriformes phylogenetic tree

The cachara (Pseudoplatystoma reticulatum) is a Neotropical freshwater catfish from family Pimelodidae (Siluriformes) native to Brazil. The species is of relative economic importance for local aquaculture production and basic biological information is under development to help boost efforts to domesticate and raise the species in commercial systems. The complete cachara mitochondrial genome was obtained by assembling Illumina RNA-seq data from pooled samples. The full mitogenome was found to be 16,576xa0bp in length, showing the same basic structure, order, and genetic organization observed in other Pimelodidae, with 13 protein-coding genes, 2 rNA genes, 22 trNAs, and a control region. Observed base composition was 24.63% T, 28.47% C, 31.45% A, and 15.44% G. With the exception of NAD6 and eight tRNAs, all of the observed mitochondrial genes were found to be coded on the H strand. A total of 107 SNPs were identified in P. reticulatum mtDNA, 67 of which were located in coding regions. Of these SNPs, 10 result in amino acid changes. Analysis of the obtained sequence with 94 publicly available full Siluriformes mitogenomes resulted in a phylogenetic tree that generally agreed with available phylogenetic proposals for the order. The first report of the complete Pseudoplatystoma reticulatum mitochondrial genome sequence revealed general gene organization, structure, content, and order similar to most vertebrates. Specific sequence and content features were observed and may have functional attributes which are now available for further investigation.

Validation of a microsatellite panel for parentage testing of locally adapted and commercial goats in Brazil

Brazilian goats are generally kept in small herds and extensive rearing systems, mainly in the northeastern region of the country. Despite production improvement in recent years, the lack of pedigree control has affected genetic progress. This study aimed to validate a panel of 16 microsatellites for parentage testing in locally adapted and commercial goats breeds raised in Brazil, as well as to compare its efficiency with the panel recommended by the Brazilian Ministry of Agriculture, Livestock and Supplies (MAPA) in 2004. The number of alleles and expected heterozygosity (He) per marker ranged from four to 18, and from 0.051 to 0.831, respectively. Using all markers, 100% of parentage cases of the validation dataset were resolved with a strict confidence level of 95%. The 16 microsatellites panel showed adequate exclusion power (99.99%) and identity accuracy (99.99%). Suggestions for improvement of the marker panel endorsed by MAPA are provided.