dgfr: an R package to assess sequence diversity of gene families.

BACKGROUND: Gene families are groups of homologous genes that often have similar biological functions. These families are formed by gene duplication events throughout evolution, resulting in multiple copies of an ancestral gene. Over time, these copies can acquire mutations and structural variations, resulting in members that may vary in size, motif ordering and sequence. Multigene families have been described in a broad range of organisms, from single-celled bacteria to complex multicellular organisms, and have been linked to an array of phenomena, such as host-pathogen interactions, immune evasion and embryonic development. Despite the importance of gene families, few approaches have been developed for estimating and graphically visualizing their diversity patterns and expression profiles in genome-wide studies. RESULTS: Here, we introduce an R package named dgfr, which estimates and enables the visualization of sequence divergence within gene families, as well as the visualization of secondary data such as gene expression. The package takes as input a multi-fasta file containing the coding sequences (CDS) or amino acid sequences from a multigene family, performs a pairwise alignment among all sequences, and estimates their distance, which is subjected to dimension reduction, optimal cluster determination, and gene assignment to each cluster. The result is a dataset that allows for the visualization of sequence divergence and expression within the gene family, an approximation of the number of clusters present in the family. CONCLUSIONS: dgfr provides a way to estimate and study the diversity of gene families, as well as visualize the dispersion and secondary profile of the sequences. The dgfr package is available at https://github.com/lailaviana/dgfr under the GPL-3 license.

Genome sequencing of captive white tigers from Bangladesh.

OBJECTIVES: The Bengal tiger Panthera tigris tigris, is an emblematic animal for Bangladesh. Despite being the apex predator in the wild, their number is decreasing due to anthropogenic activities such as hunting, urbanization, expansion of agriculture and deforestation. By contrast, captive tigers are flourishing due to practical conservation efforts. Breeding within the small captive population can produce inbreeding depression and genetic bottlenecks, which may limit the success of conservation efforts. Despite past decades of research, a comprehensive database on genetic variation in the captive and wild Bengal tigers in Bangladesh still needs to be included. Therefore, this research aimed to investigate the White Bengal tiger genome to create a resource for future studies to understand variation underlying important functional traits. DATA DESCRIPTION: Blood samples from Chattogram Zoo were collected for three white Bengal tigers. Genomic DNA for all collected samples were extracted using a commercial DNA extraction kit. Whole genome sequencing was performed using a DNBseq platform. We generated 77 Gb of whole-genome sequencing (WGS) data for three white Bengal tigers (Average 11X coverage/sample). The data we generated will establish a paradigm for tiger research in Bangladesh by providing a genomic resource for future functional studies on the Bengal white tiger.

Analysis of population structure and genetic diversity of Camellia tachangensis in Guizhou based on SNP markers.

BACKGROUND: Camellia tachangensis F. C. Zhang is a five-compartment species in the ovary of tea group plants, which represents the original germline of early differentiation of some tea group plants. METHODS AND RESULTS: In this study, we analyzed single-nucleotide polymorphisms (SNPs) at the genome level, constructed a phylogenetic tree, analyzed the genetic diversity, and further investigated the population structure of 100 C. tachangensis accessions using the genotyping-by-sequencing (GBS) method. A total of 91,959 high-quality SNPs were obtained. Population structure analysis showed that the 100 C. tachangensis accessions clustered into three groups: YQ-1 (Village Group), YQ-2 (Forest Group) and YQ-3 (Transition Group), which was further consistent with the results of phylogenetic analysis and principal component analyses (PCA). In addition, a comparative analysis of the genetic diversity among the three populations (Forest, Village, and Transition Groups) detected the highest genetic diversity in the Transition Group and the highest differentiation between Forest and Village Groups. CONCLUSIONS: C. tachangensis plants growing in the forest had different genetic backgrounds from those growing in villages. This study provides a basis for the effective protection and utilization of C. tachangensis populations and lays a foundation for future C. tachangensis breeding.

Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes.

BACKGROUND: Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes. METHODS: Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis. RESULTS: The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%. CONCLUSION: The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.

Identification of novel and known genetic variants associated with hereditary hearing loss in iranian families using whole exome sequencing.

BACKGROUND: Hearing loss (HL) is a common sensory impairment worldwide, with genetic and environmental factors contributing to its occurrence. Next Generation Sequencing (NGS) plays a crucial role in identifying the genetic factors involved in this heterogeneous disorder. METHODS AND RESULTS: In this study, a total of 9 unrelated Iranian families, each having at least one affected individual who tested negative for mutations in GJB2, underwent screening using whole exome sequencing (WES). The pathogenicity and novelty of the identified variant was checked using various databases. Co-segregation study was also performed to confirm the presence of the candidate variants in parents. Plus, The pathogenicity of the detected variant was assessed through in silico analysis using a number of mutation prediction software tools. Among the 9 investigated families, hearing loss-causing genes were identified in 6 families. the mutations were observed in USH2A, CLRN1, BSND, SLC26A4, and MITF, with two of the identified mutations being novel. CONCLUSION: Discovering additional variants and broadening the range of mutations associated with hearing impairment has the potential to enhance the diagnostic effectiveness of molecular testing in patient screening, and can also lead to improved counseling aimed at reducing the risk of affected offspring for high-risk couples.

Molecular diversity of black corals from the Saudi Arabian Red Sea: a first assessment.

Black corals occur as part of benthic assemblages from shallow to deep waters in all oceans. Despite the importance in many benthic ecosystems, where these act as biodiversity aggregators, antipatharians remain poorly studied, with 75% of the known species occurring below recreational SCUBA diving depth limits. Currently, information regarding the diversity and evolutionary history is limited, with most studies focusing on Hawaii and the South Pacific Ocean. Other regions of the world have received less attention, such as the Red Sea, where only two black coral families and four genera have been recorded. We provide the first analysis of the molecular diversity of black corals in the eastern Gulf of Aqaba and the northern and central Saudi Arabian Red Sea, based on a dataset of 161 antipatharian colonies collected down to 627 m deep. Based on specimen morphology, we ascribed our material to 11 genera belonging to 4 of the 7 known Antipatharia families, i.e. Antipathidae, Aphanipathidae, Myriopathidae and Schizopathidae. The genus level phylogeny of three intergenic mitochondrial regions, the trnW-IGR-nad2 (IgrW ), nad5-IGR-nad1 (IgrN ) and cox3-IGR-cox1 was reconstructed including previously published material. Overall, we recovered six molecular clades that included exclusively Red Sea sequences, with the highest diversity occurring at mesophotic depths. This study highlights that diversity of black corals in the Red Sea is much higher than previously known, with seven new generic records, suggesting that this basin may be a hotspot for antipatharian diversity as is known for other taxa. Our results recovered unresolved relationships within the order at the familial and generic levels. This emphasises the urgent need for an integration of genomic-wide data with a re-examination of informative morphological features necessary to revise the systematics of the order at all taxonomic levels.

Analyzing the functional effects of DNA variants with gene editing.

Continual advancements in genomics have led to an ever-widening disparity between the rate of discovery of genetic variants and our current understanding of their functions and potential roles in disease. Systematic methods for phenotyping DNA variants are required to effectively translate genomics data into improved outcomes for patients with genetic diseases. To make the biggest impact, these approaches must be scalable and accurate, faithfully reflect disease biology, and define complex disease mechanisms. We compare current methods to analyze the function of variants in their endogenous DNA context using genome editing strategies, such as saturation genome editing, base editing and prime editing. We discuss how these technologies can be linked to high-content readouts to gain deep mechanistic insights into variant effects. Finally, we highlight key challenges that need to be addressed to bridge the genotype to phenotype gap, and ultimately improve the diagnosis and treatment of genetic diseases.

Alzheimer's diseases in America, Europe, and Asian regions: a global genetic variation.

Background: Alzheimer's disease (AD) is one of the multifaceted neurodegenerative diseases influenced by many genetic and epigenetic factors. Genetic factors are merely not responsible for developing AD in the whole population. The studies of genetic variants can provide significant insights into the molecular basis of Alzheimer's disease. Our research aimed to show how genetic variants interact with environmental influences in different parts of the world. Methodology: We searched PubMed and Google Scholar for articles exploring the relationship between genetic variations and global regions such as America, Europe, and Asia. We aimed to identify common genetic variations susceptible to AD and have no significant heterogeneity. To achieve this, we analyzed 35 single-nucleotide polymorphisms (SNPs) from 17 genes (ABCA7, APOE, BIN1, CD2AP, CD33, CLU, CR1, EPHA1, TOMM40, MS4A6A, ARID5B, SORL1, APOC1, MTHFD1L, BDNF, TFAM, and PICALM) from different regions based on previous genomic studies of AD. It has been reported that rs3865444, CD33, is the most common polymorphism in the American and European populations. From TOMM40 and APOE rs2075650, rs429358, and rs6656401, CR1 is the common investigational polymorphism in the Asian population. Conclusion: The results of all the research conducted on AD have consistently shown a correlation between genetic variations and the incidence of AD in the populations of each region. This review is expected to be of immense value in future genetic research and precision medicine on AD, as it provides a comprehensive understanding of the genetic factors contributing to the development of this debilitating disease.

Meta-analysis of genomic variants in power and endurance sports to decode the impact of genomics on athletic performance and success.

Association between genomic variants and athletic performance has seen a high degree of controversy, as there is often conflicting data as far as the association of genomic variants with endurance, speed and strength is concerned. Here, findings from a thorough meta-analysis from 4228 articles exploring the association of genomic variants with athletic performance in power and endurance sports are summarized, aiming to confirm or overrule the association of genetic variants with athletic performance of all types. From the 4228 articles, only 107 were eligible for further analysis, including 37 different genes. From these, there were 21 articles for the ACE gene, 29 articles for the ACTN3 gene and 8 articles for both the ACE and ACTN3 genes, including 54,382 subjects in total, from which 11,501 were endurance and power athletes and 42,881 control subjects. These data show that there is no statistically significant association between genomic variants and athletic performance either for endurance or power sports, underlying the fact that it is highly risky and even unethical to make such genetic testing services for athletic performance available to the general public. Overall, a strict regulatory monitoring should be exercised by health and other legislative authorities to protect the public from such services from an emerging discipline that still lacks the necessary scientific evidence and subsequent regulatory approval.

Conserved DNA sequence analysis reveals the phylogeography and evolutionary events of Akebia trifoliata in the region across the eastern edge of the Tibetan Plateau and subtropical China.

BACKGROUND: The eastern edge of the Qinghai‒Tibet Plateau (QTP) and subtropical China have various regions where plant species originate and thrive, but these regions have been the focus of very few integrative studies. Here, we elucidated the phylogeographic structure of a continuous and widespread Akebia trifoliata population across these two regions. RESULTS: Sixty-one populations consisting of 391 genotypes were examined to assess population diversity and structure via network distribution analysis, maximum likelihood phylogenetic tree reconstruction, divergence time estimation, demographic history inference, and ancestral area reconstruction of both conserved internal transcribed spacer (ITS) and chloroplast (rps16) DNA sequences. The results showed that the ITS region was more variable than the rps16 region and could be suitable for studying intraspecific phylogeography. The A. trifoliata population displayed high genetic diversity, genetic differentiation and obvious phylogeographical structure, possibly originating on the eastern QTP, expanding during the last glacial-interglacial cycle, diverging in the early Pleistocene and middle Pleistocene, and extensively migrating thereafter. The migration route from west to east along rivers could be largely responsible for the long-distance dispersal of this species, while three main refuges (Qinba Mountains, Nanling Mountains and Yunnan-Guizhou Plateau) with multiple ice shelters facilitated its wide distribution. CONCLUSIONS: Our results suggested that the from west to east long migration accompanying with the minor short reciprocal migration in the south-north direction, and the three main refuges (the Qinba Mountains, Nanling Mountains and Yunnan-Guizhou Plateau) contributed to the extant geographical distribution of A. trifoliata. In addition, this finding also strongly reduced the discrepancy between glacial contraction and postglacial expansion and the in situ survival hypothesis by simultaneously considering the existence of many similar climate-related ecological niches and migration influences.

Rare genetic variation in fibronectin 1 (FN1) protects against APOEε4 in Alzheimer's disease.

The risk of developing Alzheimer's disease (AD) significantly increases in individuals carrying the APOEε4 allele. Elderly cognitively healthy individuals with APOEε4 also exist, suggesting the presence of cellular mechanisms that counteract the pathological effects of APOEε4; however, these mechanisms are unknown. We hypothesized that APOEε4 carriers without dementia might carry genetic variations that could protect them from developing APOEε4-mediated AD pathology. To test this, we leveraged whole-genome sequencing (WGS) data in the National Institute on Aging Alzheimer's Disease Family Based Study (NIA-AD FBS), Washington Heights/Inwood Columbia Aging Project (WHICAP), and Estudio Familiar de Influencia Genetica en Alzheimer (EFIGA) cohorts and identified potentially protective variants segregating exclusively among unaffected APOEε4 carriers. In homozygous unaffected carriers above 70 years old, we identified 510 rare coding variants. Pathway analysis of the genes harboring these variants showed significant enrichment in extracellular matrix (ECM)-related processes, suggesting protective effects of functional modifications in ECM proteins. We prioritized two genes that were highly represented in the ECM-related gene ontology terms, (FN1) and collagen type VI alpha 2 chain (COL6A2) and are known to be expressed at the blood-brain barrier (BBB), for postmortem validation and in vivo functional studies. An independent analysis in a large cohort of 7185 APOEε4 homozygous carriers found that rs140926439 variant in FN1 was protective of AD (OR = 0.29; 95% CI [0.11, 0.78], P = 0.014) and delayed age at onset of disease by 3.37 years (95% CI [0.42, 6.32], P = 0.025). The FN1 and COL6A2 protein levels were increased at the BBB in APOEε4 carriers with AD. Brain expression of cognitively unaffected homozygous APOEε4 carriers had significantly lower FN1 deposition and less reactive gliosis compared to homozygous APOEε4 carriers with AD, suggesting that FN1 might be a downstream driver of APOEε4-mediated AD-related pathology and cognitive decline. To validate our findings, we used zebrafish models with loss-of-function (LOF) mutations in fn1b-the ortholog for human FN1. We found that fibronectin LOF reduced gliosis, enhanced gliovascular remodeling, and potentiated the microglial response, suggesting that pathological accumulation of FN1 could impair toxic protein clearance, which is ameliorated with FN1 LOF. Our study suggests that vascular deposition of FN1 is related to the pathogenicity of APOEε4, and LOF variants in FN1 may reduce APOEε4-related AD risk, providing novel clues to potential therapeutic interventions targeting the ECM to mitigate AD risk.

Kernel density estimation of allele frequency including undetected alleles.

Whereas undetected species contribute to estimation of species diversity, undetected alleles have not been used to estimated genetic diversity. Although random sampling guarantees unbiased estimation of allele frequency and genetic diversity measures, using undetected alleles may provide biased but more precise estimators useful for conservation. We newly devised kernel density estimation (KDE) for allele frequency including undetected alleles and tested it in estimation of allele frequency and nucleotide diversity using population generated by coalescent simulation as well as well as real population data. Contrary to expectations, nucleotide diversity estimated by KDE had worse bias and accuracy. Allele frequency estimated by KDE was also worse except when the sample size was small. These might be due to finity of population and/or the curse of dimensionality. In conclusion, KDE of allele frequency does not contribute to genetic diversity estimation.

Thymus algeriensis Boiss. et Reut., a north African endemic plant species: genetic diversity and population structure as assessed by molecular markers, a pioneer step for conservation implications.

BACKGROUND: Thymus algeriensis Boiss. et Reut. is one of the most widespread North African species of the genus Thymus L. The species is subshrub growing primarily in subtropical biome of Morocco, Algeria, Tunisia and Libya. In Tunisia, the plant species is under high pressure of anthropogenic activities including over-collecting. The assessment of genetic diversity and population structure of T. algeriensis is a pioneer step to retrace its evolutionary history and to perform appropriate conservation strategies of the plant species. METHODS AND RESULTS: Seven wild populations growing, widely, in different bioclimatic zones were selected and analysed using two molecular markers systems. Fifteen Simple Sequence Repeats (SSRs) and fifteen Inter-Simple Sequence Repeats (ISSRs) were used to characterize genetically 140 different genotypes. The results showed a high molecular variation within populations and among the studied genotypes. The intra-populations genetic diversity revealed by SSRs was higher (P = 80.95%, Na = 2.143 and He = 0.364) than that based on ISSRs (P = 78.12%, Na = 1.632, He = 0.265 and I = 0.398). As demonstrated by inbreeding coefficients, a significant level of differentiation and a low level of gene flow were detected among studied populations (FST = 0.161 for SSRs and ΦST = 0.197 for ISSRs). Furthermore, the results of ISSRs marker suggest land strips as barriers in population genetic structure. While SSRs marker reflects a relatively structured bioclimatic patterns of studied populations. The Bayesian analysis showed a specific adaptation of populations to local environments. CONCLUSIONS: The used molecular markers (ISSRs and SSRs) seem to be effective in deciphering genetic polymorphism of Tunisian genotypes of T. algeriensis. Therefore, the genetic structure of the studied genotypes could constitute a starting point for further conservation, characterization and breeding programs.

Discovery of a peculiar insular race of Ravenna nivea (Nire, 1920) (Lepidoptera: Lycaenidae) endemic to Yinggeling Mountain of Hainan, suggesting heterogeneous geological history of mountain formation of the island.

A peculiar population of Ravenna nivea (Nire, 1920) was discovered from the Yinggeling Mountain Mass of central Hainan. Its wing pattern and COI barcode data show considerable distinction from other geographic populations of R. nivea, including that of Bawangling, approximately only 40 km away and also located in Hainan. The p-distance value of the COI barcode between the Yinggeling and Bawangling populations was 1.1%, considerably higher than the value (0.6%) between Bawangling population and populations in eastern China, where the subspecific name howarthi Saigusa, 1993 applies. The population is regarded as a distinct subspecies ngiunmoiae Lo & Hsu, subsp. nov. The distinctness and high degree of COI haplotype diversity of R. nivea found in Hainan and Taiwan suggest continental islands may serve as glacial refugees for the butterfly and other organisms during previous glaciations, and the presence of the relict populations of montane butterflies like R. nivea may provide useful clues towards a better understanding of the geological history of mountain formation within islands.

Applicability of SCoT markers in unraveling genetic variation and population structure among sugar beet (Beta vulgaris L.) germplasm.

BACKGROUND: Sugar beet (Beta vulgaris L.) holds significant importance as a crop globally cultivated for sugar production. The genetic diversity present in sugar beet accessions plays a crucial role in crop improvement programs. METHODS AND RESULTS: During the present study, we collected 96 sugar beet accessions from different regions and extracted DNA from their leaves. Genomic DNA was amplified using SCoT primers, and the resulting fragments were separated by gel electrophoresis. The data were analyzed using various genetic diversity indices, and constructed a population STRUCTURE, applied the unweighted pair-group method with arithmetic mean (UPGMA), and conducted Principle Coordinate Analysis (PCoA). The results revealed a high level of genetic diversity among the sugar beet accessions, with 265 bands produced by the 10 SCoT primers used. The percentage of polymorphic bands was 97.60%, indicating substantial genetic variation. The study uncovered significant genetic variation, leading to higher values for overall gene diversity (0.21), genetic distance (0.517), number of effective alleles (1.36), Shannon's information index (0.33), and polymorphism information contents (0.239). The analysis of molecular variance suggested a considerable amount of genetic variation, with 89% existing within the population. Using STRUCTURE and UPGMA analysis, the sugar beet germplasm was divided into two major populations. Structure analysis partitioned the germplasm based on the origin and domestication history of sugar beet, resulting in neighboring countries clustering together. CONCLUSION: The utilization of SCoT markers unveiled a noteworthy degree of genetic variation within the sugar beet germplasm in this study. These findings can be used in future breeding programs with the objective of enhancing both sugar beet yield and quality.

Phenotypic variation seems not to be associated with the genetic profile in Zygopetalum (Orchidaceae): a case study of a high-elevation rocky complex.

BACKGROUND: Hybridization associated with polyploidy studies is rare in the tropics. The genus Zygopetalum (Orchidaceae) was investigated here as a case study of Neotropical plants. In the rocky highlands of the Ibitipoca State Park (ISP), southeast Brazil, individuals with intermediate colors and forms between the species Z. maculatum and Z. triste were commonly identified. METHODS AND RESULTS: Chromosomal analysis and DNA quantity showed a uniform population. Regardless of the aspects related to the color and shape of floral structures, all individuals showed 2n = 96 chromosomes and an average of 14.05 pg of DNA. Irregularities in meiosis associated with chromosome number and C value suggest the occurrence of polyploidy. The genetic distance estimated using ISSR molecular markers revealed the existence of genetic variability not related to morphological clusters. Morphometric measurements of the flower pieces revealed that Z. maculatum shows higher variation than Z. triste although lacking a defined circumscription. CONCLUSION: The observed variation can be explained by the polyploid and phenotypic plasticity resulting from the interaction of the genotypes with the heterogeneous environments observed in this habitat.

Genetic divergence accessed with microsatellite markers reflects the time of Crassostrea gigas genetic breeding in Brazil.

The Pacific Oyster was introduced on Santa Catarina Island in 1987, experiencing processes of selection and genetic breeding since then. Such procedures may have led to the establishment of specific strains, given the saltier and warmer conditions of the Atlantic Ocean. This study employed microsatellite markers to compare allelic patterns of oysters cultivated in Santa Catarina, the USA, and Asia. Specific allelic patterns were revealed in the Santa Catarina samples, reflecting the time of selection/breeding of the oyster in this region. This result supports the effectiveness of the selection/breeding procedures and the demand for protection of this commercially important genetic resource.

Optimizing Bangkaew dog breed identification using DNA technology.

BACKGROUND: The Bangkaew dog is an indigenous dog breed in the Phitsanulok province of Thailand. This breed is recognized by the Fédération Cynologique Internationale (FCI), a global canine organization. The unique traits of the Bangkaew breed lead to purebred selection for breeding, while only their traits and pedigree from parental history are recorded. Determination of the risk of inbreeding depression and the origin of unknown DNA profiles is essential due to the challenges in predicting puppy characteristics, which are crucial for breed management and conservation. OBJECTIVE: This study aimed to emphasize that current allelic frequency data for the Bangkaew dog breed must be considered for precise individual identification. METHODS: Approximately 82 Bangkaew dogs from various Thai localities were studied using 15 microsatellite markers for genotypic monitoring and individual identification. Maternal genetic inheritance was assessed via mtDNA D-loop analysis. RESULTS: The results revealed high genetic diversity in the Bangkaew breed, indicating low potential for inbreeding. We also found that using a 15 loci microsatellite panel was effective for the identification of Bangkaew dogs. The optimized 10 loci microsatellite genotyping panel developed in this study presents improved identification testing efficiency, promoting both time- and cost-effectiveness. CONCLUSION: Analysis of microsatellite DNA markers in Bangkaew dogs using an optimized panel of 10 loci selected from 15 loci effectively facilitated individual identification. This approach not only enhances time and cost efficiency, but also provides accurate allelic frequency estimates, which are crucial for the realistic evaluation of DNA evidence.

Critical assessment of variant prioritization methods for rare disease diagnosis within the rare genomes project.

BACKGROUND: A major obstacle faced by families with rare diseases is obtaining a genetic diagnosis. The average "diagnostic odyssey" lasts over five years and causal variants are identified in under 50%, even when capturing variants genome-wide. To aid in the interpretation and prioritization of the vast number of variants detected, computational methods are proliferating. Knowing which tools are most effective remains unclear. To evaluate the performance of computational methods, and to encourage innovation in method development, we designed a Critical Assessment of Genome Interpretation (CAGI) community challenge to place variant prioritization models head-to-head in a real-life clinical diagnostic setting. METHODS: We utilized genome sequencing (GS) data from families sequenced in the Rare Genomes Project (RGP), a direct-to-participant research study on the utility of GS for rare disease diagnosis and gene discovery. Challenge predictors were provided with a dataset of variant calls and phenotype terms from 175 RGP individuals (65 families), including 35 solved training set families with causal variants specified, and 30 unlabeled test set families (14 solved, 16 unsolved). We tasked teams to identify causal variants in as many families as possible. Predictors submitted variant predictions with estimated probability of causal relationship (EPCR) values. Model performance was determined by two metrics, a weighted score based on the rank position of causal variants, and the maximum F-measure, based on precision and recall of causal variants across all EPCR values. RESULTS: Sixteen teams submitted predictions from 52 models, some with manual review incorporated. Top performers recalled causal variants in up to 13 of 14 solved families within the top 5 ranked variants. Newly discovered diagnostic variants were returned to two previously unsolved families following confirmatory RNA sequencing, and two novel disease gene candidates were entered into Matchmaker Exchange. In one example, RNA sequencing demonstrated aberrant splicing due to a deep intronic indel in ASNS, identified in trans with a frameshift variant in an unsolved proband with phenotypes consistent with asparagine synthetase deficiency. CONCLUSIONS: Model methodology and performance was highly variable. Models weighing call quality, allele frequency, predicted deleteriousness, segregation, and phenotype were effective in identifying causal variants, and models open to phenotype expansion and non-coding variants were able to capture more difficult diagnoses and discover new diagnoses. Overall, computational models can significantly aid variant prioritization. For use in diagnostics, detailed review and conservative assessment of prioritized variants against established criteria is needed.

Rare GPR37L1 Variants Reveal Potential Association between GPR37L1 and Disorders of Anxiety and Migraine.

GPR37L1 is an orphan receptor that couples through heterotrimeric G-proteins to regulate physiological functions. Since its role in humans is not fully defined, we used an unbiased computational approach to assess the clinical significance of rare G-protein-coupled receptor 37-like 1 (GPR37L1) genetic variants found among 51,289 whole-exome sequences from the DiscovEHR cohort. Rare GPR37L1 coding variants were binned according to predicted pathogenicity and analyzed by sequence kernel association testing to reveal significant associations with disease diagnostic codes for epilepsy and migraine, among others. Since associations do not prove causality, rare GPR37L1 variants were functionally analyzed in SK-N-MC cells to evaluate potential signaling differences and pathogenicity. Notably, receptor variants exhibited varying abilities to reduce cAMP levels, activate mitogen-activated protein kinase (MAPK) signaling, and/or upregulate receptor expression in response to the agonist prosaptide (TX14(A)), as compared with the wild-type receptor. In addition to signaling changes, knock-out (KO) of GPR37L1 or expression of certain rare variants altered cellular cholesterol levels, which were also acutely regulated by administration of the agonist TX14(A) via activation of the MAPK pathway. Finally, to simulate the impact of rare nonsense variants found in the large patient cohort, a KO mouse line lacking Gpr37l1 was generated. Although KO animals did not recapitulate an acute migraine phenotype, the loss of this receptor produced sex-specific changes in anxiety-related disorders often seen in chronic migraineurs. Collectively, these observations define the existence of rare GPR37L1 variants associated with neuropsychiatric conditions in the human population and identify the signaling changes contributing to pathological processes.