Gene-nutrient interactions that impact magnesium homeostasis increase risk for neural tube defects in mice exposed to dolutegravir.

In 2018, data from a surveillance study in Botswana evaluating adverse birth outcomes raised concerns that women on antiretroviral therapy (ART) containing dolutegravir (DTG) may be at increased risk for neural tube defects (NTDs). The mechanism of action for DTG involves chelation of Mg2+ ions in the active site of the viral integrase. Plasma Mg2+ homeostasis is maintained primarily through dietary intake and reabsorption in the kidneys. Inadequate dietary Mg2+ intake over several months results in slow depletion of plasma Mg2+ and chronic latent hypomagnesemia, a condition prevalent in women of reproductive age worldwide. Mg2+ is critical for normal embryonic development and neural tube closure. We hypothesized that DTG therapy might slowly deplete plasma Mg2+ and reduce the amount available to the embryo, and that mice with pre-existing hypomagnesemia due to genetic variation and/or dietary Mg2+ insufficiency at the time of conception and initiation of DTG treatment would be at increased risk for NTDs. We used two different approaches to test our hypothesis: 1) we selected mouse strains that had inherently different basal plasma Mg2+ levels and 2) placed mice on diets with different concentrations of Mg2+. Plasma and urine Mg2+ were determined prior to timed mating. Pregnant mice were treated daily with vehicle or DTG beginning on the day of conception and embryos examined for NTDs on gestational day 9.5. Plasma DTG was measured for pharmacokinetic analysis. Our results demonstrate that hypomagnesemia prior to conception, due to genetic variation and/or insufficient dietary Mg2+ intake, increases the risk for NTDs in mice exposed to DTG. We also analyzed whole-exome sequencing data from inbred mouse strains and identified 9 predicted deleterious missense variants in Fam111a that were unique to the LM/Bc strain. Human FAM111A variants are associated with hypomagnesemia and renal Mg2+ wasting. The LM/Bc strain exhibits this same phenotype and was the strain most susceptible to DTG-NTDs. Our results suggest that monitoring plasma Mg2+ levels in patients on ART regimens that include DTG, identifying other risk factors that impact Mg2+ homeostasis, and correcting deficiencies in this micronutrient might provide an effective strategy for mitigating NTD risk.

Metabolome-based signature of disease pathology in MS.

BACKGROUND: Diagnostic delays are common for multiple sclerosis (MS) since diagnosis typically depends on the presentation of nonspecific clinical symptoms together with radiologically-determined central nervous system (CNS) lesions. It is important to reduce diagnostic delays as earlier initiation of disease modifying therapies mitigates long-term disability. Developing a metabolomic blood-based MS biomarker is attractive, but prior efforts have largely focused on specific subsets of metabolite classes or analytical platforms. Thus, there are opportunities to interrogate metabolite profiles using more expansive and comprehensive approaches for developing MS biomarkers and for advancing our understanding of MS pathogenesis. METHODS: To identify putative blood-based MS biomarkers, we comprehensively interrogated the metabolite profiles in 12 non-Hispanic white, non-smoking, male MS cases who were drug naïve for 3 months prior to biospecimen collection and 13 non-Hispanic white, non-smoking male controls who were frequency matched to cases by age and body mass index. We performed untargeted two-dimensional gas chromatography and time-of-flight mass spectrometry (GCxGC-TOFMS) and targeted lipidomic and amino acid analysis on serum. 325 metabolites met quality control and supervised machine learning was used to identify metabolites most informative for MS status. The discrimination potential of these select metabolites were assessed using receiver operator characteristic curves based on logistic models; top candidate metabolites were defined as having area under the curves (AUC) >80%. The associations between whole-genome expression data and the top candidate metabolites were examined, followed by pathway enrichment analyses. Similar associations were examined for 175 putative MS risk variants and the top candidate metabolites. RESULTS: 12 metabolites were determined to be informative for MS status, of which 6 had AUCs >80%: pyroglutamate, laurate, acylcarnitine C14:1, N-methylmaleimide, and 2 phosphatidylcholines (PC ae 40:5, PC ae 42:5). These metabolites participate in glutathione metabolism, fatty acid metabolism/oxidation, cellular membrane composition, and transient receptor potential channel signaling. Pathway analyses based on the gene expression association for each metabolite suggested enrichment for pathways associated with apoptosis and mitochondrial dysfunction. Interestingly, the predominant MS genetic risk allele HLA-DRB1×15:01 was associated with one of the 6 top metabolites. CONCLUSION: Our analysis represents the most comprehensive description of metabolic changes associated with MS in serum, to date, with the inclusion of genomic and genetic information. We identified atypical metabolic processes that differed between MS patients and controls, which may enable the development of biological targets for diagnosis and treatment.

Follow-up examination of linkage and association to chromosome 1q43 in multiple sclerosis.

Multiple sclerosis (MS) is a debilitating neuroimmunological and neurodegenerative disease affecting >4,00,000 individuals in the United States. Population and family-based studies have suggested that there is a strong genetic component. Numerous genomic linkage screens have identified regions of interest for MS loci. Our own second-generation genome-wide linkage study identified a handful of non-major histocompatibility complex regions with suggestive linkage. Several of these regions were further examined using single-nucleotide polymorphisms (SNPs) with average spacing between SNPs of approximately 1.0 Mb in a dataset of 173 multiplex families. The results of that study provided further evidence for the involvement of the chromosome 1q43 region. This region is of particular interest given linkage evidence in studies of other autoimmune and inflammatory diseases including rheumatoid arthritis and systemic lupus erythematosus. In this follow-up study, we saturated the region with approximately 700 SNPs (average spacing of 10 kb per SNP) in search of disease-associated variation within this region. We found preliminary evidence to suggest that common variation within the RGS7 locus may be involved in disease susceptibility.

Polymorphisms of the tumor suppressor gene LSAMP are associated with left main coronary artery disease.

Previous association mapping on chromosome 3q13-21 detected evidence for association at the limbic system-associated membrane protein (LSAMP) gene in individuals with late-onset coronary artery disease (CAD). LSAMP has never been implicated in the pathogenesis of CAD. We sought to thoroughly characterize the association and the gene. Non-redundant single nucleotide polymorphisms (SNPs) across the gene were examined in an initial dataset (168 cases with late-onset CAD, 149 controls). Stratification analysis on left main CAD (N = 102) revealed stronger association, which was further validated in a validation dataset (141 cases with left main CAD, 215 controls), a third control dataset (N = 255), and a family-based dataset (N = 2954). A haplotype residing in a novel alternative transcript of the LSAMP gene was significant in all independent case-control datasets (p = 0.0001 to 0.0205) and highly significant in the joint analysis (p = 0.00004). Lower expression of the novel alternative transcript was associated with the risk haplotype (p = 0.0002) and atherosclerosis burden in human aortas (p = 0.0001). Furthermore, silencing LSAMP expression in human aortic smooth muscle cells (SMCs) substantially augmented SMC proliferation (p<0.01). Therefore, the risk conferred by the LSAMP haplotype appears to be mediated by LSAMP down-regulation, which may promote SMC proliferation in the arterial wall and progression of atherosclerosis.

1p36 is a preferential target of chromosome 1 deletions in astrocytic tumours and homozygously deleted in a subset of glioblastomas.

Astrocytic, oligodendroglial and mixed gliomas are the commonest gliomas in adults. They have distinct phenotypes and clinical courses, but as they exist as a continuous histological spectrum, differentiating them can be difficult. Co-deletions of total 1p and 19q are found in the majority of oligodendrogliomas and considered as a diagnostic marker and a prognostic indicator. The 1p status of astrocytomas has not yet been thoroughly examined. Using a chromosome 1 tile path array, we investigated 108 adult astrocytic tumours for copy number alterations. Total 1p deletions were rare (2%), however partial deletions involving 1p36 were frequently identified in anaplastic astrocytomas (22%) and glioblastomas (34%). Multivariate analysis showed that patients with total 1p deletions had significantly longer survival (P=0.005). In nine glioblastomas homozygous deletions at 1p36 were identified. No somatic mutations were found among the five genes located in the homozygously deleted region. However, the CpG island of TNFRSF9 was hypermethylated in 19% of astrocytic tumours and 87% of glioma cell lines. TNFRSF9 expression was upregulated after demethylation of glioma cell lines. Akt3 amplifications were found in four glioblastomas. Our results indicate that 1p deletions are common anaplastic astrocytomas and glioblastomas but are distinct from the 1p abnormalities in oligodendrogliomas.

Array CGH profiling of favourable histology Wilms tumours reveals novel gains and losses associated with relapse.

Despite the excellent survival of Wilms tumour patients treated with multimodality therapy, approximately 15% will suffer from tumour relapse, where response rates are markedly reduced. We have carried out microarray-based comparative genomic hybridisation on a series of 76 Wilms tumour samples, enriched for cases which recurred, to identify changes in DNA copy number associated with clinical outcome. Using 1Mb-spaced genome-wide BAC arrays, the most significantly different genomic changes between favourable histology tumours that did (n = 37), and did not (n = 39), subsequently relapse were gains on 1q, and novel deletions at 12q24 and 18q21. Further relapse-associated loci included losses at 1q32.1, 2q36.3-2q37.1, and gain at 13q31. 1q gains correlated strongly with loss of 1p and/or 16q. In 3 of 11 cases with concurrent 1p(-)/1q(+), a breakpoint was identified at 1p13. Multiple low-level sub-megabase gains along the length of 1q were identified using chromosome 1 tiling-path arrays. One such recurrent region at 1q22-q23.1 included candidate genes RAB25, NES, CRABP2, HDGF and NTRK1, which were screened for mRNA expression using quantitative RT-PCR. These data provide a high-resolution catalogue of genomic copy number changes in relapsing favourable histology Wilms tumours.

Allelic association of sequence variants in the herpes virus entry mediator-B gene (PVRL2) with the severity of multiple sclerosis.

Discrepant findings have been reported regarding an association of the apolipoprotein E (APOE) gene with the clinical course of multiple sclerosis (MS). To resolve these discrepancies, we examined common sequence variation in six candidate genes residing in a 380-kb genomic region surrounding and including the APOE locus for an association with MS severity. We genotyped at least three polymorphisms in each of six candidate genes in 1,540 Caucasian MS families (729 single-case and multiple-case families from the United States, 811 single-case families from the UK). By applying the quantitative transmission/disequilibrium test to a recently proposed MS severity score, the only statistically significant (P=0.003) association with MS severity was found for an intronic variant in the Herpes Virus Entry Mediator-B Gene PVRL2. Additional genotyping extended the association to a 16.6 kb block spanning intron 1 to intron 2 of the gene. Sequencing of PVRL2 failed to identify variants with an obvious functional role. In conclusion, the analysis of a very large data set suggests that genetic polymorphisms in PVRL2 may influence MS severity and supports the possibility that viral factors may contribute to the clinical course of MS, consistent with previous reports.

The DNA sequence and biological annotation of human chromosome 1.

The reference sequence for each human chromosome provides the framework for understanding genome function, variation and evolution. Here we report the finished sequence and biological annotation of human chromosome 1. Chromosome 1 is gene-dense, with 3,141 genes and 991 pseudogenes, and many coding sequences overlap. Rearrangements and mutations of chromosome 1 are prevalent in cancer and many other diseases. Patterns of sequence variation reveal signals of recent selection in specific genes that may contribute to human fitness, and also in regions where no function is evident. Fine-scale recombination occurs in hotspots of varying intensity along the sequence, and is enriched near genes. These and other studies of human biology and disease encoded within chromosome 1 are made possible with the highly accurate annotated sequence, as part of the completed set of chromosome sequences that comprise the reference human genome.

Examination of seven candidate regions for multiple sclerosis: strong evidence of linkage to chromosome 1q44.

Multiple sclerosis (MS) is a debilitating neuroimmunological and neurodegenerative disease with a strong genetic component. Numerous studies have failed to consistently identify genes that confer disease susceptibility except for association with HLA-DR. Seven non-HLA regions (1q, 2q, 9q, 13q, 16q, 18p and 19q) identified in a recent genomic screen were investigated by genotyping approximately 20 single-nucleotide polymorphisms (SNPs) at approximately 1 Mb intervals. Non-parametric multipoint analyses identified a peak LOD* score of 2.99 for the 1q44 region and substantially narrowed the linkage peak to approximately 7 Mb. Ordered subset analyses (OSA) identified significant LOD score increases for 2q35 and 18p11 when ranking families by HLA-DR status and identified a significant LOD score increase in region 2q35 when ranking families by linkage to chromosome 1q44. 1q44 is particularly interesting because of linkage evidence for this region in studies of both rheumatoid arthritis and systemic lupus erythematosus.

Matroshka and ectopic polymorphisms: Two new classes of DNA sequence variation identified at the Van der Woude syndrome locus on 1q32-q41.

Van der Woude syndrome (VWS) is an orofacial clefting disorder with an autosomal dominant pattern of inheritance. In our efforts to clone the VWS gene, 900 kb of genomic sequence from the VWS candidate region at chromosome 1q32-q41 was analyzed for new DNA sequence variants. We observed that in clone CTA-321i20 a 7922 bp sequence is absent relative to the sequence present in PAC clone RP4-782d21 at positions 1669-9590, suggesting the presence of a deletion/insertion (del/ins) polymorphism. Embedded in this 7922 bp region was a TTCC short tandem repeat (STR). Genotype analysis showed that both the internal STR and the (del/ins) mutation were true polymorphisms. This is a novel example of intraallelic variation, a polymorphism within a polymorphism, and we suggest that it be termed a "Matroshka" polymorphism. Further genetic and DNA sequence analysis indicated that the ancestral state of the 1669-9590 del/ins polymorphism was the insertion allele and that the original deletion mutation probably occurred only once. A second class of novel DNA sequence variation was discovered on chromosome 5 that shared a 328 bp identical sequence with this region on chromosome 1. A single nucleotide polymorphism (SNP) was detected by SSCP using a pair of primers derived from the chromosome 1 sequence. Surprisingly, these primers also amplified the identical locus on chromosome 5, and the SNP was only located on chromosome 5. Since the probe unexpectedly detected alleles from another locus, we suggest that this type of sequence variant be termed an "ectopic" polymorphism. These two novel classes of DNA sequence polymorphisms have the potential to confound genetic and DNA sequence analysis and may also contribute to variation in disease phenotypes.

An SSLP marker-anchored BAC framework map of the mouse genome.

We have constructed a BAC framework map of the mouse genome consisting of 2,808 PCR-confirmed BAC clusters, using a previously described method. Fingerprints of BACs from selected clusters confirm the accuracy of the map. Combined with BAC fingerprint data, the framework map covers 37% of the mouse genome.

Detailed molecular analysis of 1p36 in neuroblastoma.

BACKGROUND: Several lines of evidence es tablish that chromosome band 1p36 is frequently deleted in neuroblastoma primary tumors and cell lines, suggesting that a tumor suppressor gene within this region is involved in the development of this tumor. PROCEDURE: We analyzed the status of 1p36 in primary neuroblastomas and cell lines to define the region of consistent rearrangement. RESULTS: Loss of heterozygosity (LOH) studies of primary neuro blastomas identified allelic loss in 135 of 503 tumors (27%), with the smallest region of overlap (SRO) defined distal to D15214 (1p36.3). No homozygous deletions were detected at 120 loci mapping to 1p36.1-p36.3 in a panel of 46 neuroblastoma cell lines. A recently identified patient with neuroblastoma was found to have a constitutional deletion within 1p36.2-p36.3, and this deletion, when combined with the LOH results, defined a smaller SRO of one megabase within 1p36.3. We constructed a comprehensive integrated map of chromosome 1 containing 11,000 markers and large-insert clones, a high-resolution radiation hybrid (RH) map of 1p36, and a P1-artificial chromosome (PAC) contig spanning the SRO, to further characterize the region of interest. Over 768 kb (75%) of the SRO has been sequenced to completion. Further analysis of distal 1p identified 113 transcripts localizing to 1p36, 21 of which were mapped within the SRO. CONCLUSION: This analysis will identify suitable positional candidate transcripts for mutational screening and subsequent identification of the 1p36.3 neuroblastoma suppressor gene.

A physical map of the human genome.

The human genome is by far the largest genome to be sequenced, and its size and complexity present many challenges for sequence assembly. The International Human Genome Sequencing Consortium constructed a map of the whole genome to enable the selection of clones for sequencing and for the accurate assembly of the genome sequence. Here we report the construction of the whole-genome bacterial artificial chromosome (BAC) map and its integration with previous landmark maps and information from mapping efforts focused on specific chromosomal regions. We also describe the integration of sequence data with the map.

The physical maps for sequencing human chromosomes 1, 6, 9, 10, 13, 20 and X.

We constructed maps for eight chromosomes (1, 6, 9, 10, 13, 20, X and (previously) 22), representing one-third of the genome, by building landmark maps, isolating bacterial clones and assembling contigs. By this approach, we could establish the long-range organization of the maps early in the project, and all contig extension, gap closure and problem-solving was simplified by containment within local regions. The maps currently represent more than 94% of the euchromatic (gene-containing) regions of these chromosomes in 176 contigs, and contain 96% of the chromosome-specific markers in the human gene map. By measuring the remaining gaps, we can assess chromosome length and coverage in sequenced clones.

A 6-Mb high-resolution physical and transcription map encompassing the hereditary prostate cancer 1 (HPC1) region.

Several hereditary disease loci have been genetically mapped to the chromosome 1q24-q31 interval, including the hereditary prostate cancer 1 (HPC1) locus. Here, we report the construction of a 20-Mb yeast artificial chromosome contig and a high-resolution 6-Mb sequence-ready bacterial artificial chromosome (BAC)/P1-derived artificial chromosome (PAC) contig of 1q25 by sequence and computational analysis, STS content mapping, and chromosome walking. One hundred thirty-six new STSs, including 10 novel simple sequence repeat polymorphisms that are being used for genetic refinement of multiple disease loci, have been generated from this contig and are shown to map to the 1q25 interval. The integrity of the 6-Mb BAC/PAC contig has been confirmed by restriction fingerprinting, and this contig is being used as a template for human chromosome 1 genome sequencing. A transcription mapping effort has resulted in the precise localization of 18 known genes and 31 ESTs by database searching, exon trapping, direct cDNA hybridization, and sample sequencing of BACs from the 1q25 contig. An additional 11 known genes and ESTs have been placed within the larger 1q24-q31 interval. These transcription units represent candidate genes for multiple hereditary diseases, including HPC1.

A preliminary gene map for the Van der Woude syndrome critical region derived from 900 kb of genomic sequence at 1q32-q41.

Van der Woude syndrome (VWS) is a common form of syndromic cleft lip and palate and accounts for approximately 2% of all cleft lip and palate cases. Distinguishing characteristics include cleft lip with or without cleft palate, isolated cleft palate, bilateral lip pits, hypodontia, normal intelligence, and an autosomal-dominant mode of transmission with a high degree of penetrance. Previously, the VWS locus was mapped to a 1.6-cM region in 1q32-q41 between D1S491 and D1S205, and a 4.4-Mb contig of YAC clones of this region was constructed. In the current investigation, gene-based and anonymous STSs were developed from the existing physical map and were then used to construct a contig of sequence-ready bacterial clones across the entire VWS critical region. All STSs and BAC clones were shared with the Sanger Centre, which developed a contig of PAC clones over the same region. A subset of 11 clones from both contigs was selected for high-throughput sequence analysis across the approximately 1.1-Mb region; all but two of these clones have been sequenced completely. Over 900 kb of genomic sequence, including the 350-kb VWS critical region, were analyzed and revealed novel polymorphisms, including an 8-kb deletion/insertion, and revealed 4 known genes, 11 novel genes, 9 putative genes, and 3 psuedogenes. The positional candidates LAMB3, G0S2, HIRF6, and HSD11 were excluded as the VWS gene by mutation analysis. A preliminary gene map for the VWS critical region is as follows: [see text] 41-TEL. The data provided here will help lead to the identification of the VWS gene, and this study provides a model for how laboratories that have a regional interest in the human genome can contribute to the sequencing efforts of the entire human genome.

Genome mapping by fluorescent fingerprinting.

The construction of sequence-ready maps of overlapping genomic clones is central to large-scale genome sequencing. We have implemented a method for fluorescent fingerprinting of bacterial clones to assemble contig maps. The method utilizes three spectrally distinct fluorescently tagged dideoxy ATPs to specifically label the HindIII termini in HindIII and Sau3AI restriction digests of clones that are multiplexed prior to electrophoresis and data collection. There is excellent reproducibility of raw data, improved resolution of large fragments, and concordance between the results obtained using this and the equivalent radioactive protocol. This method also allows detection of smaller overlaps between clones when compared to the analysis of restriction digests on nondenaturing agarose gels.