A genome-wide survey of the prevalence and evolutionary forces acting on human nonsense SNPs

Human-specific elimination of epithelial Siglec-XII suppresses the risk of inflammation-driven colorectal cancers

Carcinomas are common in humans but rare among closely related "great apes." Plausible explanations, including human-specific genomic alterations affecting the biology of sialic acids, are proposed, but causality remains unproven. Here, an integrated evolutionary genetics-phenome-transcriptome approach studied the role of SIGLEC12 gene (encoding Siglec-XII) in epithelial transformation and cancer. Exogenous expression of the protein in cell lines and genetically engineered mice recapitulated approximately 30% of the human population in whom the protein is expressed in a form that cannot bind ligand because of a fixed, homozygous, human-universal missense mutation. Siglec-XII-null cells/mice recapitulated the remaining approximately 70% of the human population in whom an additional polymorphic frameshift mutation eliminates the entire protein. Siglec-XII expression drove several pro-oncogenic phenotypes in cell lines and increased tumor burden in mice challenged with chemical carcinogen and inflammation. Transcriptomic studies yielded a 29-gene signature of Siglec-XII-positive disease and when used as a computational tool for navigating human data sets, pinpointed with surprising precision that SIGLEC12 expression (model) recapitulates a very specific type of colorectal carcinomas (disease) that is associated with mismatch-repair defects and inflammation, disproportionately affects European Americans, and carries a favorable prognosis. They revealed a hitherto-unknown evolutionary genetic mechanism for an ethnic/environmental predisposition of carcinogenesis.

Analysis of SIGLEC12 expression, immunomodulation and prognostic value in renal cancer using multiomic databases

Siglecs belong to a family of immune regulatory receptors predominantly found on hematopoietic cells. They interact with Sia, resulting in the activation or inhibition of the immune response. Previous reports have suggested that the SIGLEC12 gene, which encodes the Siglec-XII protein, is expressed in the epithelial tissues and upregulated in carcinomas. However, studies deciphering the role of Siglec-XII in renal cancer (RC) are still unavailable, and here we provide insights on this question. We conducted expression analysis using the Human Protein Atlas and UALCAN databases. The impact of SIGLEC12 on RC prognosis was determined using the KM plotter, and an assessment of immune infiltration with SIGLEC12 was performed using the TIMER database. GSEA was conducted to identify the pathways affected by SIGLEC12. Finally, using GeneMania, we identified Siglec-XII interacting proteins. Our findings indicated that macrophages express SIGLEC12 in the kidney. Furthermore, we hypothesize that Siglec-XII expression might be involved in the increase of primary RC, but this effect may not be dependent on the age of the patient. In the tumour microenvironment, oncogenic pathways appeared to be upregulated by SIGLEC12. Similarly, our analysis suggested that SIGLEC12-related kidney renal papillary cell carcinomas may be more suitable for targeted immunotherapy, such as CTLA-4 and PD-1/PD-L1 inhibitors. These preliminary results suggested that high expression of SIGLEC12 is associated with poor prognosis for RC. Future studies to assess its clinical utility are necessitated.

rs67047829 genotypes of ERV3-1/ZNF117 are associated with lower body mass index in the Polish population

There is now substantial evidence that zinc-finger proteins are implicated in adiposity. Aims were to datamine for high-frequency (near-neutral selection) pretermination-codon (PTC) single-nucleotide polymorphisms (SNPs; n = 141) from a database with > 550,000 variants and analyze possible association with body mass index in a large Polish sample (n = 5757). BMI was regressed (males/females together or separately) against genetic models. Regression for rs67047829 uncovered an interaction-independent association with BMI with both sexes together: mean ± standard deviation, kg/m2: [G];[G], 25.4 ± 4.59 (n = 3650); [G](;)[A], 25.0 ± 4.28 (n = 731); [A];[A], 23.4 ± 3.60 (n = 44); additive model adjusted for age and sex: p = 4.08 × 10-5; beta: - 0.0458, 95% confidence interval (CI) - 0.0732 : - 0.0183; surviving Bonferroni correction; for males: [G];[G], 24.8 ± 4.94 (n = 1878); [G](;)[A], 24.2 ± 4.31 (n = 386); [A];[A], 22.4 ± 3.69 (n = 23); p = 4.20 × 10-4; beta: - 0.0573, CI - 0.0947 : - 0.0199. For average-height males the difference between [G];[G] and [A];[A] genotypes would correspond to ~ 6 kg, suggesting considerable protection against increased BMI. rs67047829 gives a pretermination codon in ERV3-1 which shares an exonic region and possibly promoter with ZNF117, previously associated with adiposity and type-2 diabetes. As this result occurs in a near-neutral Mendelian setting, a drug targetting ERV3-1/ZNF117 might potentially provide considerable benefits with minimal side-effects. This result needs to be replicated, followed by analyses of splice-variant mRNAs and protein expression.

Integrative analysis of risk factors for immune-related adverse events of checkpoint blockade therapy in cancer

Immune-related adverse events (irAEs) induced by checkpoint inhibitors involve a multitude of different risk factors. Here, to interrogate the multifaceted underlying mechanisms, we compiled germline exomes and blood transcriptomes with clinical data, before and after checkpoint inhibitor treatment, from 672 patients with cancer. Overall, irAE samples showed a substantially lower contribution of neutrophils in terms of baseline and on-therapy cell counts and gene expression markers related to neutrophil function. Allelic variation of HLA-B correlated with overall irAE risk. Analysis of germline coding variants identified a nonsense mutation in an immunoglobulin superfamily protein, TMEM162. In our cohort and the Cancer Genome Atlas (TCGA) data, TMEM162 alteration was associated with higher peripheral and tumor-infiltrating B cell counts and suppression of regulatory T cells in response to therapy. We developed machine learning models for irAE prediction, validated using additional data from 169 patients. Our results provide valuable insights into risk factors of irAE and their clinical utility.

Novel start codons introduce novel coding sequences in the human genomes

Start-gain mutations can introduce novel start codons and generate novel coding sequences that may affect the function of genes. In this study, we systematically investigated the novel start codons that were either polymorphic or fixed in the human genomes. 829 polymorphic start-gain SNVs were identified in the human populations, and the novel start codons introduced by these SNVs have significantly higher activity in translation initiation. Some of these start-gain SNVs were reported to be associated with phenotypes and diseases in previous studies. By comparative genomic analysis, we found 26 human-specific start codons that were fixed after the divergence between the human and chimpanzee, and high-level translation initiation activity was observed on them. The negative selection signal was detected in the novel coding sequences introduced by these human-specific start codons, indicating the important function of these novel coding sequences.

Discovery, classification, evolution and diversity of Siglecs

Immunoglobulin (Ig) superfamily proteins play diverse roles in vertebrates, including regulation of cellular responses by sensing endogenous or exogenous ligands. Siglecs are a family of glycan-recognizing proteins belonging to the Ig superfamily (i.e., I-type lectins). Siglecs are expressed on various leukocyte types and are involved in diverse aspects of immunity, including the regulation of inflammatory responses, leukocyte proliferation, host-microbe interaction, and cancer immunity. Sialoadhesin/Siglec-1, CD22/Siglec-2, and myelin-associated glycoprotein/Siglec-4 were among the first to be characterized as members of the Siglec family, and along with Siglec-15, they are relatively well-conserved among tetrapods. Conversely, CD33/Siglec-3-related Siglecs (CD33rSiglecs, so named as they show high sequence similarity with CD33/Siglec-3) are encoded in a gene cluster with many interspecies variations and even intraspecies variations within some lineages such as humans. The rapid evolution of CD33rSiglecs expressed on leukocytes involved in innate immunity likely reflects the selective pressure by pathogens that interact and possibly exploit these Siglecs. Human Siglecs have several additional unique and/or polymorphic properties as compared with closely related great apes, changes possibly related to the loss of the sialic acid Neu5Gc, another distinctly human event in sialobiology. Multiple changes in human CD33rSiglecs compared to great apes include many examples of human-specific expression in non-immune cells, coinciding with human-specific diseases involving such cell types. Some Siglec gene polymorphisms have dual consequences-beneficial in a situation but detrimental in another. The association of human Siglec gene polymorphisms with several infectious and non-infectious diseases likely reflects the ongoing competition between the host and microbial pathogens.

Adaptive Evolution of the OAS Gene Family Provides New Insights into the Antiviral Ability of Laurasiatherian Mammals

Many mammals risk damage from virus invasion due to frequent environmental changes. The oligoadenylate synthesis (OAS) gene family, which is an important component of the immune system, provides an essential response to the antiviral activities of interferons by regulating immune signal pathways. However, little is known about the evolutionary characteristics of OASs in Laurasiatherian mammals. Here, we examined the evolution of the OAS genes in 64 mammals to explore the accompanying molecular mechanisms of the antiviral ability of Laurasiatherian mammals living in different environments. We found that OAS2 and OAS3 were found to be pseudogenes in Odontoceti species. This may be related to the fact that they live in water. Some Antilopinae, Caprinae, and Cervidae species lacked the OASL gene, which may be related to their habitats being at higher altitudes. The OASs had a high number of positive selection sites in Cetartiodactyla, which drove the expression of strong antiviral ability. The OAS gene family evolved in Laurasiatherian mammals at different rates and was highly correlated with the species' antiviral ability. The gene evolution rate in Cetartiodactyla was significantly higher than that in the other orders. Compared to other species of the Carnivora family, the higher selection pressure on the OAS gene and the absence of positive selection sites in Canidae may be responsible for its weak resistance to rabies virus. The OAS gene family was relatively conserved during evolution. Conserved genes are able to provide better maintenance of gene function. The rate of gene evolution and the number of positively selected sites combine to influence the resistance of a species to viruses. The positive selection sites demonstrate the adaptive evolution of the OAS gene family to the environment. Adaptive evolution combined with conserved gene function improves resistance to viruses. Our findings offer insights into the molecular and functional evolution of the antiviral ability of Laurasian mammals.

Nonsense Mutations in Eukaryotes

Nonsense mutations are a type of mutations which results in a premature termination codon occurrence. In general, these mutations have been considered to be among the most harmful ones which lead to premature protein translation termination and result in shortened nonfunctional polypeptide. However, there is evidence that not all nonsense mutations are harmful as well as some molecular mechanisms exist which allow to avoid pathogenic effects of these mutations. This review addresses relevant information on nonsense mutations in eukaryotic genomes, characteristics of these mutations, and different molecular mechanisms preventing or mitigating harmful effects thereof.

Molecular profiling of individual FDA-approved clinical drugs identifies modulators of nonsense-mediated mRNA decay

Nonsense-mediated mRNA decay (NMD) degrades transcripts with premature stop codons. Given the prevalence of nonsense single nucleotide polymorphisms (SNPs) in the general population, it is urgent to catalog the effects of clinically approved drugs on NMD activity: any interference could alter the expression of nonsense SNPs, inadvertently inducing adverse effects. This risk is higher for patients with disease-causing nonsense mutations or an illness linked to dysregulated nonsense transcripts. On the other hand, hundreds of disorders are affected by cellular NMD efficiency and may benefit from NMD-modulatory drugs. Here, we profiled individual FDA-approved drugs for their impact on cellular NMD efficiency using a sensitive method that directly probes multiple endogenous NMD targets for a robust readout of NMD modulation. We found most FDA-approved drugs cause unremarkable effects on NMD, while many elicit clear transcriptional responses. Besides several potential mild NMD modulators, the anticancer drug homoharringtonine (HHT or omacetaxine mepesuccinate) consistently upregulates various endogenous NMD substrates in a dose-dependent manner in multiple cell types. We further showed translation inhibition mediates HHT's NMD effect. In summary, many FDA drugs induce transcriptional changes, and a few impact global NMD, and direct measurement of endogenous NMD substrate expression is robust to monitor cellular NMD.

Prioritising positively selected variants in whole-genome sequencing data using FineMAV

Background

In population genomics, polymorphisms that are highly differentiated between geographically separated populations are often suggestive of Darwinian positive selection. Genomic scans have highlighted several such regions in African and non-African populations, but only a handful of these have functional data that clearly associates candidate variations driving the selection process. Fine-Mapping of Adaptive Variation (FineMAV) was developed to address this in a high-throughput manner using population based whole-genome sequences generated by the 1000 Genomes Project. It pinpoints positively selected genetic variants in sequencing data by prioritizing high frequency, population-specific and functional derived alleles.

Results

We developed a stand-alone software that implements the FineMAV statistic. To graphically visualise the FineMAV scores, it outputs the statistics as bigWig files, which is a common file format supported by many genome browsers. It is available as a command-line and graphical user interface. The software was tested by replicating the FineMAV scores obtained using 1000 Genomes Project African, European, East and South Asian populations and subsequently applied to whole-genome sequencing datasets from Singapore and China to highlight population specific variants that can be subsequently modelled. The software tool is publicly available at https://github.com/fadilla-wahyudi/finemav.

Conclusions

The software tool described here determines genome-wide FineMAV scores, using low or high-coverage whole-genome sequencing datasets, that can be used to prioritize a list of population specific, highly differentiated candidate variants for in vitro or in vivo functional screens. The tool displays these scores on the human genome browsers for easy visualisation, annotation and comparison between different genomic regions in worldwide human populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04506-9.

A Positively Selected MAGEE2 LoF Allele Is Associated with Sexual Dimorphism in Human Brain Size and Shows Similar Phenotypes in Magee2 Null Mice

A nonsense allele at rs1343879 in human MAGEE2 on chromosome X has previously been reported as a strong candidate for positive selection in East Asia. This premature stop codon causing ∼80% protein truncation is characterized by a striking geographical pattern of high population differentiation: common in Asia and the Americas (up to 84% in the 1000 Genomes Project East Asians) but rare elsewhere. Here, we generated a Magee2 mouse knockout mimicking the human loss-of-function mutation to study its functional consequences. The Magee2 null mice did not exhibit gross abnormalities apart from enlarged brain structures (13% increased total brain area, P = 0.0022) in hemizygous males. The area of the granular retrosplenial cortex responsible for memory, navigation, and spatial information processing was the most severely affected, exhibiting an enlargement of 34% (P = 3.4×10-6). The brain size in homozygous females showed the opposite trend of reduced brain size, although this did not reach statistical significance. With these insights, we performed human association analyses between brain size measurements and rs1343879 genotypes in 141 Chinese volunteers with brain MRI scans, replicating the sexual dimorphism seen in the knockout mouse model. The derived stop gain allele was significantly associated with a larger volume of gray matter in males (P = 0.00094), and smaller volumes of gray (P = 0.00021) and white (P = 0.0015) matter in females. It is unclear whether or not the observed neuroanatomical phenotypes affect behavior or cognition, but it might have been the driving force underlying the positive selection in humans.

Genetic Variation and the Distribution of Variant Types in the Horse

Genetic variation is a key contributor to health and disease. Understanding the link between an individual's genotype and the corresponding phenotype is a major goal of medical genetics. Whole genome sequencing (WGS) within and across populations enables highly efficient variant discovery and elucidation of the molecular nature of virtually all genetic variation. Here, we report the largest catalog of genetic variation for the horse, a species of importance as a model for human athletic and performance related traits, using WGS of 534 horses. We show the extent of agreement between two commonly used variant callers. In data from ten target breeds that represent major breed clusters in the domestic horse, we demonstrate the distribution of variants, their allele frequencies across breeds, and identify variants that are unique to a single breed. We investigate variants with no homozygotes that may be potential embryonic lethal variants, as well as variants present in all individuals that likely represent regions of the genome with errors, poor annotation or where the reference genome carries a variant. Finally, we show regions of the genome that have higher or lower levels of genetic variation compared to the genome average. This catalog can be used for variant prioritization for important equine diseases and traits, and to provide key information about regions of the genome where the assembly and/or annotation need to be improved.

A map of copy number variations in the Tunisian population: a valuable tool for medical genomics in North Africa

Copy number variation (CNV) is considered as the most frequent type of structural variation in the human genome. Some CNVs can act on human phenotype diversity, encompassing rare Mendelian diseases and genomic disorders. The North African populations remain underrepresented in public genetic databases in terms of single-nucleotide variants as well as for larger genomic mutations. In this study, we present the first CNV map for a North African population using the Affymetrix Genome-Wide SNP (single-nucleotide polymorphism) array 6.0 array genotyping intensity data to call CNVs in 102 Tunisian healthy individuals. Two softwares, PennCNV and Birdsuite, were used to call CNVs in order to provide reliable data. Subsequent bioinformatic analyses were performed to explore their features and patterns. The CNV map of the Tunisian population includes 1083 CNVs spanning 61.443 Mb of the genome. The CNV length ranged from 1.017 kb to 2.074 Mb with an average of 56.734 kb. Deletions represent 57.43% of the identified CNVs, while duplications and the mixed loci are less represented. One hundred and three genes disrupted by CNVs are reported to cause 155 Mendelian diseases/phenotypes. Drug response genes were also reported to be affected by CNVs. Data on genes overlapped by deletions and duplications segments and the sequence properties in and around them also provided insights into the functional and health impacts of CNVs. These findings represent valuable clues to genetic diversity and personalized medicine in the Tunisian population as well as in the ethnically similar populations from North Africa.

Human-specific polymorphic pseudogenization of SIGLEC12 protects against advanced cancer progression

Compared with our closest living evolutionary cousins, humans appear unusually prone to develop carcinomas (cancers arising from epithelia). The SIGLEC12 gene, which encodes the Siglec-XII protein expressed on epithelial cells, has several uniquely human features: a fixed homozygous missense mutation inactivating its natural ligand recognition property; a polymorphic frameshift mutation eliminating full-length protein expression in ~60%-70% of worldwide human populations; and, genomic features suggesting a negative selective sweep favoring the pseudogene state. Despite the loss of canonical sialic acid binding, Siglec-XII still recruits Shp2 and accelerates tumor growth in a mouse model. We hypothesized that dysfunctional Siglec-XII facilitates human carcinoma progression, correlating with known tumorigenic signatures of Shp2-dependent cancers. Immunohistochemistry was used to detect Siglec-XII expression on tissue microarrays. PC-3 prostate cancer cells were transfected with Siglec-XII and transcription of genes enriched with Siglec-XII was determined. Genomic SIGLEC12 status was determined for four different cancer cohorts. Finally, a dot blot analysis of human urinary epithelial cells was established to determine the Siglec-XII expressors versus non-expressors. Forced expression in a SIGLEC12 null carcinoma cell line enriched transcription of genes associated with cancer progression. While Siglec-XII was detected as expected in ~30%-40% of normal epithelia, ~80% of advanced carcinomas showed strong expression. Notably, >80% of late-stage colorectal cancers had a functional SIGLEC12 allele, correlating with overall increased mortality. Thus, advanced carcinomas are much more likely to occur in individuals whose genomes have an intact SIGLEC12 gene, likely because the encoded Siglec-XII protein recruits Shp2-related oncogenic pathways. The finding has prognostic, diagnostic, and therapeutic implications.

Less Is More, Natural Loss-of-Function Mutation Is a Strategy for Adaptation

Gene gain and loss are crucial factors that shape the evolutionary success of diverse organisms. In the past two decades, more attention has been paid to the significance of gene gain through gene duplication or de novo genes. However, gene loss through natural loss-of-function (LoF) mutations, which is prevalent in the genomes of diverse organisms, has been largely ignored. With the development of sequencing techniques, many genomes have been sequenced across diverse species and can be used to study the evolutionary patterns of gene loss. In this review, we summarize recent advances in research on various aspects of LoF mutations, including their identification, evolutionary dynamics in natural populations, and functional effects. In particular, we discuss how LoF mutations can provide insights into the minimum gene set (or the essential gene set) of an organism. Furthermore, we emphasize their potential impact on adaptation. At the genome level, although most LoF mutations are neutral or deleterious, at least some of them are under positive selection and may contribute to biodiversity and adaptation. Overall, we highlight the importance of natural LoF mutations as a robust framework for understanding biological questions in general.

High-mobility group nucleosomal binding domain 2 protects against microcephaly by maintaining global chromatin accessibility during corticogenesis

The surface area of the human cerebral cortex undergoes dramatic expansion during late fetal development, leading to cortical folding, an evolutionary feature not present in rodents. Microcephaly is a neurodevelopmental disorder defined by an abnormally small brain, and many gene mutations have been found to be associated with primary microcephaly. However, mouse models generated by ablating primary microcephaly-associated genes often fail to recapitulate the severe loss of cortical surface area observed in individuals with this pathology. Here, we show that a mouse model with deficient expression of high-mobility group nucleosomal binding domain 2 (HMGN2) manifests microcephaly with reduced cortical surface area and almost normal radial corticogenesis, with a pattern of incomplete penetrance. We revealed that altered cleavage plane and mitotic delay of ventricular radial glia may explain the rising ratio of intermediate progenitor cells to radial glia and the displacement of neural progenitor cells in microcephalic mutant mice. These led to decreased self-renewal of the radial glia and reduction in lateral expansion. Furthermore, we found that HMGN2 protected corticogenesis by maintaining global chromatin accessibility mainly at promoter regions, thereby ensuring the correct regulation of the transcriptome. Our findings underscore the importance of the regulation of chromatin structure in cortical development and highlight a mouse model with critical insights into the etiology of microcephaly.

Revisiting the morbid genome of Mendelian disorders

BACKGROUND

The pathogenicity of many Mendelian variants has been challenged by large-scale sequencing efforts. However, many rare and benign "disease mutations" are difficult to analyze due to their rarity. The Saudi Arabian variome is enriched for homozygosity due to inbreeding, a key advantage that can be exploited for the critical examination of previously published variants.

RESULTS

We collated all "disease-related mutations" listed in the Human Gene Mutation Database (HGMD) and ClinVar, including "variants of uncertain significance" (VOUS). We find that the use of public databases including 1000 Genomes, ExAC, and Kaviar can reclassify many of these variants as likely benign. Our Saudi Human Genome Program (SHGP) can reclassify many variants that are rare in public databases. Furthermore, SGPD allows us to observe many previously reported variants in the homozygous state and our extensive phenotyping of participants makes it possible to demonstrate the lack of phenotype for these variants, thus challenging their pathogenicity despite their rarity. We also find that 18 VOUS BRCA1 and BRCA2 variants that are listed in BRCA Exchange are present at least once in the homozygous state in patients who lack features of Fanconi anemia. Reassuringly, we could reciprocally demonstrate that none of those labeled as "pathogenic" were observed in the homozygous statue in individuals who lack Fanconi phenotype in our database.

CONCLUSION

Our study shows the importance of revisiting disease-related databases using public resources as well as of population-specific resources to improve the specificity of the morbid genome of Mendelian diseases in humans.

RNA-Seq-Mediated Transcriptome Analysis of a Fiberless Mutant Cotton and Its Possible Origin Based on SNP Markers

As the longest known single-celled trichomes, cotton (Gossypium L.) fibers constitute a classic model system to investigate cell initiation and elongation. In this study, we used a high-throughput transcriptome sequencing technology to identify fiber-initiation-related single nucleotide polymorphism (SNP) markers and differentially expressed genes (DEGs) between the wild-type (WT) Upland cotton (G. hirsutum) Xuzhou 142 and its natural fuzzless-lintless mutant Xuzhou 142 fl. Approximately 700 million high-quality cDNA reads representing over 58 Gb of sequences were obtained, resulting in the identification of 28,610 SNPs--of which 17,479 were novel--from 13,960 expressed genes. Of these SNPs, 50% of SNPs in fl were identical to those of G. barbadense, which suggests the likely origin of the fl mutant from an interspecific hybridization between Xuzhou 142 and an unknown G. barbadense genotype. Of all detected SNPs, 15,555, 12,750, and 305 were classified as non-synonymous, synonymous, and pre-terminated ones, respectively. Moreover, 1,352 insertion/deletion polymorphisms (InDels) were also detected. A total of 865 DEGs were identified between the WT and fl in ovules at -3 and 0 days post-anthesis, with 302 candidate SNPs selected from these DEGs for validation by a high-resolution melting analysis and Sanger sequencing in seven cotton genotypes. The number of genotypic pairwise polymorphisms varied from 43 to 302, indicating that the identified SNPs are reliable. These SNPs should serve as good resources for breeding and genetic studies in cotton.

Premature termination codons in modern human genomes

The considerable range of genetic variation in human populations may partly reflect distinctive processes of adaptation to variable environmental conditions. However, the adaptive genomic signatures remain to be completely elucidated. This research explores candidate loci under selection at the population level by characterizing recently arisen premature termination codons (PTCs), some of which indicate a human knockout. From a total of 7595 participants from two population exome projects, 246 PTCs were found where natural selection has resulted in new alleles with a high frequency (from 1% to 96%) of derived alleles and various levels of population differentiation (F_ST = 0.00139–0.626). The PTC genes formed protein and regulatory networks limited to 15 biological processes or gene families, of which seven categories were previously unreported. PTC mutations have a strong tendency to be introduced into members of the same gene family, even during modern human evolution, although the exact nature of the selection is not fully known. The findings here suggest the ongoing evolutionary plasticity of modern humans at the genetic level and also partly provide insights into common human knockouts.

A catalog of hemizygous variation in 127 22q11 deletion patients

The 22q11.2 deletion syndrome is the most common microdeletion disorder, with wide phenotypic variability. To investigate variation within the non-deleted allele we performed targeted resequencing of the 22q11.2 region for 127 patients, identifying multiple deletion sizes, including two deletions with atypical breakpoints. We cataloged ~12,000 hemizygous variant positions, of which 84% were previously annotated. Within the coding regions 95 non-synonymous variants, three stop gains, and two frameshift insertions were identified, some of which we speculate could contribute to atypical phenotypes. We also catalog tolerability of 22q11 gene mutations based on related autosomal recessive disorders in man, embryonic lethality in mice, cross-species conservation and observations that some genes harbor more or less variants than expected. This extensive catalog of hemizygous variants will serve as a blueprint for future experiments to correlate 22q11DS variation with phenotype.

The distribution pattern of genetic variation in the transcript isoforms of the alternatively spliced protein-coding genes in the human genome

By enabling the transcription of multiple isoforms from the same gene locus, alternative-splicing mechanisms greatly expand the diversity of the human transcriptome and proteome. Currently, the alternatively spliced transcripts from each protein-coding gene locus in the human genome can be classified as either principal or non-principal isoforms, providing that they differ with respect to cross-species conservation or biological features. By mapping the variants from the 1000 Genomes Project onto the coding region of each isoform, an interesting pattern of the genetic variation distributions of the coding regions for these two types of transcript isoforms was revealed on a whole-genome scale: compared with the principal isoform-specific coding regions, the non-principal isoform-specific coding regions are significantly enriched in amino acid-changing variants, particularly those that have a strong impact on protein function and have higher derived allele frequencies, suggesting that non-principal isoform-specific substitutions are less likely to be related to phenotype changes or disease. The results herein can help us better understand the potential consequences of alternatively spliced products from a population perspective.

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.

Multiple loss-of-function variants of taste receptors in modern humans

Despite recent advances in the knowledge of interindividual taste differences, the underlying genetic backgrounds have remained to be fully elucidated. Much of the taste variation among different mammalian species can be explained by pseudogenization of taste receptors. Here I investigated whether the most recent disruptions of taste receptor genes segregate with their intact forms in modern humans by analyzing 14 ethnically diverse populations. The results revealed an unprecedented prevalence of 25 segregating loss-of-function (LoF) taste receptor variants, identifying one of the most pronounced cases of functional population diversity in the human genome. LoF variant frequency in taste receptors (2.10%) was considerably higher than the overall LoF frequency in human genome (0.16%). In particular, molecular evolutionary rates of candidate sour (14.7%) and bitter (1.8%) receptors were far higher in humans than those of sweet (0.02%), salty (0.05%), and umami (0.17%) receptors compared with other carnivorous mammals, although not all of the taste receptors were identified. Many LoF variants are population-specific, some of which arose even after population differentiation, not before divergence of the modern and archaic human. I conclude that modern humans might have been losing some sour and bitter receptor genes because of high-frequency LoF variants.

Selection Against Maternal microRNA Target Sites in Maternal Transcripts

In animals, before the zygotic genome is expressed, the egg already contains gene products deposited by the mother. These maternal products are crucial during the initial steps of development. In Drosophila melanogaster, a large number of maternal products are found in the oocyte, some of which are indispensable. Many of these products are RNA molecules, such as gene transcripts and ribosomal RNAs. Recently, microRNAs (small RNA gene regulators) have been detected early during development and are important in these initial steps. The presence of some microRNAs in unfertilized eggs has been reported, but whether they have a functional impact in the egg or early embryo has not being explored. I have extracted and sequenced small RNAs from Drosophila unfertilized eggs. The unfertilized egg is rich in small RNAs and contains multiple microRNA products. Maternal microRNAs often are encoded within the intron of maternal genes, suggesting that many maternal microRNAs are the product of transcriptional hitchhiking. Comparative genomics analyses suggest that maternal transcripts tend to avoid target sites for maternal microRNAs. I also developed a microRNA target mutation model to study the functional impact of polymorphisms at microRNA target sites. The analysis of Drosophila populations suggests that there is selection against maternal microRNA target sites in maternal transcripts. A potential role of the maternal microRNA mir-9c in maternal-to-zygotic transition is also discussed. In conclusion, maternal microRNAs in Drosophila have a functional impact in maternal protein-coding transcripts.

Genetic signatures of adaptation revealed from transcriptome sequencing of Arctic and red foxes

Background

The genus Vulpes (true foxes) comprises numerous species that inhabit a wide range of habitats and climatic conditions, including one species, the Arctic fox (Vulpes lagopus) which is adapted to the arctic region. A close relative to the Arctic fox, the red fox (Vulpes vulpes), occurs in subarctic to subtropical habitats. To study the genetic basis of their adaptations to different environments, transcriptome sequences from two Arctic foxes and one red fox individual were generated and analyzed for signatures of positive selection. In addition, the data allowed for a phylogenetic analysis and divergence time estimate between the two fox species.

Results

The de novo assembly of reads resulted in more than 160,000 contigs/transcripts per individual. Approximately 17,000 homologous genes were identified using human and the non-redundant databases. Positive selection analyses revealed several genes involved in various metabolic and molecular processes such as energy metabolism, cardiac gene regulation, apoptosis and blood coagulation to be under positive selection in foxes. Branch site tests identified four genes to be under positive selection in the Arctic fox transcriptome, two of which are fat metabolism genes. In the red fox transcriptome eight genes are under positive selection, including molecular process genes, notably genes involved in ATP metabolism. Analysis of the three transcriptomes and five Sanger re-sequenced genes in additional individuals identified a lower genetic variability within Arctic foxes compared to red foxes, which is consistent with distribution range differences and demographic responses to past climatic fluctuations. A phylogenomic analysis estimated that the Arctic and red fox lineages diverged about three million years ago.

Conclusions

Transcriptome data are an economic way to generate genomic resources for evolutionary studies. Despite not representing an entire genome, this transcriptome analysis identified numerous genes that are relevant to arctic adaptation in foxes. Similar to polar bears, fat metabolism seems to play a central role in adaptation of Arctic foxes to the cold climate, as has been identified in the polar bear, another arctic specialist.

Electronic supplementary material

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

Loss of gene function and evolution of human phenotypes

Humans have acquired many distinct evolutionary traits after the human-chimpanzee divergence. These phenotypes have resulted from genetic changes that occurred in the human genome and were retained by natural selection. Comparative primate genome analyses reveal that loss-of-function mutations are common in the human genome. Some of these gene inactivation events were revealed to be associated with the emergence of advantageous phenotypes and were therefore positively selected and fixed in modern humans (the "less-ismore" hypothesis). Representative cases of human gene inactivation and their functional implications are presented in this review. Functional studies of additional inactive genes will provide insight into the molecular mechanisms underlying acquisition of various human-specific traits.

Human knockout research: new horizons and opportunities

Although numerous approaches have been pursued to understand the function of human genes, Mendelian genetics has by far provided the most compelling and medically actionable dataset. Biallelic loss-of-function (LOF) mutations are observed in the majority of autosomal recessive Mendelian disorders, representing natural human knockouts and offering a unique opportunity to study the physiological and developmental context of these genes. The restriction of such context to 'disease' states is artificial, however, and the recent ability to survey entire human genomes for biallelic LOF mutations has revealed a surprising landscape of knockout events in 'healthy' individuals, sparking interest in their role in phenotypic diversity beyond disease causation. As I discuss in this review, the potentially wide implications of human knockout research warrant increased investment and multidisciplinary collaborations to overcome existing challenges and reap its benefits.

High-throughput screening of mouse gene knockouts identifies established and novel skeletal phenotypes

Screening gene function in vivo is a powerful approach to discover novel drug targets. We present high-throughput screening (HTS) data for 3 762 distinct global gene knockout (KO) mouse lines with viable adult homozygous mice generated using either gene-trap or homologous recombination technologies. Bone mass was determined from DEXA scans of male and female mice at 14 weeks of age and by microCT analyses of bones from male mice at 16 weeks of age. Wild-type (WT) cagemates/littermates were examined for each gene KO. Lethality was observed in an additional 850 KO lines. Since primary HTS are susceptible to false positive findings, additional cohorts of mice from KO lines with intriguing HTS bone data were examined. Aging, ovariectomy, histomorphometry and bone strength studies were performed and possible non-skeletal phenotypes were explored. Together, these screens identified multiple genes affecting bone mass: 23 previously reported genes (Calcr, Cebpb, Crtap, Dcstamp, Dkk1, Duoxa2, Enpp1, Fgf23, Kiss1/Kiss1r, Kl (Klotho), Lrp5, Mstn, Neo1, Npr2, Ostm1, Postn, Sfrp4, Slc30a5, Slc39a13, Sost, Sumf1, Src, Wnt10b), five novel genes extensively characterized (Cldn18, Fam20c, Lrrk1, Sgpl1, Wnt16), five novel genes with preliminary characterization (Agpat2, Rassf5, Slc10a7, Slc26a7, Slc30a10) and three novel undisclosed genes coding for potential osteoporosis drug targets.

Selection on a variant associated with improved viral clearance drives local, adaptive pseudogenization of interferon lambda 4 (IFNL4)

Interferon lambda 4 gene (IFNL4) encodes IFN-λ4, a new member of the IFN-λ family with antiviral activity. In humans IFNL4 open reading frame is truncated by a polymorphic frame-shift insertion that eliminates IFN-λ4 and turns IFNL4 into a polymorphic pseudogene. Functional IFN-λ4 has antiviral activity but the elimination of IFN-λ4 through pseudogenization is strongly associated with improved clearance of hepatitis C virus (HCV) infection. We show that functional IFN-λ4 is conserved and evolutionarily constrained in mammals and thus functionally relevant. However, the pseudogene has reached moderately high frequency in Africa, America, and Europe, and near fixation in East Asia. In fact, the pseudogenizing variant is among the 0.8% most differentiated SNPs between Africa and East Asia genome-wide. Its raise in frequency is associated with additional evidence of positive selection, which is strongest in East Asia, where this variant falls in the 0.5% tail of SNPs with strongest signatures of recent positive selection genome-wide. Using a new Approximate Bayesian Computation (ABC) approach we infer that the pseudogenizing allele appeared just before the out-of-Africa migration and was immediately targeted by moderate positive selection; selection subsequently strengthened in European and Asian populations resulting in the high frequency observed today. This provides evidence for a changing adaptive process that, by favoring IFN-λ4 inactivation, has shaped present-day phenotypic diversity and susceptibility to disease.

The genetic association with clearance of Hepatitis C virus (HCV) is one of the strongest and most elusive known associations with disease. The genetic variant more strongly associated with improved HCV clearance inactivates the recently discovered IFNL4 gene, which encodes for antiviral IFN-λ4 protein, and turns it into a polymorphic pseudogene. We show that functional IFN-λ4 is conserved and functionally important in mammals. In humans though the inactivating mutation appeared in Africa just before the out-of-Africa migration and quickly became advantageous, with the strength of selection (the degree of advantage) varying across human groups. In particular, selection became stronger out of Africa and was strongest in East Asia, raising the frequency of the pseudogene and resulting in the virtual loss of functional IFN-λ4 protein in several Asian populations. Although the environmental force driving selection is unknown, this process resulted in variable clearance of HCV in modern human populations. The complex selective history of IFNL4-inactivating allele has thus shaped present-day heterogeneity across populations not only in genetic variation, but also in relevant phenotypes and susceptibility to disease.

Analysis of stop-gain and frameshift variants in human innate immunity genes

Loss-of-function variants in innate immunity genes are associated with Mendelian disorders in the form of primary immunodeficiencies. Recent resequencing projects report that stop-gains and frameshifts are collectively prevalent in humans and could be responsible for some of the inter-individual variability in innate immune response. Current computational approaches evaluating loss-of-function in genes carrying these variants rely on gene-level characteristics such as evolutionary conservation and functional redundancy across the genome. However, innate immunity genes represent a particular case because they are more likely to be under positive selection and duplicated. To create a ranking of severity that would be applicable to innate immunity genes we evaluated 17,764 stop-gain and 13,915 frameshift variants from the NHLBI Exome Sequencing Project and 1,000 Genomes Project. Sequence-based features such as loss of functional domains, isoform-specific truncation and nonsense-mediated decay were found to correlate with variant allele frequency and validated with gene expression data. We integrated these features in a Bayesian classification scheme and benchmarked its use in predicting pathogenic variants against Online Mendelian Inheritance in Man (OMIM) disease stop-gains and frameshifts. The classification scheme was applied in the assessment of 335 stop-gains and 236 frameshifts affecting 227 interferon-stimulated genes. The sequence-based score ranks variants in innate immunity genes according to their potential to cause disease, and complements existing gene-based pathogenicity scores. Specifically, the sequence-based score improves measurement of functional gene impairment, discriminates across different variants in a given gene and appears particularly useful for analysis of less conserved genes.

There are well-characterized severe immunodeficiencies associated with loss-of-function variants in innate immunity genes. Genome sequencing projects identify rare stop-gain and frameshift variants in innate immunity genes whose phenotype is uncharacterized. Current methods to estimate the severity of rare stop-gains and frameshifts are based on evolutionary conservation of the gene, the likelihood for redundancy in its function or mutational burden. These parameters are not always applicable to innate immunity genes. We evaluated sequence-level characteristics of more than 30'000 stop-gains and frameshifts and prioritized variants according to their predicted functional consequences. Our scoring approach complements existing tools in the prediction of innate immunity OMIM disease variants and associates with functional readouts such as gene expression. In this framework, we show that many individuals do carry highly pathogenic variants in genes participating in antiviral defense. The clinical assessment of these variants is of significant interest.

Effect of duplicate genes on mouse genetic robustness: an update

In contrast to S. cerevisiae and C. elegans, analyses based on the current knockout (KO) mouse phenotypes led to the conclusion that duplicate genes had almost no role in mouse genetic robustness. It has been suggested that the bias of mouse KO database toward ancient duplicates may possibly cause this knockout duplicate puzzle, that is, a very similar proportion of essential genes (P_E) between duplicate genes and singletons. In this paper, we conducted an extensive and careful analysis for the mouse KO phenotype data and corroborated a strong effect of duplicate genes on mouse genetics robustness. Moreover, the effect of duplicate genes on mouse genetic robustness is duplication-age dependent, which holds after ruling out the potential confounding effect from coding-sequence conservation, protein-protein connectivity, functional bias, or the bias of duplicates generated by whole genome duplication (WGD). Our findings suggest that two factors, the sampling bias toward ancient duplicates and very ancient duplicates with a proportion of essential genes higher than that of singletons, have caused the mouse knockout duplicate puzzle; meanwhile, the effect of genetic buffering may be correlated with sequence conservation as well as protein-protein interactivity.

Associations of genetic polymorphisms of Siglecs with human diseases

Genetic polymorphism studies in humans provide unique opportunities to understand human biology and the mechanisms of diseases. Correlations between polymorphisms in the genes encoding human Siglecs and various diseases have been reported. Leading examples, such as the CD33 polymorphism associated with late-onset Alzheimer's disease, are well supported by genetic replication and mechanistic studies, while some others (such as SIGLEC8 polymorphism associated with bronchial asthma and SIGLEC14 polymorphism associated with exacerbation of chronic obstructive pulmonary disease) may benefit reinforcement by independent genetic replication or mechanistic studies. In a few cases, such as MAG polymorphism associated with psychological disorder and CD22 polymorphism associated with autoimmune disease, the phenotype associated with a genetic polymorphism of a Siglec gene and that of an enzyme gene involved in the biosynthesis of Siglec ligand show some overlap, providing indirect support for the observed genotype-phenotype association. Although studies using engineered mutant mice have provided invaluable insights into the biological functions and mechanisms of diseases, it is not always possible to develop appropriate mouse model to replicate human situations because of significant species-to-species differences, which can be a major obstacle in understanding the biology of some of human CD33/Siglec-3-related Siglecs. Further studies in genetic polymorphisms of human Siglecs, combined with appropriate functional studies, may reveal unexpected biological roles of human Siglecs, and identify possible targets for prevention and/or treatment of certain diseases.

Evaluation of MAGE-D4 expression in hepatocellular carcinoma in Japanese patients

BACKGROUND AND OBJECTIVES

Though Melanoma-associated antigen (MAGE) family genes have received lots of attention as cancer-related genes and targets for immunotherapy, MAGE-D4 expression in hepatocellular carcinoma (HCC) has not yet been evaluated.

METHODS

MAGE-D4 mRNA expression was assayed in nine HCC cell lines and 94 HCC surgical specimens obtained from Japanese patients by quantitative real-time reverse transcription polymerase chain reaction, and the correlations between MAGE-D4 mRNA expression and clinicopathological factors were evaluated. The expression and distribution of MAGE-D4b protein were evaluated immunohistochemically.

RESULTS

MAGE-D4 mRNA was overexpressed in five of nine HCC cell lines and 34 of 94 primary HCCs (36.2%). Median overall survival (14.8 vs. 118 months, P < 0.001) and relapse-free survival (2.7 vs. 18.3 months, P < 0.001) were significantly shorter in patients with high than with low-moderate MAGE-D4 expression. Multivariate analysis for overall survival showed that MAGE-D4 overexpression was independently prognostic for survival (hazard ratio 2.88, P = 0.009) and significantly associated with high alpha-fetoprotein concentration (P < 0.001), poor tumor differentiation (P = 0.003) and vascular invasion (P = 0.021). MAGE-D4b protein expression patterns were consistent with those of MAGE-D4 mRNA.

CONCLUSIONS

Overexpression of MAGE-D4 may be a predictive marker of early recurrence and mortality in patients with HCC.

Genome-wide discovery and information resource development of DNA polymorphisms in cassava

Cassava (Manihot esculenta Crantz) is an important crop that provides food security and income generation in many tropical countries, and is known for its adaptability to various environmental conditions. Its draft genome sequence and many expressed sequence tags are now publicly available, allowing the development of cassava polymorphism information. Here, we describe the genome-wide discovery of cassava DNA polymorphisms. Using the alignment of predicted transcribed sequences from the cassava draft genome sequence and ESTs from GenBank, we discovered 10,546 single-nucleotide polymorphisms and 647 insertions and deletions. To facilitate molecular marker development for cassava, we designed 9,316 PCR primer pairs to amplify the genomic region around each DNA polymorphism. Of the discovered SNPs, 62.7% occurred in protein-coding regions. Disease-resistance genes were found to have a significantly higher ratio of nonsynonymous-to-synonymous substitutions. We identified 24 read-through (changes of a stop codon to a coding codon) and 38 premature stop (changes of a coding codon to a stop codon) single-nucleotide polymorphisms, and found that the 5 gene ontology terms in biological process were significantly different in genes with read-through single-nucleotide polymorphisms compared with all cassava genes. All data on the discovered DNA polymorphisms were organized into the Cassava Online Archive database, which is available at http://cassava.psc.riken.jp/.

Deleterious- and disease-allele prevalence in healthy individuals: insights from current predictions, mutation databases, and population-scale resequencing

We have assessed the numbers of potentially deleterious variants in the genomes of apparently healthy humans by using (1) low-coverage whole-genome sequence data from 179 individuals in the 1000 Genomes Pilot Project and (2) current predictions and databases of deleterious variants. Each individual carried 281-515 missense substitutions, 40-85 of which were homozygous, predicted to be highly damaging. They also carried 40-110 variants classified by the Human Gene Mutation Database (HGMD) as disease-causing mutations (DMs), 3-24 variants in the homozygous state, and many polymorphisms putatively associated with disease. Whereas many of these DMs are likely to represent disease-allele-annotation errors, between 0 and 8 DMs (0-1 homozygous) per individual are predicted to be highly damaging, and some of them provide information of medical relevance. These analyses emphasize the need for improved annotation of disease alleles both in mutation databases and in the primary literature; some HGMD mutation data have been recategorized on the basis of the present findings, an iterative process that is both necessary and ongoing. Our estimates of deleterious-allele numbers are likely to be subject to both overcounting and undercounting. However, our current best mean estimates of ~400 damaging variants and ~2 bona fide disease mutations per individual are likely to increase rather than decrease as sequencing studies ascertain rare variants more effectively and as additional disease alleles are discovered.

Effects of premature termination codon polymorphisms in the Drosophila pseudoobscura subclade

Premature termination codon (PTC) mutations can have dramatic effects--both adaptive and deleterious--on gene expression and function. Here, we examine the number and selective effects of PTC mutations within the Drosophila pseudoobscura subclade using 18 resequenced genomes aligned to the reference genome. We located and characterized 1,679 PTC mutations in 605 genes across each of these genomes relative to the D. pseudoobscura reference genome, and use RT-PCR to confirm transcription of a subset of these genes containing PTC mutations. We confirm previous findings that genes containing PTC mutations are less selectively constrained and less broadly expressed than non-PTC-containing genes, suggesting that the most of these mutations are at least mildly deleterious. Further, we find highly significant codon usage bias in regions downstream of the PTC in 38 of these PTC-containing genes, suggesting that some of these PTC mutations--if not alternatively spliced out of the transcript--have neutral effects. Ultimately, these analyzes support the view that the PTC mutations are mostly detrimental, but are nonetheless common enough in genomes that a subset could be effectively neutral.

Sequencing and analysis of a South Asian-Indian personal genome

BACKGROUND

With over 1.3 billion people, India is estimated to contain three times more genetic diversity than does Europe. Next-generation sequencing technologies have facilitated the understanding of diversity by enabling whole genome sequencing at greater speed and lower cost. While genomes from people of European and Asian descent have been sequenced, only recently has a single male genome from the Indian subcontinent been published at sufficient depth and coverage. In this study we have sequenced and analyzed the genome of a South Asian Indian female (SAIF) from the Indian state of Kerala.

RESULTS

We identified over 3.4 million SNPs in this genome including over 89,873 private variations. Comparison of the SAIF genome with several published personal genomes revealed that this individual shared ~50% of the SNPs with each of these genomes. Analysis of the SAIF mitochondrial genome showed that it was closely related to the U1 haplogroup which has been previously observed in Kerala. We assessed the SAIF genome for SNPs with health and disease consequences and found that the individual was at a higher risk for multiple sclerosis and a few other diseases. In analyzing SNPs that modulate drug response, we found a variation that predicts a favorable response to metformin, a drug used to treat diabetes. SNPs predictive of adverse reaction to warfarin indicated that the SAIF individual is not at risk for bleeding if treated with typical doses of warfarin. In addition, we report the presence of several additional SNPs of medical relevance.

CONCLUSIONS

This is the first study to report the complete whole genome sequence of a female from the state of Kerala in India. The availability of this complete genome and variants will further aid studies aimed at understanding genetic diversity, identifying clinically relevant changes and assessing disease burden in the Indian population.

Personal receptor repertoires: olfaction as a model

Background

Information on nucleotide diversity along completely sequenced human genomes has increased tremendously over the last few years. This makes it possible to reassess the diversity status of distinct receptor proteins in different human individuals. To this end, we focused on the complete inventory of human olfactory receptor coding regions as a model for personal receptor repertoires.

Results

By performing data-mining from public and private sources we scored genetic variations in 413 intact OR loci, for which one or more individuals had an intact open reading frame. Using 1000 Genomes Project haplotypes, we identified a total of 4069 full-length polypeptide variants encoded by these OR loci, average of ~10 per locus, constituting a lower limit for the effective human OR repertoire. Each individual is found to harbor as many as 600 OR allelic variants, ~50% higher than the locus count. Because OR neuronal expression is allelically excluded, this has direct effect on smell perception diversity of the species. We further identified 244 OR segregating pseudogenes (SPGs), loci showing both intact and pseudogene forms in the population, twenty-six of which are annotatively “resurrected” from a pseudogene status in the reference genome. Using a custom SNP microarray we validated 150 SPGs in a cohort of 468 individuals, with every individual genome averaging 36 disrupted sequence variations, 15 in homozygote form. Finally, we generated a multi-source compendium of 63 OR loci harboring deletion Copy Number Variations (CNVs). Our combined data suggest that 271 of the 413 intact OR loci (66%) are affected by nonfunctional SNPs/indels and/or CNVs.

Conclusions

These results portray a case of unusually high genetic diversity, and suggest that individual humans have a highly personalized inventory of functional olfactory receptors, a conclusion that might apply to other receptor multigene families.

Presence of the functional CASPASE-12 allele in Indian subpopulations

Most humans lack a functional CASP12 gene, with the nonfunctional variant (CASP12p1), found in 100% of the Caucasian and east Asian population, and in approximately 80% of people of African descent. However, 20% of Africans carry an intact allele of CASP12, which produces a full-length pro-enzyme and increases the risk of sepsis. We examined CASP12 allele distribution in persons from central and southern Asia and found that CASP12 was significantly present in members of the Dravidian language group, particularly in persons from the Indian state of Tamil Nadu.

A systematic survey of loss-of-function variants in human protein-coding genes

Genome-sequencing studies indicate that all humans carry many genetic variants predicted to cause loss of function (LoF) of protein-coding genes, suggesting unexpected redundancy in the human genome. Here we apply stringent filters to 2951 putative LoF variants obtained from 185 human genomes to determine their true prevalence and properties. We estimate that human genomes typically contain ~100 genuine LoF variants with ~20 genes completely inactivated. We identify rare and likely deleterious LoF alleles, including 26 known and 21 predicted severe disease-causing variants, as well as common LoF variants in nonessential genes. We describe functional and evolutionary differences between LoF-tolerant and recessive disease genes and a method for using these differences to prioritize candidate genes found in clinical sequencing studies.

Widespread polymorphism in the positions of stop codons in Drosophila melanogaster

The mechanisms underlying evolutionary changes in protein length are poorly understood. Protein domains are lost and gained between species and must have arisen first as within-species polymorphisms. Here, we use Drosophila melanogaster population genomic data combined with between species divergence information to understand the evolutionary forces that generate and maintain polymorphisms causing changes in protein length in D. melanogaster. Specifically, we looked for protein length variations resulting from premature termination codons (PTCs) and stop codon losses (SCLs). We discovered that 438 genes contained polymorphisms resulting in truncation of the translated region (PTCs) and 119 genes contained polymorphisms predicted to lengthen the translated region (SCLs). Stop codon polymorphisms (SCPs) (especially PTCs) appear to be more deleterious than other polymorphisms, including protein amino acid changes. Genes harboring SCPs are in general less selectively constrained, more narrowly expressed, and enriched for dispensable biological functions. However, we also observed exceptional cases such as genes that have multiple independent SCPs, alleles that are shared between D. melanogaster and Drosophila simulans, and high-frequency alleles that cause extreme changes in gene length. SCPs likely have an important role in the evolution of these genes.

SIGLEC12, a human-specific segregating (pseudo)gene, encodes a signaling molecule expressed in prostate carcinomas

The primate SIGLEC12 gene encodes one of the CD33-related Siglec family of signaling molecules in immune cells. We had previously reported that this gene harbors a human-specific missense mutation of the codon for an Arg residue required for sialic acid recognition. Here we show that this R122C mutation of the Siglec-XII protein is fixed in the human population, i.e. it occurred prior to the origin of modern humans. Additional mutations have since completely inactivated the SIGLEC12 gene in some but not all humans. The most common inactivating mutation with a global allele frequency of 58% is a single nucleotide frameshift that markedly shortens the open reading frame. Unlike other CD33-related Siglecs that are primarily found on immune cells, we found that Siglec-XII protein is expressed not only on some macrophages but also on various epithelial cell surfaces in humans and chimpanzees. We also found expression on certain human prostate epithelial carcinomas and carcinoma cell lines. This expression correlates with the presence of the nonframeshifted, intact SIGLEC12 allele. Although SIGLEC12 allele status did not predict prostate carcinoma incidence, restoration of expression in a prostate carcinoma cell line homozygous for the frameshift mutation induced altered regulation of several genes associated with carcinoma progression. These stably transfected Siglec-XII-expressing prostate cancer cells also showed enhanced growth in nude mice. Finally, monoclonal antibodies against the protein were internalized by Siglec-XII-expressing prostate carcinoma cells, allowing targeting of a toxin to such cells. Polymorphic expression of Siglec-XII in humans thus has implications for prostate cancer biology and therapeutics.

The impact of caspase-12 on susceptibility to candidemia

Candida is one of the leading causes of sepsis, and an effective host immune response to Candida critically depends on the cytokines IL-1β and IL-18, which need caspase-1 cleavage to become bioactive. Caspase-12 has been suggested to inhibit caspase-1 activation and has been implicated as a susceptibility factor for bacterial sepsis. In populations of African descent, CASPASE-12 is either functional or non-functional. Here, we have assessed the frequencies of both CASPASE-12 alleles in an African-American Candida sepsis patients cohort compared to uninfected patients with similar predisposing factors. African-American Candida sepsis patients (n = 93) and non-infected African-American patients (n = 88) were genotyped for the CASPASE-12 genotype. Serum cytokine concentrations of IL-6, IL-8, and IFNγ were measured in the serum of infected patients. Statistical comparisons were performed in order to assess the effect of the CASPASE-12 genotype on susceptibility to candidemia and on serum cytokine concentrations. Our findings demonstrate that CASPASE-12 does not influence the susceptibility to Candida sepsis, nor has any effect on the serum cytokine concentrations in Candida sepsis patients during the course of infection. Although the functional CASPASE-12 allele has been suggested to increase susceptibility to bacterial sepsis, this could not be confirmed in our larger cohort of fungal sepsis patients.

Mutations of GIPC3 cause nonsyndromic hearing loss DFNB72 but not DFNB81 that also maps to chromosome 19p

A missense mutation of Gipc3 was previously reported to cause age-related hearing loss in mice. Point mutations of human GIPC3 were found in two small families, but association with hearing loss was not statistically significant. Here, we describe one frameshift and six missense mutations in GIPC3 cosegregating with DFNB72 hearing loss in six large families that support statistically significant evidence for genetic linkage. However, GIPC3 is not the only nonsyndromic hearing impairment gene in this region; no GIPC3 mutations were found in a family cosegregating hearing loss with markers of chromosome 19p. Haplotype analysis excluded GIPC3 from the obligate linkage interval in this family and defined a novel locus spanning 4.08 Mb and 104 genes. This closely linked but distinct nonsyndromic hearing loss locus was designated DFNB81.

Human TLRs and IL-1Rs in host defense: natural insights from evolutionary, epidemiological, and clinical genetics

Toll-like receptors (TLRs) and interleukin-1 receptors (IL-1Rs) have TIR intracellular domains that engage two main signaling pathways, via the TIR-containing adaptors MyD88 (which is not used by TLR3) and TRIF (which is used only by TLR3 and TLR4). Extensive studies in inbred mice in various experimental settings have attributed key roles in immunity to TLR- and IL-1R-mediated responses, but what contribution do human TLRs and IL-1Rs actually make to host defense in the natural setting? Evolutionary genetic studies have shown that human intracellular TLRs have evolved under stronger purifying selection than surface-expressed TLRs, for which the frequency of missense and nonsense alleles is high in the general population. Epidemiological genetic studies have yet to provide convincing evidence of a major contribution of common variants of human TLRs, IL-1Rs, or their adaptors to host defense. Clinical genetic studies have revealed that rare mutations affecting the TLR3-TRIF pathway underlie herpes simplex virus encephalitis, whereas mutations in the TIR-MyD88 pathway underlie pyogenic bacterial diseases in childhood. A careful reconsideration of the contributions of TLRs and IL-1Rs to host defense in natura is required.

The role played by natural selection on Mendelian traits in humans

Detecting whether and how natural selection has targeted regions of the human genome represents a complementary strategy for identifying functionally important loci and variants involved in disease resistance and adaptation to the environment. In contrast with most complex diseases or traits, the genetic architecture of most Mendelian traits is relatively well established. Most mutations associated with Mendelian disease-related traits are highly penetrant and kept at low population frequencies because of the effects of purifying selection. However, this is not always the case. Here, we review several examples of Mendelian mutations-associated with various disease conditions or other traits of anthropological interest-that have increased in frequency in the human population as a result of past positive selection. These examples clearly illustrate the value of a population genetics approach to unravel the biological mechanisms that have been central to our past and present survival against the selective pressures imposed by diseases and other environmental factors.