Marcos De Donato

Genotyping-by-Sequencing (GBS): A Novel, Efficient and Cost-Effective Genotyping Method for Cattle Using Next-Generation Sequencing

High-throughput genotyping methods have increased the analytical power to study complex traits but high cost has remained a barrier for large scale use in animal improvement. We have adapted genotyping-by-sequencing (GBS) used in plants for genotyping 47 animals representing 7 taurine and indicine breeds of cattle from the US and Africa. Genomic DNA was digested with different enzymes, ligated to adapters containing one of 48 unique bar codes and sequenced by the Illumina HiSeq 2000. PstI was the best enzyme producing 1.4 million unique reads per animal and initially identifying a total of 63,697 SNPs. After removal of SNPs with call rates of less than 70%, 51,414 SNPs were detected throughout all autosomes with an average distance of 48.1 kb, and 1,143 SNPs on the X chromosome at an average distance of 130.3 kb, as well as 191 on unmapped contigs. If we consider only the SNPs with call rates of 90% and over, we identified 39,751 on autosomes, 850 on the X chromosome and 124 on unmapped contigs. Of these SNPs, 28,843 were not tightly linked to other SNPs. Average marker density per autosome was highly correlated with chromosome size (coefficient of correlation  = −0.798, r2  = 0.637) with higher density in smaller chromosomes. Average SNP call rate was 86.5% for all loci, with 53.0% of the loci having call rates >90% and the average minor allele frequency being 0.212. Average observed heterozygosity ranged from 0.046–0.294 among individuals, and from 0.064–0.197 among breeds, with Brangus showing the highest diversity as expected. GBS technique is novel, flexible, sufficiently high-throughput, and capable of providing acceptable marker density for genomic selection or genome-wide association studies at roughly one third of the cost of currently available genotyping technologies.

New ribosomal RNA gene locations in Gossypium hirsutum mapped by meiotic FISH

Abstract. In this study we have mapped newly identified rDNA loci in Gossypium hirsutum. Four new minor 18S-26S rDNA loci, in addition to the sites previously identified, were mapped using fluorescence in situ hybridization (FISH) to heterozygous translocation (NT) quadrivalents (IVs). The newly detected 18S-26S rDNA loci were mapped to the right arms of chromosomes 8, 9, 15, 17, 19, 20, and 23 and the left arms of chromosomes 5, 11, 12, and 14. Using the rDNA loci as common reference points, we detected several erroneous arm assignments in the previously published map of NT breakpoints. The data are summarized in the form of an integrated map for all 17 known rDNA loci, relative to centromeres, telomeres, and NT breakpoints. This information will facilitate future locus-specific research on rRNA gene evolution and function.

Morphological and microsatellite DNA diversity of Nigerian indigenous sheep

BackgroundSheep is important in the socio-economic lives of people around the world. It is estimated that more than half of our once common livestock breeds are now endangered. Since genetic characterization of Nigerian sheep is still lacking, we analyzed ten morphological traits on 402 animals and 15 microsatellite DNA markers in 384 animals of the 4 Nigerian sheep breeds to better understand genetic diversity for breeding management and germplasm conservation.ResultsMorphological traits of Uda and Balami were significantly (P < 0.05) higher than Yankasa, which were both higher than West African Dwarf (WAD) sheep. Stepwise discriminant analysis showed tail length, rump height, chest girth, ear length and chest depth as the most discriminating variables for classification. Mahalanobis distances show the least differentiation between Uda and Balami and the largest between WAD and Balami sheep. While 93.3% of WAD sheep were correctly assigned to their source genetic group, 63.9% of Yankasa, 61.2% of Balami and 45.2% of Uda were classified correctly by nearest neighbour discriminant analysis. The overall high Polymorphism Information Content (PIC) of all microsatellite markers ranged from 0.751 to 0.927 supporting their use in genetic characterization. Expected heterozygosity was high for all loci (0.783 to 0.93). Mean heterozygote deficiency across all populations (0.171 to 0.534) possibly indicate significant inbreeding (P < 0.05). Mean values for FST, FIT and FIS statistics across all loci were 0.088, 0.394 and 0.336 respectively. Yankasa and Balami are the most closely related breeds (DA = 0.184) while WAD and Balami are the farthest apart breeds (DA = 0.665), which is coincident with distance based on morphological analysis and population structure assessed by STRUCTURE.ConclusionsThese results suggest that within-breed genetic variation in Nigerian sheep is higher than between-breeds and may be a valuable tool for genetic improvement and conservation. The higher genetic variability in Yankasa suggests the presence of unique ancestral alleles reflecting the presence of certain functional genes which may result in better adaptability in more agro-ecological zones of Nigeria. These genetic characteristics are potentially useful in planning improvement and conservation strategies in Nigerian indigenous sheep.

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.

Molecular characterization of Salmonella strains in individuals with acute diarrhea syndrome in the State of Sucre, Venezuela

INTRODUCTION In Venezuela, acute diarrheic syndrome (ADS) is a primary cause of morbi-mortality, often involving the Salmonella genus. Salmonella infections are associated with acute gastroenteritis, one of the most common alimentary intoxications, and caused by the consumption of contaminated water and food, especially meat. METHODS Conventional and molecular methods were used to detect Salmonella strains from 330 fecal samples from individuals of different ages and both sexes with ADS. Polymerase chain reaction (PCR) was used for the molecular characterization of Salmonella, using invA, sefA, and fliC genes for the identification of this genus and the serotypes Enteritidis and Typhimurium, respectively. RESULTS The highest frequency of individuals with ADS was found in children 0-2 years old (39.4%), and the overall frequency of positive coprocultures was 76.9%. A total of 14 (4.2%) strains were biochemically and immunologically identified as Salmonella enterica subsp. enterica, of which 7 were classified as belonging to the Enteritidis serotype, 4 to the Typhimurium serotype, and 3 to other serotypes. The S. enterica strains were distributed more frequently in the age groups 3-4 and 9-10 years old. CONCLUSIONS The molecular characterization method used proved to be highly specific for the typing of S. enterica strains using DNA extracted from both the isolated colonies and selective enrichment broths directly inoculated with fecal samples, thus representing a complementary tool for the detection and identification of ADS-causing bacteria.

Association of SNP variants of MHC Class II DRB gene with thermo-physiological traits in tropical goats

Host defense in vertebrates depend on many secreted regulatory proteins such as major histocompatibility complex (MHC) class II which provide important regulatory and effector functions of T cells. Gene polymorphism in the second exon of Capra-DRB gene in three major Nigerian goat breeds [West African Dwarf (WAD), Red Sokoto (RS), and Sahel (SH)] was analyzed by restriction fragment length polymorphisms (RFLP). Four restriction enzymes, BsaHI, AluI, HaeIII, and SacII, were utilized. The association between the polymorphic sites and some heat tolerance traits were also investigated in a total of 70 WAD, 90 RS, and 50 SH goats. Fourteen different types of alleles identified in the Nigerian goats, four of which were found in the peptide coding region (A57G, Q89R, G104D, and T112I), indicate a high degree of polymorphism at the DRB locus in this species. An obvious excess (P < 0.01) of non-synonymous substitutions than synonymous (dN/dS) in this locus is a reflection of adaptive evolution and positive selection. The phylogenetic trees revealed largely species-wise clustering in DRB gene. BsaHI, AluI, HaeIII, and SacII genotype frequencies were in Hardy-Weinberg equilibrium (P > 0.05), except AluI in RS goats and HaeIII in WAD goats (P < 0.05). The expected heterozygosity (H), which is a measure of gene diversity in the goat populations, ranged from 0.16 to 0.50. Genotypes AA (BsaHI), GG, GC and CC (AluI) and GG, GA, AA (HaeIII) appeared better in terms of heat tolerance. The heat-tolerant ability of SH and RS goats to the hot and humid tropical environment of Nigeria seemed better than that of the WAD goats. Sex effect (P < 0.05) was mainly on pulse rate and heat stress index, while there were varying interaction effects on heat tolerance. Variation at the DRB locus may prove to be important in possible selection and breeding for genetic resistance to heat stress in the tropics.

Molecular evolution of type II MAGE genes from ancestral MAGED2 gene and their phylogenetic resolution of basal mammalian clades

Type II melanoma-associated antigens (MAGE) are a subgroup of about a dozen proteins found in various locations in the genome and expressed in normal tissues, thus are not related to cancer as the type I MAGE genes. This gene family exists as a single copy in non-mammals and monotremata, but found as two copies in metatherians and occur as a diverse group in all eutherians. Our studies suggest MAGED2 as the ancestor of this subfamily and the most likely evolutionary history of eutherian type II MAGE genes is hereby proposed based on synteny conservation, phylogenetic relations, genome location, homology conservation, and the protein and gene structures. Type II genes can be divided into two: those with 13 exons (MAGED1, MAGED2, TRO, and MAGED4) and those with only one exon (MAGEE1, MAGEE2, MAGEF1, NSMCE3, MAGEH1, MAGEL2, and NDN) with different evolutionary patterns. Our results suggest a need to change the gene nomenclature to MAGE1 (the ancestral gene), currently designated as LOC103095671 and LOC100935086, in opossum and Tasmanian devil, respectively, and MAGE2 (the duplicated one), currently designated as LOC100617402 and NDNL2, respectively, to avoid confusion. We reconstructed the phylogenetic relationships among 23 mammalian species using the combined sequences of MAGED1, MAGED2, MAGEL2, and NDN, because of their high divergence, and found high levels of support, being able to resolve the phylogenetic relationships among Euarchontoglires, Laurasiatheria, Afrotheria, and Xenarthra, as an example that small, but phylogenetically informative sequences, can be very useful for resolving basal mammalian clades.

Genetic Diversity of Bovine Major Histocompatibility Complex Class II DRB3 locus in cattle breeds from Asia compared to those from Africa and America

Genetic polymorphisms and diversity of BoLA-DRB3.2 are essential because of DRB3 genes function in innate immunity and its association with infectious diseases resistance or tolerance in cattle. The present study was aimed at assessing the level of genetic diversity of DRB3 in the exon 2 (BoLA-DRB3.2) region in African, American and Asian cattle breeds. Amplification of exon 2 in 174 cattle revealed 15 haplotypes. The breeds with the highest number of haplotypes were Brangus (10), Sokoto Gudali (10) and Dajal (9), while the lowest number of haplotypes were found in Holstein and Sahiwal with 4 haplotypes each. Medium Joining network obtained from haplotypic data showed that all haplotypes condensed around a centric area and each sequence (except in H-3, H-51 and H-106) representing almost a specific haplotype. The BoLA-DRB3.2 sequence analyses revealed a non-significant higher rate of non-synonymous (dN) compared to synonymous substitutions (dS). The ratio of dN/dS substitution across the breeds were observed to be greater than one suggesting that variation at the antigen-binding sites is under positive selection; thus increasing the chances of these breeds to respond to wide array of pathogenic attacks. An analysis of molecular variance revealed that 94.01 and 5.99% of the genetic variation was attributable to differences within and among populations, respectively. Generally, results obtained suggest that within breed genetic variation across breeds is higher than between breeds. This genetic information will be important for investigating the relationship between BoLADRB3.2 and diseases in various cattle breeds studied with attendant implication on designing breeding programs that will aim at selecting individual cattle that carry resistant alleles.

Genetic variation in N- and C-terminal regions of bovine DNAJA1 heat shock protein gene in African, Asian and American cattle

DNAJA1 or heat shock protein 40 (Hsp40) is associated with heat adaptation in various organisms. We amplified and sequenced a total of 1,142 bp of bovine Hsp40 gene representing the critical N-terminal (NTR) and C-terminal (CTR) regions in representative samples of African, Asian and American cattle breeds. Eleven and 9 different haplotypes were observed in the NTR in Asian and African breeds respectively while in American Brangus, only two mutations were observed resulting in two haplotypes. The CTR appears to be highly conserved between cattle and yak. In-silico functional analysis with PANTHER predicted putative deleterious functional impact of c.161 T>A; p. V54Q while alignment of bovine and human NTR-J domains revealed that p.Q19H, p.E20Q and p. E21X mutations occurred in helix 2 and p.V54Q missense mutation occurred in helix 3 respectively. The 124 bp insertion found in the yak DNAJA1 ortholog may have significant functional relevance warranting further investigation. Our results suggest that these genetic differences may be concomitant with population genetic history and possible functional consequences for climate adaptation in bovidae.

Computational genome-wide identification of heat shock protein genes in the bovine genome

Background: Heat shock proteins (HSPs) are molecular chaperones known to bind and sequester client proteins under stress. Methods: To identify and better understand some of these proteins, we carried out a computational genome-wide survey of the bovine genome. For this, HSP sequences from each subfamily (sHSP, HSP40, HSP70 and HSP90) were used to search the Pfam (Protein family) database, for identifying exact HSP domain sequences based on the hidden Markov model. ProtParam tool was used to compute potential physico-chemical parameters detectable from a protein sequence. Evolutionary trace (ET) method was used to extract evolutionarily functional residues of a homologous protein family. Results: We computationally identified 67 genes made up of 10, 43, 10 and 4 genes belonging to small HSP, HSP40, HSP70 and HSP90 families respectively. These genes were widely dispersed across the bovine genome, except in chromosomes 24, 26 and 27, which lack bovine HSP genes. We found an uncharacterized outer dense fiber ( ODF1) gene in cattle with an intact alpha crystallin domain, like other small HSPs. Physico-chemical characteristic of aliphatic index was higher in HSP70 and HSP90 gene families, compared to small HSP and HSP40. Grand average hydropathy showed that small HSP (sHSP), HSP40, HSP70 and HSP90 genes had negative values except for DNAJC22, a member of HSP40 gene family. The uniqueness of DNAJA3 and DNAJB13 among HSP40 members, based on multiple sequence alignment, evolutionary trace analysis and sequence identity dendrograms, suggests evolutionary distinct structural and functional features, with unique roles in substrate recognition and chaperone functions. The monophyletic pattern of the sequence identity dendrograms of cattle, human and mouse HSP sequences suggests functional similarities. Conclusions: Our computational results demonstrate the first-pass in-silico identification of heat shock proteins and calls for further investigation to better understand their functional roles and mechanisms in Bovidae.