Aberrant low expression level of bovine beta-lactoglobulin is associated with a C to A transversion in the BLG promoter region

Beta-lactoglobulin (beta-LG) is the major whey protein in cow's milk. It is well established that the predominant 2 genetic variants, beta-LG A and B, are differentially expressed. Extensive investigation of the genetic variation in the promoter region of the BLG gene revealed the existence of specific haplotypes associated with the A and B variants, respectively. However, the genetic basis for the differential expression of BLG A and B alleles is still elusive. We have previously reported a quantitative beta-LG B variant, characterized by a very low beta-LG protein expression level. Here, we report that the corresponding BLG allele (BLG B*) shows a correspondingly low mRNA expression level. Comparative DNA sequencing of 7,670 bp of the BLG B* allele and the established BLG B allele revealed a unique difference of a C to A transversion at position 215 bp upstream of the translation initiation site (g.-215C>A). This mutation segregated perfectly with the differential phenotypic expression in a paternal half-sib family and could be confirmed in 2 independent cases. The sequence of the BLG B allele in the region of the mutation is highly conserved among 4 related ruminant species. The site of the mutation corresponds to a putative consensus-binding sequence for the transcription factors c-Rel and Elk-1 as predicted by searching the TRANSFAC database. The beta-LG B* site might be relevant in the natural production of milk of low beta-LG content.

Neighboring base composition and transversion/transition bias in a comparison of rice and maize chloroplast noncoding regions.

The correspondence between the transversion/transition ratio and the neighboring base composition in chloroplast DNA is examined. For 18 noncoding regions of the chloroplast genome, alignments between rice (Oryza sativa) and maize (Zea mays) were generated by two different methods. Difficulties of aligning noncoding DNA are discussed, and the alignments are analyzed in a manner that reduces alignment artifacts. Sequence divergence is < 10%, so multiple substitutions at a site are assumed to be rare. Observed substitutions were analyzed with respect to the A+T content of the two immediately flanking bases. It is shown that as this content increases, the proportion of transversions also increases. When both the 5'- and 3'-flanking nucleotides are G or C (A+T content of 0), only 25% of the observed substitutions are transversions. However, when both the 5'- and 3'-flanking nucleotides are A or T (A+T content of 2), 57% of the observed substitutions are transversions. Therefore, the influence of flanking base composition on substitutions, previously reported for a single noncoding region, is a general feature of the chloroplast genome.

Effect of a thrombin receptor (protease-activated receptor 1, PAR-1) gene polymorphism in chronic hepatitis C liver fibrosis

Background and Aim: Tissue injury leads to activation of coagulation and generation of thrombin. Inhibition of thrombin receptor protease-activated receptor 1 (PAR-1) has been shown to reduce liver fibrosis in animals. This study aimed to evaluate the effect of PAR-1 gene polymorphism on rate of liver fibrosis (RF) in chronic hepatitis C. Methods: Polymorphisms studied: C > T transition 1426 bp upstream of translation start site (-1426C/T), 13 bp repeat of preceding -506 5`-CGGCCGCGGGAAG-3` sequence (-506I/D), and A > T transversion in intervening sequence (IVS) 14 bp upstream of exon-2 start site (IVS-14A/T). A total of 287 European and 90 Brazilian patients were studied. Results: 1426C/T polymorphism: There was a trend to higher RF in patients with the TT genotype (P = 0.06) and an association between genotype CC and slow fibrosis (P = 0.03) in Europeans. In males, RF was significantly higher in those with the TT genotype compared to CT (P = 0.003) and CC (P = 0.007). There was a significant association between TT and fast fibrosis (P = 0.04). This was confirmed in an independent cohort of Brazilians where RF was higher in TT than in CC (P = 0.03). Analysis of -506I/D showed no difference in RF and distribution of slow/fast fibrosis among different genotypes in both populations. Analysis of IVS-14A/T showed no difference between genotypes. Conclusion: In conclusion, these findings suggest that PAR-1 receptor polymorphisms influence the progression of liver fibrosis.

Are there two species of Pygmy Shrews (Sorex)? Revisiting the question using DNA sequence data

Pygmy Shrews in North America have variously been considered to be one species (Sorex hoyi) or two species (S. hoyi and S. thompsoni). Currently, only S. hoyi is recognized. In this study, we examine mitochondrial DNA sequence data for the cytochrome b gene to evaluate the level of differentiation and phylogeographic relationships among eleven samples of Pygmy Shrews from across Canada. Pygmy Shrews from eastern Canada (i.e., Ontario, Quebec, New Brunswick, Nova Scotia, and Prince Edward Island) are distinct from Pygmy Shrews from western Canada (Alberta, Yukon) and Alaska. The average level of sequence divergence between these clades (3.3%) falls within the range of values for other recognized pairs of sister species of shrews. A molecular clock based on third position transversion substitutions suggests that these two lineages diverged between 0.44 and 1.67 million years ago. These molecular phylogenetic data. combined with a reinterpretation of previously published morphological data, are suggestive of separate species status for S. hoyi and S. thompsoni as has been previously argued by others. Further analysis of specimens from geographically intermediate areas (e.g., Manitoba. northern Ontario) is required to determine if there is secondary contact and/or introgression between these two putative species.

Dyssegmental dysplasia, Silverman-Handmaker type: prenatal ultrasound findings and molecular analysis.

OBJECTIVES: The objective of this study is to describe the prenatal sonographic features and the results of DNA analysis on three fetuses with dyssegmental dysplasia, Silverman-Handmaker type (DD-SH). METHODS: A retrospective review of three fetuses with confirmed DD-SH was conducted. The fetal ultrasound findings, the radiological characteristics, and the results of the mutation analysis of the heparan sulphate perlecan gene 2 (HSPG2) were reviewed. RESULTS: There were three cases in two families with DD-SH diagnosed prenatally. The main prenatal ultrasound and the radiological features of DD-SH were severe limb shortening and vertebral segmentation and fusion defects (anisospondyly). The DNA analysis of the HSPG2 gene showed that the two affected fetuses in a nonconsanguineous family had a compound heterozygote for the c.646G > T transversion in exon 7 and a c.5788C > T transition in exon 46. The fetus born to the consanguineous couple had a homozygous mutation c.1356-27_1507 + 59del. CONCLUSION: DD-SH can be diagnosed prenatally using fetal ultrasound as early as 13 weeks. Xrays and DNA analysis of the HSPG2 gene are important for the confirmation of the diagnosis and for the preimplantation and prenatal diagnosis in pregnancies at risk. © 2013 John Wiley & Sons, Ltd.

The maritime shrew, Sorex maritimensis (Insectivora : Soricidae): a newly recognized Canadian endemic.

Previous morphological and cytological analyses have suggested that the arctic shrew (Sorex arcticus) as currently recognized may be two distinct species. Specifically, those studies demonstrated considerable differentiation between the putative subspecies S. a. maritimensis and one or both of the other two subspecies, S. a. arcticus and S. a. laricorum. Phylogenetic analysis of 546 base pairs of cytochrome b sequence data from 10 arctic shrews from across Canada indicates that maritimensis is the sister-group to arcticus + laricorum. Furthermore, there is considerable genetic divergence between maritimensis and the other two putative subspecies (similar to8-9%; Kimura's two-parameter distance). Given that maritimensis and arcticus + laricorum appear to be reciprocally monophyletic clades with considerable genetic divergence (i.e., greater than that between other recognized pairs of sister-species within the S. araneus-arcticus group), we propose that S. maritimensis be recognized as a distinct species. The proportion of third-position transversion substitutions between S. arcticus and S. maritimensis suggests that these two species shared a common ancestor approximately 2.4 million years ago.

Prenatal presentation and postnatal evolution of a patient with Jansen metaphyseal dysplasia with a novel missense mutation in PTH1R.

Wave-shaped ribs were detected at prenatal ultrasound in a 20(+1) week female fetus. At birth, skeletal radiographs showed marked hypomineralization and suggested hypophosphatasia. However, elevated blood calcium and alkaline phosphatase excluded hypophosphatasia and raised the possibility of Jansen metaphyseal dysplasia. Molecular analysis of the PTH/PTHrP receptor gene (PTH1R) showed heterozygosity for a previously undescribed transversion variant (c.1373T>A), which predicts p.Ile458Lys. In vitro evaluation of wild type and mutant PTH/PTHrP receptors supported the pathogenic role of the p.Ile458Lys substitution, and confirmed the diagnosis of Jansen metaphyseal dysplasia. This disorder may present prenatally with wavy ribs and in the newborn with hypomineralization, and may therefore be confused with hypophosphatasia. The mottled metaphyseal lesions typically associated with this disease appear only in childhood.

Development and Validation of Single Nucleotide Polymorphisms (SNPs) Markers from Two Transcriptome 454-Runs of Turbot (Scophthalmus maximus) Using High-Throughput Genotyping

The turbot (Scophthalmus maximus) is a commercially valuable flatfish and one of the most promising aquaculture species in Europe. Two transcriptome 454-pyrosequencing runs were used in order to detect Single Nucleotide Polymorphisms (SNPs) in genesrelated to immune response and gonad differentiation. A total of 866 true SNPs were detected in 140 different contigs representing 262,093 bp as a whole. Only one true SNP was analyzed in each contig. One hundred and thirteen SNPs out of the 140 analyzed were feasible (genotyped), while Ш were polymorphic in a wild population. Transition/transversion ratio (1.354) was similar to that observed in other fish studies. Unbiased gene diversity (He) estimates ranged from 0.060 to 0.510 (mean = 0.351), minimum allele frequency (MAF) from 0.030 to 0.500 (mean = 0.259) and all loci were in Hardy-Weinberg equilibrium after Bonferroni correction. A large number of SNPs (49) were located in the coding region, 33 representing synonymous and 16 non-synonymous changes. Most SNP-containing genes were related to immune response and gonad differentiation processes, and could be candidates for functional changes leading to phenotypic changes. These markers will be useful for population screening to look for adaptive variation in wild and domestic turbot

KCM, a general framework for codon-based model of evolution

With the advancement of high-throughput sequencing and dramatic increase of available genetic data, statistical modeling has become an essential part in the field of molecular evolution. Statistical modeling results in many interesting discoveries in the field, from detection of highly conserved or diverse regions in a genome to phylogenetic inference of species evolutionary history Among different types of genome sequences, protein coding regions are particularly interesting due to their impact on proteins. The building blocks of proteins, i.e. amino acids, are coded by triples of nucleotides, known as codons. Accordingly, studying the evolution of codons leads to fundamental understanding of how proteins function and evolve. The current codon models can be classified into three principal groups: mechanistic codon models, empirical codon models and hybrid ones. The mechanistic models grasp particular attention due to clarity of their underlying biological assumptions and parameters. However, they suffer from simplified assumptions that are required to overcome the burden of computational complexity. The main assumptions applied to the current mechanistic codon models are (a) double and triple substitutions of nucleotides within codons are negligible, (b) there is no mutation variation among nucleotides of a single codon and (c) assuming HKY nucleotide model is sufficient to capture essence of transition- transversion rates at nucleotide level. In this thesis, I develop a framework of mechanistic codon models, named KCM-based model family framework, based on holding or relaxing the mentioned assumptions. Accordingly, eight different models are proposed from eight combinations of holding or relaxing the assumptions from the simplest one that holds all the assumptions to the most general one that relaxes all of them. The models derived from the proposed framework allow me to investigate the biological plausibility of the three simplified assumptions on real data sets as well as finding the best model that is aligned with the underlying characteristics of the data sets. -- Avec l'avancement de séquençage à haut débit et l'augmentation dramatique des données géné¬tiques disponibles, la modélisation statistique est devenue un élément essentiel dans le domaine dé l'évolution moléculaire. Les résultats de la modélisation statistique dans de nombreuses découvertes intéressantes dans le domaine de la détection, de régions hautement conservées ou diverses dans un génome de l'inférence phylogénétique des espèces histoire évolutive. Parmi les différents types de séquences du génome, les régions codantes de protéines sont particulièrement intéressants en raison de leur impact sur les protéines. Les blocs de construction des protéines, à savoir les acides aminés, sont codés par des triplets de nucléotides, appelés codons. Par conséquent, l'étude de l'évolution des codons mène à la compréhension fondamentale de la façon dont les protéines fonctionnent et évoluent. Les modèles de codons actuels peuvent être classés en trois groupes principaux : les modèles de codons mécanistes, les modèles de codons empiriques et les hybrides. Les modèles mécanistes saisir une attention particulière en raison de la clarté de leurs hypothèses et les paramètres biologiques sous-jacents. Cependant, ils souffrent d'hypothèses simplificatrices qui permettent de surmonter le fardeau de la complexité des calculs. Les principales hypothèses retenues pour les modèles actuels de codons mécanistes sont : a) substitutions doubles et triples de nucleotides dans les codons sont négligeables, b) il n'y a pas de variation de la mutation chez les nucléotides d'un codon unique, et c) en supposant modèle nucléotidique HKY est suffisant pour capturer l'essence de taux de transition transversion au niveau nucléotidique. Dans cette thèse, je poursuis deux objectifs principaux. Le premier objectif est de développer un cadre de modèles de codons mécanistes, nommé cadre KCM-based model family, sur la base de la détention ou de l'assouplissement des hypothèses mentionnées. En conséquence, huit modèles différents sont proposés à partir de huit combinaisons de la détention ou l'assouplissement des hypothèses de la plus simple qui détient toutes les hypothèses à la plus générale qui détend tous. Les modèles dérivés du cadre proposé nous permettent d'enquêter sur la plausibilité biologique des trois hypothèses simplificatrices sur des données réelles ainsi que de trouver le meilleur modèle qui est aligné avec les caractéristiques sous-jacentes des jeux de données. Nos expériences montrent que, dans aucun des jeux de données réelles, tenant les trois hypothèses mentionnées est réaliste. Cela signifie en utilisant des modèles simples qui détiennent ces hypothèses peuvent être trompeuses et les résultats de l'estimation inexacte des paramètres. Le deuxième objectif est de développer un modèle mécaniste de codon généralisée qui détend les trois hypothèses simplificatrices, tandis que d'informatique efficace, en utilisant une opération de matrice appelée produit de Kronecker. Nos expériences montrent que sur un jeux de données choisis au hasard, le modèle proposé de codon mécaniste généralisée surpasse autre modèle de codon par rapport à AICc métrique dans environ la moitié des ensembles de données. En outre, je montre à travers plusieurs expériences que le modèle général proposé est biologiquement plausible.

MMP13 mutations are the cause of recessive metaphyseal dysplasia, Spahr type.

Metaphyseal dysplasia, Spahr type (MDST; OMIM 250400) was described in 1961 based on the observation of four children in one family who had rickets-like metaphyseal changes but normal blood chemistry and moderate short stature. Its molecular basis and nosologic status remained unknown. We followed up on those individuals and diagnosed the disorder in an additional member of the family. We used exome sequencing to ascertain the underlying mutation and explored its consequences on three-dimensional models of the affected protein. The MDST phenotype is associated with moderate short stature and knee pain in adults, while extra-skeletal complications are not observed. The sequencing showed that MDST segregated with a c.619T>G single nucleotide transversion in MMP13. The predicted non-conservative amino acid substitution, p.Trp207Gly, disrupts a crucial hydrogen bond in the calcium-binding region of the catalytic domain of the matrix metalloproteinase, MMP13. The MDST phenotype is associated with recessive MMP13 mutations, confirming the importance of this metalloproteinase in the metaphyseal growth plate. Dominant MMP13 mutations have been associated with metaphyseal anadysplasia (OMIM 602111), while a single child homozygous for a MMP13 mutation had been previously diagnosed as "recessive metaphyseal anadysplasia," that we conclude is the same nosologic entity as MDST. Molecular confirmation of MDST allows distinction of it from dominant conditions (e.g., metaphyseal dysplasia, Schmid type; OMIM # 156500) and from more severe multi-system conditions (such as cartilage-hair hypoplasia; OMIM # 250250) and to give precise recurrence risks and prognosis. © 2014 Wiley Periodicals, Inc.

A genotypic mutation system measuring mutations in restriction recognition sequences.

The RFLP/PCR approach (restriction fragment length polymorphism/polymerase chain reaction) to genotypic mutation analysis described here measures mutations in restriction recognition sequences. Wild-type DNA is restricted before the resistant, mutated sequences are amplified by PCR and cloned. We tested the capacity of this experimental design to isolate a few copies of a mutated sequence of the human c-Ha-ras1 gene from a large excess of wild-type DNA. For this purpose we constructed a 272 bp fragment with 2 mutations in the PvuII recognition sequence 1727-1732 and studied the rescue by RFLP/PCR of a few copies of this 'PvuII mutant standard'. Following amplification with Taq-polymerase and cloning into lambda gt10, plaques containing wild-type sequence, PvuII mutant standard or Taq-polymerase induced bp changes were quantitated by hybridization with specific oligonucleotide probes. Our results indicate that 10 PvuII mutant standard copies can be rescued from 10(8) to 10(9) wild-type sequences. Taq polymerase errors originating from unrestricted, residual wild-type DNA were sequence dependent and consisted mostly of transversions originating at G.C bp. In contrast to a doubly mutated 'standard' the capacity to rescue single bp mutations by RFLP/PCR is limited by Taq-polymerase errors. Therefore, we assessed the capacity of our protocol to isolate a G to T transversion mutation at base pair 1698 of the MspI-site 1695-1698 of the c-Ha-ras1 gene from excess wild-type ras1 DNA. We found that 100 copies of the mutated ras1 fragment could be readily rescued from 10(8) copies of wild-type DNA.

Molecular mechanism involved in the mutagenicity induced by Aflatoxin B1

The aflatoxin B1 (AFB1) is a mycotoxin that has been identified as the most potent hepatocarcinogen. The metabolite resulting from detoxification process of AFB1 in liver, has the ability to react with the genomic DNA, generating AFB1-DNA adducts; during DNA replication process, this adduct induced the G:C→T:A transversion. Polymorphism in genes encoding for enzymes involved in the activation and detoxification of AFB1 and DNA repair enzymes has been associated with the risk of hepatocellular carcinoma (HCC) development. Additionally, in populations of high exposure to aflatoxin and high prevalence of hepatitis B virus (HBV) infection, has been demonstrated a synergism between these two risk factors for the development of HCC.

Patterns of nucleotide substitution in angiosperm cpDNA trnL (UAA)-trnF (GAA) regions

Patterns of substitution in chloroplast encoded trnL_F regions were compared between species of Actaea (Ranunculales), Digitalis (Scrophulariales), Drosera (Caryophyllales), Panicoideae (Poales), the small chromosome species clade of Pelargonium (Geraniales), each representing a different order of flowering plants, and Huperzia (Lycopodiales). In total, the study included 265 taxa, each with > 900-bp sequences, totaling 0.24 Mb. Both pairwise and phylogeny-based comparisons were used to assess nucleotide substitution patterns. In all six groups, we found that transition/transversion ratios, as estimated by maximum likelihood on most-parsimonious trees, ranged between 0.8 and 1.0 for ingroups. These values occurred both at low sequence divergences, where substitutional saturation, i.e., multiple substitutions having occurred at the same (homologous) nucleotide position, was not expected, and at higher levels of divergence. This suggests that the angiosperm trnL-F regions evolve in a pattern different from that generally observed for nuclear and animal mtDNA (transitional/transversion ratio > or = 2). Transition/transversion ratios in the intron and the spacer region differed in all alignments compared, yet base compositions between the regions were highly similar in all six groups. A>-C transversions were significantly less frequent than the other four substitution types. This correlates with results from studies on fidelity mechanisms in DNA replication that predict A<->T and G<->C transversions to be least likely to occur. It therefore strengthens confidence in the link between mutation bias at the polymerase level and the actual fixation of substitutions as recorded on evolutionary trees, and concomitantly, in the neutrality of nucleotide substitutions as phylogenetic markers.

Polymorphism analysis of the hsp70 stress gene in Broiler chickens (Gallus gallus) of different breeds

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)