The Genetic Architecture of a Congenital Heart Defect Is Related to Its Fitness Cost.

In newborns, severe congenital heart defects are rarer than mild ones. This epidemiological relationship between heart defect severity and incidence lacks explanation. Here, an analysis of ~10,000 Nkx2-5+/- mice from two inbred strain crosses illustrates the fundamental role of epistasis. Modifier genes raise or lower the risk of specific defects via pairwise (G×GNkx) and higher-order (G×G×GNkx) interactions with Nkx2-5. Their effect sizes correlate with the severity of a defect. The risk loci for mild, atrial septal defects exert predominantly small G×GNkx effects, while the loci for severe, atrioventricular septal defects exert large G×GNkx and G×G×GNkx effects. The loci for moderately severe ventricular septal defects have intermediate effects. Interestingly, G×G×GNkx effects are three times more likely to suppress risk when the genotypes at the first two loci are from the same rather than different parental inbred strains. This suggests the genetic coadaptation of interacting G×G×GNkx loci, a phenomenon that Dobzhansky first described in Drosophila. Thus, epistasis plays dual roles in the pathogenesis of congenital heart disease and the robustness of cardiac development. The empirical results suggest a relationship between the fitness cost and genetic architecture of a disease phenotype and a means for phenotypic robustness to have evolved.

Assessment of Metastatic Colorectal Cancer (CRC) Tissues for Interpreting Genetic Data in Forensic Science by Applying 16 STR Loci among Saudi Patients.

BACKGROUND: In forensic science, there are cases when the only available provider of biological data is samples of malignant tissues. It can be useful in identification and/or paternity tests. Still, such samples have ambiguities because of microsatellite instability (MSI) and loss of heterozygosity (LOH) effects, being often related to neoplasia. METHODS: This research evaluates 16 autosomal short tandem repeat (STR) loci (traditional in forensic investigations) to get genetic data. MSI and LOH were estimated in DNA patterns derived from 73 Saudi respondents (30 healthy individuals and 43 persons with diagnosed colorectal cancer (CRC). Upon deriving DNA from blood, CRC specimens were obtained in both groups, along with the adjoining normal non-cancerous tissues (N-CRC). All specimens and 16 loci (15 STR loci and Amelogenin) were evaluated. Moreover, both colorectal samples were histologically analyzed utilizing HandE staining. RESULTS: Findings revealed non-essential variability in genetic information because of MSI and/or LOH. In CRC, mutations rates were 0.42% (MSI) and 1.62% (LOH). In N-CRC, mutation rates were 0.00% (MSI) and 0.59% (LOH). Further, LOH-related deviations were recorded in 5 loci out of 16. MSI-related deviations were recorded in 4 out of 16 loci, being present in CRC samples only. Genetic deviations within the marker loci might inform about false homozygosity/heterozygosity. Similarly, false gender might come from improper interpretation of DNA profiles. Finally, histopathological trials showed considerable histopathological alterations contrasted to N-CRC. CONCLUSION: This study is unique in demonstrating the application of 16 autosomal STRs from CRC samples and their comparison with the adjoining N-CRCs in Saudi participants, contributing to the field of forensic science. The experiment revealed no considerable distinctions, while showing that cancer tissues might display MSI and LOH effects that might challenge data interpretation, if STRs are to be applied in the forensic investigation.

Genetic Diversity Assessment of MPOB-Senegal Oil Palm Germplasm Using Microsatellite Markers.

Molecular characterization of oil palm germplasm is crucial in utilizing and conserving germplasm with promising traits. This study was conducted to evaluate the genetic diversity structures and relationships among 26 families of MPOB-Senegal oil palm germplasm using thirty-five microsatellite markers. High level of polymorphism (P = 96.26%), number of effective allele (N e = 2.653), observed heterozygosity (H o = 0.584), expected heterozygosity (H e = 0.550), total heterozygosity (H T = 0.666), and rare alleles (54) were observed which indicates that MPOB-Senegal germplasm has a broad genetic variation. Among the SSR markers, sMo00053 and sMg00133 were the most informative markers for discrimination among the MPOB-Senegal oil palm germplasm for having the highest private alleles and the rare alleles. For selection and conservation, oil palm populations with high rare alleles and Nei's gene diversity index should be considered as these populations may possess unique genes for further exploitation.

Multiple Immunoglobulin κ Gene Rearrangements within a Single Clone Unraveled by Next-Generation Sequencing-Based Clonality Assessment.

Clonality assessment of the Ig heavy- and light-chain genes (IGH and IGK) using GeneScan analysis is an important supplemental assay in diagnostic testing for lymphoma. Occasionally cases with an IGK rearrangement pattern that cannot readily be assigned to a monoclonal lymphoma are encountered, whereas the occurrence of biclonal lymphomas is rare, and the result of the IGH locus of these cases is in line with monoclonality. Three such ambiguous cases were assessed for clonality using next-generation sequencing. Information on the sequences of the rearrangements, combined with knowledge of the complex organization of the IGK locus, pointed to two explanations that can attribute seemingly biclonal IGK rearrangements to a single clone. In two cases, this explanation involved inversion rearrangements on the IGK locus, whereas in the third case, the cross-reactivity of primers generated an additional clonal product. In conclusion, next-generation sequencing-based clonality assessment allows for the detection of both inversion rearrangements and the cross-reactivity of primers, and can therefore facilitate the interpretation of cases of lymphoma with complex IGK rearrangement patterns.

Identification of putative causal loci in whole-genome sequencing data via knockoff statistics.

The analysis of whole-genome sequencing studies is challenging due to the large number of rare variants in noncoding regions and the lack of natural units for testing. We propose a statistical method to detect and localize rare and common risk variants in whole-genome sequencing studies based on a recently developed knockoff framework. It can (1) prioritize causal variants over associations due to linkage disequilibrium thereby improving interpretability; (2) help distinguish the signal due to rare variants from shadow effects of significant common variants nearby; (3) integrate multiple knockoffs for improved power, stability, and reproducibility; and (4) flexibly incorporate state-of-the-art and future association tests to achieve the benefits proposed here. In applications to whole-genome sequencing data from the Alzheimer's Disease Sequencing Project (ADSP) and COPDGene samples from NHLBI Trans-Omics for Precision Medicine (TOPMed) Program we show that our method compared with conventional association tests can lead to substantially more discoveries.

Novel risk factors for craniofacial microsomia and assessment of their utility in clinic diagnosis.

Craniofacial microsomia (CFM, OMIM%164 210) is one of the most common congenital facial abnormalities worldwide, but it's genetic risk factors and environmental threats are poorly investigated, as well as their interaction, making the diagnosis and prenatal screening of CFM impossible. We perform a comprehensive association study on the largest CFM cohort of 6074 samples. We identify 15 significant (P < 5 × 10-8) associated genomic loci (including eight previously reported) and decipher 107 candidates based on multi-omics data. Gene Ontology term enrichment found that these candidates are mainly enriched in neural crest cell (NCC) development and hypoxic environment. Single-cell RNA-seq data of mouse embryo demonstrate that nine of them show dramatic expression change during early cranial NCC development whose dysplasia is involved in pathogeny of CFM. Furthermore, we construct a well-performed CFM risk-predicting model based on polygenic risk score (PRS) method and estimate seven environmental risk factors that interacting with PRS. Single-nucleotide polymorphism-based PRS is significantly associated with CFM [P = 7.22 × 10-58, odds ratio = 3.15, 95% confidence interval (CI) 2.74-3.63], and the top fifth percentile has a 6.8-fold CFM risk comparing with the 10th percentile. Father's smoking increases CFM risk as evidenced by interaction parameter of -0.324 (95% CI -0.578 to -0.070, P = 0.011) with PRS. In conclusion, the newly identified risk loci will significantly improve our understandings of genetics contribution to CFM. The risk prediction model is promising for CFM prediction, and father's smoking is a key environmental risk factor for CFM through interacting with genetic factors.

Genetic evidence of shared ancestry among diverse ethno-linguistic human populations of Himachal Pradesh.

In Indian montane system, human populations often exhibit an unparallel social organization where inter-caste marriages are still not common. This attribute affects the demography and population genetic structure of the resident populations. Further, human populations residing in the mountains in India are poorly studied for their genetic make-up and allele distribution patterns. In the present study, we genotyped 594 unrelated individuals using PowerPlex® 21 System (Promega, USA) from eight different populations belonging to 12 districts of Himachal Pradesh which differed in ethnicity, language, geography and social organization. Altogether, we obtained 1415 alleles with a mean of 8.84 ± 0.26 alleles per locus and 0.80 mean observed heterozygosity. Locus Penta E showed the highest combined power of discrimination and was found most informative for forensic purposes. Interestingly, phylogenetic analysis grouped the populations of Rajputs, Scheduled castes and Brahmins into one cluster, which indicated a deep genetic admixture in the ancestral populations. This study documents the first-ever report on the population genetic assignment of various castes in Himachal Pradesh and unveils the facts of cryptic gene flow among the diverse castes in the northern hilly state of India. Our results showed a genetic relationship among the various ethno-linguistically diverse populations of India.

Low Hydrophobic Mismatch Scores Calculated for HLA-A/B/DR/DQ Loci Improve Kidney Allograft Survival.

We evaluated the impact of human leukocyte antigen (HLA) disparity (immunogenicity; IM) on long-term kidney allograft survival. The IM was quantified based on physicochemical properties of the polymorphic linear donor/recipient HLA amino acids (the Cambridge algorithm) as a hydrophobic, electrostatic, amino acid mismatch scores (HMS\AMS\EMS) or eplet mismatch (EpMM) load. High-resolution HLA-A/B/DRB1/DQB1 types were imputed to calculate HMS for primary/re-transplant recipients of deceased donor transplants. The multiple Cox regression showed the association of HMS with graft survival and other confounders. The HMS integer 0-10 scale showed the most survival benefit between HMS 0 and 3. The Kaplan-Meier analysis showed that: the HMS=0 group had 18.1-year median graft survival, a 5-year benefit over HMS>0 group; HMS ≤ 3.0 had 16.7-year graft survival, a 3.8-year better than HMS>3.0 group; and, HMS ≤ 7.8 had 14.3-year grafts survival, a 1.8-year improvement over HMS>7.8 group. Stratification based on EMS, AMS or EpMM produced similar results. Additionally, the importance of HLA-DR with/without -DQ IM for graft survival was shown. In our simulation of 1,000 random donor/recipient pairs, 75% with HMS>3.0 were re-matched into HMS ≤ 3.0 and the remaining 25% into HMS≥7.8: after re-matching, the 13.5 years graft survival would increase to 16.3 years. This approach matches donors to recipients with low/medium IM donors thus preventing transplants with high IM donors.

Molecular Typing of Different Salmonella Serotypes by Multiple-Locus Analysis of Tandem Repeats.

Salmonella enterica subspecies enterica causes salmonellosis in humans and animals and is an important cause of food infections worldwide. In recent years, the multiple-locus variable-number of tandem repeat (VNTR) analysis (MLVA), a fast molecular typing method with strong epidemiological discrimination, has facilitated the effective control of diverse infections. This study aimed at the typing of 28 human origined Salmonella enteritidis, Salmonella infantis, and Salmonella typhimurium strains by using a single MLVA protocol. Previously these strains have been identified by pulsed field gel electrophoresis (PFGE) method and it has been shown that each strain produced a distinct PFGE banding profile. One MLVA protocol was tested on 3 serotypes simultaneously and it produced three banding patterns specific to each of the three common Salmonella serotypes. MLVA also constitute a relatively more cost-effective and faster method than PFGE.

A high-density SNP-based genetic map and several economic traits-related loci in Pelteobagrus vachelli.

BACKGROUND: A high-density genetic linkage map is essential for QTL fine mapping, comparative genome analysis, identification of candidate genes and marker-assisted selection in aquaculture species. Pelteobagrus vachelli is a very popular commercial species in Asia. However, some specific characters hindered achievement of the traditional selective breeding based on phenotypes, such as lack of large-scale genomic resource and short of markers tightly associated with growth, sex determination and hypoxia tolerance related traits. RESULTS: By making use of 5059 ddRAD markers in P. vachelli, a high-resolution genetic linkage map was successfully constructed. The map' length was 4047.01 cM by using an interval of 0.11 cm, which is an average marker standard. Comparative genome mapping revealed that a high proportion (83.2%) of markers with a one-to-one correspondence were observed between P. vachelli and P. fulvidraco. Based on the genetic map, 8 significant genome-wide QTLs for 4 weight, 1 body proportion, 2 sex determination, and 1 hypoxia tolerance related traits were detected on 4 LGs. Some SNPs from these significant genome-wide QTLs were observably associated with these phenotypic traits in other individuals by Kompetitive Allele Specific PCR. In addition, two candidate genes for weight, Sipa1 and HSD11B2, were differentially expressed between fast-, medium- and slow-growing P. vachelli. Sema7a, associated with hypoxia tolerance, was induced after hypoxia exposure and reoxygenation. CONCLUSIONS: We mapped a set of suggestive and significant QTLs as well as candidate genes for 12 growth, 1 sex determination and 1 hypoxia tolerance related traits based on a high-density genetic linkage map by making use of SNP markers for P. fulvidraco. Our results have offered a valuable method about the much more efficient production of all-male, fast growth and hypoxia tolerance P. vachelli for the aquaculture industry.

DNA methylation study of Huntington's disease and motor progression in patients and in animal models.

Although Huntington's disease (HD) is a well studied Mendelian genetic disorder, less is known about its associated epigenetic changes. Here, we characterize DNA methylation levels in six different tissues from 3 species: a mouse huntingtin (Htt) gene knock-in model, a transgenic HTT sheep model, and humans. Our epigenome-wide association study (EWAS) of human blood reveals that HD mutation status is significantly (p < 10-7) associated with 33 CpG sites, including the HTT gene (p = 6.5 × 10-26). These Htt/HTT associations were replicated in the Q175 Htt knock-in mouse model (p = 6.0 × 10-8) and in the transgenic sheep model (p = 2.4 × 10-88). We define a measure of HD motor score progression among manifest HD cases based on multiple clinical assessments. EWAS of motor progression in manifest HD cases exhibits significant (p < 10-7) associations with methylation levels at three loci: near PEX14 (p = 9.3 × 10-9), GRIK4 (p = 3.0 × 10-8), and COX4I2 (p = 6.5 × 10-8). We conclude that HD is accompanied by profound changes of DNA methylation levels in three mammalian species.

Detecting and Quantifying Natural Selection at Two Linked Loci from Time Series Data of Allele Frequencies with Forward-in-Time Simulations.

Recent advances in DNA sequencing techniques have made it possible to monitor genomes in great detail over time. This improvement provides an opportunity for us to study natural selection based on time serial samples of genomes while accounting for genetic recombination effect and local linkage information. Such time series genomic data allow for more accurate estimation of population genetic parameters and hypothesis testing on the recent action of natural selection. In this work, we develop a novel Bayesian statistical framework for inferring natural selection at a pair of linked loci by capitalising on the temporal aspect of DNA data with the additional flexibility of modeling the sampled chromosomes that contain unknown alleles. Our approach is built on a hidden Markov model where the underlying process is a two-locus Wright-Fisher diffusion with selection, which enables us to explicitly model genetic recombination and local linkage. The posterior probability distribution for selection coefficients is computed by applying the particle marginal Metropolis-Hastings algorithm, which allows us to efficiently calculate the likelihood. We evaluate the performance of our Bayesian inference procedure through extensive simulations, showing that our approach can deliver accurate estimates of selection coefficients, and the addition of genetic recombination and local linkage brings about significant improvement in the inference of natural selection. We also illustrate the utility of our method on real data with an application to ancient DNA data associated with white spotting patterns in horses.

The assessment of epigenetic diversity, differentiation, and structure in the 'Fuji' mutation line implicates roles of epigenetic modification in the occurrence of different mutant groups as well as spontaneous mutants.

The 'Fuji' line includes many varieties with a similar genetic background and consistent inducement factors with epigenetic occurrence, thus it may be considered an ideal candidate for epigenetic research. In this study, 91 bud mutations of 'Fuji' apple were used as the test materials. Using the genetic variation within 'Fuji' as the control, the characteristics of epigenetic variation at different levels in both varieties and mutant groups were examined. The results showed that: (1) the global genomic DNA methylation level of the 91 bud mutants of 'Fuji' ranged from 29.120%-45.084%, with an average of 35.910%. Internal cytosine methylation was the main DNA methylation pattern. Regarding the variation of methylation patterns of 'Fuji' mutants, the vast majority of loci maintained the original methylation pattern existed in 'Fuji'. CHG methylation variation was the main type of variation; (2) the variation in methylation patterns between the mutant groups was greater than that of methylation levels. Among these patterns, the variation in CHG methylation patterns (including CHG hypermethylation and CHG demethylation) was expected to be dominant. The observed variation in methylation levels was more important in the Color mutant group; however, the variation in methylation patterns was more obvious in both the early maturation and Spur mutant groups. Moreover, the range of variation in the Early-maturation group was much wider than that in the Spur mutant group; (3) epigenetic diversity and genetic diversity were both low between the mutant groups. In the 'Fuji' mutant groups, there was few correlation between genetic and epigenetic variation, and epigenetic differentiation resulted in more loci with moderate or greater differentiation; (4) the purifying selection seemed to play a major role in the differentiation of different groups of 'Fuji' mutants (65.618%), but epigenetic diversity selection still occurred at nearly 35% of loci. Sixteen epigenetic outlier loci were detected.

Assessment of candidate variants causative of inborn metabolic diseases in SUDI cases in South Africa, and a case report.

Sudden unexpected death in infants (SUDI) is a devastating event, and unfortunately is still a burden in many parts of the world, including in South Africa. Due to the absence of routine testing for inborn metabolic diseases in newborns and in a post-mortem context, little is known about the presence of metabolic diseases in local SUDI cases. The aim of this study was to genotype five candidate variants previously associated with metabolic disorders in a cohort of SUDI cases (n = 169) from Salt River Mortuary, Cape Town. DNA was isolated from blood, and SNaPshot® PCR and Sanger sequencing were used to genotype the following variants: ACADM: c.583G > A, ACADM: c.985A > G, GCDH: c.877G > A/T, GALT: c.404C > G/T and GALT: c.563A > G. Four carriers of GCDH: c.877G > A/T were identified, while one infant was homozygous for the founder mutation GALT: c.404C > G/T; the latter which is causative of galactosaemia and was previously undiagnosed. During the follow-up with the family, it emerged that the affected infant's identical twin had subsequently demised. The findings in this study highlight possible new candidate variants to assess in South African SUDI cases, and these results directly contribute to the development of a molecular autopsy which is locally relevant. It is evident that until newborn screening becomes routine and accessible in South Africa, molecular autopsies should include testing for inherited metabolic disorders, as it holds potential to save lives.

Long-Molecule Sequencing: A New Approach for Identification of Clinically Significant DNA Variants in α-Thalassemia and ß-Thalassemia Carriers.

Multiple molecular tests are currently needed for accurate carrier testing for thalassemia. Therefore, long-molecule sequencing (LMS) was evaluated as an alternate on the PacBio Sequel platform for genotyping carriers of α-thalassemia or ß-thalassemia. Multiplex long PCR was used to generate representative amplicons for the α (HBA1/2) and ß (HBB) gene loci. Following LMS, circular consensus sequencing reads were aligned to the hg19 reference genome and variants called using FreeBayes software version 1.2.0. In a blinded study of 64 known carrier samples, all HBA1/2 and HBB variants detected by LMS were concordant with those independently assigned by targeted PCR assays. For HBA1/2 carrier samples, LMS accurately detected the common South East Asian, -α3.7, and -α4.2 deletions and four different rare single-nucleotide variants (SNVs). For HBB carrier samples, LMS accurately detected the most common Chinese insertion and deletion variant c.126_129delCTTT and 14 different SNVs/insertions and deletions and could discriminate compound heterozygous SNVs (trans configuration) and identify variants linked to benign SNPs (cis configuration). Overall, LMS displayed the hallmarks of a scalable, accurate, and cost-effective genotyping method. With further test coverage to additionally include detection of other clinically significant HBA1/2 copy number variations, such as the Thai, Mediterranean, and Filipino deletions, LMS may eventually serve as a comprehensive method for large-scale thalassemia carrier screening.

Systematic meta-analyses, field synopsis and global assessment of the evidence of genetic association studies in colorectal cancer.

OBJECTIVE: To provide an understanding of the role of common genetic variations in colorectal cancer (CRC) risk, we report an updated field synopsis and comprehensive assessment of evidence to catalogue all genetic markers for CRC (CRCgene2). DESIGN: We included 869 publications after parallel literature review and extracted data for 1063 polymorphisms in 303 different genes. Meta-analyses were performed for 308 single nucleotide polymorphisms (SNPs) in 158 different genes with at least three independent studies available for analysis. Scottish, Canadian and Spanish data from genome-wide association studies (GWASs) were incorporated for the meta-analyses of 132 SNPs. To assess and classify the credibility of the associations, we applied the Venice criteria and Bayesian False-Discovery Probability (BFDP). Genetic associations classified as 'positive' and 'less-credible positive' were further validated in three large GWAS consortia conducted in populations of European origin. RESULTS: We initially identified 18 independent variants at 16 loci that were classified as 'positive' polymorphisms for their highly credible associations with CRC risk and 59 variants at 49 loci that were classified as 'less-credible positive' SNPs; 72.2% of the 'positive' SNPs were successfully replicated in three large GWASs and the ones that were not replicated were downgraded to 'less-credible' positive (reducing the 'positive' variants to 14 at 11 loci). For the remaining 231 variants, which were previously reported, our meta-analyses found no evidence to support their associations with CRC risk. CONCLUSION: The CRCgene2 database provides an updated list of genetic variants related to CRC risk by using harmonised methods to assess their credibility.

Assessment of three essential tremor genetic loci in sporadic Parkinson's disease in Eastern China.

OBJECTIVE: The aim of this study was to investigate potential genetic overlap between essential tremor and Parkinson's disease in a cohort of 825 subjects from an Eastern Chinese population. METHODS: A total of 441 Parkinson's disease patients and 384 healthy controls were recruited. The MassARRAY System was used to detect three essential tremor-related single nucleotide polymorphisms. Odds ratio (OR) and 95% confidential interval (CI) were calculated to assess the relationship between polymorphisms and Parkinson's disease susceptibility. RESULTS: Our results demonstrated that the odds ratios of rs3794087 of SLC1A2, rs9652490 of LINGO1, and rs17590046 of PPARGC1A were 0.71 (95% CI = 0.55-0.91), 0.99 (95% CI = 0.78-1.26), and 0.88 (95% CI = 0.62-1.25), respectively. CONCLUSION: An essential tremor SNP (rs3794087 of SLC1A2) is associated with a decreased risk of PD in the Eastern Han Chinese population, while rs9652490 (LINGO1) and rs17590046 (PPARGC1A) do not show an association.

Assessment of genetic polymorphism at 20 autosomal STR loci in scheduled caste population of Himachal Pradesh.

In this study, we evaluated the genetic diversity indices and forensic parameters of scheduled caste population of Himachal Pradesh, India, at 20 autosomal STR loci. Altogether, 233 alleles were observed with an average of 11.65 ± 0.88 alleles per locus. In exception of locus D13S317, no loci deviated from Hardy Weinberg equilibrium. Locus Penta E was found to be the most polymorphic and discriminative loci. The combined power of discrimination and the combined power of exclusion were 1 and 0.999. Further, we established the phylogenetic relationship between the scheduled caste population of Himachal Pradesh and 14 other populations of India. We found 20 autosomal STR loci used in the present study were polymorphic and can be used in population genetic studies and forensic related case works.

Cost-Effective Profiling of Mutator Transposon Insertions in Maize by Next-Generation Sequencing.

Transposable elements can be highly mutagenic because when they transpose they can insert into genes and disrupt their function, a propensity which has been exploited in many organisms to generate tagged mutant alleles. The Mutator (Mu) family transposon is a family of DNA-type transposons in maize with a particularly high duplication frequency, which results in large numbers of new mutations in lineages that carry active Mu elements. Here we describe a rapid and cost-effective Miseq-based Mu transposon profiling pipeline. This method can also be used for identifying flanking sequences of other types of long insertions such as T-DNAs.