The zebrafish lens opaque (lop) mutant was previously isolated in a genetic screen and shown to lack rod and cone photoreceptors and exhibit lens opacity, or cataract, at 7 days post-fertilization (dpf). In this manuscript, we provide four different lines of evidence demonstrating that the lop phenotype results from a defect in the cdipt (phosphatidylinositol (PI) synthase; CDP-diacylglycerol-inositol 3-phosphatidyltransferase) gene. First, DNA sequence analysis revealed that the lop mutant contained a missense mutation in the lop open reading frame, which yields a nonconservative amino acid substitution (Ser-111-Cys) within the PI synthase catalytic domain. Second, morpholino-mediated knockdown of the cdipt-encoded PI synthase protein phenocopied the cdipt(lop/lop) mutant, with abnormal lens epithelial and secondary fiber cell morphologies and reduced numbers of photoreceptors. Third, microinjection of in vitro transcribed, wild-type cdipt mRNA into 1-4 cell stage cdipt(lop/lop) embryos significantly reduced the percentage of larvae displaying lens opacity at 7 dpf. Fourth, a cdipt retroviral-insertion allele, cdipt(hi559), exhibited similar lens and retinal abnormalities and failed to complement the cdipt(lop) mutant phenotype. To determine the initial cellular defects associated with the cdipt mutant, we examined homozygous cdipt(hi559/hi559) mutants prior to gross lens opacification at 6 dpf. The cdipt(hi559/hi559) mutants first exhibited photoreceptor layer disruption and photoreceptor cell death at 3 and 4 dpf, respectively, followed by lens dismorphogenesis by 5 dpf. RT-PCR revealed that the cdipt gene is maternally expressed and continues to be transcribed throughout development and into adulthood, in a wide variety of tissues. Using an anti-zebrafish PI synthase polyclonal antiserum, we localized the protein throughout the developing eye, including the photoreceptor layer and lens cortical secondary fiber cells. As expected, the polyclonal antiserum revealed that the PI synthase protein was reduced in amount in both the cdipt(lop/lop) and cdipt(hi559/hi559) mutants. Furthermore, we used a heterologous yeast phenotypic complementation assay to confirm that the wild-type zebrafish cdipt allele encodes functional PI synthase activity. Taken together, the cdipt-encoded PI synthase is required for survival of photoreceptor cells and lens epithelial and secondary cortical fiber cells. These zebrafish cdipt alleles represent excellent in vivo genetic tools to study the role of phosphatidylinositol and its phosphorylated derivatives in lens and photoreceptor development and maintenance.

Drug exposure during critical periods of development is known to have lasting effects, increasing one's risk for developing mental health disorders. Emerging evidence has also indicated the possibility for drug exposure to even impact subsequent generations. Our previous work demonstrated that adolescent exposure to Δ(9)-tetrahydrocannabinol (THC), the main psychoactive component of marijuana (Cannabis sativa), in a Long-Evans rat model affects reward-related behavior and gene regulation in the subsequent (F1) generation unexposed to the drug. Questions, however, remained regarding potential epigenetic consequences. In the current study, using the same rat model, we employed Enhanced Reduced Representation Bisulfite Sequencing to interrogate the epigenome of the nucleus accumbens, a key brain area involved in reward processing. This analysis compared 16 animals with parental THC exposure and 16 without to characterize relevant systems-level changes in DNA methylation. We identified 1027 differentially methylated regions (DMRs) associated with parental THC exposure in F1 adults, each represented by multiple CpGs. These DMRs fell predominantly within introns, exons, and intergenic intervals, while showing a significant depletion in gene promoters. From these, we identified a network of DMR-associated genes involved in glutamatergic synaptic regulation, which also exhibited altered mRNA expression in the nucleus accumbens. These data provide novel insight into drug-related cross-generational epigenetic effects, and serve as a useful resource for investigators to explore novel neurobiological systems underlying drug abuse vulnerability.

Certain human traits such as neurodevelopmental disorders (NDs) and congenital anomalies (CAs) are believed to be primarily genetic in origin. However, even after whole-genome sequencing (WGS), a substantial fraction of such disorders remain unexplained. We hypothesize that some cases of ND-CA are caused by aberrant DNA methylation leading to dysregulated genome function. Comparing DNA methylation profiles from 489 individuals with ND-CAs against 1534 controls, we identify epivariations as a frequent occurrence in the human genome. De novo epivariations are significantly enriched in cases, while RNAseq analysis shows that epivariations often have an impact on gene expression comparable to loss-of-function mutations. Additionally, we detect and replicate an enrichment of rare sequence mutations overlapping CTCF binding sites close to epivariations, providing a rationale for interpreting non-coding variation. We propose that epivariations contribute to the pathogenesis of some patients with unexplained ND-CAs, and as such likely have diagnostic relevance.

Gut microbiota composition is associated with human and rodent Plasmodium infections, yet the mechanism by which gut microbiota affects the severity of malaria remains unknown. Humoral immunity is critical in mediating the clearance of Plasmodium blood stage infections, prompting the hypothesis that mice with gut microbiota-dependent decreases in parasite burden exhibit better germinal center (GC) responses. In support of this hypothesis, mice with a low parasite burden exhibit increases in GC B cell numbers and parasite-specific antibody titers, as well as better maintenance of GC structures and a more targeted, qualitatively different antibody response. This enhanced humoral immunity affects memory, as mice with a low parasite burden exhibit robust protection against challenge with a heterologous, lethal Plasmodium species. These results demonstrate that gut microbiota composition influences the biology of spleen GCs as well as the titer and repertoire of parasite-specific antibodies, identifying potential approaches to develop optimal treatments for malaria.

High-throughput sequencing of the adaptive immune receptor repertoire (AIRR-seq) is providing unprecedented insights into the immune response to disease and into the development of immune disorders. The accurate interpretation of AIRR-seq data depends on the existence of comprehensive germline gene reference sets. Current sets are known to be incomplete and unrepresentative of the degree of polymorphism and diversity in human and animal populations. A key issue is the complexity of the genomic regions in which they lie, which, because of the presence of multiple repeats, insertions and deletions, have not proved tractable with short-read whole genome sequencing. Recently, tools and methods for inferring such gene sequences from AIRR-seq datasets have become available, and a community approach has been developed for the expert review and publication of such inferences. Here, we present OGRDB, the Open Germline Receptor Database (https://ogrdb.airr-community.org), a public resource for the submission, review and publication of previously unknown receptor germline sequences together with supporting evidence.

An incomplete ascertainment of genetic variation within the highly polymorphic immunoglobulin heavy chain locus (IGH) has hindered our ability to define genetic factors that influence antibody-mediated processes. Due to locus complexity, standard high-throughput approaches have failed to accurately and comprehensively capture IGH polymorphism. As a result, the locus has only been fully characterized two times, severely limiting our knowledge of human IGH diversity. Here, we combine targeted long-read sequencing with a novel bioinformatics tool, IGenotyper, to fully characterize IGH variation in a haplotype-specific manner. We apply this approach to eight human samples, including a haploid cell line and two mother-father-child trios, and demonstrate the ability to generate high-quality assemblies (>98% complete and >99% accurate), genotypes, and gene annotations, identifying 2 novel structural variants and 15 novel IGH alleles. We show multiplexing allows for scaling of the approach without impacting data quality, and that our genotype call sets are more accurate than short-read (>35% increase in true positives and >97% decrease in false-positives) and array/imputation-based datasets. This framework establishes a desperately needed foundation for leveraging IG genomic data to study population-level variation in antibody-mediated immunity, critical for bettering our understanding of disease risk, and responses to vaccines and therapeutics.

Despite the remarkable success of immunotherapy in many types of cancer, pancreatic ductal adenocarcinoma has yet to benefit. Innate immune cells are critical to anti-tumor immunosurveillance and recent studies have revealed that these populations possess a form of memory, termed trained innate immunity, which occurs through transcriptomic, epigenetic, and metabolic reprograming. Here we demonstrate that yeast-derived particulate β-glucan, an inducer of trained immunity, traffics to the pancreas, which causes a CCR2-dependent influx of monocytes/macrophages to the pancreas that display features of trained immunity. These cells can be activated upon exposure to tumor cells and tumor-derived factors, and show enhanced cytotoxicity against pancreatic tumor cells. In orthotopic models of pancreatic ductal adenocarcinoma, β-glucan treated mice show significantly reduced tumor burden and prolonged survival, which is further enhanced when combined with immunotherapy. These findings characterize the dynamic mechanisms and localization of peripheral trained immunity and identify an application of trained immunity to cancer.

Conventional immunization strategies will likely be insufficient for the development of a broadly neutralizing antibody (bnAb) vaccine for HIV or other difficult pathogens because of the immunological hurdles posed, including B cell immunodominance and germinal center (GC) quantity and quality. We found that two independent methods of slow delivery immunization of rhesus monkeys (RMs) resulted in more robust T follicular helper (TFH) cell responses and GC B cells with improved Env-binding, tracked by longitudinal fine needle aspirates. Improved GCs correlated with the development of >20-fold higher titers of autologous nAbs. Using a new RM genomic immunoglobulin locus reference, we identified differential IgV gene use between immunization modalities. Ab mapping demonstrated targeting of immunodominant non-neutralizing epitopes by conventional bolus-immunized animals, whereas slow delivery-immunized animals targeted a more diverse set of epitopes. Thus, alternative immunization strategies can enhance nAb development by altering GCs and modulating the immunodominance of non-neutralizing epitopes.

Alzheimer's disease affects ~13% of people in the United States 65 years and older, making it the most common neurodegenerative disorder. Recent work has identified roles for environmental, genetic, and epigenetic factors in Alzheimer's disease risk.

Many models of evolution by sexual selection predict a coevolution of sensory systems and mate preferences, but the genomic architecture (number and arrangement of contributing loci) underlying these characters could constrain this coevolution. Here, we examine how the genomic organization and evolution of the opsin genes (responsible for tuning color vision) can influence the evolutionary trajectory of sexually selected traits across 15 species in the family Poeciliidae, which includes classic systems for studies of color-mediated sexual selection such as guppies, swordtails, and mollies. Although male coloration patterns and the importance of this coloration in female mate choice vary widely within and among genera, sequencing revealed low variability at amino acid sites that tune Long Wavelength-Sensitive (LWS) opsins in this speciose family. Although most opsin genes in these species appear to have evolved along traditional mutation-selection dynamics, we identified high rates of gene conversion between two of the LWS loci (LWS-1 and LWS-3), likely due to the inverted tandem repeat nature of these genes. Yet members of the subgenus Lebistes appear to resist LWS gene conversion. The LWS opsins are responsible for detecting and discriminating red and orange coloration-a key sexually selected trait in members of the subgenus Lebistes. Taken together these results suggest selection is acting against the homogenizing effects of gene conversion to maintain LWS-1/LWS-3 differences within this subgenus.

No abstract available

Inference of germline polymorphisms in immunoglobulin genes from B cell receptor repertoires is complicated by somatic hypermutations, sequencing/PCR errors, and by varying length of reference alleles. The light chain inference is particularly challenging owing to large gene duplications and absence of D genes. We analyzed the light chain cDNA sequences from naïve B cell receptor repertoires from 100 individuals. We optimized light chain allele inference by tweaking parameters of the TIgGER functions, extending the germline reference sequences, and establishing mismatch frequency patterns at polymorphic positions to filter out false-positive candidates. We identified 48 previously unreported variants of light chain variable genes. We selected 14 variants for validation and successfully validated 11 by Sanger sequencing. Clustering of light chain 5'UTR, L-PART1, and L-PART2 revealed partial intron retention in 11 kappa and 9 lambda V alleles. Our results provide insight into germline variation in human light chain immunoglobulin loci.

In adaptive immune receptor repertoire analysis, determining the germline variable (V) allele associated with each T- and B-cell receptor sequence is a crucial step. This process is highly impacted by allele annotations. Aligning sequences, assigning them to specific germline alleles, and inferring individual genotypes are challenging when the repertoire is highly mutated, or sequence reads do not cover the whole V region. Here, we propose an alternative naming scheme for the V alleles, as well as a novel method to infer individual genotypes. We demonstrate the strengths of the two by comparing their outcomes to other genotype inference methods. We validate the genotype approach with independent genomic long-read data. The naming scheme is compatible with current annotation tools and pipelines. Analysis results can be converted from the proposed naming scheme to the nomenclature determined by the International Union of Immunological Societies (IUIS). Both the naming scheme and the genotype procedure are implemented in a freely available R package (PIgLET https://bitbucket.org/yaarilab/piglet). To allow researchers to further explore the approach on real data and to adapt it for their uses, we also created an interactive website (https://yaarilab.github.io/IGHV_reference_book).

GATA binding protein 3 (GATA3) is a regulator of mammary luminal cell differentiation, and an estrogen receptor (ER) associated marker in breast cancer. Tumor suppressor functions of GATA3 have been demonstrated primarily in basal-like breast cancers. Here, we focused on its function in luminal breast cancer, where GATA3 is frequently mutated, and its levels are significantly elevated.

Multiple sclerosis (MS) is a complex disease with underlying genetic and environmental factors. Although the contribution of alleles within the major histocompatibility complex (MHC) are known to exert strong effects on MS risk, much remains to be learned about the contributions of loci with more modest effects identified by genome-wide association studies (GWASs), as well as loci that remain undiscovered. We use a recently developed method to estimate the proportion of variance in disease liability explained by 475,806 single nucleotide polymorphisms (SNPs) genotyped in 1,854 MS cases and 5,164 controls. We reveal that ~30% of MS genetic liability is explained by SNPs in this dataset, the majority of which is accounted for by common variants. These results suggest that the unaccounted for proportion could be explained by variants that are in imperfect linkage disequilibrium with common GWAS SNPs, highlighting the potential importance of rare variants in the susceptibility to MS.

Long noncoding RNAs (lncRNAs) are major regulators of many cellular processes including cell cycle progression and tumorigenesis. In this study, we identify a novel lncRNA, MA-linc1, and reveal its effects on cell cycle progression and cancer growth. Inhibition of MA-linc1 expression alters cell cycle distribution, leading to a decrease in the number of G1 cells and a concomitant increase in all other stages of the cell cycle, and in particular G2/M, suggesting its involvement in the regulation of M phase. Accordingly, knock down of MA-linc1 inhibits M phase exit upon release from a mitotic block. We further demonstrate that MA-linc1 predominantly functions in cis to repress expression of its neighboring gene, Purα, which is often deleted in human cancers and whose ectopic expression inhibits cell cycle progression. Knock down of Purα partially rescues the MA-linc1 dependent inhibition of M phase exit. In agreement with its suggested role in M phase, inhibition of MA-linc1 enhances apoptotic cell death induced by the antimitotic drug, Paclitaxel and this enhancement of apoptosis is rescued by Purα knockdown. Furthermore, high levels of MA-linc1 are associated with reduced survival in human breast and lung cancer patients.Taken together, our data identify MA-linc1 as a novel lncRNA regulator of cell cycle and demonstrate its potential role in cancer progression and treatment.

The immunoglobulin heavy-chain locus (IGH) encodes variable (IGHV), diversity (IGHD), joining (IGHJ), and constant (IGHC) genes and is responsible for antibody heavy-chain biosynthesis, which is vital to the adaptive immune response. Programmed V-(D)-J somatic rearrangement and the complex duplicated nature of the locus have impeded attempts to reconcile its genomic organization based on traditional B-lymphocyte derived genetic material. As a result, sequence descriptions of germline variation within IGHV are lacking, haplotype inference using traditional linkage disequilibrium methods has been difficult, and the human genome reference assembly is missing several expressed IGHV genes. By using a hydatidiform mole BAC clone resource, we present the most complete haplotype of IGHV, IGHD, and IGHJ gene regions derived from a single chromosome, representing an alternate assembly of ∼1 Mbp of high-quality finished sequence. From this we add 101 kbp of previously uncharacterized sequence, including functional IGHV genes, and characterize four large germline copy-number variants (CNVs). In addition to this germline reference, we identify and characterize eight CNV-containing haplotypes from a panel of nine diploid genomes of diverse ethnic origin, discovering previously unmapped IGHV genes and an additional 121 kbp of insertion sequence. We genotype four of these CNVs by using PCR in 425 individuals from nine human populations. We find that all four are highly polymorphic and show considerable evidence of stratification (Fst = 0.3-0.5), with the greatest differences observed between African and Asian populations. These CNVs exhibit weak linkage disequilibrium with SNPs from two commercial arrays in most of the populations tested.

Both genetic and environmental factors contribute to the aetiology of multiple sclerosis (MS). More than 50 genomic regions have been associated with MS susceptibility and vitamin D status also influences the risk of this complex disease. However, how these factors interact in disease causation is unclear. We aimed to investigate the relationship between vitamin D receptor (VDR) binding in lymphoblastoid cell lines (LCLs), chromatin states in LCLs and MS-associated genomic regions. Using the Genomic Hyperbrowser, we found that VDR-binding regions overlapped with active regulatory regions [active promoter (AP) and strong enhancer (SE)] in LCLs more than expected by chance [45.3-fold enrichment for SE (P < 2.0e-05) and 63.41-fold enrichment for AP (P < 2.0e-05)]. Approximately 77% of VDR regions were covered by either AP or SE elements. The overlap between VDR binding and regulatory elements was significantly greater in LCLs than in non-immune cells (P < 2.0e-05). VDR binding also occurred within MS regions more than expected by chance (3.7-fold enrichment, P < 2.0e-05). Furthermore, regions of joint overlap SE-VDR and AP-VDR were even more enriched within MS regions and near to several disease-associated genes. These findings provide relevant insights into how vitamin D influences the immune system and the risk of MS through VDR interactions with the chromatin state inside MS regions. Furthermore, the data provide additional evidence for an important role played by B cells in MS. Further analyses in other immune cell types and functional studies are warranted to fully elucidate the role of vitamin D in the immune system.

Germline variation at immunoglobulin (IG) loci is critical for pathogen-mediated immunity, but establishing complete haplotype sequences in these regions has been problematic because of complex sequence architecture and diploid source DNA. We sequenced BAC clones from the effectively haploid human hydatidiform mole cell line, CHM1htert, across the light chain IG loci, kappa (IGK) and lambda (IGL), creating single haplotype representations of these regions. The IGL haplotype generated here is 1.25 Mb of contiguous sequence, including four novel IGLV alleles, one novel IGLC allele, and an 11.9-kb insertion. The CH17 IGK haplotype consists of two 644 kb proximal and 466 kb distal contigs separated by a large gap of unknown size; these assemblies added 49 kb of unique sequence extending into this gap. Our analysis also resulted in the characterization of seven novel IGKV alleles and a 16.7-kb region exhibiting signatures of interlocus sequence exchange between distal and proximal IGKV gene clusters. Genetic diversity in IGK/IGL was compared with that of the IG heavy chain (IGH) locus within the same haploid genome, revealing threefold (IGK) and sixfold (IGL) higher diversity in the IGH locus, potentially associated with increased levels of segmental duplication and the telomeric location of IGH.

An important challenge in vaccine development is to figure out why a vaccine succeeds in some individuals and fails in others. Although antibody repertoires hold the key to answering this question, there have been very few personalized immunogenomics studies so far aimed at revealing how variations in immunoglobulin genes affect a vaccine response. We conducted an immunosequencing study of 204 calves vaccinated against bovine respiratory disease (BRD) with the goal to reveal variations in immunoglobulin genes and somatic hypermutations that impact the efficacy of vaccine response. Our study represents the largest longitudinal personalized immunogenomics study reported to date across all species, including humans. To analyze the generated data set, we developed an algorithm for identifying variations of the immunoglobulin genes (as well as frequent somatic hypermutations) that affect various features of the antibody repertoire and titers of neutralizing antibodies. In contrast to relatively short human antibodies, cattle have a large fraction of ultralong antibodies that have opened new therapeutic opportunities. Our study reveals that ultralong antibodies are a key component of the immune response against the costliest disease of beef cattle in North America. The detected variants of the cattle immunoglobulin genes, which are implicated in the success/failure of the BRD vaccine, have the potential to direct the selection of individual cattle for ongoing breeding programs.