(Table S2) Relative abundance of archaeal 16S rRNA genes in sediment cores

Deep drilling into the marine sea floor has uncovered a vast sedimentary ecosystem of microbial cells (Parkes et al., 1994, doi:10.1038/371410a0; D'Hondt et al., 2004, doi:10.1126/science.1101155). Extrapolation of direct counts of stained microbial cells to the total volume of habitable marine subsurface sediments suggests that between 56 Pg (Parkes et al., 1994, doi:10.1038/371410a0) and 303 Pg (Whitman et al., 1998) of cellular carbon could be stored in this largely unexplored habitat. From recent studies using various culture-independent techniques, no clear picture has yet emerged as to whether Archaea or Bacteria are more abundant in this extensive ecosystem (Schippers et al., doi:10.1038/nature03302; Inagaki et al., doi:10.1073/pnas.0511033103 ; Mauclaire et al., doi:10.1111/j.1472-4677.2004.00035.x; Biddle et al., doi:10.1073/pnas.0600035103). Here we show that in subsurface sediments buried deeper than 1 m in a wide range of oceanographic settings at least 87% of intact polar membrane lipids, biomarkers for the presence of live cells (Biddle et al., doi:10.1073/pnas.0600035103; Sturt et al., 2004, doi:10.1002/rcm.1378), are attributable to archaeal membranes, suggesting that Archaea constitute a major fraction of the biomass. Results obtained from modified quantitative polymerase chain reaction and slot-blot hybridization protocols support the lipid-based evidence and indicate that these techniques have previously underestimated archaeal biomass. The lipid concentrations are proportional to those of total organic carbon. On the basis of this relationship, we derived an independent estimate of amounts of cellular carbon in the global marine subsurface biosphere. Our estimate of 90 Pg of cellular carbon is consistent, within an order of magnitude, with previous estimates, and underscores the importance of marine subsurface habitats for global biomass budgets.

Development of molecular approaches to assess the diversity and abundance of marine nematodes assemblages for assessment of Good Environmental Status of the estuarine and marine ecosystems

In Europe, the concerns with the status of marine ecosystems have increased, and the Marine Directive has as main goal the achievement of Good Environmental Status (GES) of EU marine waters by 2020. Molecular tools are seen as promising and emerging approaches to improve ecosystem monitoring, and have led ecology into a new era, representing perhaps the most source of innovation in marine monitoring techniques. Benthic nematodes are considered ideal organisms to be used as biological indicator of natural and anthropogenic disturbances in aquatic ecosystems underpinning monitoring programmes on the ecological quality of marine ecosystems, very useful to assess the GES of the marine environment. dT-RFLP (directed Terminal-Restriction Fragment Length Polymorphism) allows to assess the diversity of nematode communities, but also allows studying the functioning of the ecosystem, and combined with relative real-time PCR (qPCR), provides a high-throughput semi-quantitative characterization of nematode communities. These characteristics make the two molecular tools good descriptors for the good environmental status assessment. The main aim of this study is to develop and optimize the dT-RFLP and qPCR in Mira estuary (SW coast, Portugal). A molecular phylogenetic analysis of marine and estuarine nematodes is being performed combining morphological and molecular analysis to evaluate the diversity of free-living marine nematodes in Mira estuary. After morphological identification, barcoding of 18S rDNA and COI genes are being determined for each nematode species morphologically identified. So far we generated 40 new sequences belonging to 32 different genus and 17 families, and the study has shown a good degree of concordance between traditional morphology-based identification and DNA sequences. These results will improve the assessment of marine nematode diversity and contribute to a more robust nematode taxonomy. The DNA sequences are being used to develop the dT-RFLP with the ability to easily process large sample numbers (hundreds and thousands), rather than typical of classical taxonomic or low throughput molecular analyses. A preliminary study showed that the digest enzymes used in dT-RFLP for terrestrial assemblages separated poorly the marine nematodes at taxonomic level for functional group analysis. A new digest combination was designed using the software tool DRAT (Directed Terminal Restriction Analysis Tool) to distinguished marine nematode taxa. Several solutions were provided by DRAT and tested empirically to select the solution that cuts most efficiently. A combination of three enzymes and a single digest showed to be the best solution to separate the different clusters. Parallel to this, another tool is being developed to estimate the population size (qPCR). An improvement in qPCR estimation of gene copy number using an artificial reference is being performed for marine nematodes communities to quantify the abundance. Once developed, it is proposed to validate both methodologies by determining the spatial and temporal variability of benthic nematodes assemblages across different environments. The application of these high-throughput molecular approaches for benthic nematodes will improve sample throughput and their implementation more efficient and faster as indicator of ecological status of marine ecosystems.

Cuff-induced vascular intima thickening is influenced by titration of the Ace gene in mice

Lacchini S, Heimann AS, Evangelista FS, Cardoso L, Silva GJ, Krieger JE. Cuff-induced vascular intima thickening is influenced by titration of the Ace gene in mice. Physiol Genomics 37: 225-230, 2009. First published March 3, 2009; doi:10.1152/physiolgenomics.90288.2008.-We tested the hypothesis that small changes in angiotensin I-converting enzyme (ACE) expression can alter the vascular response to injury. Male mice containing one, two, three, and four copies of the Ace gene with no detectable vascular abnormality or changes in blood pressure were submitted to cuff-induced femoral artery injury. Femoral thickening was higher in 3- and 4-copy mice (42.4 +/- 4.3% and 45.7 +/- 6.5%, respectively) compared with 1- and 2-copy mice (8.3 +/- 1.3% and 8.5 +/- 0.9%, respectively). Femoral ACE levels from control and injured vessels were assessed in 1- and 3-copy Ace mice, which represent the extremes of the observed response. ACE vascular activity was higher in 3- vs. 1-copy Ace mice (2.4-fold, P < 0.05) in the control uninjured vessel. Upon injury, ACE activity significantly increased in both groups [2.41-fold and 2.14-fold (P < 0.05) for 1- and 3-copy groups, respectively] but reached higher levels in 3- vs. 1-copy Ace mice (P < 0.05). Pharmacological interventions were then used as a counterproof and to indirectly assess the role of angiotensin II (ANG II) on this response. Interestingly, ACE inhibition (enalapril) and ANG II AT(1) receptor blocker (losartan) reduced intima thickening in 3-copy mice to 1-copy mouse values (P < 0.05) while ANG II treatment significantly increased intima thickening in 1-copy mice to 3-copy mouse levels (P < 0.05). Together, these data indicate that small physiologically relevant changes in ACE, not associated with basal vascular abnormalities or blood pressure levels, do influence the magnitude of cuff-induced neointima thickening in mice.

Steroidogenic Factor 1 Overexpression and Gene Amplification Are More Frequent in Adrenocortical Tumors from Children than from Adults

Background: Steroidogenic factor 1 (SF-1) is a key determinant of endocrine development and function of adrenal cortex. SF-1 overexpression and gene amplification were previously demonstrated in a small group of pediatric adrenocortical tumors. Objective: Our objective was to determine the frequency of SF-1 protein expression and gene amplification in a large cohort of pediatric and adult adrenocortical tumors. Patients: SF-1 protein expression was assessed in a cohort of 103 adrenocortical tumors from 36 children and 67 adults, whereas gene amplification was studied in 38 adrenocortical tumors ( 17 from children). Methods: Tissue microarray, multiplex ligation-dependent probe amplification, and quantitative real-time PCR were used. Results: Astrong nuclear SF-1 expression was detected by tissue microarray in 56% (20 of 36) and 19% (13 of 67) of the pediatric and adult adrenocortical tumors, respectively (P = 0.0004). Increased SF-1 copy number was identified in 47% (eight of 17) and 10% (two of 21) of the pediatric and adult adrenocortical tumors, respectively (P = 0.02). All adrenocortical tumors with SF-1 gene amplification showed a strong SF-1 staining, whereas most of the tumors (61%) without SF-1 amplification displayed a weak or negative staining (P = 0.0008). Interestingly, a strong SF-1 staining was identified in five (29%) pediatric adrenocortical tumors without SF-1 amplification. The frequency of SF-1 overexpression and gene amplification was similar in adrenocortical adenomas and carcinomas. Conclusion: We demonstrated a higher frequency of SF-1 overexpression and gene amplification in pediatric than in adult adrenocortical tumors, suggesting an important role of SF-1 in pediatric adrenocortical tumorigenesis. (J Clin Endocrinol Metab 95: 1458-1462, 2010)

High-Resolution Genomic Mapping Reveals Consistent Amplification of the Fibroblast Growth Factor Receptor Substrate 2 Gene in Well-Differentiated and Dedifferentiated Liposarcoma

Well-differentiated liposarcoma (WDLS) is one of the most common malignant mesenchymal tumors and dedifferentiated liposarcoma (DDLS) is a malignant tumor consisting of both WDLS and a transformed nonlipogenic sarcomatous component. Cytogenetically, WDLS is characterized by the presence of ring or giant rod chromosomes containing several amplified genes, including MDM2, TSPAN31 CDK4, and others mainly derived from chromosome bands 12q13-15. However, the 12q13-15 amplicon is large and discontinuous. The focus of this study was to identify novel critical genes that are consistently amplified in primary (nonrecurrent) WDLS and with potential relevance for future targeted therapy. Using a high-resolution (5.0 kb) ""single nucleotide polymorphism""/copy number variation microarray to screen the whole genome in a series of primary WDLS, two consistently amplified areas were found on chromosome 12: one region containing the MDM2 and CPM genes, and another region containing the FRS2 gene. Based on these findings, we further validated FRS2 amplification in both WDLS and DDLS. Fluorescence in situ hybridization confirmed FRS2 amplification in all WDLS and DDLS tested (n = 57). Real time PCR showed FRS2 mRNA transcriptional upregulation in WDLS (n = 19) and DDLS (n = 13) but not in lipoma (n = 5) and normal fat (n = 9). Immunoblotting revealed high expression levels of phospho-FRS2 at 1436 and slightly overexpression of total FRS2 protein in liposarcoma but not in normal fat or preadipocytes. Considering the critical role of FRS2 in mediating fibroblast growth factor receptor signaling, our findings indicate that FRS2 signaling should be further investigated as a potential therapeutic target for liposarcoma. (C) 2011 Wiley-Liss, Inc.

Analysis of PAX2 gene alterations in renal cell tumors

Dissertation for applying to a Master’s Degree in Molecular Genetics and Biomedicine submitted to the Sciences and Technology Faculty of New University of Lisbon

Organization of lin genes and IS6100 among different strains of hexachlorocyclohexane-degrading Sphingomonas paucimobilis: evidence for horizontal gene transfer.

The organization of lin genes and IS6100 was studied in three strains of Sphingomonas paucimobilis (B90A, Sp+, and UT26) which degraded hexachlorocyclohexane (HCH) isomers but which had been isolated at different geographical locations. DNA-DNA hybridization data revealed that most of the lin genes in these strains were associated with IS6100, an insertion sequence classified in the IS6 family and initially found in Mycobacterium fortuitum. Eleven, six, and five copies of IS6100 were detected in B90A, Sp+, and UT26, respectively. IS6100 elements in B90A were sequenced from five, one, and one regions of the genomes of B90A, Sp+, and UT26, respectively, and were found to be identical. DNA-DNA hybridization and DNA sequencing of cosmid clones also revealed that S. paucimobilis B90A contains three and two copies of linX and linA, respectively, compared to only one copy of these genes in strains Sp+ and UT26. Although the copy number and the sequence of the remaining genes of the HCH degradative pathway (linB, linC, linD, and linE) were nearly the same in all strains, there were striking differences in the organization of the linA genes as a result of replacement of portions of DNA sequences by IS6100, which gave them a strange mosaic configuration. Spontaneous deletion of linD and linE from B90A and of linA from Sp+ occurred and was associated either with deletion of a copy of IS6100 or changes in IS6100 profiles. The evidence gathered in this study, coupled with the observation that the G+C contents of the linA genes are lower than that of the remaining DNA sequence of S. paucimobilis, strongly suggests that all these strains acquired the linA gene through horizontal gene transfer mediated by IS6100. The association of IS6100 with the rest of the lin genes further suggests that IS6100 played a role in shaping the current lin gene organization.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5′-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

Whole genome sequencing in patients with retinitis pigmentosa reveals pathogenic DNA structural changes and NEK2 as a new disease gene.

We performed whole genome sequencing in 16 unrelated patients with autosomal recessive retinitis pigmentosa (ARRP), a disease characterized by progressive retinal degeneration and caused by mutations in over 50 genes, in search of pathogenic DNA variants. Eight patients were from North America, whereas eight were Japanese, a population for which ARRP seems to have different genetic drivers. Using a specific workflow, we assessed both the coding and noncoding regions of the human genome, including the evaluation of highly polymorphic SNPs, structural and copy number variations, as well as 69 control genomes sequenced by the same procedures. We detected homozygous or compound heterozygous mutations in 7 genes associated with ARRP (USH2A, RDH12, CNGB1, EYS, PDE6B, DFNB31, and CERKL) in eight patients, three Japanese and five Americans. Fourteen of the 16 mutant alleles identified were previously unknown. Among these, there was a 2.3-kb deletion in USH2A and an inverted duplication of ∼446 kb in EYS, which would have likely escaped conventional screening techniques or exome sequencing. Moreover, in another Japanese patient, we identified a homozygous frameshift (p.L206fs), absent in more than 2,500 chromosomes from ethnically matched controls, in the ciliary gene NEK2, encoding a serine/threonine-protein kinase. Inactivation of this gene in zebrafish induced retinal photoreceptor defects that were rescued by human NEK2 mRNA. In addition to identifying a previously undescribed ARRP gene, our study highlights the importance of rare structural DNA variations in Mendelian diseases and advocates the need for screening approaches that transcend the analysis of the coding sequences of the human genome.

Recurrent deletions and reciprocal duplications of 10q11.21q11.23 including CHAT and SLC18A3 are likely mediated by complex low-copy repeats.

We report 24 unrelated individuals with deletions and 17 additional cases with duplications at 10q11.21q21.1 identified by chromosomal microarray analysis. The rearrangements range in size from 0.3 to 12 Mb. Nineteen of the deletions and eight duplications are flanked by large, directly oriented segmental duplications of >98% sequence identity, suggesting that nonallelic homologous recombination (NAHR) caused these genomic rearrangements. Nine individuals with deletions and five with duplications have additional copy number changes. Detailed clinical evaluation of 20 patients with deletions revealed variable clinical features, with developmental delay (DD) and/or intellectual disability (ID) as the only features common to a majority of individuals. We suggest that some of the other features present in more than one patient with deletion, including hypotonia, sleep apnea, chronic constipation, gastroesophageal and vesicoureteral refluxes, epilepsy, ataxia, dysphagia, nystagmus, and ptosis may result from deletion of the CHAT gene, encoding choline acetyltransferase, and the SLC18A3 gene, mapping in the first intron of CHAT and encoding vesicular acetylcholine transporter. The phenotypic diversity and presence of the deletion in apparently normal carrier parents suggest that subjects carrying 10q11.21q11.23 deletions may exhibit variable phenotypic expressivity and incomplete penetrance influenced by additional genetic and nongenetic modifiers.

Tandem repeat sequence variation as causative cis-eQTLs for protein-coding gene expression variation: the case of CSTB.

Association studies have revealed expression quantitative trait loci (eQTLs) for a large number of genes. However, the causative variants that regulate gene expression levels are generally unknown. We hypothesized that copy-number variation of sequence repeats contribute to the expression variation of some genes. Our laboratory has previously identified that the rare expansion of a repeat c.-174CGGGGCGGGGCG in the promoter region of the CSTB gene causes a silencing of the gene, resulting in progressive myoclonus epilepsy. Here, we genotyped the repeat length and quantified CSTB expression by quantitative real-time polymerase chain reaction in 173 lymphoblastoid cell lines (LCLs) and fibroblast samples from the GenCord collection. The majority of alleles contain either two or three copies of this repeat. Independent analysis revealed that the c.-174CGGGGCGGGGCG repeat length is strongly associated with CSTB expression (P = 3.14 × 10(-11)) in LCLs only. Examination of both genotyped and imputed single-nucleotide polymorphisms (SNPs) within 2 Mb of CSTB revealed that the dodecamer repeat represents the strongest cis-eQTL for CSTB in LCLs. We conclude that the common two or three copy variation is likely the causative cis-eQTL for CSTB expression variation. More broadly, we propose that polymorphic tandem repeats may represent the causative variation of a fraction of cis-eQTLs in the genome.

Mirror extreme BMI phenotypes associated with gene dosage at the chromosome 16p11.2 locus.

Both obesity and being underweight have been associated with increased mortality. Underweight, defined as a body mass index (BMI) ≤ 18.5 kg per m(2) in adults and ≤ -2 standard deviations from the mean in children, is the main sign of a series of heterogeneous clinical conditions including failure to thrive, feeding and eating disorder and/or anorexia nervosa. In contrast to obesity, few genetic variants underlying these clinical conditions have been reported. We previously showed that hemizygosity of a ∼600-kilobase (kb) region on the short arm of chromosome 16 causes a highly penetrant form of obesity that is often associated with hyperphagia and intellectual disabilities. Here we show that the corresponding reciprocal duplication is associated with being underweight. We identified 138 duplication carriers (including 132 novel cases and 108 unrelated carriers) from individuals clinically referred for developmental or intellectual disabilities (DD/ID) or psychiatric disorders, or recruited from population-based cohorts. These carriers show significantly reduced postnatal weight and BMI. Half of the boys younger than five years are underweight with a probable diagnosis of failure to thrive, whereas adult duplication carriers have an 8.3-fold increased risk of being clinically underweight. We observe a trend towards increased severity in males, as well as a depletion of male carriers among non-medically ascertained cases. These features are associated with an unusually high frequency of selective and restrictive eating behaviours and a significant reduction in head circumference. Each of the observed phenotypes is the converse of one reported in carriers of deletions at this locus. The phenotypes correlate with changes in transcript levels for genes mapping within the duplication but not in flanking regions. The reciprocal impact of these 16p11.2 copy-number variants indicates that severe obesity and being underweight could have mirror aetiologies, possibly through contrasting effects on energy balance.

Understanding mutational and selective processes that govern gene duplication in mammals

Abstract : Gene duplication is an essential source of material for the origin of genetic novelties. The reverse transcription of source gene mRNA followed by the genomic insertion of the resulting cDNA - retroposition - has provided the human genome with at least ~3600 detectable retrocopies. We find that ~30% of these retrocopies are transcribed, generally in testes. Their transcription often relies on preexisting regulatory elements (or open chromatin) close to their insertion site, which is illustrated by mRNA molecules containing retrocopies fused to their neighboring genes. Retrocopies appear to have been profoundly shaped by selection. Consistently, human retrocopies with an intact open reading (ORF) are more often transcribed than retropseudogenes, which leads to a minimal estimate of 120 functional retrogenes present in our genome. We also performed an analysis of Ka/Ks for human retrocopies. This analysis demonstrates that several intact retrocopies evolved under purifying selection and yields an estimated formation rate of ~1 retrogene per million year in the primate lineage. Using DNA sequencing and evolutionary simulations, we have identified 7 such primate-specific retrogenes that emerged on the lineage leading to humans In therian genomes, we found an excess of retrogenes with X-linked parents. Expression analyses support the idea that this "out of X" movement was driven by natural selection to produce autosomal functional counterparts for X-linked genes, which are silenced during male meiosis. Phylogenetic dating of this "out of X" movement suggests that our sex chromosomes arose about 180 MYA ago and are thus much younger than previously thought. Finally, we have also analyzed young gene duplications (and deletions) that arose by non allelic-homologous recombination and are not fixed in species. Using wild-caught and laboratory animals, we detected thousands of DNA segments that are polymorphic in copy number in mice. These copy number variants were found to profoundly alter the transcriptome of several mouse tissues. Strikingly, their influence on gene expression is not limited to the gene they contain but seems to extend to genes located up to 1.5 million bases away.

Validation of real-time methylation-specific PCR to determine O6-methylguanine-DNA methyltransferase gene promoter methylation in glioma.

Epigenetic silencing of the DNA repair protein O(6)-methylguanine-DNA methyltransferase (MGMT) by promoter methylation predicts successful alkylating agent therapy, such as with temozolomide, in glioblastoma patients. Stratified therapy assignment of patients in prospective clinical trials according to tumor MGMT status requires a standardized diagnostic test, suitable for high-throughput analysis of small amounts of formalin-fixed, paraffin-embedded tumor tissue. A direct, real-time methylation-specific PCR (MSP) assay was developed to determine methylation status of the MGMT gene promoter. Assay specificity was obtained by selective amplification of methylated DNA sequences of sodium bisulfite-modified DNA. The copy number of the methylated MGMT promoter, normalized to the beta-actin gene, provides a quantitative test result. We analyzed 134 clinical glioma samples, comparing the new test with the previously validated nested gel-based MSP assay, which yields a binary readout. A cut-off value for the MGMT methylation status was suggested by fitting a bimodal normal mixture model to the real-time results, supporting the hypothesis that there are two distinct populations within the test samples. Comparison of the tests showed high concordance of the results (82/91 [90%]; Cohen's kappa = 0.80; 95% confidence interval, 0.82-0.95). The direct, real-time MSP assay was highly reproducible (Pearson correlation 0.996) and showed valid test results for 93% (125/134) of samples compared with 75% (94/125) for the nested, gel-based MSP assay. This high-throughput test provides an important pharmacogenomic tool for individualized management of alkylating agent chemotherapy.

A functional microsatellite of the macrophage migration inhibitory factor gene associated with meningococcal disease.

Macrophage migration inhibitory factor (MIF) is an abundantly expressed proinflammatory cytokine playing a critical role in innate immunity and sepsis and other inflammatory diseases. We examined whether functional MIF gene polymorphisms (-794 CATT(5-8) microsatellite and -173 G/C SNP) were associated with the occurrence and outcome of meningococcal disease in children. The CATT(5) allele was associated with the probability of death predicted by the Pediatric Index of Mortality 2 (P=0.001), which increased in correlation with the CATT(5) copy number (P=0.04). The CATT(5) allele, but not the -173 G/C alleles, was also associated with the actual mortality from meningoccal sepsis [OR 2.72 (1.2-6.4), P=0.02]. A family-based association test (i.e., transmission disequilibrium test) performed in 240 trios with 1 afflicted offspring indicated that CATT(5) was a protective allele (P=0.02) for the occurrence of meningococcal disease. At baseline and after stimulation with Neisseria meningitidis in THP-1 monocytic cells or in a whole-blood assay, CATT(5) was found to be a low-expression MIF allele (P=0.005 and P=0.04 for transcriptional activity; P=0.09 and P=0.09 for MIF production). Taken together, these data suggest that polymorphisms of the MIF gene affecting MIF expression are associated with the occurrence, severity, and outcome of meningococcal disease in children.