Genetics of schizophrenia in the South African Xhosa.

Africa, the ancestral home of all modern humans, is the most informative continent for understanding the human genome and its contribution to complex disease. To better understand the genetics of schizophrenia, we studied the illness in the Xhosa population of South Africa, recruiting 909 cases and 917 age-, gender-, and residence-matched controls. Individuals with schizophrenia were significantly more likely than controls to harbor private, severely damaging mutations in genes that are critical to synaptic function, including neural circuitry mediated by the neurotransmitters glutamine, γ-aminobutyric acid, and dopamine. Schizophrenia is genetically highly heterogeneous, involving severe ultrarare mutations in genes that are critical to synaptic plasticity. The depth of genetic variation in Africa revealed this relationship with a moderate sample size and informed our understanding of the genetics of schizophrenia worldwide.

Challenges of Francisella classification exemplified by an atypical clinical isolate.

The accumulation of sequenced Francisella strains has made it increasingly apparent that the 16S rRNA gene alone is not enough to stratify the Francisella genus into precise and clinically useful classifications. Continued whole-genome sequencing of isolates will provide a larger base of knowledge for targeted approaches with broad applicability. Additionally, examination of genomic information on a case-by-case basis will help resolve outstanding questions regarding strain stratification. We report the complete genome sequence of a clinical isolate, designated here as F. novicida-like strain TCH2015, acquired from the lymph node of a 6-year-old male. Two features were atypical for F. novicida: exhibition of functional oxidase activity and additional gene content, including proposed virulence determinants. These differences, which could potentially impact virulence and clinical diagnosis, emphasize the need for more comprehensive methods to profile Francisella isolates. This study highlights the value of whole-genome sequencing, which will lead to a more robust database of environmental and clinical genomes and inform strategies to improve detection and classification of Francisella strains.

Whole-Exome Sequencing Identifies Novel Variants for Tooth Agenesis.

Tooth agenesis is a common craniofacial abnormality in humans and represents failure to develop 1 or more permanent teeth. Tooth agenesis is complex, and variations in about a dozen genes have been reported as contributing to the etiology. Here, we combined whole-exome sequencing, array-based genotyping, and linkage analysis to identify putative pathogenic variants in candidate disease genes for tooth agenesis in 10 multiplex Turkish families. Novel homozygous and heterozygous variants in LRP6, DKK1, LAMA3, and COL17A1 genes, as well as known variants in WNT10A, were identified as likely pathogenic in isolated tooth agenesis. Novel variants in KREMEN1 were identified as likely pathogenic in 2 families with suspected syndromic tooth agenesis. Variants in more than 1 gene were identified segregating with tooth agenesis in 2 families, suggesting oligogenic inheritance. Structural modeling of missense variants suggests deleterious effects to the encoded proteins. Functional analysis of an indel variant (c.3607+3_6del) in LRP6 suggested that the predicted resulting mRNA is subject to nonsense-mediated decay. Our results support a major role for WNT pathways genes in the etiology of tooth agenesis while revealing new candidate genes. Moreover, oligogenic cosegregation was suggestive for complex inheritance and potentially complex gene product interactions during development, contributing to improved understanding of the genetic etiology of familial tooth agenesis.

Comparative genomic study of arachnid immune systems indicates loss of beta-1,3-glucanase-related proteins and the immune deficiency pathway.

Analyses of arthropod genomes have shown that the genes in the different innate humoral immune responses are conserved. These genes encode proteins that are involved in immune signalling pathways that recognize pathogens and activate immune responses. These immune responses include phagocytosis, encapsulation of the pathogen and production of effector molecules for pathogen elimination. So far, most studies have focused on insects leaving other major arthropod groups largely unexplored. Here, we annotate the immune-related genes of six arachnid genomes and present evidence for a conserved pattern of some immune genes, but also evolutionary changes in the arachnid immune system. Specifically, our results suggest that the family of recognition molecules of beta-1,3-glucanase-related proteins (ßGRPs) and the genes from the immune deficiency (IMD) signalling pathway have been lost in a common ancestor of arachnids. These findings are consistent with previous work suggesting that the humoral immune effector proteins are constitutively produced in arachnids in contrast to insects, where these have to be induced. Further functional studies are needed to verify this. We further show that the full haemolymph clotting cascade found in the horseshoe crab is retrieved in most arachnid genomes. Tetranychus lacks at least one major component, although it is possible that this cascade could still function through recruitment of a different protein. The gel-forming protein in horseshoe crabs, coagulogen, was not recovered in any of the arachnid genomes; however, it is possible that the arachnid clot consists of a related protein, spätzle, that is present in all of the genomes.

High-Quality Draft Genome Sequence of Francisella tularensis subsp. holarctica Strain OR96-0246.

The bacterial pathogen Francisella tularensis was recently renewed as a tier-one select agent. F. tularensis subsp. tularensis (type A) and holarctica (type B) are of clinical relevance. Here, we report the complete genome of a virulent F. tularensis type B strain and describe its usefulness in comparative genomics.

Whole-exome sequencing points to considerable genetic heterogeneity of cerebral palsy.

Cerebral palsy (CP) is a common, clinically heterogeneous group of disorders affecting movement and posture. Its prevalence has changed little in 50 years and the causes remain largely unknown. The genetic contribution to CP causation has been predicted to be ~2%. We performed whole-exome sequencing of 183 cases with CP including both parents (98 cases) or one parent (67 cases) and 18 singleton cases (no parental DNA). We identified and validated 61 de novo protein-altering variants in 43 out of 98 (44%) case-parent trios. Initial prioritization of variants for causality was by mutation type, whether they were known or predicted to be deleterious and whether they occurred in known disease genes whose clinical spectrum overlaps CP. Further, prioritization used two multidimensional frameworks-the Residual Variation Intolerance Score and the Combined Annotation-dependent Depletion score. Ten de novo mutations in three previously identified disease genes (TUBA1A (n=2), SCN8A (n=1) and KDM5C (n=1)) and in six novel candidate CP genes (AGAP1, JHDM1D, MAST1, NAA35, RFX2 and WIPI2) were predicted to be potentially pathogenic for CP. In addition, we identified four predicted pathogenic, hemizygous variants on chromosome X in two known disease genes, L1CAM and PAK3, and in two novel candidate CP genes, CD99L2 and TENM1. In total, 14% of CP cases, by strict criteria, had a potentially disease-causing gene variant. Half were in novel genes. The genetic heterogeneity highlights the complexity of the genetic contribution to CP. Function and pathway studies are required to establish the causative role of these putative pathogenic CP genes.

Vaccine-associated varicella and rubella infections in severe combined immunodeficiency with isolated CD4 lymphocytopenia and mutations in IL7R detected by tandem whole exome sequencing and chromosomal microarray.

In areas without newborn screening for severe combined immunodeficiency (SCID), disease-defining infections may lead to diagnosis, and in some cases, may not be identified prior to the first year of life. We describe a female infant who presented with disseminated vaccine-acquired varicella (VZV) and vaccine-acquired rubella infections at 13 months of age. Immunological evaluations demonstrated neutropenia, isolated CD4 lymphocytopenia, the presence of CD8(+) T cells, poor lymphocyte proliferation, hypergammaglobulinaemia and poor specific antibody production to VZV infection and routine immunizations. A combination of whole exome sequencing and custom-designed chromosomal microarray with exon coverage of primary immunodeficiency genes detected compound heterozygous mutations (one single nucleotide variant and one intragenic copy number variant involving one exon) within the IL7R gene. Mosaicism for wild-type allele (20-30%) was detected in pretransplant blood and buccal DNA and maternal engraftment (5-10%) demonstrated in pretransplant blood DNA. This may be responsible for the patient's unusual immunological phenotype compared to classical interleukin (IL)-7Rα deficiency. Disseminated VZV was controlled with anti-viral and immune-based therapy, and umbilical cord blood stem cell transplantation was successful. Retrospectively performed T cell receptor excision circle (TREC) analyses completed on neonatal Guthrie cards identified absent TREC. This case emphasizes the danger of live viral vaccination in severe combined immunodeficiency (SCID) patients and the importance of newborn screening to identify patients prior to high-risk exposures. It also illustrates the value of aggressive pathogen identification and treatment, the influence newborn screening can have on morbidity and mortality and the significant impact of newer genomic diagnostic tools in identifying the underlying genetic aetiology for SCID patients.

Next-generation sequencing study finds an excess of rare, coding single-nucleotide variants of ADAMTS13 in patients with deep vein thrombosis.

BACKGROUND: The considerable genetic predisposition to deep vein thrombosis (DVT) is only partially accounted for by known genetic risk variants. Rare single-nucleotide variants (SNVs) of the coding areas of hemostatic genes may explain part of this missing heritability. The ADAMTS13 and VWF genes encode two interconnected proteins with fundamental hemostatic functions, the disruption of which may result in thrombosis. OBJECTIVES: To study the distribution and burden of rare coding SNVs of ADAMTS13 and VWF found by sequencing in cases and controls of DVT. PATIENTS/METHODS: The protein-coding areas of 186 hemostatic/proinflammatory genes were sequenced by next-generation technology in 94 thrombophilia-negative patients with DVT and 98 controls. Gene-specific information on ADAMTS13 and VWF was used to study the association between DVT and rare coding SNVs of the two genes. RESULTS: More than 70 billion base pairs of raw sequence data were produced to sequence the 700-kb target area with a median redundancy of × 45 in 192 individuals. Most of the 4366 SNVs identified were rare and non-synonymous, indicating pathogenetic potential. Rare (frequency of < 1%) and low-frequency (< 5%) coding SNVs of ADAMTS13 were associated with DVT (prevalence 17% vs. 4%; odds ratio [OR] 4.8 and 95% confidence interval [CI] 1.6-15.0 for rare coding; prevalence 36% vs. 23%, OR 1.9 and 95% CI 1.0-3.5 for low-frequency coding). Patients with rare coding SNVs of ADAMTS13 had lower plasma levels of ADAMTS-13 activity than patients without them. SNVs of VWF were not associated with DVT. CONCLUSIONS: We found an excess of rare coding SNVs of the ADAMTS13 gene in patients with DVT.

Next-generation sequencing for disorders of low and high bone mineral density.

UNLABELLED: To achieve an efficient molecular diagnosis of osteogenesis imperfecta (OI), Ehlers-Danlos syndrome (EDS), and osteopetrosis (OPT), we designed a next-generation sequencing (NGS) platform to sequence 34 genes. We validated this platform on known cases and have successfully identified the causative mutation in most patients without a prior molecular diagnosis. INTRODUCTION: Osteogenesis imperfecta, Ehlers-Danlos syndrome, and osteopetrosis are collectively common inherited skeletal diseases. Evaluation of subjects with these conditions often includes molecular testing which has important counseling and therapeutic and sometimes legal implications. Since several different genes have been implicated in these conditions, Sanger sequencing of each gene can be a prohibitively expensive and time-consuming way to reach a molecular diagnosis. METHODS: In order to circumvent these problems, we have designed and tested a NGS platform that would allow simultaneous sequencing on a single diagnostic platform of different genes implicated in OI, OPT, EDS, and other inherited conditions, leading to low or high bone mineral density. We used a liquid-phase probe library that captures 602 exons (~100 kb) of 34 selected genes and have applied it to test clinical samples from patients with bone disorders. RESULTS: NGS of the captured exons by Illumina HiSeq 2000 resulted in an average coverage of over 900X. The platform was successfully validated by identifying mutations in six patients with known mutations. Moreover, in four patients with OI or OPT without a prior molecular diagnosis, the assay was able to detect the causative mutations. CONCLUSIONS: In conclusion, our NGS panel provides a fast and accurate method to arrive at a molecular diagnosis in most patients with inherited high or low bone mineral density disorders.

Characterizing polymorphisms and allelic diversity of von Willebrand factor gene in the 1000 Genomes.

BACKGROUND: The von Willebrand factor (VWF) gene is highly polymorphic, with variants correlated with VWF antigen levels, adhesion activity, clearance and factor VIII binding. VWF mutations are detected in patients with von Willebrand disease (VWD), whereas polymorphic variants could be associated with thrombosis. However, information on the ethnic diversity of VWF variants and their association with diseases is limited. OBJECTIVES: To characterize novel VWF variants from different ethnicities in the general population. PATIENTS/METHODS: We analyzed samples from 1092 subjects of 14 ethnicities available in the 1000 Genomes database for VWF variants and their potential functional impacts. RESULTS: We identified 2728 SNPs and 91 insertions and deletions that had a high level of ethnic diversity, with Africans having the highest number of variants. The highest level of diversity was found in the D' and D2 domains. Among 94 non-synonymous variants, 31 were predicted to be deleterious, including 19 that were previously associated with VWD. Most of these 'VWD variants' had allele frequencies consistent with disease incidence in European subjects, but some had a significantly higher frequency in other ethnicities. The mutations R2185Q, H817Q and M740I associated with type 1 and type 2N VWD were present in more than 13% of African subjects. CONCLUSIONS: These results highlight the complexity of VWF variations in different ethnic groups and emphasize the importance of interrogating variations on multiple ethnic backgrounds for associations with bleeding and thrombosis.

CRHR1 genotypes, neural circuits and the diathesis for anxiety and depression.

The corticotrophin-releasing hormone (CRH) system integrates the stress response and is associated with stress-related psychopathology. Previous reports have identified interactions between childhood trauma and sequence variation in the CRH receptor 1 gene (CRHR1) that increase risk for affective disorders. However, the underlying mechanisms that connect variation in CRHR1 to psychopathology are unknown. To explore potential mechanisms, we used a validated rhesus macaque model to investigate association between genetic variation in CRHR1, anxious temperament (AT) and brain metabolic activity. In young rhesus monkeys, AT is analogous to the childhood risk phenotype that predicts the development of human anxiety and depressive disorders. Regional brain metabolism was assessed with (18)F-labeled fluoro-2-deoxyglucose (FDG) positron emission tomography in 236 young, normally reared macaques that were also characterized for AT. We show that single nucleotide polymorphisms (SNPs) affecting exon 6 of CRHR1 influence both AT and metabolic activity in the anterior hippocampus and amygdala, components of the neural circuit underlying AT. We also find evidence for association between SNPs in CRHR1 and metabolism in the intraparietal sulcus and precuneus. These translational data suggest that genetic variation in CRHR1 affects the risk for affective disorders by influencing the function of the neural circuit underlying AT and that differences in gene expression or the protein sequence involving exon 6 may be important. These results suggest that variation in CRHR1 may influence brain function before any childhood adversity and may be a diathesis for the interaction between CRHR1 genotypes and childhood trauma reported to affect human psychopathology.

An outbreak of Cyclospora infection on a cruise ship.

In 2010, an outbreak of cyclosporiasis affected passengers and crew on two successive voyages of a cruise ship that departed from and returned to Fremantle, Australia. There were 73 laboratory-confirmed and 241 suspected cases of Cyclospora infection reported in passengers and crew from the combined cruises. A case-control study performed in crew members found that illness was associated with eating items of fresh produce served onboard the ship, but the study was unable conclusively to identify the responsible food(s). It is likely that one or more of the fresh produce items taken onboard at a south-east Asian port during the first cruise was contaminated. If fresh produce supplied to cruise ships is sourced from countries or regions where Cyclospora is endemic, robust standards of food production and hygiene should be applied to the supply chain.

Genomic sequencing of key genes in mouse pancreatic cancer cells.

Pancreatic cancer is a multiple genetic disorder with many mutations identified during the progression. Two mouse pancreatic cancer cell lines were established which showed different phenotype in vivo: a non-metastatic cell line, Panc02, and a highly metastatic cell line, Panc02-H7, a derivative of Panc02. In order to investigate whether the genetic mutations of key genes in pancreatic cancer such as KRAS, TP53 (p53), CDKN2A (p16), SMAD4, ZIP4, and PDX-1 contribute to the phenotypic difference of these two mouse pancreatic cancer cells, we sequenced the exonic regions of these key genes in both cell lines and in the normal syngeneic mouse pancreas and compared them with the reference mouse genome sequence. The exons of KRAS, SMAD4, CDKN2A (p16), TP53 (p53), ZIP4, and PDX-1 genes were amplified and the genotype of these genes was determined by Sanger sequencing. The sequences were analyzed with Sequencher software. A mutation in SMAD4 was identified in both cell lines. This homozygote G to T mutation in the first position of codon 174 (GAA) generated a stop codon resulting in the translation of a truncated protein. Further functional analysis indicates that different TGF-ß/SMAD signaling pathways were involved in those two mouse cell lines, which may explain the phonotypic difference between the two cells. A single nucleotide polymorphism (SNP) in KRAS gene (TAT to TAC at codon 32) was also identified in the normal pancreas DNA of the syngenic mouse and in both derived tumoral Panc02 and Panc02-H7 cells. No mutation or SNP was found in CDKN2A (p16), TP53 (p53), ZIP4, and PDX-1 genes in these two cell lines. The absence of mutations in genes such as KRAS, TP53, and CDKN2A, which are considered as key genes in the development of human pancreatic cancer suggests that SMAD4 might play a central and decisive role in mouse pancreatic cancer. These results also suggest that other mechanisms are involved in the substantial phenotypic difference between these two mouse pancreatic cancer cell lines. Further studies are warranted to elucidate the molecular pathways that lead to the aggressive metastatic potential of Panc02-H7.

A full-length cDNA resource for the pea aphid, Acyrthosiphon pisum.

Large collections of full-length cDNAs are important resources for genome annotation and functional genomics. We report the creation of a collection of 50 599 full-length cDNA clones from the pea aphid, Acyrthosiphon pisum. Sequencing from 5' and 3' ends of the clones generated 97 828 high-quality expressed sequence tags, representing approximately 9000 genes. These sequences were imported to AphidBase and are shown to play crucial roles in both automatic gene prediction and manual annotation. Our detailed analyses demonstrated that the full-length cDNAs can further improve gene models and can even identify novel genes that are not included in the current version of the official gene set. This full-length cDNA collection can be utilized for a wide variety of functional studies, serving as a community resource for the study of the functional genomics of the pea aphid.

Comparison of algal and fish sources on the oxidative stability of poultry meat and its enrichment with omega-3 polyunsaturated fatty acids.

Human consumption of long-chain n-3 polyunsaturated fatty acids (LC n-3 PUFA) is below recommendations, and enriching chicken meat (by incorporating LC n-3 PUFA into broiler diets) is a viable means of increasing consumption. Fish oil is the most common LC n-3 PUFA supplement used but is unsustainable and reduces the oxidative stability of the mat. The objective of this experiment was to compare fresh fish oil (FFO) with fish oil encapsulated (EFO) in a gelatin matrix (to maintain its oxidative stability) and algal biomass at a low (LAG, 11), medium (MAG, 22), or high (HAG, 33 g/kg of diet) level of inclusion. The C22:6n-3 contents of the FFO, EFO, and MAG diets were equal. A control (CON) diet using blended vegetable oil was also made. As-hatched 1-d-old Ross 308 broilers (144) were reared (21 d) on a common starter diet then allocated to treatment pens (4 pens per treatment, 6 birds per pen) and fed treatment diets for 21 d before being slaughtered. Breast and leg meat was analyzed (per pen) for fatty acids, and cooked samples (2 pens per treatment) were analyzed for volatile aldehydes. Concentrations (mg/100 g of meat) of C20:5n-3, C22:5n-3, and C22:6n-3 were (respectively) CON: 4, 15, 24; FFO: 31, 46, 129; EFO: 18, 27, 122; LAG: 9, 19, 111; MAG: 6, 16, 147; and HAG: 9, 14, 187 (SEM: 2.4, 3.6, 13.1) in breast meat and CON: 4, 12, 9; FFO: 58, 56, 132; EFO: 63, 49, 153; LAG: 13, 14, 101; MAG: 11, 15, 102; HAG: 37, 37, 203 (SEM: 7.8, 6.7, 14.4) in leg meat. Cooked EFO and HAG leg meat was more oxidized (5.2 mg of hexanal/kg of meat) than the other meats (mean 2.2 mg/kg, SEM 0.63). It is concluded that algal biomass is as effective as fish oil at enriching broiler diets with C22:6 LC n-3 PUFA, and at equal C22:6n-3 contents, there is no significant difference between these 2 supplements on the oxidative stability of the meat that is produced.

Rats in the genomic era.

The rat genome project and the resources that it has generated are transforming the translation of rat biology to human medicine. The rat genome was sequenced to a high quality "draft," the structure and location of the genes were predicted, and a global assessment was published (Gibbs RA et al., Nature 428: 493-521, 2004). Since that time, researchers have made use of the genome sequence and annotations and related resources. We take this opportunity to review the currently available rat genome resources and to discuss the progress and future plans for the rat genome.