An Integrated Bioinformatics Approach to Identify Network-Derived Hub Genes in Starving Zebrafish.

The present study was aimed at identifying causative hub genes within modules formed by co-expression and protein-protein interaction (PPI) networks, followed by Bayesian network (BN) construction in the liver transcriptome of starved zebrafish. To this end, the GSE11107 and GSE112272 datasets from the GEO databases were downloaded and meta-analyzed using the MetaDE package, an add-on R package. Differentially expressed genes (DEGs) were identified based upon expression intensity N(µ = 0.2, σ2 = 0.4). Reconstruction of BNs was performed by the bnlearn R package on genes within modules using STRINGdb and CEMiTool. ndufs5 (shared among PPI, BN and COEX), rps26, rpl10, sdhc (shared between PPI and BN), ndufa6, ndufa10, ndufb8 (shared between PPI and COEX), skp1, atp5h, ndufb10, rpl5b, zgc:193613, zgc:123327, zgc:123178, wu:fc58f10, zgc:111986, wu:fc37b12, taldo1, wu:fb62f08, zgc:64133 and acp5a (shared between COEX and BN) were identified as causative hub genes affecting gene expression in the liver of starving zebrafish. Future work will shed light on using integrative analyses of miRNA and DNA microarrays simultaneously, and performing in silico and experimental validation of these hub-causative (CST) genes affecting starvation in zebrafish.

Constraint-Based, Score-Based and Hybrid Algorithms to Construct Bayesian Gene Networks in the Bovine Transcriptome.

Bayesian gene networks are powerful for modelling causal relationships and incorporating prior knowledge for making inferences about relationships. We used three algorithms to construct Bayesian gene networks around genes expressed in the bovine uterus and compared the efficacies of the algorithms. Dataset GSE33030 from the Gene Expression Omnibus (GEO) repository was analyzed using different algorithms for hub gene expression due to the effect of progesterone on bovine endometrial tissue following conception. Six different algorithms (grow-shrink, max-min parent children, tabu search, hill-climbing, max-min hill-climbing and restricted maximum) were compared in three higher categories, including constraint-based, score-based and hybrid algorithms. Gene network parameters were estimated using the bnlearn bundle, which is a Bayesian network structure learning toolbox implemented in R. The results obtained indicated the tabu search algorithm identified the highest degree between genes (390), Markov blankets (25.64), neighborhood sizes (8.76) and branching factors (4.38). The results showed that the highest number of shared hub genes (e.g., proline dehydrogenase 1 (PRODH), Sam-pointed domain containing Ets transcription factor (SPDEF), monocyte-to-macrophage differentiation associated 2 (MMD2), semaphorin 3E (SEMA3E), solute carrier family 27 member 6 (SLC27A6) and actin gamma 2 (ACTG2)) was seen between the hybrid and the constraint-based algorithms, and these genes could be recommended as central to the GSE33030 data series. Functional annotation of the hub genes in uterine tissue during progesterone treatment in the pregnancy period showed that the predicted hub genes were involved in extracellular pathways, lipid and protein metabolism, protein structure and post-translational processes. The identified hub genes obtained by the score-based algorithms had a role in 2-arachidonoylglycerol and enzyme modulation. In conclusion, different algorithms and subsequent topological parameters were used to identify hub genes to better illuminate pathways acting in response to progesterone treatment in the bovine uterus, which should help with our understanding of gene regulatory networks in complex trait expression.

Evolutionary Pattern of Interferon Alpha Genes in Bovidae and Genetic Diversity of IFNAA in the Bovine Genome.

Interferons are secretory proteins induced in response to specific extracellular stimuli which stimulate intra- and intercellular networks for regulating innate and acquired immunity, resistance to viral infections, and normal and tumor cell survival and death. Type 1 interferons plays a major role in the CD8 T-cell response to viral infection. The genomic analysis carried out here for type I interferons within Bovidae family shows that cattle, bison, water buffalo, goat, and sheep (all Bovidae), have different number of genes of the different subtypes, with a large increase in the numbers, compared to human and mouse genomes. A phylogenetic analysis of the interferon alpha (IFNA) proteins in this group shows that the genes do not follow the evolutionary pattern of the species, but rather a cycle of duplications and deletions in the different species. In this study we also studied the genetic diversity of the bovine interferon alpha A (IFNAA), as an example of the IFNA genes in cattle, sequencing a fragment of the coding sequence in 18 breeds of cattle from Pakistan, Nigeria and USA. Similarity analysis allowed the allocation of sequences into 22 haplotypes. Bhagnari, Brangus, Sokoto Gudali, and White Fulani, had the highest number of haplotypes, while Angus, Hereford and Nari Master had the least. However, when analyzed by the average haplotype count, Angus, Bhagnari, Hereford, Holstein, Muturu showed the highest values, while Cholistani, Lohani, and Nari Master showed the lowest values. Haplotype 4 was found in the highest number of individuals (74), and in 15 breeds. Sequences for yak, bison, and water buffalo, were included within the bovine haplotypes. Medium Joining network showed that the sequences could be divided into 4 groups: one with highly similar haplotypes containing mostly Asian and African breeds, one with almost all of the Bos taurus American breeds, one mid-diverse group with mostly Asian and African sequences, and one group with highly divergent haplotypes with five N'Dama sequences and one from each of White Fulani, Dhanni, Tharparkar, and Bhagnari. The large genetic diversity found in IFNAA could be a very good indication of the genetic variation among the different genes of IFNA and could be an adaptation for these species in response to viral challenges they face.

RNA-seq profiling of skin in temperate and tropical cattle.

Skin is a major thermoregulatory organ in the body controlling homeothermy, a critical function for climate adaptation. We compared genes expressed between tropical- and temperate-adapted cattle to better understand genes involved in climate adaptation and hence thermoregulation. We profiled the skin of representative tropical and temperate cattle using RNA-seq. A total of 214,754,759 reads were generated and assembled into 72,993,478 reads and were mapped to unique regions in the bovine genome. Gene coverage of unique regions of the reference genome showed that of 24,616 genes, only 13,130 genes (53.34%) displayed more than one count per million reads for at least two libraries and were considered suitable for downstream analyses. Our results revealed that of 255 genes expressed differentially, 98 genes were upregulated in tropically-adapted White Fulani (WF; Bos indicus) and 157 genes were down regulated in WF compared to Angus, AG (Bos taurus). Fifteen pathways were identified from the differential gene sets through gene ontology and pathway analyses. These include the significantly enriched melanin metabolic process, proteinaceous extracellular matrix, inflammatory response, defense response, calcium ion binding and response to wounding. Quantitative PCR was used to validate six representative genes which are associated with skin thermoregulation and epithelia dysfunction (mean correlation 0.92; p < 0.001). Our results contribute to identifying genes and understanding molecular mechanisms of skin thermoregulation that may influence strategic genomic selection in cattle to withstand climate adaptation, microbial invasion and mechanical damage.

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.

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 gene's 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.

Sequence variation of necdin gene in Bovidae.

BACKGROUND: Necdin (NDN), a member of the melanoma antigen family showing imprinted pattern of expression, has been implicated as causing Prader-Willi symptoms, and known to participate in cellular growth, cellular migration and differentiation. The region where NDN is located has been associated to QTLs affecting reproduction and early growth in cattle, but location and functional analysis of the molecular mechanisms have not been established. METHODS: Here we report the sequence variation of the entire coding sequence from 72 samples of cattle, yak, buffalo, goat and sheep, and discuss its variation in Bovidae. Median-joining network analysis was used to analyze the variation found in the species. Synonymous and non-synonymous substitution rates were determined for the analysis of all the polymorphic sites. Phylogenetic analysis were carried out among the species of Bovidae to reconstruct their relationships. RESULTS: From the phylogenetic analysis with the consensus sequences of the studied Bovidae species, we found that only 11 of the 26 nucleotide changes that differentiate them produced amino acid changes. All the SNPs found in the cattle breeds were novel and showed similar percentages of nucleotides with non-synonymous substitutions at the N-terminal, MHD and C-terminal (12.3, 12.8 and 12.5%, respectively), and were much higher than the percentage of synonymous substitutions (2.5, 2.6 and 4.9%, respectively). Three mutations in cattle and one in sheep, detected in heterozygous individuals were predicted to be deleterious. Additionally, the analysis of the biochemical characteristics in the most common form of the proteins in each species show very little difference in molecular weight, pI, net charge, instability index, aliphatic index and GRAVY (Table 4) in the Bovidae species, except for sheep, which had a higher molecular weight, instability index and GRAVY. CONCLUSIONS: There is sufficient variation in this gene within and among the studied species, and because NDN carry key functions in the organism, it can have effects in economically important traits in the production of these species. NDN sequence is phylogenetically informative in this group, thus we propose this gene as a phylogenetic marker to study the evolution and conservation in Bovidae.

Conservation of Repeats at the Mammalian KCNQ1OT1-CDKN1C Region Suggests a Role in Genomic Imprinting.

KCNQ1OT1 is located in the region with the highest number of genes showing genomic imprinting, but the mechanisms controlling the genes under its influence have not been fully elucidated. Therefore, we conducted a comparative analysis of the KCNQ1/KCNQ1OT1-CDKN1C region to study its conservation across the best assembled eutherian mammalian genomes sequenced to date and analyzed potential elements that may be implicated in the control of genomic imprinting in this region. The genomic features in these regions from human, mouse, cattle, and dog show a higher number of genes and CpG islands (detected using cpgplot from EMBOSS), but lower number of repetitive elements (including short interspersed nuclear elements and long interspersed nuclear elements), compared with their whole chromosomes (detected by RepeatMasker). The KCNQ1OT1-CDKN1C region contains the highest number of conserved noncoding sequences (CNS) among mammals, where we found 16 regions containing about 38 different highly conserved repetitive elements (using mVista), such as LINE1 elements: L1M4, L1MB7, HAL1, L1M4a, L1Med, and an LTR element: MLT1H. From these elements, we found 74 CNS showing high sequence identity (>70%) between human, cattle, and mouse, from which we identified 13 motifs (using Multiple Em for Motif Elicitation/Motif Alignment and Search Tool) with a significant probability of occurrence, 3 of which were the most frequent and were used to find transcription factor-binding sites. We detected several transcription factors (using JASPAR suite) from the families SOX, FOX, and GATA. A phylogenetic analysis of these CNS from human, marmoset, mouse, rat, cattle, dog, horse, and elephant shows branches with high levels of support and very similar phylogenetic relationships among these groups, confirming previous reports. Our results suggest that functional DNA elements identified by comparative genomics in a region densely populated with imprinted mammalian genes may be related to the regulation of imprinted gene expression.

Genetic diversity among Trypanosoma vivax strains detected in naturally infected cattle in Nigeria based on ITS1 of rDNA and diagnostic antigen gene sequences.

Trypanosoma vivax (sub-genus Duttonella) is largely responsible for non profitable livestock production in sub-Sahara Africa. In Nigeria, no study has addressed the molecular characteristic of T. vivax except Y486. Hence, we characterized and assessed the genetic diversity among T. vivax detected in naturally infected cattle in Nigeria using internal transcribed spacer 1 (ITS1) of ribosoma DNA (rDNA) and diagnostic antigen gene (DAG) sequences. The length of ITS1 and DAG sequences range from 215-220 to 257-338 bp, respectively and the mean G-C contents were 60 and 61.5 %. Homology search revealed 93-99 and 95-100 % homologies to T. vivax DAG and ITS1 sequences from GenBank. Aligned sequences revealed both ITS1 rDNA and DAG to be less polymorphic but DAG sequences of the Y486 strain and its clone showed marked variation from autochthonous strains. Phylogenetic analysis yielded tree that grouped T. vivax ITS1rDNA gene and DAG sequences into two main clades each. Considering the ITI1 rDNA sequences, clade A contained autochthonous T. vivax within which the South American sequences clustered, clade B contained the sequences of T. vivax from East Africa. Analysis of DAG revealed that the clade A contains autochthonous T. vivax sequences but clade B contained the Y486 and its clones. In conclusion, the diagnostic antigen gene sequences of the T. vivax detected in this study may have undergone considerable gene recombination through time and suggests that more than one strain of T. vivax exist among cattle population in Nigeria.

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.

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.

Genetic Variation of Goat Interferon Regulatory Factor 3 Gene and Its Implication in Goat Evolution.

The immune systems are fundamentally vital for evolution and survival of species; as such, selection patterns in innate immune loci are of special interest in molecular evolutionary research. The interferon regulatory factor (IRF) gene family control many different aspects of the innate and adaptive immune responses in vertebrates. Among these, IRF3 is known to take active part in very many biological processes. We assembled and evaluated 1356 base pairs of the IRF3 gene coding region in domesticated goats from Africa (Nigeria, Ethiopia and South Africa) and Asia (Iran and China) and the wild goat (Capra aegagrus). Five segregating sites with θ value of 0.0009 for this gene demonstrated a low diversity across the goats' populations. Fu and Li tests were significantly positive but Tajima's D test was significantly negative, suggesting its deviation from neutrality. Neighbor joining tree of IRF3 gene in domesticated goats, wild goat and sheep showed that all domesticated goats have a closer relationship than with the wild goat and sheep. Maximum likelihood tree of the gene showed that different domesticated goats share a common ancestor and suggest single origin. Four unique haplotypes were observed across all the sequences, of which, one was particularly common to African goats (MOCH-K14-0425, Poitou and WAD). In assessing the evolution mode of the gene, we found that the codon model dN/dS ratio for all goats was greater than one. Phylogenetic Analysis by Maximum Likelihood (PAML) gave a ω0 (dN/dS) value of 0.067 with LnL value of -6900.3 for the first Model (M1) while ω2 = 1.667 in model M2 with LnL value of -6900.3 with positive selection inferred in 3 codon sites. Mechanistic empirical combination (MEC) model for evaluating adaptive selection pressure on particular codons also confirmed adaptive selection pressure in three codons (207, 358 and 408) in IRF3 gene. Positive diversifying selection inferred with recent evolutionary changes in domesticated goat IRF3 led us to conclude that the gene evolution may have been influenced by domestication processes in goats.

High density genome wide genotyping-by-sequencing and association identifies common and low frequency SNPs, and novel candidate genes influencing cow milk traits.

High-throughput sequencing technologies have increased the ability to detect sequence variations for complex trait improvement. A high throughput genome wide genotyping-by-sequencing (GBS) method was used to generate 515,787 single nucleotide polymorphisms (SNPs), from which 76,355 SNPs with call rates >85% and minor allele frequency ≥1.5% were used in genome wide association study (GWAS) of 44 milk traits in 1,246 Canadian Holstein cows. GWAS was accomplished with a mixed linear model procedure implementing the additive and dominant models. A strong signal within the centromeric region of bovine chromosome 14 was associated with test day fat percentage. Several SNPs were associated with eicosapentaenoic acid, docosapentaenoic acid, arachidonic acid, CLA:9c11t and gamma linolenic acid. Most of the significant SNPs for 44 traits studied are novel and located in intergenic regions or introns of genes. Novel potential candidate genes for milk traits or mammary gland functions include ERCC6, TONSL, NPAS2, ACER3, ITGB4, GGT6, ACOX3, MECR, ADAM12, ACHE, LRRC14, FUK, NPRL3, EVL, SLCO3A1, PSMA4, FTO, ADCK5, PP1R16A and TEP1. Our study further demonstrates the utility of the GBS approach for identifying population-specific SNPs for use in improvement of complex dairy traits.

Phylogeny of Trypanosoma brucei and Trypanosoma evansi in naturally infected cattle in Nigeria by analysis of repetitive and ribosomal DNA sequences.

In continuing efforts to better understand the genetics of bovine trypanosomosis, we assessed genetic diversity of Trypanosoma brucei and Trypanosoma evansi in naturally infected Nigerian cattle using repetitive DNA and internal transcribed spacer 1 of rDNA sequences and compared these sequences to species from other countries. The length of repetitive DNA sequences in both species ranged from 161 to 244 bp and 239 to 240 bp for T. brucei and T. evansi, respectively, while the ITS1 rDNA sequences length range from 299 to 364 bp. The mean GC content of ITS1 rDNA sequences was 33.57 %, and that of repetitive sequences were 39.9 and 31.1 % for T. brucei and T. evansi, respectively. Result from sequence alignment revealed both T. brucei and T. evansi repetitive DNA sequences to be more polymorphic than ITS1 rDNA sequences, with moderate points of deletion and insertions. T. brucei separated into two clades when subjected to phylogenetic analysis. T. evansi repetitive DNA sequences clustered tightly within the T. brucei clade while the ITS1 rDNA sequences of T. brucei were clearly separated from T. theileri and T. vivax individually used as outgroups. This study suggest that ITS1 rDNA sequences may not be suitable for phylogenetic differentiation of the Trypanozoon group and also suggest that T. evansi may be a phenotypic variant of T. brucei which may have potential implications in designing prevention and therapeutic strategies.

Extensive genomic variability of knops blood group polymorphisms is associated with sickle cell disease in Africa.

Sickle cell disease (SCD) is a multisystem disorder characterized by chronic hemolytic anemia, vaso-occlusive crises, and marked variability in disease severity. Patients require transfusions to manage disease complications, with complements, directed by complement regulatory genes (CR1) and its polymorphisms, implicated in the development of alloantibodies. We hypothesize that CR1 polymorphisms affect complement regulation and function, leading to adverse outcome in SCD. To this end, we determined the genomic diversity of complement regulatory genes by examining single nucleotide polymorphisms associated with Knops blood group antigens. Genomic DNA samples from 130 SCD cases and 356 control Africans, 331 SCD cases and 497 control African Americans, and 254 Caucasians were obtained and analyzed, utilizing a PCR-RFLP (polymerase chain reaction-restriction fragment length polymorphism) assay. Analyzing for ethnic diversity, we found significant differences in the genotypic and allelic frequencies of Sl1/Sl2 (rs17047661) and McCa/b (rs17047660) polymorphisms between Africans, African Americans, and Caucasians (P < 0.05). The homozygote mutant variants had significantly higher frequencies in Africans and African Americans but were insignificant in Caucasians (80.2% and 59.6% vs 5.9% for Sl1/2; and 36% and 24% vs 1.8% for McCa/b). With SCD, we did not detect any difference among cases and controls either in Africa or in the United States. However, we found significant difference in genotypic (P < 0.0001) and allelic frequencies (P < 0.0001) of Sl1/Sl2 (rs17047661) and McCa/b (rs17047660) polymorphisms between SCD groups from Africa and the United States. There was no difference in haplotype frequencies of these polymorphisms among or between groups. The higher frequency of CR1 homozygote mutant variants in Africa but not United States indicates a potential pathogenic role, possibly associated with complicated disease pathophysiology in the former and potentially protective in the latter. The difference between sickle cell groups suggests potential genetic drift or founder effect imposed on the disease in the United States, but not in Africa, and a possible confirmation of the ancestral susceptibility hypothesis. The lower haplotype frequencies among sickle cell and control populations in the United States may be due to the admixture and the dilution of African genetic ancestry in the African American population.

Interethnic diversity of the CD209 (rs4804803) gene promoter polymorphism in African but not American sickle cell disease.

Elucidating the genomic diversity of CD209 gene promoter polymorphism could assist in clarifying disease pathophysiology as well as contribution to co-morbidities. CD209 gene promoter polymorphism has been shown to be associated with susceptibility to infection. We hypothesize that CD209 mutant variants occur at a higher frequency among Africans and in sickle cell disease. We analyzed the frequency of the CD209 gene (rs4804803) in healthy control and sickle cell disease (SCD) populations and determined association with disease. Genomic DNA was extracted from blood samples collected from 145 SCD and 231 control Africans (from Mali), 331 SCD and 379 control African Americans and 159 Caucasians. Comparative analysis among and between groups was carried out by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Per ethnic diversification, we found significant disparity in genotypic (23.4% versus 16.9% versus 3.2%) and allelic frequencies (48.7% versus 42.1% versus 19.8%) of the homozygote mutant variant of the CD209 (snp 309A/G) gene promoter between Africans, African Americans and Caucasians respectively. Comparative evaluation between disease and control groups reveal a significant difference in genotypic (10.4% versus 23.4%; p = 0.002) and allelic frequencies (39.7% versus 48.7%; p = 0.02) of the homozygote mutant variant in African SCD and healthy controls respectively, an observation that is completely absent among Americans. Comparing disease groups, we found no difference in the genotypic (p = 0.19) or allelic (p = 0.72) frequencies of CD209 homozygote mutant variant between Africans and Americans with sickle cell disease. The higher frequency of CD209 homozygote mutant variants in the African control group reveals a potential impairment of the capacity to mount an immune response to infectious diseases, and possibly delineate susceptibility to or severity of infectious co-morbidities within and between groups.

Endothelin-1 but not Endothelial Nitric Oxide Synthase Gene Polymorphism is Associated with Sickle Cell Disease in Africa.

Sickle cell disease shows marked variability in severity and pathophysiology among individuals, probably linked to differential expression of various adhesion molecules. In this study, we investigated the differential distribution, genomic diversity and haplotype frequency of endothelial nitric oxide synthase (eNOS) and endothelin-1 (ET-1) polymorphisms, recently implicated as important in modification of disease severity. One hundred and forty five sickle cell disease patients (HbSS) and 244 adult and pediatric controls, without sickle cell disease (HbAA), were recruited from Mali. Genotypic analysis of the functionally significant eNOS variants (T786C, G894T and intron 4) and endothelin-1 (G5665T) was carried out with a polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay. Our results show that the wild type alleles are the most frequent for all eNOS variants between cases and controls. Allelic and genotypic frequencies of eNOS polymorphic groups are not significantly different between cases and controls (P > 0.05). In addition, there is no association between eNOS variants and sickle cell disease, contrary to published reports. On the other hand, we report that endothelin-1 (G5665T) mutant variant had the lowest allelic frequency, and is significantly associated with sickle cell disease in Africa (P < 0.05). Similarly, haplotype frequencies were the same between cases and controls, except for the haplotype combining all mutant variants (T, C, 4a; P = 0.01). eNOS polymorphic variants are less frequent, with no significance with sickle cell disease in Africa. On the other hand, endothelin-1 is associated with sickle cell disease, and has the capacity to redefine pathophysiology and possibly serve as modulator of disease phenotype.

A novel TaqI polymorphism in the coding region of the ovine TNXB gene in the MHC class III region: morphostructural and physiological influences.

The tenascin-XB (TNXB) gene has antiadhesive effects, functions in matrix maturation in connective tissues, and localizes to the major histocompatibility complex class III region. We hypothesized that it may influence adaptive physiological response through an effect on blood vessel function. We identified a novel g.1324 A→G polymorphism at a TaqI recognition site in a 454 bp fragment of ovine TNXB and genotyped it in 150 Nigerian sheep using PCR-RFLP. The missense mutation changes glutamic acid (GAA) to glycine (GGA). Among SNP genotypes, significant differences (P < 0.05) were observed in body weight and fore cannon bone length. Interaction effects of breed, SNP genotype, and geographic location had a significant effect (P < 0.05) on chest girth. The SNP genotype was significantly (P < 0.05) associated with physiological traits of pulse rate and skin temperature. The observed effect of this novel polymorphism may be mediated through its role in connective tissue biology, requiring further association and functional studies.