Identifying protein-coding genes in genomic sequences

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

16q24.1 microdeletion in a premature newborn: Usefulness of array-based comparative genomic hybridization in persistent pulmonary hypertension of the newborn.

OBJECTIVE:: Report of a 16q24.1 deletion in a premature newborn, demonstrating the usefulness of array-based comparative genomic hybridization in persistent pulmonary hypertension of the newborn and multiple congenital malformations. DESIGN:: Descriptive case report. SETTING:: Genetic department and neonatal intensive care unit of a tertiary care children's hospital. INTERVENTIONS:: None. PATIENT:: We report the case of a preterm male infant, born at 26 wks of gestation. A cardiac malformation and bilateral hydronephrosis were diagnosed at 19 wks of gestation. Karyotype analysis was normal, and a 22q11.2 microdeletion was excluded by fluorescence in situ hybridization analysis. A cesarean section was performed due to fetal distress. The patient developed persistent pulmonary hypertension unresponsive to mechanical ventilation and nitric oxide treatment and expired at 16 hrs of life. MEASUREMENTS AND MAIN RESULTS:: An autopsy revealed partial atrioventricular canal malformation and showed bilateral dilation of the renal pelvocaliceal system with bilateral ureteral stenosis and annular pancreas. Array-based comparative genomic hybridization analysis (Agilent oligoNT 44K, Agilent Technologies, Santa Clara, CA) showed an interstitial microdeletion encompassing the forkhead box gene cluster in 16q24.1. Review of the pulmonary microscopic examination showed the characteristic features of alveolar capillary dysplasia with misalignment of pulmonary veins. Some features were less prominent due to the gestational age. CONCLUSIONS:: Our review of the literature shows that alveolar capillary dysplasia with misalignment of pulmonary veins is rare but probably underreported. Prematurity is not a usual presentation, and histologic features are difficult to interpret. In our case, array-based comparative genomic hybridization revealed a 16q24.1 deletion, leading to the final diagnosis of alveolar capillary dysplasia with misalignment of pulmonary veins. It emphasizes the usefulness of array-based comparative genomic hybridization analysis as a diagnostic tool with implications for both prognosis and management decisions in newborns with refractory persistent pulmonary hypertension and multiple congenital malformations.

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition.

Isolation and characterization of highly polymorphic microsatellite loci in the bladder campion, Silene vulgaris (Caryophyllaceae)

This study reports the isolation and characterization of seven highly polymorphic microsatellite loci in Silene vulgaris (Caryophyllaceae). The loci were isolated from two libraries constructed from genomic DNA enriched for CA and GA repeats. These markers yielded nine to 40 alleles per locus (mean 22.1) in a survey of 45 individuals from a single population located in the western Swiss Alps. Average observed heterozygosity ranged from 16.2 to 77.4%. These microsatellite loci should be valuable tools for studying fine-scale genetic structure.

Hundreds of variants clustered in genomic loci and biological pathways affect human height.

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

Hybridization Capture Using RAD Probes (hyRAD), a New Tool for Performing Genomic Analyses on Collection Specimens.

In the recent years, many protocols aimed at reproducibly sequencing reduced-genome subsets in non-model organisms have been published. Among them, RAD-sequencing is one of the most widely used. It relies on digesting DNA with specific restriction enzymes and performing size selection on the resulting fragments. Despite its acknowledged utility, this method is of limited use with degraded DNA samples, such as those isolated from museum specimens, as these samples are less likely to harbor fragments long enough to comprise two restriction sites making possible ligation of the adapter sequences (in the case of double-digest RAD) or performing size selection of the resulting fragments (in the case of single-digest RAD). Here, we address these limitations by presenting a novel method called hybridization RAD (hyRAD). In this approach, biotinylated RAD fragments, covering a random fraction of the genome, are used as baits for capturing homologous fragments from genomic shotgun sequencing libraries. This simple and cost-effective approach allows sequencing of orthologous loci even from highly degraded DNA samples, opening new avenues of research in the field of museum genomics. Not relying on the restriction site presence, it improves among-sample loci coverage. In a trial study, hyRAD allowed us to obtain a large set of orthologous loci from fresh and museum samples from a non-model butterfly species, with a high proportion of single nucleotide polymorphisms present in all eight analyzed specimens, including 58-year-old museum samples. The utility of the method was further validated using 49 museum and fresh samples of a Palearctic grasshopper species for which the spatial genetic structure was previously assessed using mtDNA amplicons. The application of the method is eventually discussed in a wider context. As it does not rely on the restriction site presence, it is therefore not sensitive to among-sample loci polymorphisms in the restriction sites that usually causes loci dropout. This should enable the application of hyRAD to analyses at broader evolutionary scales.

Genomic analysis of the blood attributed to Louis XVI (1754-1793), king of France.

A pyrographically decorated gourd, dated to the French Revolution period, has been alleged to contain a handkerchief dipped into the blood of the French king Louis XVI (1754-1793) after his beheading but recent analyses of living males from two Bourbon branches cast doubts on its authenticity. We sequenced the complete genome of the DNA contained in the gourd at low coverage (similar to 2.5x) with coding sequences enriched at a higher similar to 7.3x coverage. We found that the ancestry of the gourd's genome does not seem compatible with Louis XVI's known ancestry. From a functional perspective, we did not find an excess of alleles contributing to height despite being described as the tallest person in Court. In addition, the eye colour prediction supported brown eyes, while Louis XVI had blue eyes. This is the first draft genome generated from a person who lived in a recent historical period; however, our results suggest that this sample may not correspond to the alleged king.

Partial characterization of nif genes from the bacterium Azospirillum amazonense

Azospirillum amazonense revealed genomic organization patterns of the nitrogen fixation genes similar to those of the distantly related species A. brasilense. Our work suggests that A. brasilense nifHDK, nifENX, fixABC operons and nifA and glnB genes may be structurally homologous to the counterpart genes of A. amazonense. This is the first analysis revealing homology between A. brasilense nif genes and the A. amazonense genome. Sequence analysis of PCR amplification products revealed similarities between the amino acid sequences of the highly conserved nifD and glnB genes of A. amazonense and related genes of A. brasilense and other bacteria. However, the A. amazonense non-coding regions (the upstream activator sequence region and the region between the nifH and nifD genes) differed from related regions of A. brasilense even in nitrogenase structural genes which are highly conserved among diazotrophic bacteria. The feasibility of the 16S ribosomal RNA gene-based PCR system for specific detection of A. amazonense was shown. Our results indicate that the PCR primers for 16S rDNA defined in this article are highly specific to A. amazonense and can distinguish this species from A. brasilense.

Mutations in SRY and WT1 genes required for gonadal development are not responsible for XY partial gonadal dysgenesis

The WT1 transcription factor regulates SRY expression during the initial steps of the sex determination process in humans, activating a gene cascade leading to testis differentiation. In addition to causing Wilms' tumor, mutations in WT1 are often responsible for urogenital defects in men, while SRY mutations are mainly related to 46,XY pure gonadal dysgenesis. In order to evaluate their role in abnormal testicular organogenesis, we screened for SRY and WT1 gene mutations in 10 children with XY partial gonadal dysgenesis, 2 of whom with a history of Wilms' tumor. The open reading frame and 360 bp of the 5' flanking sequence of the SRY gene, and the ten exons and intron boundaries of the WT1 gene were amplified by PCR of genomic DNA. Single-strand conformation polymorphism was initially used for WT1 mutation screening. Since shifts in fragment migration were only observed for intron/exon 4, the ten WT1 exons from all patients were sequenced manually. No mutations were detected in the SRY 5' untranslated region or within SRY open-reading frame sequences. WT1 sequencing revealed one missense mutation (D396N) in the ninth exon of a patient who also had Wilms' tumor. In addition, two silent point mutations were found in the first exon including one described here for the first time. Some non-coding sequence variations were detected, representing one new (IVS4+85A>G) and two already described (-7ATG T>G, IVS9-49 T>C) single nucleotide polymorphisms. Therefore, mutations in two major genes required for gonadal development, SRY and WT1, are not responsible for XY partial gonadal dysgenesis.

Protective effects of organoselenium compounds against methylmercury-induced oxidative stress in mouse brain mitochondrial-enriched fractions

We evaluated the potential neuroprotective effect of 1-100 µM of four organoselenium compounds: diphenyl diselenide, 3’3-ditri-fluoromethyldiphenyl diselenide, p-methoxy-diphenyl diselenide, and p-chloro-diphenyl diselenide, against methylmercury-induced mitochondrial dysfunction and oxidative stress in mitochondrial-enriched fractions from adult Swiss mouse brain. Methylmercury (10-100 µM) significantly decreased mitochondrial activity, assessed by MTT reduction assay, in a dose-dependent manner, which occurred in parallel with increased glutathione oxidation, hydroperoxide formation (xylenol orange assay) and lipid peroxidation end-products (thiobarbituric acid reactive substances, TBARS). The co-incubation with diphenyl diselenide (100 µM) completely prevented the disruption of mitochondrial activity as well as the increase in TBARS levels caused by methylmercury. The compound 3’3-ditrifluoromethyldiphenyl diselenide provided a partial but significant protection against methylmercury-induced mitochondrial dysfunction (45.4 ± 5.8% inhibition of the methylmercury effect). Diphenyl diselenide showed a higher thiol peroxidase activity compared to the other three compounds. Catalase blocked methylmercury-induced TBARS, pointing to hydrogen peroxide as a vector during methylmercury toxicity in this model. This result also suggests that thiol peroxidase activity of organoselenium compounds accounts for their protective actions against methylmercury-induced oxidative stress. Our results show that diphenyl diselenide and potentially other organoselenium compounds may represent important molecules in the search for an improved therapy against the deleterious effects of methylmercury as well as other mercury compounds.

Partial chemical and functional characterization of milk whey products obtained by different processes

Whey protein samples (S-1 to S-5) were tested in vivo and in vitro for nutritional properties and selected bioactivities. Weanling male Wistar rats fed modified AIN-93G (12 g protein.100 g-1) diets for 21 days were used the in vivo studies. The nutritional parameters did not differ among the protein diets tested. Erythrocyte glutathione content was considered high and was higher for S-3, but liver glutathione was the same for all dietary groups. For S-3, cytokine secretion (IL-10 and TNF-α) by human peripheral blood mononuclear cells (in RPMI-1640 medium) was higher in the absence of antigen than in the presence of BCG antigen. Interleukin-4 secretion was repressed in all treatments. The IC50, whey protein concentration required to inhibit 50% of the melanoma cell proliferation, was 2.68 mg.mL-1 of culture medium for the S-3 sample and 3.66 mg.mL-1 for the S-2 sample. Based on these results, it was concluded that S-3 (whey protein concentrate enriched with TGF-β and lactoferrin) produced better nutritional and immunological responses than the other products tested.

Genomic architecture of sickle cell disease clinical variation in children from West Africa : a case-control study design

Contexte : L’anémie falciforme ou drépanocytose est un problème de santé important, particulièrement pour les patients d’origine africaine. La variation phénotypique de l’anémie falciforme est problématique pour le suivi et le traitement des patients. L’architecture génomique responsable de cette variabilité est peu connue. Principe : Mieux saisir la contribution génétique de la variation clinique de cette maladie facilitera l’identification des patients à risque de développer des phénotypes sévères, ainsi que l’adaptation des soins. Objectifs : L’objectif général de cette thèse est de combler les lacunes relatives aux connaissances sur l’épidémiologie génomique de l’anémie falciforme à l’aide d’une cohorte issue au Bénin. Les objectifs spécifiques sont les suivants : 1) caractériser les profils d’expressions génomiques associés à la sévérité de l’anémie falciforme ; 2) identifier des biomarqueurs de la sévérité de l’anémie falciforme ; 3) identifier la régulation génétique des variations transcriptionelles ; 4) identifier des interactions statistiques entre le génotype et le niveau de sévérité associé à l’expression ; 5) identifier des cibles de médicaments pour améliorer l’état des patients atteints d’anémie falciforme. Méthode : Une étude cas-témoins de 250 patients et 61 frères et soeurs non-atteints a été menée au Centre de Prise en charge Médical Intégré du Nourrisson et de la Femme Enceinte atteints de Drépanocytose, au Bénin entre février et décembre 2010. Résultats : Notre analyse a montré que des profils d’expressions sont associés avec la sévérité de l’anémie falciforme. Ces profils sont enrichis de génes des voies biologiques qui contribuent à la progression de la maladie : l’activation plaquettaire, les lymphocytes B, le stress, l’inflammation et la prolifération cellulaire. Des biomarqueurs transcriptionnels ont permis de distinguer les patients ayant des niveaux de sévérité clinique différents. La régulation génétique de la variation de l’expression des gènes a été démontrée et des interactions ont été identifiées. Sur la base de ces résultats génétiques, des cibles de médicaments sont proposées. Conclusion: Ce travail de thèse permet de mieux comprendre l’impact de la génomique sur la sévérité de l’anémie falciforme et ouvre des perspectives de développement de traitements ciblés pour améliorer les soins offerts aux patients.

Genomic architecture of adaptive color pattern divergence and convergence in Heliconius butterflies

Identifying the genetic changes driving adaptive variation in natural populations is key to understanding the origins of biodiversity. The mosaic of mimetic wing patterns in Heliconius butterflies makes an excellent system for exploring adaptive variation using next-generation sequencing. In this study, we use a combination of techniques to annotate the genomic interval modulating red color pattern variation, identify a narrow region responsible for adaptive divergence and convergence in Heliconius wing color patterns, and explore the evolutionary history of these adaptive alleles. We use whole genome resequencing from four hybrid zones between divergent color pattern races of Heliconius erato and two hybrid zones of the co-mimic Heliconius melpomene to examine genetic variation across 2.2 Mb of a partial reference sequence. In the intergenic region near optix, the gene previously shown to be responsible for the complex red pattern variation in Heliconius, population genetic analyses identify a shared 65-kb region of divergence that includes several sites perfectly associated with phenotype within each species. This region likely contains multiple cis-regulatory elements that control discrete expression domains of optix. The parallel signatures of genetic differentiation in H. erato and H. melpomene support a shared genetic architecture between the two distantly related co-mimics; however, phylogenetic analysis suggests mimetic patterns in each species evolved independently. Using a combination of next-generation sequencing analyses, we have refined our understanding of the genetic architecture of wing pattern variation in Heliconius and gained important insights into the evolution of novel adaptive phenotypes in natural populations.

Molecular genetic differentiation in earthworms inhabiting a heterogeneous Pb-polluted landscape

A Pb-mine site situated on acidic soil, but comprising of Ca-enriched islands around derelict buildings was used to study the spatial pattern of genetic diversity in Lumbricus rubellus. Two distinct genetic lineages ('A' and 'B'), differentiated at both the mitochondrial (mtDNA COII) and nuclear level (AFLPs) were revealed with a mean inter-lineage mtDNA sequence divergence of approximately 13%, indicative of a cryptic species complex. AFLP analysis indicates that lineage A individuals within one central 'ecological island' site are uniquely clustered, with little genetic overlap with lineage A individuals at the two peripheral sites. FTIR microspectroscopy of Pb-sequestering chloragocytes revealed different phosphate profiles in residents of adjacent acidic and calcareous islands. Bioinformatics found over-representation of Ca pathway genes in ESTPb libraries. Subsequent sequencing of a Ca-transport gene, SERCA, revealed mutations in the protein's cytosolic domain. We recommend the mandatory genotyping of all individuals prior to field-based ecotoxicological assays, particularly those using discriminating genomic technologies.

Biodiversity of human faecal bacteria isolated from phytic acid enriched chemostat fermenters

Background: Myo-inositol hexaphosphate (IP6) or phytic acid is found mostly in cereals and legumes and is thought to possess anti-carcinogenic properties. Aim: To isolate and identify faecal bacteria capable of phytic acid metabolism and to assess the effectiveness of prebiotics (dietary oligosaccharides, metabolised by selective colonic bacteria) in preserving the integrity of phytic acid. Methods: Faecal samples from three volunteers were used in continuous culture experiments under varying conditions of pH, substrate concentration and dilution rates, seventy three different isolates cultured at steady state were then screened for phytic acid metabolism and identified through partial sequencing of their 16S rRNA genes (16S ribosomal ribonucleic acid). Utilisation of phytic acid was also assessed in a continuous culture system enriched with prebiotic fructooligosaccharides (FOS). Results: Bacteroides spp., Clostridium spp. and facultatively anaerobic bacteria generally appeared to maintain viable counts in the presence of phytic acid. Bifidobacterium spp. and Lactobacillus spp. appeared less able to maintain viable counts in the presence of phytic acid. These results were confirmed by an increase in viable counts of Bacteroides spp., Clostridium spp. and a decrease in viable counts of Bifidobacterium spp. and Lactobacillus spp. once phytic acid was introduced to a FOS enriched continuous culture. Conclusions: The phytate metabolising biodiversity from the human large intestine does not appear to encompass major bacterial genera associated with beneficial or benign health effects (e.g. Lactobacillus spp. and Bifidobacterium spp).