A preliminary study of copy number variation in Tibetans

Cerebrospinal fluid oligoclonal bands in Chinese patients with multiple sclerosis: the prevalence and its association with clinical features

Background

Cerebrospinal fluid oligoclonal band (CSF-OCB) is an established biomarker in diagnosing multiple sclerosis (MS), however, there are no nationwide data on CSF-OCB prevalence and its diagnostic performance in Chinese MS patients, especially in the virtue of common standard operation procedure (SOP).

Methods

With a consensus SOP and the same isoelectric focusing system, we conducted a nationwide multi-center study on OCB status in consecutively, and recruited 483 MS patients and 880 non-MS patients, including neuro-inflammatory diseases (NID, n = 595) and non-inflammatory neurological diseases (NIND, n=285). Using a standardized case report form (CRF) to collect the clinical, radiological, immunological, and CSF data, we explored the association of CSF-OCB positivity with patient characters and the diagnostic performance of CSF-OCB in Chinese MS patients. Prospective source data collection, and retrospective data acquisition and statistical data analysis were used.

Findings

369 (76.4%) MS patients were OCB-positive, while 109 NID patients (18.3%) and 6 NIND patients (2.1%) were OCB-positive, respectively. Time from symptom onset to diagnosis was significantly shorter in OCB-positive than that in OCB-negative MS patients (13.2 vs 23.7 months, P=0.020). The prevalence of CSF-OCB in Chinese MS patients was significantly higher in high-latitude regions (41°-50°N)(P=0.016), and at high altitudes (>1000m)(P=0.025). The diagnostic performance of CSF-OCB differentiating MS from non-MS patients yielded a sensitivity of 76%, a specificity of 87%.

Interpretation

The nationwide prevalence of CSF-OCB was 76.4% in Chinese MS patients, and demonstrated a good diagnostic performance in differentiating MS from other CNS diseases. The CSF-OCB prevalence showed a correlation with high latitude and altitude in Chinese MS patients.

Why is the beautyberry so colourful? Evolution, biogeography, and diversification of fruit colours in Callicarpa (Lamiaceae)

Fruit colour is essential to seed dispersal, speciation, and biological diversity in global ecosystems. The relationship between fruit-colour variation and species diversification has long been of interest in evolutionary biology, but remains poorly understood at the genus level. Here, we used Callicarpa, a typical representative of pantropical angiosperm, to analyse whether fruit colours are correlated with biogeographic distribution, dispersal events, and diversification rate. We estimated a time-calibrated phylogeny for Callicarpa and reconstructed ancestral fruit colour. Utilizing phylogenetic methods, we estimated the major dispersal events across the phylogenetic tree and the most likely fruit colours related to each dispersal event, and tested whether the dispersal frequencies and distances of the four fruit colours between major biogeographical areas were equal. We then tested whether fruit colours are correlated with latitude, elevation, and diversification rate. Biogeographical reconstructions showed that Callicarpa originated in the East Asia and Southeast Asia during the Eocene (∼35.53 Ma) and diverse species diverged mainly in the Miocene and lasted into the Pleistocene. Large-scale dispersal events were significantly associated with violet-fruited lineages. Furthermore, different fruit colours were markedly correlated with different latitudes and elevations (e.g., violet fruits were correlated with higher latitudes and elevations; red fruits and black fruits with lower latitudes; white fruits with higher elevations). Notably, violet fruits were statistically associated with highest diversification rates, driving fruit colour variation among different regions globally. Our results contribute to further understanding why fruit colour is so variable at the genus level of angiosperms in different areas around the world.

A study of normal copy number variations in Israeli population

The population of Israel is ethnically diverse, and individuals from different ethnic groups share specific genetic variations. These variations, which have been passed on from common ancestors, are usually reported in public databases as rare variants. Here, we aimed to identify ethnicity-based benign copy number variants (CNVs) and generate the first Israeli CNV database. We applied a data-mining approach to the results of 10,193 chromosomal microarray tests, of which 2150 tests were from individuals of 13 common ethnic backgrounds (n ≥ 10). We found 165 CNV regions (> 50 kbp) that are unique to specific ethnic groups (uCNVRs). The frequency of more than 19% of these uCNVRs is between 1 and 20% of the common ethnic origin, while their frequency in the overall cohort is between 0.5 and 1.6%. Of these 165 uCNVRs, 98 are reported as variants of unknown significance or as not available in dbVar; we postulate that these uCNVRs should be annotated as either "likely benign" or "benign". The ethnic-specific CNVs extracted in this study will allow geneticists to distinguish between relevant pathogenic genomic aberrations and benign ethnicity-related variations, thus preventing variant misinterpretation that may lead to unnecessary pregnancy terminations.

Neurodevelopmental trajectory and modifiers of 16p11.2 microdeletion: A follow-up study of four Chinese children carriers

Background

Neurodevelopmental disorders (NDDs) are a group of disorders with high genetic and phenotypic heterogeneities. The 16p11.2 microdeletion has been implicated as an important genetic risk factor for NDDs.

Methods

Multiple genetic tests were used to detect the 16p11.2 microdeletion from 918 Chinese children with NDDs. Targeted sequencing of genes in the 16p11.2 interval was performed in all carriers of the 16p11.2 microdeletion, and whole‐genome expression profiling analysis was performed for the patient carriers and normal carriers in their intra‐family.

Results

Three patients carrying the 16p11.2 microdeletion were screened out, indicating a frequency of 0.33% for the 16p11.2 microdeletion in this cohort. We reviewed the neurodevelopmental trajectories of the 16p11.2 microdeletion carriers from childhood to puberty and confirmed that this microdeletion was associated with abnormal neurodevelopment, with varied neurodevelopmental phenotypes. A differential PRRT2 genotype (rs10204, T>C) was identified between patients and normal carriers of the 16p11.2 microdeletion. Moreover, the determination of differential whole‐genome expression profiling demonstrated the destruction of the top‐ranked network in neurogenesis and accounted for observation of abnormal neurodevelopmental phenotypes in the 16p11.2 microdeletion carriers.

Conclusions

We have provided the frequency of the 16p11.2 microdeletion in a Chinese pediatric NDD cohort with a variable NDD phenotype from childhood to puberty, which is useful for Chinese geneticists/pediatricians to conduct the 16p11.2 microdeletion testing in children with NDDs.

Population Structure, and Selection Signatures Underlying High-Altitude Adaptation Inferred From Genome-Wide Copy Number Variations in Chinese Indigenous Cattle

Copy number variations (CNVs) have been demonstrated as crucial substrates for evolution, adaptation and breed formation. Chinese indigenous cattle breeds exhibit a broad geographical distribution and diverse environmental adaptability. Here, we analyzed the population structure and adaptation to high altitude of Chinese indigenous cattle based on genome-wide CNVs derived from the high-density BovineHD SNP array. We successfully detected the genome-wide CNVs of 318 individuals from 24 Chinese indigenous cattle breeds and 37 yaks as outgroups. A total of 5,818 autosomal CNV regions (683 bp-4,477,860 bp in size), covering ~14.34% of the bovine genome (UMD3.1), were identified, showing abundant CNV resources. Neighbor-joining clustering, principal component analysis (PCA), and population admixture analysis based on these CNVs support that most Chinese cattle breeds are hybrids of Bos taurus taurus (hereinafter to be referred as Bos taurus) and Bos taurus indicus (Bos indicus). The distribution patterns of the CNVs could to some extent be related to the geographical backgrounds of the habitat of the breeds, and admixture among cattle breeds from different districts. We analyzed the selective signatures of CNVs positively involved in high-altitude adaptation using pairwise Fst analysis within breeds with a strong Bos taurus background (taurine-type breeds) and within Bos taurus×Bos indicus hybrids, respectively. CNV-overlapping genes with strong selection signatures (at top 0.5% of Fst value), including LETM1 (Fst = 0.490), TXNRD2 (Fst = 0.440), and STUB1 (Fst = 0.420) within taurine-type breeds, and NOXA1 (Fst = 0.233), RUVBL1 (Fst = 0.222), and SLC4A3 (Fst=0.154) within hybrids, were potentially involved in the adaptation to hypoxia. Thus, we provide a new profile of population structure from the CNV aspects of Chinese indigenous cattle and new insights into high-altitude adaptation in cattle.

What has GWAS done for HLA and disease associations?

The major histocompatibility complex (MHC) is located in chromosome 6p21 and contains crucial regulators of immune response, including human leucocyte antigen (HLA) genes, alongside other genes with nonimmunological roles. More recently, a repertoire of noncoding RNA genes, including expressed pseudogenes, has also been identified. The MHC is the most gene dense and most polymorphic part of the human genome. The region exhibits haplotype-specific linkage disequilibrium patterns, contains the strongest cis- and trans-eQTLs/meQTLs in the genome and is known as a hot spot for disease associations. Another layer of complexity is provided to the region by the extreme structural variation and copy number variations. While the HLA-B gene has the highest number of alleles, the HLA-DR/DQ subregion is structurally most variable and shows the highest number of disease associations. Reliance on a single reference sequence has complicated the design, execution and analysis of GWAS for the MHC region and not infrequently, the MHC region has even been excluded from the analysis of GWAS data. Here, we contrast features of the MHC region with the rest of the genome and highlight its complexities, including its functional polymorphisms beyond those determined by single nucleotide polymorphisms or single amino acid residues. One of the several issues with customary GWAS analysis is that it does not address this additional layer of polymorphisms unique to the MHC region. We highlight alternative approaches that may assist with the analysis of GWAS data from the MHC region and unravel associations with all functional polymorphisms beyond single SNPs. We suggest that despite already showing the highest number of disease associations, the true extent of the involvement of the MHC region in disease genetics may not have been uncovered.

High-altitude adaptation in humans: from genomics to integrative physiology

About 1.2 to 33% of high-altitude populations suffer from Monge's disease or chronic mountain sickness (CMS). Number of factors such as age, sex, and population of origin (older, male, Andean) contribute to the percentage reported from a variety of samples. It is estimated that there are around 83 million people who live at altitudes > 2500 m worldwide and are at risk for CMS. In this review, we focus on a human "experiment in nature" in various high-altitude locations in the world-namely, Andean, Tibetan, and Ethiopian populations that have lived under chronic hypoxia conditions for thousands of years. We discuss the adaptive as well as mal-adaptive changes at the genomic and physiological levels. Although different genes seem to be involved in adaptation in the three populations, we can observe convergence at genetic and signaling, as well as physiological levels. What is important is that we and others have shown that lessons learned from the genes mined at high altitude can be helpful in better understanding and treating diseases that occur at sea level. We discuss two such examples: EDNRB and SENP1 and their role in cardiac tolerance and in the polycythemic response, respectively.

Identifying the risk of producing aneuploids using meiotic recombination genes as biomarkers: A copy number variation approach

Background & objectives:

Aneuploids are the most common chromosomal abnormality in liveborns and are usually the result of non-disjunction (NDJ) in meiosis. Copy number variations (CNVs) are large structural variations affecting the human genome. CNVs influence critical genes involved in causing NDJ by altering their copy number which affects the clinical outcome. In this study influence of CNVs on critical meiotic recombination was examined using new computational technologies to assess their role in causing aneuploidy.

Methods:

This investigation was based on the analysis of 12 random normal populations consisting of 1714 individuals for aneuploid causing genes under CNV effect. To examine the effect of CNVs on genes causing aneuploidy, meiotic recombination genes were analyzed using EnrichR, WebGestalt and Ingenuity Pathway Analysis (IPA).

Results:

Forty three NDJ genes were found under CNV burden; IPA (Ingenuity Pathway Analysis) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis of CNV in meiotic recombination genes revealed a significant role of breast cancer gene 1, amyloid protein precursor, mitogen-activated protein kinase and nerve growth factor as key molecular players involved in causing aneuploidy. Interaction between these genes with other CNV-overlapping genes involved in cell cycle, recombination and meiosis might lead to increased incidences of aneuploidy.

Interpretation & conclusions:

The findings of this study implied that the effect of CNVs on normal genome contributed in amplifying the occurrences of chromosomal aneuploidies. The normal individuals consisting of variations in the susceptible genes causing aneuploids in the population remain undetected until the disorder genes express in the succeeding generations.

Genome-wide identification of genes probably relevant to the adaptation of schizothoracins (Teleostei: Cypriniformes) to the uplift of the Qinghai-Tibet Plateau

BACKGROUND

Molecular adaptation to the severe environments present during the uplift of the Qinghai-Tibet Plateau has attracted the attention of researchers. The divergence of the three specialization groups of schizothoracins (Primitive, Specialized and Highly Specialized) may correspond to the three phases of plateau uplift. Based on the transcripts of representative species of the three specialized groups and an outgroup, genes in schizothoracins that may have played important roles during the adaptation to new environments were investigated.

RESULTS

The contigs of Gymnodiptychus dybowskii and Schizothorax pseudaksaiensis were compared with those of Gymnocypris przewalskii ganzihonensis and the outgroup Sinocyclocheilus angustiporus, and 5,894 ortholog groups with an alignment length longer than 90 nt after deleting gaps were retained. Evolutionary analyses indicated that the average evolutionary rate of the branch leading to the Specialized group was faster than that of the branch leading to the Highly Specialized group. Moreover, the numbers of gene categories in which more than half of the genes evolved faster than the average values of the genome were 117 and 15 along the branches leading to the Specialized and Highly Specialized groups, respectively. A total of 40, 36, and 55 genes were likely subject to positive selection along the branches leading to the Primitive, Specialized and Highly Specialized groups, respectively, and many of these genes are likely relevant to adaptation to the cold temperatures, low oxygen concentrations, and strong ultraviolet radiation that result from elevation.

CONCLUSIONS

By selecting representative species of the three groups of schizothoracins and applying next-generation sequencing technology, several candidate genes corresponding to adaptation to the three phases of plateau uplift were identified. Some of the genes identified in this report that were likely subject to positive selection are good candidates for subsequent evolutionary and functional analyses of adaptation to high altitude.

Copy number variations in the genome of the Qatari population

Background

The populations of the Arabian Peninsula remain the least represented in public genetic databases, both in terms of single nucleotide variants and of larger genomic mutations. We present the first high-resolution copy number variation (CNV) map for a Gulf Arab population, using a hybrid approach that integrates array genotyping intensity data and next-generation sequencing reads to call CNVs in the Qatari population.

Methods

CNVs were detected in 97 unrelated Qatari individuals by running two calling algorithms on each of two primary datasets: high-resolution genotyping (Illumina Omni 2.5M) and high depth whole-genome sequencing (Illumina PE 100bp). The four call-sets were integrated to identify high confidence CNV regions, which were subsequently annotated for putative functional effect and compared to public databases of CNVs in other populations. The availability of genome sequence was leveraged to identify tagging SNPs in high LD with common deletions in this population, enabling their imputation from genotyping experiments in the future.

Results

Genotyping intensities and genome sequencing data from 97 Qataris were analyzed with four different algorithms and integrated to discover 16,660 high confidence CNV regions (CNVRs) in the total population, affecting ~28 Mb in the median Qatari genome. Up to 40 % of all CNVs affected genes, including novel CNVs affecting Mendelian disease genes, segregating at different frequencies in the 3 major Qatari subpopulations, including those with Bedouin, Persian/South Asian, and African ancestry. Consistent with high consanguinity levels in the Bedouin subpopulation, we found an increased burden for homozygous deletions in this group. In comparison to known CNVs in the comprehensive Database of Genomic Variants, we found that 5 % of all CNVRs in Qataris were completely novel, with an enrichment of CNVs affecting several known chromosomal disorder loci and genes known to regulate sugar metabolism and type 2 diabetes in the Qatari cohort. Finally, we leveraged the availability of genome sequence to find suitable tagging SNPs for common deletions in this population.

Conclusion

We combine four independently generated datasets from 97 individuals to study CNVs for the first time at high-resolution in a Gulf Arab population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1991-5) contains supplementary material, which is available to authorized users.

Global patterns of large copy number variations in the human genome reveal complexity in chromosome organization

Global patterns of copy number variations (CNVs) in chromosomes are required to understand the dynamics of genome organization and complexity. For this study, analysis was performed using the Affymetrix Genome-Wide Human SNP Array 6.0 chip and CytoScan High-Density arrays. We identified a total of 44 109 CNVs from 1715 genomes with a mean of 25 CNVs in an individual, which established the first drafts of population-specific CNV maps providing a rationale for prioritizing chromosomal regions. About 19 905 ancient CNVs were identified across all chromosomes and populations at varying frequencies. CNV count, and sometimes CNV size, contributed to the bulk CNV size of the chromosome. Population specific lengthening and shortening of chromosomal length was observed. Sex bias for CNV presence was largely dependent on ethnicity. Lower CNV inheritance rate was observed for India, compared to YRI and CEU. A total of 33 candidate CNV hotspots from 5382 copy number (CN) variable region (CNVR) clusters were identified. Population specific CNV distribution patterns in p and q arms disturbed the assumption that CNV counts in the p arm are less common compared to long arms, and the CNV occurrence and distribution in chromosomes is length independent. This study unraveled the force of independent evolutionary dynamics on genome organization and complexity across chromosomes and populations.

Fifteen years of genomewide scans for selection: trends, lessons and unaddressed genetic sources of complication

Genomewide scans for natural selection (GWSS) have become increasingly common over the last 15 years due to increased availability of genome-scale genetic data. Here, we report a representative survey of GWSS from 1999 to present and find that (i) between 1999 and 2009, 35 of 49 (71%) GWSS focused on human, while from 2010 to present, only 38 of 83 (46%) of GWSS focused on human, indicating increased focus on nonmodel organisms; (ii) the large majority of GWSS incorporate interpopulation or interspecific comparisons using, for example F(ST), cross-population extended haplotype homozygosity or the ratio of nonsynonymous to synonymous substitutions; (iii) most GWSS focus on detection of directional selection rather than other modes such as balancing selection; and (iv) in human GWSS, there is a clear shift after 2004 from microsatellite markers to dense SNP data. A survey of GWSS meant to identify loci positively selected in response to severe hypoxic conditions support an approach to GWSS in which a list of a priori candidate genes based on potential selective pressures are used to filter the list of significant hits a posteriori. We also discuss four frequently ignored determinants of genomic heterogeneity that complicate GWSS: mutation, recombination, selection and the genetic architecture of adaptive traits. We recommend that GWSS methodology should better incorporate aspects of genomewide heterogeneity using empirical estimates of relevant parameters and/or realistic, whole-chromosome simulations to improve interpretation of GWSS results. Finally, we argue that knowledge of potential selective agents improves interpretation of GWSS results and that new methods focused on correlations between environmental variables and genetic variation can help automate this approach.

Ladakh, India: the land of high passes and genetic heterogeneity reveals a confluence of migrations

Owing to its geographic location near the longitudinal center of Asia, Ladakh, the land of high passes, has witnessed numerous demographic movements during the past millenniums of occupation. In an effort to view Ladakh's multicultural history from a paternal genetic perspective, we performed a high-resolution Y-chromosomal survey of Ladakh, within the context of Y haplogroup and haplotype distributions of 41 Asian reference populations. The results of this investigation highlight the rich ethnic and genetic diversity of Ladkah which includes genetic contributions from disparate regions of the continent including, West, East, South and Central Asia. The phylogenetic signals from Ladakh are consistent with the Indo-Aryans' occupation during the Neolithic age and its historic connection with Tibet, as well as the East-West gene flow associated with the Silk Road.

Global spectrum of copy number variations reveals genome organizational plasticity and proposes new migration routes

Global spectrum of CNVs is required to catalog variations to provide a high-resolution on the dynamics of genome-organization and human migration. In this study, we performed genome-wide genotyping using high-resolution arrays and identified 44,109 CNVs from 1,715 genomes across 12 populations. The study unraveled the force of independent evolutionary dynamics on genome-organizational plasticity across populations. We demonstrated the use of CNV tool to study human migration and identified a second major settlement establishing new migration routes in addition to existing ones.

Impact of copy number variations burden on coding genome in humans using integrated high resolution arrays

Copy number variations (CNVs) alter the transcriptional and translational levels of genes by disrupting the coding structure and this burden of CNVs seems to be a significant contributor to phenotypic variations. Therefore it was necessary to assess the complexities of CNV burden on the coding genome. A total of 1715 individuals from 12 populations were used for CNV analysis in the present investigation. Analysis was performed using Affymetrix Genome-Wide Human SNP Array 6·0 chip and CytoScan High-Density arrays. CNVs were more frequently observed in the coding region than in the non-coding region. CNVs were observed vastly more frequently in the coding region than the non-coding region. CNVs were found to be enriched in the regions containing functional genes (83-96%) compared with the regions containing pseudogenes (4-17%). CNVs across the genome of an individual showed multiple hits across many genes, whose proteins interact physically and function under the same pathway. We identified varying numbers of proteins and degrees of interactions within protein complexes of single individual genomes. This study represents the first draft of a population-specific CNV genes map as well as a cross-populational map. The complex relationship of CNVs on genes and their physically interacting partners unravels many complexities involved in phenotype expression. This study identifies four mechanisms contributing to the complexities caused by the presence of multiple CNVs across many genes in the coding part of the genome.

Two functional loci in the promoter of EPAS1 gene involved in high-altitude adaptation of Tibetans

EPAS1 involves in the hypoxic response and is suggested to be responsible for the genetic adaptation of high-altitude hypoxia in Tibetans. However, the detailed molecular mechanism remains unknown. In this study, a single nucleotide polymorphism rs56721780:G>C and an insertion/deletion (indel) polymorphism -742 indel in the promoter region showed divergence between Tibetans and non-Tibetan lowlanders. rs56721780:G>C regulated the transcription of EPAS1 by IKAROS family zinc finger 1 (IKZF1), which was identified as a new transcriptional repressor for EPAS1 gene. It demonstrated that the C allele of rs56721780:G>C decreased the binding of IKZF1, leading to the attenuated transcriptional repression of EPAS1 gene. The insertion at -742 indel provided a new binding site for Sp1 and was related to the activation of EPAS1 promoter. Further functional analysis revealed that lysyl oxidase (LOX) gene, which was reported to be responsible for extracellular matrix protein cross-linking of amnion previously, was a direct target of EPAS1. The CC genotype at rs56721780:G>C and the insertion genotype at -742 indel were found associated with higher EPAS1 and LOX expression levels in amnion, as well as higher birth weight of Tibetan newborns, suggesting that EPAS1 gene might play important roles in the development of amnion, fetus growth and high-altitude adaptation of Tibetans.

Insertion-deletions burden in copy number polymorphisms of the Tibetan population

BACKGROUND:

Many studies have been conducted to identify either insertions-deletions (inDels) or copy number variations (CNVs) in humans, but few studies have been conducted to identify both of these forms coexisting in the same region.

AIMS AND OBJECTIVES:

To map the functionally significant sites within human genes that are likely to influence human traits and diseases.

MATERIALS AND METHODS:

In this report, we describe an inDel map in the 1051 Tibetan CNV regions obtained through CNV genotyping using Affymetrix Genome-wide single nucleotide polymorphism 6.0 chip. InDel polymorphisms in these copy number polymorphism regions were identified with a computational approach using the 2500 deoxyribonucleic acid sequences obtained from the 1000 Genome Project.

RESULTS:

The study identified a total of 95935 inDels that range from 1 bp to several bps in length which were found scattered across regulatory regions, exons and in introns of genes underlying the CNVs. A study on the distribution of inDels revealed that the majority of inDels were found in coding regions of the genome than the noncoding, while within the genes, inDels in intron regions were more followed by exonic regions and finally the regulatory regions.

CONCLUSION:

Study of inDels in CNV regions contribute to the enhanced understanding of the role played by the two variations and their collective influence on the genome. Further, a collection of these inDel genetic markers will aid in genetic mapping, further understanding of the phenotypic variability, identification of disease genes and in detecting novel CNVs.

Identification of independent primary tumors and intrapulmonary metastases using DNA rearrangements in non-small-cell lung cancer

PURPOSE

Distinguishing independent primary tumors from intrapulmonary metastases in non-small-cell carcinoma remains a clinical dilemma with significant clinical implications. Using next-generation DNA sequencing, we developed a chromosomal rearrangement-based approach to differentiate multiple primary tumors from metastasis.

METHODS

Tumor specimens from patients with known independent primary tumors and metastatic lesions were used for lineage test development, which was then applied to multifocal tumors. Laser capture microdissection was performed separately for each tumor. Genomic DNA was isolated using direct in situ whole-genome amplification methodology, and next-generation sequencing was performed using an Illumina mate-pair library protocol. Sequence reads were mapped to the human genome, and primers spanning the fusion junctions were used for validation polymerase chain reaction.

RESULTS

A total of 41 tumor samples were sequenced (33 adenocarcinomas [ADs] and eight squamous cell carcinomas [SQCCs]), with a range of three to 276 breakpoints per tumor identified. Lung tumors predicted to be independent primary tumors based on different histologic subtype did not share any genomic rearrangements. In patients with lung primary tumors and paired distant metastases, shared rearrangements were identified in all tumor pairs, emphasizing the patient specificity of identified breakpoints. Multifocal AD and SQCC samples were reviewed independently by two pulmonary pathologists. Concordance between histology and genomic data occurred in the majority of samples. Discrepant tumor samples were resolved by genome sequencing.

CONCLUSION

A diagnostic lineage test based on genomic rearrangements from mate-pair sequencing demonstrates promise for distinguishing independent primary from metastatic disease in lung cancer.

Genome-wide linkage study suggests a susceptibility locus for isolated bilateral microtia on 4p15.32-4p16.2

Microtia is a congenital deformity where the external ear is underdeveloped. Genetic investigations have identified many susceptibility genes of microtia-related syndromes. However, no causal genes were reported for isolated microtia, the main form of microtia. We conducted a genome-wide linkage analysis on a 5-generation Chinese pedigree with isolated bilateral microtia. We identified a suggestive linkage locus on 4p15.32-4p16.2 with parametric LOD score of 2.70 and nonparametric linkage score (Zmean) of 12.28 (simulated occurrence per genome scan equal to 0.46 and 0.47, respectively). Haplotype reconstruction analysis of the 4p15.32-4p16.2 region further confined the linkage signal to a 10-Mb segment located between rs12505562 and rs12649803 (9.65-30.24 cM; 5.54-15.58 Mb). Various human organ developmental genes reside in this 10-Mb susceptibility region, such as EVC, EVC2, SLC2A9, NKX3-2, and HMX1. The coding regions of three genes, EVC known for cartilage development and NKX3-2, HMX1 involved in microtia, were selected for sequencing with 5 individuals from the pedigree. Of the 38 identified sequence variants, none segregates along with the disease phenotype. Other genes or DNA sequences of the 10-Mb region warrant for further investigation. In conclusion, we report a susceptibility locus of isolated microtia, and this finding will encourage future studies on the genetic basis of ear deformity.

Novel population specific autosomal copy number variation and its functional analysis amongst Negritos from Peninsular Malaysia

Copy number variation (CNV) has been recognized as a major contributor to human genome diversity. It plays an important role in determining phenotypes and has been associated with a number of common and complex diseases. However CNV data from diverse populations is still limited. Here we report the first investigation of CNV in the indigenous populations from Peninsular Malaysia. We genotyped 34 Negrito genomes from Peninsular Malaysia using the Affymetrix SNP 6.0 microarray and identified 48 putative novel CNVs, consisting of 24 gains and 24 losses, of which 5 were identified in at least 2 unrelated samples. These CNVs appear unique to the Negrito population and were absent in the DGV, HapMap3 and Singapore Genome Variation Project (SGVP) datasets. Analysis of gene ontology revealed that genes within these CNVs were enriched in the immune system (GO:0002376), response to stimulus mechanisms (GO:0050896), the metabolic pathways (GO:0001852), as well as regulation of transcription (GO:0006355). Copy number gains in CNV regions (CNVRs) enriched with genes were significantly higher than the losses (P value <0.001). In view of the small population size, relative isolation and semi-nomadic lifestyles of this community, we speculate that these CNVs may be attributed to recent local adaptation of Negritos from Peninsular Malaysia.

Compilation of copy number variants identified in phenotypically normal and parous Japanese women

With increasing public concern about infertility and the frequent involvement of chromosomal anomalies in miscarriage, analyses of copy number variations (CNVs) have been used to identify the genomic regions responsible for each process of childbearing. Although associations between CNVs and diseases have been reported, many CNVs have also been identified in healthy individuals. Like other types of mutations, phenotypically indefinite CNVs may have been retained and accumulated during anthropogenesis. Therefore to distinguish causative variants from other variants is a formidable task. Furthermore, because previous studies have predominantly focused on European and African populations, comprehensive detection of common Asian CNVs is eagerly awaited. Here, using a high-resolution genotyping array and samples from 411 Japanese women with normal parity without significant complications, we have compiled 1043 copy number variable regions. In total, the collected regions cover 164 Mb, or up to 0.5% of the genome. The copy number differences in these regions may be irrelevant not only to infertility but also to a wide range of diseases. The utility of this resource in reducing the candidate pathogenetic variants, especially in Japanese subjects, is also demonstrated.

‘Ome’ on the range: update on high-altitude acclimatization/adaptation and disease

The main physiological challenge in high-altitude plateau environments is hypoxia.