Validation and genotyping of multiple human polymorphic inversions mediated by inverted repeats reveals a high degree of recurrence.

In recent years different types of structural variants (SVs) have been discovered in the human genome and their functional impact has become increasingly clear. Inversions, however, are poorly characterized and more difficult to study, especially those mediated by inverted repeats or segmental duplications. Here, we describe the results of a simple and fast inverse PCR (iPCR) protocol for high-throughput genotyping of a wide variety of inversions using a small amount of DNA. In particular, we analyzed 22 inversions predicted in humans ranging from 5.1 kb to 226 kb and mediated by inverted repeat sequences of 1.6-24 kb. First, we validated 17 of the 22 inversions in a panel of nine HapMap individuals from different populations, and we genotyped them in 68 additional individuals of European origin, with correct genetic transmission in ∼ 12 mother-father-child trios. Global inversion minor allele frequency varied between 1% and 49% and inversion genotypes were consistent with Hardy-Weinberg equilibrium. By analyzing the nucleotide variation and the haplotypes in these regions, we found that only four inversions have linked tag-SNPs and that in many cases there are multiple shared SNPs between standard and inverted chromosomes, suggesting an unexpected high degree of inversion recurrence during human evolution. iPCR was also used to check 16 of these inversions in four chimpanzees and two gorillas, and 10 showed both orientations either within or between species, providing additional support for their multiple origin. Finally, we have identified several inversions that include genes in the inverted or breakpoint regions, and at least one disrupts a potential coding gene. Thus, these results represent a significant advance in our understanding of inversion polymorphism in human populations and challenge the common view of a single origin of inversions, with important implications for inversion analysis in SNP-based studies.

A common 56-kilobase deletion in a primate-specific segmental duplication creates a novel butyrophilin-like protein.

BACKGROUND: The Butyrophilin-like (BTNL) proteins are likely to play an important role in inflammation and immune response. Like the B7 protein family, many human and murine BTNL members have been shown to control T lymphocytes response, and polymorphisms in human BTNL2 have been linked to several inflammatory diseases, such as pulmonary sarcoidosis, inflammatory bowel disease and neonatal lupus. RESULTS: In this study we provide a comprehensive population, genomic and transcriptomic analysis of a 56-kb deletion copy number variant (CNV), located within two segmental duplications of two genes belonging to the BTNL family, namely BTNL8 and BTNL3. We confirm the presence of a novel BTNL8*3 fusion-protein product, and show an influence of the deletion variant on the expression level of several genes involved in immune function, including BTNL9, another member of the same family. Moreover, by genotyping HapMap and human diversity panel (HGDP) samples, we demonstrate a clear difference in the stratification of the BTNL8_BTNL3-del allele frequency between major continental human populations. CONCLUSION: Despite tremendous progress in the field of structural variation, rather few CNVs have been functionally characterized so far. Here, we show clear functional consequences of a new deletion CNV (BTNL8_BTNL3-del) with potentially important implication in the human immune system and in inflammatory and proliferative disorders. In addition, the marked population differences found of BTNL8_BTNL3-del frequencies suggest that this deletion CNV might have evolved under positive selection due to environmental conditions in some populations, with potential phenotypic consequences.

MLPAstats: an R GUI package for the integrated analysis of copy number alterations using MLPA data.

BACKGROUND: Multiplex-Dependent Probe Amplification (MLPA) is a cost-effective experimental method for candidate gene studies, aimed at the identification of copy number alterations. The analysis of such genetic variants, from electropherogram peak intensities, involves two main stages. First, peak normalization for each probe is required to remove the contribution of probe size to peak intensity. Second, the statistical significance of peak alteration between case and control samples is estimated. A number of methods have been proposed in each step with varying levels of complexity and precision. However, there is no single framework from which the results of each method and possible combinations at each step can be assessed. RESULTS: We present MLPAstats, an R package designed to integrate the methods for exploring different analysis scenarios in a reliable way. A GUI has been developed to allow researchers to find their optimal analysis strategy. CONCLUSIONS: MLPAstats is an analysis tool that promotes the use of cost-effective MLPA suitable for candidate gene studies. Its R implementation allows future methods to be easily incorporated, while its GUI will facilitate its use by non-expert analysts. A vignette describing a set-by-step tutorial is also available with the package.

Novel UBE3A mutations causing Angelman syndrome: different parental origin for single nucleotide changes and multiple nucleotide deletions or insertions.

Angelman syndrome (AS) is a genetic disorder caused by a deficiency of UBE3A imprinted gene expression from the maternal chromosome 15. In 10% of AS cases the genetic cause is a mutation affecting the maternal copy of the UBE3A gene. In two large Spanish series of clinically stringently selected and nonstringently selected patients, we have identified 11 pathological mutations--eight of them novel mutations--and 14 sequence changes considered polymorphic variants. Remarkably, single nucleotide substitutions are more likely to be inherited, while multiple nucleotide deletions or insertions are less frequently inherited, thus indicating that single nucleotide substitutions are more likely to originate from the paternal germline. Additionally, there seems to be a different distribution of nucleotide changes and multiple nucleotide deletions or insertions along the UBE3A gene sequence.

[Prader Willi syndrome patients: study of 77 patients]. / Síndrome de Prader Willi: estudio de 77 pacientes.

BACKGROUND: The Prader-Willi syndrome (PWS) is a disease of genetic origin. It is characterized by neonatal hypotonia, hypogonadism, hiperfagia leading to obesity, low stature, developmental delay, moderate mental retardation, abnormal behavior and characteristic facial appearance. It is caused by the loss or the inactivation of paternal genes of the imprinted region 15q11-13. There are different genetic causes: paternal 15q11-q13 deletion in 70% of patients, maternal uniparental disomy in the 20-25% and less than 5% have an imprinting defect. We present the results obtained in the transverse clinical - genetic study of 77 PWS patients. PATIENTS AND METHODS: There has been realized the study of 374 suspected PWS patients. Cytogenetics studies of bands G and hybridization in situ fluorescent (FISH) and molecular genetics analysis of microsatellites, Southern blot, MS-PCR and sequenciation were carried out. Holm's criteria use for the correlation phenotype - genotype in 48 patients. RESULTS: PWS was confirmed in 77 patients, 46 deletion, 16 uniparental disomy, two imprinting defect and 13 only PWS methylation pattern. Significant differences do not observe in the correlation phenotype - genotype. CONCLUSIONS: The frequencies of the molecular alterations, 71.87 % deletion, 25 % UPD and 3.12 % DI, they are similar to described in the literature. It presents the algorithm of diagnosis used with the MS-PCR as rapid technology to confirm PWS.

Evolutionary and functional impact of common polymorphic inversions in the human genome.

Inversions are one type of structural variants linked to phenotypic differences and adaptation in multiple organisms. However, there is still very little information about polymorphic inversions in the human genome due to the difficulty of their detection. Here, we develop a new high-throughput genotyping method based on probe hybridization and amplification, and we perform a complete study of 45 common human inversions of 0.1-415 kb. Most inversions promoted by homologous recombination occur recurrently in humans and great apes and they are not tagged by SNPs. Furthermore, there is an enrichment of inversions showing signatures of positive or balancing selection, diverse functional effects, such as gene disruption and gene-expression changes, or association with phenotypic traits. Therefore, our results indicate that the genome is more dynamic than previously thought and that human inversions have important functional and evolutionary consequences, making possible to determine for the first time their contribution to complex traits.

Prospective analysis of liquid biopsies of advanced non-small cell lung cancer patients after progression to targeted therapies using GeneReader NGS platform.

BACKGROUND: In a significant percentage of advanced non-small cell lung cancer (NSCLC) patients, tumor tissue is unavailable or insufficient for genetic analyses at time to progression. We prospectively analyzed the appearance of genetic alterations associated with resistance in liquid biopsies of advanced NSCLC patients progressing to targeted therapies using the NGS platform. METHODS: A total of 24 NSCLC patients were included in the study, 22 progressing to tyrosine kinase inhibitors and two to other treatments. Liquid biopsies samples were obtained and analyzed using the GeneReadTM QIAact Lung DNA UMI Panel, designed to enrich specific target regions and containing 550 variant positions in 19 selected genes frequently altered in lung cancer tumors. Previously, a retrospective validation of the panel was performed in clinical samples. RESULTS: Of the 21 patients progressing to tyrosine kinase inhibitors with valid results in liquid biopsy, NGS analysis identified a potential mechanism of resistance in 12 (57%). The most common were acquired mutations in ALK and EGFR, which appeared in 8/21 patients (38%), followed by amplifications in 5/21 patients (24%), and KRAS mutations in one patient (5%). Loss of the p.T790M was also identified in two patients progressing to osimertinib. Three of the 21 (14%) patients presented two or more concomitant alterations associated with resistance. Finally, an EGFR amplification was found in the only patient progressing to immunotherapy included in the study. CONCLUSIONS: NGS analysis in liquid biopsies of patients progressing to targeted therapies using the GeneReader platform is feasible and can help the oncologist to make treatment decisions.

Probe-specific mixed-model approach to detect copy number differences using multiplex ligation-dependent probe amplification (MLPA).

BACKGROUND: MLPA method is a potentially useful semi-quantitative method to detect copy number alterations in targeted regions. In this paper, we propose a method for the normalization procedure based on a non-linear mixed-model, as well as a new approach for determining the statistical significance of altered probes based on linear mixed-model. This method establishes a threshold by using different tolerance intervals that accommodates the specific random error variability observed in each test sample. RESULTS: Through simulation studies we have shown that our proposed method outperforms two existing methods that are based on simple threshold rules or iterative regression. We have illustrated the method using a controlled MLPA assay in which targeted regions are variable in copy number in individuals suffering from different disorders such as Prader-Willi, DiGeorge or Autism showing the best performace. CONCLUSION: Using the proposed mixed-model, we are able to determine thresholds to decide whether a region is altered. These threholds are specific for each individual, incorporating experimental variability, resulting in improved sensitivity and specificity as the examples with real data have revealed.

ProSeeK: a web server for MLPA probe design.

BACKGROUND: The technological evolution of platforms for detecting genome-wide copy number imbalances has allowed the discovery of an unexpected amount of human sequence that is variable in copy number among individuals. This type of human variation can make an important contribution to human diversity and disease susceptibility. Multiplex Ligation-dependent Probe Amplification (MLPA) is a targeted method to assess copy number differences for up to 40 genomic loci in one single experiment. Although specific MLPA assays can be ordered from MRC-Holland (the proprietary company of the MLPA technology), custom designs are also developed in many laboratories worldwide. After our own experience, an important drawback of custom MLPA assays is the time spent during the design of the specific oligonucleotides that are used as probes. Due to the large number of probes included in a single assay, a number of restrictions need to be met in order to maximize specificity and to increase success likelihood. RESULTS: We have developed a web tool for facilitating and optimising custom probe design for MLPA experiments. The algorithm only requires the target sequence in FASTA format and a set of parameters, that are provided by the user according to each specific MLPA assay, to identify the best probes inside the given region. CONCLUSION: To our knowledge, this is the first available tool for optimizing custom probe design of MLPA assays. The ease-of-use and speed of the algorithm dramatically reduces the turn around time of probe design. ProSeeK will become a useful tool for all laboratories that are currently using MLPA in their research projects for CNV studies.

Imprinting center analysis in Prader-Willi and Angelman syndrome patients with typical and atypical phenotypes.

Prader-Willi syndrome (PWS) and Angelman syndrome (AS) are genetic disorders caused by a deficiency of imprinted gene expression from the paternal or maternal chromosome 15, respectively. This deficiency is due to the deletion of the 15q11-q13 region, parental uniparental disomy of the chromosome 15, or imprinting defect (ID). Mutation of the UBE3A gene causes approximately 10% of AS cases. In this present study, we describe the molecular analysis and phenotypes of two PWS patients and four AS patients with ID. One of the PWS patients has a non-familial imprinting center (IC) deletion and displayed a severe phenotype with an atypical PWS appearance, hyperactivity and psychiatric vulnerability. The other PWS and AS patients did not present genetic abnormalities in the IC, suggesting an epimutation as the genetic cause. The methylation pattern of two AS patients showed a faint maternal band corresponding to a mosaic ID. One of these mosaic patients displayed a mild AS phenotype while the other displayed a PWS-like phenotype.

KRAS mutations in the circulating free DNA (cfDNA) of non-small cell lung cancer (NSCLC) patients.

Circulating free DNA (cfDNA) is obtained from serum or plasma by non-invasive methods such as a simple blood draw, a technique known as "liquid biopsy". Genetic analyses of driver alterations in cfDNA have proved very effective to predict survival and treatment response of cancer patients according to tumoral cfDNA burden in blood. Non-small cell lung cancer (NSCLC) patients with higher concentration of tumoral cfDNA in blood have, on average, shorter progression-free survival (PFS) and overall survival (OS). Regarding specific genetic alterations, KRAS proto-oncogene, GTPase (KRAS) is one of the main genes involved in NSCLC and several studies have been performed to determine its value as a predictive and prognostic biomarker in liquid biopsy. Unfortunately, to date no strong conclusions can be drawn since they have yielded contradictory results. Therefore, further investigations are necessary to establish the value of KRAS testing in liquid biopsy as prognostic or predictive factor in NSCLC. Herein, we review the current knowledge on the importance of KRAS as prognostic and predictive biomarker using non-invasive approaches and the scientific data available regarding its application in clinical practice for treatment of NSCLC.

Liquid Biopsy in Non-Small Cell Lung Cancer.

Liquid biopsy analyses are already incorporated in the routine clinical practice in many hospitals and oncology departments worldwide, improving the selection of treatments and monitoring of lung cancer patients. Although they have not yet reached its full potential, liquid biopsy-based tests will soon be as widespread as "standard" biopsies and imaging techniques, offering invaluable diagnostic, prognostic, and predictive information. This review summarizes the techniques available for the isolation and analysis of circulating free DNA and RNA, exosomes, tumor-educated platelets, and circulating tumor cells from the blood of cancer patients, presents the methodological challenges associated with each of these materials, and discusses the clinical applications of liquid biopsy testing in lung cancer.

Human inversions and their functional consequences.

Polymorphic inversions are a type of structural variants that are difficult to analyze owing to their balanced nature and the location of breakpoints within complex repeated regions. So far, only a handful of inversions have been studied in detail in humans and current knowledge about their possible functional effects is still limited. However, inversions have been related to phenotypic changes and adaptation in multiple species. In this review, we summarize the evidences of the functional impact of inversions in the human genome. First, given that inversions have been shown to inhibit recombination in heterokaryotes, chromosomes displaying different orientation are expected to evolve independently and this may lead to distinct gene-expression patterns. Second, inversions have a role as disease-causing mutations both by directly affecting gene structure or regulation in different ways, and by predisposing to other secondary arrangements in the offspring of inversion carriers. Finally, several inversions show signals of being selected during human evolution. These findings illustrate the potential of inversions to have phenotypic consequences also in humans and emphasize the importance of their inclusion in genome-wide association studies.

Genome-wide analysis of single nucleotide polymorphisms and copy number variants in fibromyalgia suggest a role for the central nervous system.

Fibromyalgia (FM) is a highly disabling syndrome defined by a low pain threshold and a permanent state of pain. The mechanisms explaining this complex disorder remain unclear, and its genetic factors have not yet been identified. With the aim of elucidating FM genetic susceptibility factors, we selected 313 FM cases having low comorbidities, and we genotyped them on the Illumina 1 million duo array. Genotypic data from 220 control women (Illumina 610k array) was obtained for genome-wide association scan (GWAS) analysis. Copy number variants in FM susceptibility were analyzed by array comparative genomic hybridization (aCGH) experiments on pooled samples using the Agilent 2×400K platform. No single nucleotide polymorphism (SNP) reached GWAS association threshold, but 21 of the most associated SNPs were chosen for replication in 952 cases and 644 controls. Four of the SNPs selected for replication showed a nominal association in the joint analysis, and rs11127292 (MYT1L) was found to be associated to FM with low comorbidities (P=4.28×10(-5), odds ratio [95% confidence interval]=0.58 [0.44-0.75]). aCGH detected 5 differentially hybridized regions. They were followed up, and an intronic deletion in NRXN3 was demonstrated to be associated to female cases of FM with low levels of comorbidities (P=.021, odds ratio [95% confidence interval]=1.46 [1.05-2.04]). Both GWAS and aCGH results point to a role for the central nervous system in FM genetic susceptibility. If the proposed FM candidate genes were further validated in replication studies, this would highlight a neurocognitive involvement in agreement with latest reports.

Identification of copy number variants defining genomic differences among major human groups.

BACKGROUND: Understanding the genetic contribution to phenotype variation of human groups is necessary to elucidate differences in disease predisposition and response to pharmaceutical treatments in different human populations. METHODOLOGY/PRINCIPAL FINDINGS: We have investigated the genome-wide profile of structural variation on pooled samples from the three populations studied in the HapMap project by comparative genome hybridization (CGH) in different array platforms. We have identified and experimentally validated 33 genomic loci that show significant copy number differences from one population to the other. Interestingly, we found an enrichment of genes related to environment adaptation (immune response, lipid metabolism and extracellular space) within these regions and the study of expression data revealed that more than half of the copy number variants (CNVs) translate into gene-expression differences among populations, suggesting that they could have functional consequences. In addition, the identification of single nucleotide polymorphisms (SNPs) that are in linkage disequilibrium with the copy number alleles allowed us to detect evidences of population differentiation and recent selection at the nucleotide variation level. CONCLUSIONS: Overall, our results provide a comprehensive view of relevant copy number changes that might play a role in phenotypic differences among major human populations, and generate a list of interesting candidates for future studies.