Identification of single nucleotide variants in the Moroccan population by whole-genome sequencing.

BACKGROUND: Large-scale human sequencing projects have described around a hundred-million single nucleotide variants (SNVs). These studies have predominately involved individuals with European ancestry despite the fact that genetic diversity is expected to be highest in Africa where Homo sapiens evolved and has maintained a large population for the longest time. The African Genome Variation Project examined several African populations but these were all located south of the Sahara. Morocco is on the northwest coast of Africa and mostly lies north of the Sahara, which makes it very attractive for studying genetic diversity. The ancestry of present-day Moroccans is unknown and may be substantially different from Africans found South of the Sahara desert, Recent genomic data of Taforalt individuals in Eastern Morocco revealed 15,000-year-old modern humans and suggested that North African individuals may be genetically distinct from previously studied African populations. RESULTS: We present SNVs discovered by whole genome sequencing (WGS) of three Moroccans. From a total of 5.9 million SNVs detected, over 200,000 were not identified by 1000G and were not in the extensive gnomAD database. We summarise the SNVs by genomic position, type of sequence gene context and effect on proteins encoded by the sequence. Analysis of the overall genomic information of the Moroccan individuals to individuals from 1000G supports the Moroccan population being distinct from both sub-Saharan African and European populations. CONCLUSIONS: We conclude that Moroccan samples are genetically distinct and lie in the middle of the previously observed cline between populations of European and African ancestry. WGS of Moroccan individuals can identify a large number of novel SNVs and aid in functional characterisation of the genome.

Mass spectrometry imaging of endogenous metabolites in response to doxorubicin in a novel 3D osteosarcoma cell culture model.

Three-dimensional (3D) cell culture is a rapidly emerging field, which mimics some of the physiological conditions of human tissues. In cancer biology, it is considered a useful tool in predicting in vivo chemotherapy responses, compared with conventional two-dimensional (2D) cell culture. We have developed a novel 3D cell culture model of osteosarcoma composed of aggregated proliferative tumour spheroids, which shows regions of tumour heterogeneity formed by aggregated spheroids of polyclonal tumour cells. Aggregated spheroids show local necrotic and apoptotic regions and have sizes suitable for the study of spatial distribution of metabolites by mass spectrometry imaging (MSI). We have used this model to perform a proof-of-principle study showing a heterogeneous distribution of endogenous metabolites that colocalise with the necrotic core and apoptotic regions in this model. Cytotoxic chemotherapy (doxorubicin) responses were significantly attenuated in our 3D cell culture model compared with those of standard cell culture, as determined by resazurin assay, despite sufficient doxorubicin diffusion demonstrated by localisation throughout the 3D constructs. Finally, changes to the distribution of endogenous metabolites in response to doxorubicin were readily detected by MSI. Principal component analysis identified 50 metabolites which differed most in their abundance between treatment groups, and of these, 10 were identified by both in-software t test and mixed-effects analysis of variance (ANOVA). Subsequent independent MSIs of identified species were consistent with principle component analysis findings. This proof-of-principle study shows for the first time that chemotherapy-induced changes in metabolite abundance and distribution may be determined in 3D cell culture by MSI, highlighting this method as a potentially useful tool in the elucidation of chemotherapy responses as an alternative to in vivo testing.

Consequences of Epistasis on Growth in an Erhualian × White Duroc Pig Cross.

Epistasis describes an interaction between the effects of loci. We included epistasis in quantitative trait locus (QTL) mapping of growth at a series of ages in a cross of a Chinese pig breed, Erhualian, with a commercial line, White Duroc. Erhualian pigs have much lower growth rates than White Duroc. We improved a method for genomewide testing of epistasis and present a clear analysis workflow. We also suggest a new approach for interpreting epistasis results where significant additive and dominance effects of a locus in specific backgrounds are determined. In total, seventeen QTL were found and eleven showed epistasis. Loci on chromosomes 2, 3, 4 and 7 were highlighted as affecting growth at more than one age or forming an interaction network. Epistasis resulted in both the QTL on chromosomes 3 and 7 having effects in opposite directions. We believe it is the first time for the chromosome 7 locus that an allele from a Chinese breed has been found to decrease growth. The consequences of epistasis were diverse. Results were impacted by using growth rather than body weight as the phenotype and by correcting for an effect of mother. Epistasis made a considerable contribution to growth in this population and modelling epistasis was important for accurately determining QTL effects.

Rare loss-of-function variants in SETD1A are associated with schizophrenia and developmental disorders.

By analyzing the whole-exome sequences of 4,264 schizophrenia cases, 9,343 controls and 1,077 trios, we identified a genome-wide significant association between rare loss-of-function (LoF) variants in SETD1A and risk for schizophrenia (P = 3.3 × 10(-9)). We found only two heterozygous LoF variants in 45,376 exomes from individuals without a neuropsychiatric diagnosis, indicating that SETD1A is substantially depleted of LoF variants in the general population. Seven of the ten individuals with schizophrenia carrying SETD1A LoF variants also had learning difficulties. We further identified four SETD1A LoF carriers among 4,281 children with severe developmental disorders and two more carriers in an independent sample of 5,720 Finnish exomes, both with notable neuropsychiatric phenotypes. Together, our observations indicate that LoF variants in SETD1A cause a range of neurodevelopmental disorders, including schizophrenia. Combining these data with previous common variant evidence, we suggest that epigenetic dysregulation, specifically in the histone H3K4 methylation pathway, is an important mechanism in the pathogenesis of schizophrenia.

The UK10K project identifies rare variants in health and disease.

The contribution of rare and low-frequency variants to human traits is largely unexplored. Here we describe insights from sequencing whole genomes (low read depth, 7×) or exomes (high read depth, 80×) of nearly 10,000 individuals from population-based and disease collections. In extensively phenotyped cohorts we characterize over 24 million novel sequence variants, generate a highly accurate imputation reference panel and identify novel alleles associated with levels of triglycerides (APOB), adiponectin (ADIPOQ) and low-density lipoprotein cholesterol (LDLR and RGAG1) from single-marker and rare variant aggregation tests. We describe population structure and functional annotation of rare and low-frequency variants, use the data to estimate the benefits of sequencing for association studies, and summarize lessons from disease-specific collections. Finally, we make available an extensive resource, including individual-level genetic and phenotypic data and web-based tools to facilitate the exploration of association results.

MAPfastR: quantitative trait loci mapping in outbred line crosses.

MAPfastR is a software package developed to analyze quantitative trait loci data from inbred and outbred line-crosses. The package includes a number of modules for fast and accurate quantitative trait loci analyses. It has been developed in the R language for fast and comprehensive analyses of large datasets. MAPfastR is freely available at: http://www.computationalgenetics.se/?page_id=7.

Identification of null alleles and deletions from SNP genotypes for an intercross between domestic and wild chickens.

We analyzed genotypes from ~10K single-nucleotide polymorphisms (SNPs) in two families of an F2 intercross between Red Junglefowl and White Leghorn chickens. Possible null alleles were found by patterns of incompatible and missing genotypes. We estimated that 2.6% of SNPs had null alleles compared with 2.3% with genotyping errors and that 40% of SNPs in which a parent and offspring were genotyped as different homozygotes had null alleles. Putative deletions were identified by null alleles at adjacent markers. We found two candidate deletions that were supported by fluorescence intensity data from a 60K SNP chip. One of the candidate deletions was from the Red Junglefowl, and one was present in both the Red Junglefowl and White Leghorn. Both candidate deletions spanned protein-coding regions and were close to a previously detected quantitative trait locus affecting body weight in this population. This study demonstrates that the ~50K SNP genotyping arrays now available for several agricultural species can be used to identify null alleles and deletions in data from large families. We suggest that our approach could be a useful complement to linkage analysis in experimental crosses.

Mapping QTL affecting a systemic sclerosis-like disorder in a cross between UCD-200 and red jungle fowl chickens.

Systemic sclerosis (SSc) or scleroderma is a rare, autoimmune, multi-factorial disease characterized by early microvascular alterations, inflammation, and fibrosis. Chickens from the UCD-200 line develop a hereditary SSc-like disease, showing all the hallmarks of the human disorder, which makes this line a promising model to study genetic factors underlying the disease. A backcross was generated between UCD-200 chickens and its wild ancestor - the red jungle fowl and a genome-scan was performed to identify loci affecting early (21 days of age) and late (175 days of age) ischemic lesions of the comb. A significant difference in frequency of disease was observed between sexes in the BC population, where the homogametic males were more affected than females, and there was evidence for a protective W chromosome effect. Three suggestive disease predisposing loci were mapped to chromosomes 2, 12 and 14. Three orthologues of genes implicated in human SSc are located in the QTL region on chromosome 2, TGFRB1, EXOC2-IRF4 and COL1A2, as well as CCR8, which is more generally related to immune function. IGFBP3 is also located within the QTL on chromosome 2 and earlier studies have showed increased IGFBP3 serum levels in SSc patients. To our knowledge, this study is the first to reveal a potential genetic association between IGFBP3 and SSc. Another gene with an immunological function, SOCS1, is located in the QTL region on chromosome 14. These results illustrate the usefulness of the UCD-200 chicken as a model of human SSc and motivate further in-depth functional studies of the implicated candidate genes.

An Improved Method for Estimating Chromosomal Line Origin in QTL Analysis of Crosses Between Outbred Lines.

Estimating the line origin of chromosomal sections from marker genotypes is a vital step in quantitative trait loci analyses of outbred line crosses. The original, and most commonly used, algorithm can only handle moderate numbers of partially informative markers. The advent of high-density genotyping with SNP chips motivates a new method because the generic sets of markers on SNP chips typically result in long stretches of partially informative markers. We validated a new method for inferring line origin, triM (tracing inheritance with Markov models), with simulated data. A realistic pattern of marker information was achieved by replicating the linkage disequilibrium from an existing chicken intercross. There were approximately 1500 SNP markers and 800 F(2) individuals. The performance of triM was compared to GridQTL, which uses a variant of the original algorithm but modified for larger datasets. triM estimated the line origin with an average error of 2%, was 10% more accurate than GridQTL, considerably faster, and better at inferring positions of recombination. GridQTL could not analyze all simulated replicates and did not estimate line origin for around a third of individuals at many positions. The study shows that triM has computational benefits and improved estimation over available algorithms and is valuable for analyzing the large datasets that will be standard in future.

Comparison of analyses of the QTLMAS XII common dataset. II: genome-wide association and fine mapping.

As part of the QTLMAS XII workshop, a simulated dataset was distributed and participants were invited to submit analyses of the data based on genome-wide association, fine mapping and genomic selection. We have evaluated the findings from the groups that reported fine mapping and genome-wide association (GWA) efforts to map quantitative trait loci (QTL). Generally the power to detect QTL was high and the Type 1 error was low. Estimates of QTL locations were generally very accurate. Some methods were much better than others at estimating QTL effects, and with some the accuracy depended on simulated effect size or minor allele frequency. There were also indications of bias in the effect estimates. No epistasis was simulated, but the two studies that included searches for epistasis reported several interacting loci, indicating a problem with controlling the Type I error rate in these analyses. Although this study is based on a single dataset, it indicates that there is a need to improve fine mapping and GWA methods with respect to estimation of genetic effects, appropriate choice of significance thresholds and analysis of epistasis.