Search

Simulation is a key tool in population genetics for both methods development and empirical research, but producing simulations that recapitulate the main features of genomic datasets remains a major obstacle. Today, more realistic simulations are possible thanks to large increases in the quantity and quality of available genetic data, and the sophistication of inference and simulation software. However, implementing these simulations still requires substantial time and specialized knowledge. These challenges are especially pronounced for simulating genomes for species that are not well-studied, since it is not always clear what information is required to produce simulations with a level of realism sufficient to confidently answer a given question. The community-developed framework stdpopsim seeks to lower this barrier by facilitating the simulation of complex population genetic models using up-to-date information. The initial version of stdpopsim focused on establishing this framework using six well-characterized model species (Adrion et al., 2020). Here, we report on major improvements made in the new release of stdpopsim (version 0.2), which includes a significant expansion of the species catalog and substantial additions to simulation capabilities. Features added to improve the realism of the simulated genomes include non-crossover recombination and provision of species-specific genomic annotations. Through community-driven efforts, we expanded the number of species in the catalog more than threefold and broadened coverage across the tree of life. During the process of expanding the catalog, we have identified common sticking points and developed the best practices for setting up genome-scale simulations. We describe the input data required for generating a realistic simulation, suggest good practices for obtaining the relevant information from the literature, and discuss common pitfalls and major considerations. These improvements to stdpopsim aim to further promote the use of realistic whole-genome population genetic simulations, especially in non-model organisms, making them available, transparent, and accessible to everyone.

Mutations in interferon regulatory factor 6 (IRF6) account for ∼70% of cases of Van der Woude syndrome (VWS), the most common syndromic form of cleft lip and palate. In 8 of 45 VWS-affected families lacking a mutation in IRF6, we found coding mutations in grainyhead-like 3 (GRHL3). According to a zebrafish-based assay, the disease-associated GRHL3 mutations abrogated periderm development and were consistent with a dominant-negative effect, in contrast to haploinsufficiency seen in most VWS cases caused by IRF6 mutations. In mouse, all embryos lacking Grhl3 exhibited abnormal oral periderm and 17% developed a cleft palate. Analysis of the oral phenotype of double heterozygote (Irf6(+/-);Grhl3(+/-)) murine embryos failed to detect epistasis between the two genes, suggesting that they function in separate but convergent pathways during palatogenesis. Taken together, our data demonstrated that mutations in two genes, IRF6 and GRHL3, can lead to nearly identical phenotypes of orofacial cleft. They supported the hypotheses that both genes are essential for the presence of a functional oral periderm and that failure of this process contributes to VWS.

Sex chromosomes have evolved repeatedly across the tree of life and often exhibit extreme size dimorphism due to genetic degeneration of the sex-limited chromosome (e.g. the W chromosome of some birds and Y chromosome of mammals). However, in some lineages, ancient sex-limited chromosomes have escaped degeneration. Here, we study the evolutionary maintenance of sex chromosomes in the ostrich (Struthio camelus), where the W remains 65% the size of the Z chromosome, despite being more than 100 million years old. Using genome-wide resequencing data, we show that the population scaled recombination rate of the pseudoautosomal region (PAR) is higher than similar sized autosomes and is correlated with pedigree-based recombination rate in the heterogametic females, but not homogametic males. Genetic variation within the sex-linked region (SLR) (π = 0.001) was significantly lower than in the PAR, consistent with recombination cessation. Conversely, genetic variation across the PAR (π = 0.0016) was similar to that of autosomes and dependent on local recombination rates, GC content and to a lesser extent, gene density. In particular, the region close to the SLR was as genetically diverse as autosomes, likely due to high recombination rates around the PAR boundary restricting genetic linkage with the SLR to only ~50Kb. The potential for alleles with antagonistic fitness effects in males and females to drive chromosome degeneration is therefore limited. While some regions of the PAR had divergent male-female allele frequencies, suggestive of sexually antagonistic alleles, coalescent simulations showed this was broadly consistent with neutral genetic processes. Our results indicate that the degeneration of the large and ancient sex chromosomes of the ostrich may have been slowed by high recombination in the female PAR, reducing the scope for the accumulation of sexually antagonistic variation to generate selection for recombination cessation., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Yazdi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Simulation is a key tool in population genetics for both methods development and empirical research, but producing simulations that recapitulate the main features of genomic datasets remains a major obstacle. Today, more realistic simulations are possible thanks to large increases in the quantity and quality of available genetic data, and the sophistication of inference and simulation software. However, implementing these simulations still requires substantial time and specialized knowledge. These challenges are especially pronounced for simulating genomes for species that are not well-studied, since it is not always clear what information is required to produce simulations with a level of realism sufficient to confidently answer a given question. The community-developed framework stdpopsim seeks to lower this barrier by facilitating the simulation of complex population genetic models using up-to-date information. The initial version of stdpopsim focused on establishing this framework using six well-characterized model species (Adrion et al., 2020). Here, we report on major improvements made in the new release of stdpopsim (version 0.2), which includes a significant expansion of the species catalog and substantial additions to simulation capabilities. Features added to improve the realism of the simulated genomes include non-crossover recombination and provision of species-specific genomic annotations. Through community-driven efforts, we expanded the number of species in the catalog more than threefold and broadened coverage across the tree of life. During the process of expanding the catalog, we have identified common sticking points and developed the best practices for setting up genome-scale simulations. We describe the input data required for generating a realistic simulation, suggest good practices for obtaining the relevant information from the literature, and discuss common pitfalls and major considerations. These improvements to stdpopsim aim to further promote the use of realistic whole-genome population genetic simulations, especially in non-model organisms, making them available, transparent, and accessible to everyone., Competing Interests: ML, MC, GG, NP, GT, SB, VC, JC, IE, BH, AH, XH, LI, EN, JO, VP, AP, DP, MP, MR, CS, JS, AT, ST, PU, JV, RW, AW, YW, FB, RC, GG, RG, JK, AK, AR, PR, DS, IG No competing interests declared, AG is an employee of Embark Veterinary, Inc. The author declares that no other competing interests exist, JA is an employee of Ancestry DNA. The author declares that no other competing interests exist, AB is an employee of 54Gene, Inc. The author declares that no other competing interests exist, (© 2023, Lauterbur et al.)

Background: Germline genetic variants are an important cause of dilated cardiomyopathy (DCM). However, recent sequencing studies have revealed rare variants in DCM-associated genes also in individuals without known heart disease. In this study, we investigate variant prevalence and genotype-phenotype correlations in Swedish DCM patients, and compare their genetic variants to those detected in reference cohorts., Methods and Results: We sequenced the coding regions of 41 DCM-associated genes in 176 unrelated patients with idiopathic DCM and found 102 protein-altering variants with an allele frequency of <0.04% in reference cohorts; the majority were missense variants not previously described in DCM. Fifty-five (31%) patients had one variant, and 24 (14%) patients had two or more variants in the analysed genes. Detection of genetic variants in any gene, and in LMNA, MYH7 or TTN alone, was associated with early onset disease and reduced transplant-free survival. As expected, nonsense and frameshift variants were more common in DCM patients than in healthy individuals of the reference cohort 1000 Genomes Europeans. Surprisingly however, the prevalence, conservation and pathogenicity scores, and localization of missense variants were similar in DCM patients and healthy reference individuals., Conclusion: To our knowledge, this is the first study to identify correlations between genotype and prognosis when sequencing a large number of genes in unselected DCM patients. The similar distribution of missense variants in DCM patients and healthy reference individuals questions the pathogenic role of many variants, and suggests that results from genetic testing of DCM patients should be interpreted with caution., (Copyright © 2016 Elsevier Ireland Ltd. All rights reserved.)

Helminth infections and allergic diseases are associated with IgE hyperresponsiveness but the genetics of this phenotype remain to be defined. Susceptibility to Ascaris lumbricoides infection and antibody levels to this helminth are associated with polymorphisms in locus 13q33-34. We aimed to explore this and other genomic regions to identify genetic variants associated with the IgE responsiveness in humans. Forty-eight subjects from Cartagena, Colombia, with extreme values of specific IgE to Ascaris and ABA-1, a resistance marker of this nematode, were selected for targeted resequencing. Burden analyses were done comparing extreme groups for IgE values. One-hundred one SNPs were genotyped in 1258 individuals of two well-characterized populations from Colombia and Sweden. Two low-frequency coding variants in the gene encoding the Acidic Mammalian Chitinase (CHIA rs79500525, rs139812869, tagged by rs10494133) were found enriched in high IgE responders to ABA-1 and confirmed by genetic association analyses. The SNP rs4950928 in the Chitinase 3 Like 1 gene (CHI3L1) was associated with high IgE to ABA-1 in Colombians and with high IgE to Bet v 1 in the Swedish population. CHIA rs10494133 and ABDH13 rs3783118 were associated with IgE responses to Ascaris. SNPs in the Tumor Necrosis Factor Superfamily Member 13b gene (TNFSF13B) encoding the cytokine B cell activating Factor were associated with high levels of total IgE in both populations. This is the first report on the association between low-frequency and common variants in the chitinases-related genes CHIA and CHI3L1 with the intensity of specific IgE to ABA-1 in a population naturally exposed to Ascaris and with Bet v 1 in a Swedish population. Our results add new information about the genetic influences of human IgE responsiveness; since the genes encode for enzymes involved in the immune response to parasitic infections, they could be helpful for understanding helminth immunity and allergic responses. We also confirmed that TNFSF13B has an important and conserved role in the regulation of total IgE levels, which supports potential evolutionary links between helminth immunity and allergic response., Competing Interests: RV has received research grants from the European Union, by Biomay AG, Vienna, Austria, Thermofisher, Uppsala, Sweden and Fresenius Medical Care, Bad Homburg, Germany. He serves as a consultant for Biomay AG, Vienna, Austria, Thermofisher, Uppsala, Sweden and Fresenius Medical Care, Bad Homburg, Germany. There are no patents, products in development or marketed products to declare. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials.

We applied a targeted sequencing approach to identify germline mutations conferring a moderately to highly increased risk of cutaneous and uveal melanoma. Ninety-two high-risk melanoma patients were screened for inherited variation in 120 melanoma candidate genes. Observed gene variants were filtered based on frequency in reference populations, cosegregation with melanoma in families and predicted functional effect. Several novel or rare genetic variants in genes involved in DNA damage response, cell-cycle regulation and transcriptional control were identified in melanoma patients. Among identified genetic alterations was an extremely rare variant (minor allele frequency of 0.00008) in the BRIP1 gene that was found to cosegregate with the melanoma phenotype. We also found a rare nonsense variant in the BRCA2 gene (rs11571833), previously associated with cancer susceptibility but not with melanoma, which showed weak association with melanoma susceptibility in the Swedish population. Our results add to the growing knowledge about genetic factors associated with melanoma susceptibility and also emphasize the role of DNA damage response as an important factor in melanoma etiology. © 2016 Wiley Periodicals, Inc., (© 2016 Wiley Periodicals, Inc.)

Refractory anaemia with ring sideroblasts (RARS) is distinguished by hyperplastic inefficient erythropoiesis, aberrant mitochondrial ferritin accumulation and anaemia. Heterozygous mutations in the spliceosome gene SF3B1 are found in a majority of RARS cases. To explore the link between SF3B1 mutations and anaemia, we studied mutated RARS CD34(+) marrow cells with regard to transcriptome sequencing, splice patterns and mutational allele burden during erythroid differentiation. Transcriptome profiling during early erythroid differentiation revealed a marked up-regulation of genes involved in haemoglobin synthesis and in the oxidative phosphorylation process, and down-regulation of mitochondrial ABC transporters compared to normal bone marrow. Moreover, mis-splicing of genes involved in transcription regulation, particularly haemoglobin synthesis, was confirmed, indicating a compromised haemoglobinization during RARS erythropoiesis. In order to define the phase during which erythroid maturation of SF3B1 mutated cells is most affected, we assessed allele burden during erythroid differentiation in vitro and in vivo and found that SF3B1 mutated erythroblasts showed stable expansion until late erythroblast stage but that terminal maturation to reticulocytes was significantly reduced. In conclusion, SF3B1 mutated RARS progenitors display impaired splicing with potential downstream consequences for genes of key importance for haemoglobin synthesis and terminal erythroid differentiation., (© 2015 The Authors. British Journal of Haematology published by John Wiley & Sons Ltd.)

Infertility is a worldwide concern that can be treated with in vitro fertilization (IVF). Improvements in IVF and infertility treatment depend largely on better understanding of the molecular mechanisms for human preimplantation development. Several large-scale studies have been conducted to identify gene expression patterns for the first five days of human development, and many functional studies utilize mouse as a model system. We have identified genes of possible importance for this time period by analyzing human microarray data and available data from online databases. We selected 70 candidate genes for human preimplantation development and investigated their expression in the early mouse development from oocyte to the 8-cell stage. Maternally loaded genes expectedly decreased in expression during development both in human and mouse. We discovered that 25 significantly upregulated genes after fertilization in human included 13 genes whose orthologs in mouse behaved differently and mimicked the expression profile of maternally expressed genes. Our findings highlight many significant differences in gene expression patterns during mouse and human preimplantation development. We also describe four cancer-testis antigen families that are also highly expressed in human embryos: PRAME, SSX, GAGE and MAGEA.

Background: The transition from fertilized egg to embryo is accompanied by a multitude of changes in gene expression, and the transcriptional events that underlie these processes have not yet been fully characterized. In this study RNA-Seq is used to compare the transcription profiles of four early developmental stages in zebrafish (Danio rerio) on a global scale., Results: An average of 79 M total reads were detected from the different stages. Out of the total number of reads 65% - 73% reads were successfully mapped and 36% - 44% out of those were uniquely mapped. The total number of detected unique gene transcripts was 11187, of which 10096 were present at 1-cell stage. The largest number of common transcripts was observed between 1-cell stage and 16-cell stage. An enrichment of gene transcripts with molecular functions of DNA binding, protein folding and processing as well as metal ion binding was observed with progression of development. The sequence data (accession number ERP000635) is available at the European Nucleotide Archive., Conclusion: Clustering of expression profiles shows that a majority of the detected gene transcripts are present at steady levels, and thus a minority of the gene transcripts clusters as increasing or decreasing in expression over the four investigated developmental stages. The three earliest developmental stages were similar when comparing highly expressed genes, whereas the 50% epiboly stage differed from the other three stages in the identity of highly expressed genes, number of uniquely expressed genes and enrichment of GO molecular functions. Taken together, these observations indicate a major transition in gene regulation and transcriptional activity taking place between the 512-cell and 50% epiboly stages, in accordance with previous studies.

Recent studies on chromosome conformation show that chromosomes colocalize in the nucleus, bringing together active genes in transcription factories. This spatial proximity of actively transcribing genes could provide a means for RNA interaction at the transcript level. We have screened public databases for chimeric EST and mRNA sequences with the intent of mapping transcription-induced interchromosomal interactions. We suggest that chimeric transcripts may be the result of close encounters of active genes, either as functional products or "noise" in the transcription process, and that they could be used as probes for chromosome interactions. We have found a total of 5,614 chimeric ESTs and 587 chimeric mRNAs that meet our selection criteria. Due to their higher quality, the mRNA findings are of particular interest and we hope that they may serve as food for thought for specialists in diverse areas of molecular biology.

Trees present a life form of paramount importance for terrestrial ecosystems and human societies because of their ecological structure and physiological function and provision of energy and industrial materials. The genus Populus is the internationally accepted model for molecular tree biology. We have analyzed 102,019 Populus ESTs that clustered into 11,885 clusters and 12,759 singletons. We also provide >4,000 assembled full clone sequences to serve as a basis for the upcoming annotation of the Populus genome sequence. A public web-based EST database (POPULUSDB) provides digital expression profiles for 18 tissues that comprise the majority of differentiated organs. The coding content of Populus and Arabidopsis genomes shows very high similarity, indicating that differences between these annual and perennial angiosperm life forms result primarily from differences in gene regulation. The high similarity between Populus and Arabidopsis will allow studies of Populus to directly benefit from the detailed functional genomic information generated for Arabidopsis, enabling detailed insights into tree development and adaptation. These data will also valuable for functional genomic efforts in Arabidopsis.

The human genome sequence defines our inherent biological potential; the realization of the biology encoded therein requires knowledge of the function of each gene. Currently, our knowledge in this area is still limited. Several lines of investigation have been used to elucidate the structure and function of the genes in the human genome. Even so, gene prediction remains a difficult task, as the varieties of transcripts of a gene may vary to a great extent. We thus performed an exhaustive integrative characterization of 41,118 full-length cDNAs that capture the gene transcripts as complete functional cassettes, providing an unequivocal report of structural and functional diversity at the gene level. Our international collaboration has validated 21,037 human gene candidates by analysis of high-quality full-length cDNA clones through curation using unified criteria. This led to the identification of 5,155 new gene candidates. It also manifested the most reliable way to control the quality of the cDNA clones. We have developed a human gene database, called the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/). It provides the following: integrative annotation of human genes, description of gene structures, details of novel alternative splicing isoforms, non-protein-coding RNAs, functional domains, subcellular localizations, metabolic pathways, predictions of protein three-dimensional structure, mapping of known single nucleotide polymorphisms (SNPs), identification of polymorphic microsatellite repeats within human genes, and comparative results with mouse full-length cDNAs. The H-InvDB analysis has shown that up to 4% of the human genome sequence (National Center for Biotechnology Information build 34 assembly) may contain misassembled or missing regions. We found that 6.5% of the human gene candidates (1,377 loci) did not have a good protein-coding open reading frame, of which 296 loci are strong candidates for non-protein-coding RNA genes. In addition, among 72,027 uniquely mapped SNPs and insertions/deletions localized within human genes, 13,215 nonsynonymous SNPs, 315 nonsense SNPs, and 452 indels occurred in coding regions. Together with 25 polymorphic microsatellite repeats present in coding regions, they may alter protein structure, causing phenotypic effects or resulting in disease. The H-InvDB platform represents a substantial contribution to resources needed for the exploration of human biology and pathology., Competing Interests: The authors have declared that no conflicts of interest exist.