Robust SNP-based prediction of rheumatoid arthritis through machine-learning-optimized polygenic risk score.

BACKGROUND: The popular statistics-based Genome-wide association studies (GWAS) have provided deep insights into the field of complex disorder genetics. However, its clinical applicability to predict disease/trait outcomes remains unclear as statistical models are not designed to make predictions. This study employs statistics-free machine-learning (ML)-optimized polygenic risk score (PRS) to complement existing GWAS and bring the prediction of disease/trait outcomes closer to clinical application. Rheumatoid Arthritis (RA) was selected as a model disease to demonstrate the robustness of ML in disease prediction as RA is a prevalent chronic inflammatory joint disease with high mortality rates, affecting adults at the economic prime. Early identification of at-risk individuals may facilitate measures to mitigate the effects of the disease. METHODS: This study employs a robust ML feature selection algorithm to identify single nucleotide polymorphisms (SNPs) that can predict RA from a set of training data comprising RA patients and population control samples. Thereafter, selected SNPs were evaluated for their predictive performances across 3 independent, unseen test datasets. The selected SNPs were subsequently used to generate PRS which was also evaluated for its predictive capacity as a sole feature. RESULTS: Through robust ML feature selection, 9 SNPs were found to be the minimum number of features for excellent predictive performance (AUC > 0.9) in 3 independent, unseen test datasets. PRS based on these 9 SNPs was significantly associated with (P < 1 × 10-16) and predictive (AUC > 0.9) of RA in the 3 unseen datasets. A RA ML-PRS calculator of these 9 SNPs was developed ( https://xistance.shinyapps.io/prs-ra/ ) to facilitate individualized clinical applicability. The majority of the predictive SNPs are protective, reside in non-coding regions, and are either predicted to be potentially functional SNPs (pfSNPs) or in high linkage disequilibrium (r2 > 0.8) with un-interrogated pfSNPs. CONCLUSIONS: These findings highlight the promise of this ML strategy to identify useful genetic features that can robustly predict disease and amenable to translation for clinical application.

Third-Generation Single-Molecule Sequencing for Preimplantation Genetic Testing of Aneuploidy and Segmental Imbalances.

BACKGROUND: Current strategies for preimplantation genetic testing for aneuploidy or structural rearrangements (PGT-A/SR) rely mainly on next-generation sequencing (NGS) and microarray platforms, which are robust but require expensive instrumentation. We explored the suitability of third-generation single-molecule sequencing as a PGT-A/SR screening platform for both aneuploidy and segmental imbalance. METHODS: Single-cell and multicell replicates from aneuploid or segmentally unbalanced cell lines (n = 208) were SurePlex-amplified, randomized, and subjected to (a) Nanopore-based single-molecule sequencing (Oxford Nanopore Technologies) and (b) NGS using a leading commercial PGT-A solution (Illumina VeriSeq PGS). Archival SurePlex-amplified trophectoderm biopsy samples (n = 96) previously analyzed using the commercial kit were blinded and reanalyzed using Nanopore. RESULTS: Nanopore-based PGT-A identified the specific aberration in 95.45% (84/88) and 97.78% (88/90) of single-/multicells with an aneuploidy or segmental imbalance (10-30.5 Mb), respectively. Comparison against the commercial kit's results revealed concordances of 98.86% (87/88) and 98.89% (89/90) for the aneuploid and segmentally unbalanced (10-30.5 Mb aberration) samples, respectively. Detection sensitivity for smaller segmental imbalances (5-5.8 Mb aberration, n = 30) decreased markedly on both platforms. Nanopore-based PGT-A reanalysis of trophectoderm biopsy samples was 97.92% (94/96) concordant with the commercial kit results. CONCLUSION: Up to 24 SurePlex-amplified single-cell, multicell, or trophectoderm samples could be sequenced in a single MinION flow-cell for subsequent preimplantation genetic testing for aneuploidy or structural rearrangements (PGT-A/SR) analysis, with results obtainable in ≤3 days and at per-sample costs that are competitive with commercial offerings. Nanopore's third-generation single-molecule sequencing represents a viable alternative to current commercial NGS-based PGT-A solutions for aneuploidy and segmental imbalance (≥10 Mb) screening of single-/multicell or trophectoderm biopsy samples.

Robust preimplantation genetic testing of the common F8 Inv22 pathogenic variant of severe hemophilia A using a highly polymorphic multi-marker panel encompassing the paracentric inversion.

BACKGROUND: Hemophilia A (HEMA) is an X-linked bleeding disorder caused by reduced/absent coagulation factor VIII expression, as a result of pathogenic variants in the F8 gene. Preimplantation prevention of HEMA should ideally include direct pathogenic F8 variant detection, complemented by linkage analysis of flanking markers to identify the high-risk F8 allele. Linkage analysis is particularly indispensable when the pathogenic variant cannot be detected directly or identified. This study evaluated the suitability of a panel of F8 intragenic and extragenic short tandem repeat markers for standalone linkage-based preimplantation genetic testing for monogenic disorder (PGT-M) of the Inv22 pathogenic variant, an almost 600 kb paracentric inversion responsible for almost half of all severe HEMA globally, for which direct detection is challenging. METHODS: Thirteen markers spanning 1 Mb and encompassing both F8 and the Inv22 inversion interval were genotyped in 153 unrelated females of Viet Kinh ethnicity. RESULTS: All individuals were heterozygous for ≥ 1 marker, ~ 90% were heterozygous for ≥ 1 of the five F8 intragenic markers, and almost 98% were heterozygous for ≥ 1 upstream (telomeric) and ≥ 1 downstream (centromeric) markers. A prospective PGT-M couple at risk of transmitting F8 Inv22 were fully informative at four marker loci (2 intra-inversion, 1 centromeric, 1 telomeric) and partially informative at another five (2 intra-inversion, 3 centromeric), allowing robust phasing of low- and high-risk haplotypes. In vitro fertilization produced three embryos, all of which clearly inherited the low-risk maternal allele, enabling reliable unaffected diagnoses. A single embryo transfer produced a clinical pregnancy, which was confirmed as unaffected by amniocentesis and long-range PCR, and a healthy baby girl was delivered at term. CONCLUSION: Robust and reliable PGT-M of HEMA, including the common F8 Inv22 pathogenic variant, can be achieved with sufficient informative intragenic and flanking markers.

Single-Tube Screen for Rapid Detection of Repeat Expansions in Seven Common Spinocerebellar Ataxias.

BACKGROUND: The autosomal dominantly inherited and genetically heterogeneous spinocerebellar ataxias (SCAs) exhibit highly similar clinical presentations. Many are caused by repeat expansions, of which at least 8 involve CAG repeats. Repeat expansion detection is the only method to confirm disease status in symptomatic individuals. We present a novel strategy to simultaneously screen for the presence of CAG repeat expansion in the genes responsible for SCAs 1, 2, 3, 6, 7, 12, and dentatorubral-pallidoluysian atrophy using a simplified single-tube assay. METHODS: The method employs differentially labeled locus-specific primers and a common triplet-primed primer. Amplified products from each locus are distinguished by a combination of the product size and the fluorophore tag. The upper size limit of the normal allele range was used as the cutoff for distinguishing normal from potentially affected samples, with repeat expansion detected by presence of electrophoretic peaks extending beyond the cutoff. RESULTS: Blinded evaluation of the assay on 60 genotype-known DNA samples correctly detected repeat expansion in the expected SCA repeat locus for all 31 DNA samples. CONCLUSIONS: In principle, this strategy can be applied to the simultaneous screening of any group of disease genes sharing the same repetitive units and/or their reverse complement.

Machine learning using genetic and clinical data identifies a signature that robustly predicts methotrexate response in rheumatoid arthritis.

OBJECTIVE: To develop a hypothesis-free model that best predicts response to MTX drug in RA patients utilizing biologically meaningful genetic feature selection of potentially functional single nucleotide polymorphisms (pfSNPs) through robust machine learning (ML) feature selection methods. METHODS: MTX-treated RA patients with known response were divided in a 4:1 ratio into training and test sets. From the patients' exomes, potential features for classifier prediction were identified from pfSNPs and non-genetic factors through ML using recursive feature elimination with cross-validation incorporating the random forest classifier. Feature selection was repeated on random subsets of the training cohort, and consensus features were assembled into the final feature set. This feature set was evaluated for predictive potential using six ML classifiers, first by cross-validation within the training set, and finally by analysing its performance with the unseen test set. RESULTS: The final feature set contains 56 pfSNPs and five non-genetic factors. The majority of these pfSNPs are located in pathways related to RA pathogenesis or MTX action and are predicted to modulate gene expression. When used for training in six ML classifiers, performance was good in both the training set (area under the curve: 0.855-0.916; sensitivity: 0.715-0.892; and specificity: 0.733-0.862) and the unseen test set (area under the curve: 0.751-0.826; sensitivity: 0.581-0.839; and specificity: 0.641-0.923). CONCLUSION: Sensitive and specific predictors of MTX response in RA patients were identified in this study through a novel strategy combining biologically meaningful and machine learning feature selection and training. These predictors may facilitate better treatment decision-making in RA management.

Towards precision medicine: interrogating the human genome to identify drug pathways associated with potentially functional, population-differentiated polymorphisms.

Drug response variations amongst different individuals/populations are influenced by several factors including allele frequency differences of single nucleotide polymorphisms (SNPs) that functionally affect drug-response genes. Here, we aim to identify drugs that potentially exhibit population differences in response using SNP data mining and analytics. Ninety-one pairwise-comparisons of >22,000,000 SNPs from the 1000 Genomes Project, across 14 different populations, were performed to identify 'population-differentiated' SNPs (pdSNPs). Potentially-functional pdSNPs (pf-pdSNPs) were then selected, mapped into genes, and integrated with drug-gene databases to identify 'population-differentiated' drugs enriched with genes carrying pf-pdSNPs. 1191 clinically-approved drugs were found to be significantly enriched (Z > 2.58) with genes carrying SNPs that were differentiated in one or more population-pair comparisons. Thirteen drugs were found to be enriched with such differentiated genes across all 91 population-pairs. Notably, 82% of drugs, which were previously reported in the literature to exhibit population differences in response were also found by this method to contain a significant enrichment of population specific differentiated SNPs. Furthermore, drugs with genetic testing labels, or those suspected to cause adverse reactions, contained a significantly larger number (P < 0.01) of population-pairs with enriched pf-pdSNPs compared with those without these labels. This pioneering effort at harnessing big-data pharmacogenomics to identify 'population differentiated' drugs could help to facilitate data-driven decision-making for a more personalized medicine.

Comprehensive analysis of transcriptome profiles in hepatocellular carcinoma.

BACKGROUND: Hepatocellular carcinoma is the second most deadly cancer with late presentation and limited treatment options, highlighting an urgent need to better understand HCC to facilitate the identification of early-stage biomarkers and uncover therapeutic targets for the development of novel therapies for HCC. METHODS: Deep transcriptome sequencing of tumor and paired non-tumor liver tissues was performed to comprehensively evaluate the profiles of both the host and HBV transcripts in HCC patients. Differential gene expression patterns and the dys-regulated genes associated with clinical outcomes were analyzed. Somatic mutations were identified from the sequencing data and the deleterious mutations were predicted. Lastly, human-HBV chimeric transcripts were identified, and their distribution, potential function and expression association were analyzed. RESULTS: Expression profiling identified the significantly upregulated TP73 as a nodal molecule modulating expression of apoptotic genes. Approximately 2.5% of dysregulated genes significantly correlated with HCC clinical characteristics. Of the 110 identified genes, those involved in post-translational modification, cell division and/or transcriptional regulation were upregulated, while those involved in redox reactions were downregulated in tumors of patients with poor prognosis. Mutation signature analysis identified that somatic mutations in HCC tumors were mainly non-synonymous, frequently affecting genes in the micro-environment and cancer pathways. Recurrent mutations occur mainly in ribosomal genes. The most frequently mutated genes were generally associated with a poorer clinical prognosis. Lastly, transcriptome sequencing suggest that HBV replication in the tumors of HCC patients is rare. HBV-human fusion transcripts are a common observation, with favored HBV and host insertion sites being the HBx C-terminus and gene introns (in tumors) and introns/intergenic-regions (in non-tumors), respectively. HBV-fused genes in tumors were mainly involved in RNA binding while those in non-tumors tissues varied widely. These observations suggest that while HBV may integrate randomly during chronic infection, selective expression of functional chimeric transcripts may occur during tumorigenesis. CONCLUSIONS: Transcriptome sequencing of HCC patients reveals key cancer molecules and clinically relevant pathways deregulated/mutated in HCC patients and suggests that while HBV may integrate randomly during chronic infection, selective expression of functional chimeric transcripts likely occur during the process of tumorigenesis.

Architecture of polymorphisms in the human genome reveals functionally important and positively selected variants in immune response and drug transporter genes.

BACKGROUND: Genetic polymorphisms can contribute to phenotypic differences amongst individuals, including disease risk and drug response. Characterization of genetic polymorphisms that modulate gene expression and/or protein function may facilitate the identification of the causal variants. Here, we present the architecture of genetic polymorphisms in the human genome focusing on those predicted to be potentially functional/under natural selection and the pathways that they reside. RESULTS: In the human genome, polymorphisms that directly affect protein sequences and potentially affect function are the most constrained variants with the lowest single-nucleotide variant (SNV) density, least population differentiation and most significant enrichment of rare alleles. SNVs which potentially alter various regulatory sites, e.g. splicing regulatory elements, are also generally under negative selection. Interestingly, genes that regulate the expression of transcription/splicing factors and histones are conserved as a higher proportion of these genes is non-polymorphic, contain ultra-conserved elements (UCEs) and/or has no non-synonymous SNVs (nsSNVs)/coding INDELs. On the other hand, major histocompatibility complex (MHC) genes are the most polymorphic with SNVs potentially affecting the binding of transcription/splicing factors and microRNAs (miRNA) exhibiting recent positive selection (RPS). The drug transporter genes carry the most number of potentially deleterious nsSNVs and exhibit signatures of RPS and/or population differentiation. These observations suggest that genes that interact with the environment are highly polymorphic and targeted by RPS. CONCLUSIONS: In conclusion, selective constraints are observed in coding regions, master regulator genes, and potentially functional SNVs. In contrast, genes that modulate response to the environment are highly polymorphic and under positive selection.

FMR1 CGG repeat expansion mutation detection and linked haplotype analysis for reliable and accurate preimplantation genetic diagnosis of fragile X syndrome.

Fragile X mental retardation 1 (FMR1) full-mutation expansion causes fragile X syndrome. Trans-generational fragile X syndrome transmission can be avoided by preimplantation genetic diagnosis (PGD). We describe a robust PGD strategy that can be applied to virtually any couple at risk of transmitting fragile X syndrome. This novel strategy utilises whole-genome amplification, followed by triplet-primed polymerase chain reaction (TP-PCR) for robust detection of expanded FMR1 alleles, in parallel with linked multi-marker haplotype analysis of 13 highly polymorphic microsatellite markers located within 1 Mb of the FMR1 CGG repeat, and the AMELX/Y dimorphism for gender identification. The assay was optimised and validated on single lymphoblasts isolated from fragile X reference cell lines, and applied to a simulated PGD case and a clinical in vitro fertilisation (IVF)-PGD case. In the simulated PGD case, definitive diagnosis of the expected results was achieved for all 'embryos'. In the clinical IVF-PGD case, delivery of a healthy baby girl was achieved after transfer of an expansion-negative blastocyst. FMR1 TP-PCR reliably detects presence of expansion mutations and obviates reliance on informative normal alleles for determining expansion status in female embryos. Together with multi-marker haplotyping and gender determination, misdiagnosis and diagnostic ambiguity due to allele dropout is minimised, and couple-specific assay customisation can be avoided.

Single-Tube Dodecaplex PCR Panel of Polymorphic Microsatellite Markers Closely Linked to the DMPK CTG Repeat for Preimplantation Genetic Diagnosis of Myotonic Dystrophy Type 1.

BACKGROUND: Preimplantation genetic diagnosis (PGD) of myotonic dystrophy type 1 (DM1) currently uses conventional PCR to detect nonexpanded dystrophia myotonica protein kinase (DMPK) alleles or triplet-primed PCR to detect the CTG-expanded alleles, coupled with analysis of linked microsatellite markers to increase diagnostic accuracy. We aimed to simplify the process of identification and selection of informative linked markers for application to DM1 PGD. METHODS: An in silico search was performed to identify all markers within 1-1.5 Mb flanking the DMPK gene. Five previously known (D19S559, APOC2, D19S543, D19S112, and BV209569) and 7 novel (DM45050, DM45178, DM45209, DM45958, DM46513, DM46892, and DM47004.1) markers with potentially high heterozygosity values and polymorphism information content were selected and optimized in a single-tube multiplex PCR panel. RESULTS: Analysis of 184 DNA samples of Chinese and Caucasian individuals (91 from unrelated, anonymized cord blood of Chinese babies born at the National University Hospital, Singapore, and 93 Caucasian DNA samples from the Human Variation Panel HD100CAU) confirmed the high polymorphism indices of all markers (polymorphism information content >0.5), with observed heterozygosity values ranging from 0.62-0.93. All individuals were heterozygous for at least 6 markers, with 99.5% of individuals heterozygous for at least 2 markers on either side of the DMPK CTG repeat. The dodecaplex marker assay was successfully validated on 42 single cells and 12 whole genome amplified single cells. CONCLUSIONS: The DM1 multiplex PCR panel is suitable for use in DM1 PGD either as a standalone linkage-based assay or as a complement to DMPK CTG repeat expansion-mutation detection.

Disruption of FAT10-MAD2 binding inhibits tumor progression.

FAT10 (HLA-F-adjacent transcript 10) is a ubiquitin-like modifier that is commonly overexpressed in various tumors. It was found to play a role in mitotic regulation through its interaction with mitotic arrest-deficient 2 (MAD2). Overexpression of FAT10 promotes tumor growth and malignancy. Here, we identified the MAD2-binding interface of FAT10 to be located on its first ubiquitin-like domain whose NMR structure thus was determined. We further proceeded to demonstrate that disruption of the FAT10-MAD2 interaction through mutation of specific MAD2-binding residues did not interfere with the interaction of FAT10 with its other known interacting partners. Significantly, ablation of the FAT10-MAD2 interaction dramatically limited the promalignant capacity of FAT10, including promoting tumor growth in vivo and inducing aneuploidy, proliferation, migration, invasion, and resistance to apoptosis in vitro. Our results strongly suggest that the interaction of FAT10 with MAD2 is a key mechanism underlying the promalignant property of FAT10 and offer prospects for the development of anticancer strategies.

Identification of microsatellite markers <1 Mb from the FMR1 CGG repeat and development of a single-tube tetradecaplex PCR panel of highly polymorphic markers for preimplantation genetic diagnosis of fragile X syndrome.

PURPOSE: To develop a single-tube polymerase chain reaction (PCR) panel of highly polymorphic markers for preimplantation genetic diagnosis (PGD) of fragile X syndrome (FXS). METHODS: An in silico search was performed to identify all markers within 1 Mb flanking the FMR1 gene. Selected markers were optimized into a single-tube PCR panel and their polymorphism indices were determined from 272 female samples from three populations. The single-tube assay was also validated on 30 single cells to evaluate its applicability to FXS PGD. RESULTS: Thirteen markers with potentially high polymorphism information content (PIC) and heterozygosity values were selected and optimized into a single-tube PCR panel together with AMELX/Y for gender determination. Analysis of 272 female samples confirmed the high polymorphism (PIC > 0.5) of most markers, with expected and observed heterozygosities ranging from 0.31 to 0.87. More than 99% of individuals were heterozygous for at least three markers, with 95.8% of individuals heterozygous for at least two markers on either side of the FMR1 CGG repeat. CONCLUSION: The tetradecaplex marker assay can be performed directly on single cells or after whole-genome amplification, thus supporting its use in FXS PGD either as a standalone linkage-based assay or as a complement to FMR1 mutation detection.Genet Med 18 9, 869-875.

Identification of Novel Microsatellite Markers <1 Mb from the HTT CAG Repeat and Development of a Single-Tube Tridecaplex PCR Panel of Highly Polymorphic Markers for Preimplantation Genetic Diagnosis of Huntington Disease.

BACKGROUND: Preimplantation genetic diagnosis (PGD) of Huntington disease (HD) generally employs linkage analysis of flanking microsatellite markers to complement direct mutation testing, as well as for exclusion testing. Thus far, only 10 linked markers have been developed for use in HD PGD, with a maximum of 3 markers coamplified successfully. We aimed to develop a single-tube multiplex PCR panel of highly polymorphic markers to simplify HD PGD. METHODS: An in silico search was performed to identify all markers within 1 Mb flanking the huntingtin (HTT) gene. Selected markers were optimized in a single-tube PCR panel, and their polymorphism indices were determined in 2 populations. The panel was tested on 63 single cells to validate its utility in PGD. RESULTS: We identified 102 markers in silico, of which 56 satisfied the selection criteria. After initial testing, 12 markers with potentially high heterozygosity were optimized into a single-tube PCR panel together with a 13th more distally located marker. Analysis of DNA from 183 Chinese and Caucasian individuals revealed high polymorphism indices for all markers (polymorphism information content >0.5), with observed heterozygosities ranging from 0.5-0.92. All individuals were heterozygous for at least 5 markers, with 99.5% of individuals heterozygous for at least 2 markers upstream and downstream of the HTT CAG repeat. CONCLUSIONS: The tridecaplex marker assay amplified reliably from single cells either directly or after whole genome amplification, thus validating its standalone use in HD exclusion PGD or as a complement to HTT CAG repeat expansion-mutation detection.

Enhanced Detection and Sizing of the HTT CAG Repeat Expansion in Huntington Disease Using an Improved Triplet-Primed PCR Assay.

BACKGROUND: Accurate determination of the CAG repeat number is crucial for diagnostic and predictive testing for Huntington disease (HD). Current PCR-based assays can accurately size up to ∼110 HTT CAG repeats. OBJECTIVE: To develop an improved assay capable of detecting larger CAG repeat expansions. METHODS: A triplet-primed PCR (TP-PCR) assay was optimized and validated on 14 HD reference DNAs, including a sample carrying a large expansion of ∼180 CAG repeats. RESULTS: All 14 HD reference samples showed 100% concordance with the previously verified allele sizes. For alleles under 45 CAGs, identical repeat sizes were obtained, while alleles larger than 46 CAGs were sized to within ±1 CAG. The improved TP-PCR assay successfully detected the ∼180 CAG repeat allele in an affected sample. CONCLUSION: This method extends the detection limit of large expanded alleles to at least ∼175-180 CAG repeats, thus reducing the likelihood of requiring Southern blot analysis for any HD-affected sample.

A Patient with ß-Thalassemia Intermedia Secondary to Homozygosity for a Polyadenylation Signal Mutation (AATAAA > AATAGA) (HBB: C.*112A > G) on the ß-Globin Gene.

We describe the clinical presentation and laboratory findings of a Malay man with ß-thalassemia intermedia (ß-TI), secondary to homozygosity for a polyadenylation (polyA) signal mutation (AATAAA > AATAGA) (HBB: c.*112A > G) on the ß-globin gene, and give a brief review of the literature. This is the first report of a homozygous case of this polyA mutation, and highlights the importance of molecular analysis of the globin genes in the diagnosis of thalassemia.

Simplified strategy for rapid first-line screening of fragile X syndrome: closed-tube triplet-primed PCR and amplicon melt peak analysis.

Premutation and full-mutation hyperexpansion of CGG-triplets in the X-linked Fragile X Mental Retardation 1 (FMR1) gene have been implicated in fragile X-associated tremor/ataxia syndrome, fragile X-associated primary ovarian insufficiency, and fragile X syndrome (FXS), respectively. The currently available molecular diagnostic tests are either costly or labour-intensive, which prohibits their application as a first-line FMR1 test in large-scale population-based screening programs. In this study, we demonstrate the utility of a simplified closed-tube strategy for rapid first-line screening of FXS based on melt peak temperature (Tm) analysis of direct triplet-primed polymerase chain reaction amplicons (dTP-PCR MCA). In addition, we also evaluated the correlation between Tm and CGG-repeat size based on capillary electrophoresis (CE) of dTP-PCR amplicons. The assays were initially tested on 29 FMR1 reference DNA samples, followed by a blinded validation on 107 previously characterised patient DNA samples. The dTP-PCR MCA produced distinct melt profiles of higher Tm for samples carrying an expanded allele. Among the samples tested, we also observed a good correlation between Tm and CGG-repeat size. In the blinded validation study, dTP-PCR MCA accurately classified all normal and expansion carriers, and the FMR1 genotypic classification of all samples was completely concordant with the previously determined genotypes as well as the dTP-PCR CE results. This simple and cost-effective MCA-based assay may be useful as a first-line FXS screening tool that could rapidly screen out the large majority of unaffected individuals, thus minimising the number of samples that need to be analysed by Southern blot analysis.

Identification of novel microsatellite markers <1 Mb from the HBB gene and development of a single-tube pentadecaplex PCR panel of highly polymorphic markers for preimplantation genetic diagnosis of beta-thalassemia.

Beta (ß)-thalassemia is one of the most common monogenic diseases worldwide. Affected pregnancies can be avoided through preimplantation genetic diagnosis (PGD), which commonly involves customized assays to detect the different combinations of ß-globin (HBB) gene mutations present in couples, in conjunction with linkage analysis of flanking microsatellite markers. Currently, the limited number of reported closely linked markers hampers their utility in indirect linkage-based PGD for this disorder. To increase the available markers closely flanking the HBB gene, an in silico search was performed to identify all markers within 1 Mb flanking the HBB gene. Fifteen markers with potentially high polymorphism information content (PIC) and heterozygosity values were selected and optimized into a single-tube pentadecaplex PCR panel. Allele frequencies and polymorphism and heterozygosity indices of each marker were assessed in five populations. A total of 238 alleles were observed from the 15 markers. PIC was >0.7 for all markers, with expected heterozygosity and observed heterozygosity values ranging from 0.74 to 0.90 and 0.72 to 0.88, respectively. Greater than 99% of individuals were heterozygous for at least seven markers, with at least two heterozygous markers on either side of the HBB gene. The pentadecaplex marker assay also performed reliably on single cells either directly or after whole genome amplification, thus validating its use in standalone linkage-based ß-thalassemia PGD or in conjunction with HBB mutation detection.

Novel evidence of association with nonsyndromic cleft lip with or without cleft palate was shown for single nucleotide polymorphisms in FOXF2 gene in an Asian population.

BACKGROUND: The forkhead box F2 gene (FOXF2) located in chromosome 6p25.3 has been shown to play a crucial role in palatal development in mouse and rat models. To date, no evidence of linkage or association has been reported for this gene in humans with oral clefts. METHODS: Allelic transmission disequilibrium tests were used to robustly assess evidence of linkage and association with nonsyndromic cleft lip with or without cleft palate for nine single nucleotide polymorphisms (SNPs) in and around FOXF2 in both Asian and European trios using PLINK. RESULTS: Statistically significant evidence of linkage and association was shown for two SNPs (rs1711968 and rs732835) in 216 Asian trios where the empiric P values with permutation tests were 0.0016 and 0.005, respectively. The corresponding estimated odds ratios for carrying the minor allele at these SNPs were 2.05 (95% confidence interval = 1.41, 2.98) and 1.77 (95% confidence interval = 1.26, 2.49), respectively. CONCLUSION: Our results provided statistical evidence of linkage and association between FOXF2 and nonsyndromic cleft lip with or without cleft palate.

Single-tube nonaplex microsatellite PCR panel for preimplantation genetic diagnosis of Hb Bart's hydrops fetalis syndrome.

OBJECTIVE: To develop a single-tube multi-marker assay for improved preimplantation genetic diagnosis (PGD) of deletional and/or non-deletional Hb Bart's hydrops fetalis syndrome, providing haplotype confirmation of deletional status, and maximization of linkage informativity. METHODS: We performed in silico mining to identify novel microsatellites within 1 Mb flanking the alpha-globin gene cluster, and optimized a single-tube assay combining detection of α(0) -thalassemia deletions with multi-marker linkage analysis. We performed validation on 100 single cells prior to clinical PGD application. RESULTS: Of 42 markers encompassing the α-globin gene cluster that were identified in silico, 9 were highly polymorphic (0.68 ≤ polymorphism information content ≤ 0.92; 0.66 ≤ Ho ≤ 0.90; 10 ≤ alleles ≤ 35) and optimized to co-amplify directly from a single cell. A validation analysis of 100 single lymphoblasts yielded 100% amplification success for all markers, and individual marker allele drop-out (ADO) rates of 0-5%. Clinical application of the assay in PGD for Hb Bart's (2 cases/cycles) resulted in a twin pregnancy and healthy live birth of two baby girls. CONCLUSIONS: This single-tube nonaplex microsatellite PCR panel can be applied directly to PGD of most deletional Hb Bart's without the need for deletion-specific customization, and to linkage-based PGD of non-deletional Hb Bart's.