Rank-invariant estimation of inbreeding coefficients

The two alleles an individual carries at a locus are identical by descent (ibd) if they have descended from a single ancestral allele in a reference population, and the probability of such identity is the inbreeding coefficient of the individual. Inbreeding coefficients can be predicted from pedigrees with founders constituting the reference population, but estimation from genetic data is not possible without data from the reference population. Most inbreeding estimators that make explicit use of sample allele frequencies as estimates of allele probabilities in the reference population are confounded by average kinships with other individuals. This means that the ranking of those estimates depends on the scope of the study sample and we show the variation in rankings for common estimators applied to different subdivisions of 1000 Genomes data. Allele-sharing estimators of within-population inbreeding relative to average kinship in a study sample, however, do have invariant rankings across all studies including those individuals. They are unbiased with a large number of SNPs. We discuss how allele sharing estimates are the relevant quantities for a range of empirical applications.

AFA: Computationally efficient Ancestral Frequency estimation in Admixed populations: the Hispanic Community Health Study/Study of Latinos

We developed a computationally efficient method, Ancestral Frequency estimation in Admixed populations (AFA), to estimate the frequencies of bi-allelic variants in admixed populations with an unlimited number of ancestries. AFA uses maximum likelihood estimation by modeling the conditional probability of having an allele given proportions of genetic ancestries. It can be applied using either global or local proportions of genetic ancestries. Simulations mimicking admixture demonstrated the high accuracy of the method. We implemented the method on data from the Hispanic Community Health Study/Study of Latinos (HCHS/SOL), an admixed population with three predominant continental ancestries: Amerindian, European, and African. Comparison of the European and African estimated frequencies to the respective gnomAD frequencies demonstrated high correlations, with Pearson R2=0.97-0.99. We provide a genome-wide dataset of the estimated three ancestral allele frequencies in HCHS/SOL for all available variants with allele frequency between 5%-95% in at least one of the three ancestral populations.

CYP3A5 and UGT1A9 Polymorphisms Influence Immunosuppressive Therapy in Pediatric Kidney Transplant Recipients

Background: Tacrolimus (TAC) and mycophenolic acid (MPA) are the main immunosuppressive drugs used in pediatric kidney transplantation. Single nucleotide polymorphisms (SNPs) in metabolizing enzymes and transporters might influence plasma levels of these drugs. Herein, we sought to determine the influence of SNPs on CYP3A5, MRP2 and UGT1A9 genes in Chilean pediatric kidney recipients using TAC and MPA.Patients and Methods: A prospective study was performed on 104 pediatric kidney recipients that used TAC and MPA for immunosuppression. The median age at the time of transplantation was 8.1 years [Q1–Q3 4.5–11.6 years] and the main clinical diagnosis was a structural anomaly. In a subgroup of patients, a complete steroid withdrawal was made at day 7. The CYP3A5 polymorphism (ancestral allele *1; variant allele *3) was determined in the entire cohort, while MRP2 -24G > A, UGT1A9 -275T > A, and UGT1A9 -2152C > T polymorphisms were determined in 53 patients. Genotypes were associated with trough drug concentrations (C0), dose requirements normalized by weight (TAC-D mg/kg) or body surface (MPA-D mg/m2), trough levels normalized by dose requirements (C0/D), and area under the curve in 12 h normalized by dose requirements (AUC0–12h/D).Results: The frequencies of the variant alleles CYP3A5*3, MRP2-24A, UGT1A9-275A, and UGT1A9-2152T were 76.9, 22.1, 6.6, and 2.9%, respectively. AUC0–12h/TAC-D were 1.6-fold higher in CYP3A5*3/*3 patients than in CYP3A5*1 carriers (CYP3A5*1/*3 and CYP3A5*1/*1). When analyzing patients with steroid withdrawal, CYP3A5*3/*3 patients had 1.7-fold higher AUC0–12h/TAC-D than the other genotypes. Patients carrying the CYP3A5*3/*3 genotype had higher TAC-C0, lower TAC-D and higher TAC-C0/D, consistently in a 6-months follow-up. Creatinine clearance was stable during the follow-up, regardless of the genotype. No significant differences between MRP2 and UGT1A9 genotypes were observed in MPA-C0, MPA-D or MPA-C0/D. However, patients carrying the UGT1A9-275A allele had lower AUC0–12h/MPA-D than those carrying the UGT1A9-275T ancestral allele.Conclusions: These results support that CYP3A5 and UGT1A9 genotyping in pediatric recipients might be useful and advisable to guide TAC and MPA dosing and monitoring in children that undergo kidney transplantation.

The reversion variant (p.Arg90Leu) at the evolutionarily adaptive p.Arg90 site in CELA3B predisposes to chronic pancreatitis

ABSTRACTA gain-of-function missense variant in the CELA3B gene, p.Arg90Cys (c.268C>T), has recently been reported to cause pancreatitis in an extended pedigree. Herein, we sequenced the CELA3B gene in 644 genetically unexplained French chronic pancreatitis (CP) patients (all unrelated) and 566 controls. No predicted loss-of-function variants were identified. None of the six low frequency or common missense variants detected showed significant association with CP. Nor did the aggregate rare/very rare missense variants (n=14) show any significant association with CP. However, p.Arg90Leu (c.269G>T), which was found in 4 patients but no controls and affects the same amino acid as p.Arg90Cys, serves to revert p.Arg90 to the human elastase ancestral allele. Since p.Arg90Leu has previously been shown to exert a similar functional effect to p.Arg90Cys, our findings not only confirm the involvement of CELA3B in the etiology of CP but also pinpoint a new evolutionarily adaptive site in the human genome.

Wxlv, the Ancestral Allele of Rice Waxy Gene

In rice endosperms, the Waxy (Wx) gene is important for amylose synthesis, and various Wx alleles control the amylose content and affect the taste of cooked rice. Herein, we report the cloning of the ancestral allele Wxlv of the Wx locus, which affects the mouthfeel of rice grains by modulating the size of amylose molecules. Using evolutionary analysis, we demonstrated that Wxlv originated directly from wild rice, and the three major Wx alleles in cultivated rice (Wxb, Wxa, and Wxin) differentiated after the substitution of one base pair at the functional sites. These data indicate that the Wxlv allele played an important role in artificial selection and domestication. The findings also shed light on the evolution of various Wx alleles, which have greatly contributed to improving the eating and cooking quality of rice.

APOL1 variant-expressing endothelial cells exhibit autophagic dysfunction and mitochondrial stress

Apolipoprotein L1 (APOL1) gene risk variants (RV) associate with renal and cardiovascular disease particularly in SLE. We hypothesized that in RV-carrying human umbilical vein endothelial cells (HUVECs) cytokine-induced APOL1 expression compromises mitochondrial respiration, lysosome integrity, and autophagic flux. HUVEC cultures of each APOL1 genotype were generated. APOL1 was expressed using IFNɣ; HUVEC mitochondrial function, lysosome integrity, and autophagic flux were measured. IFNɣ increased APOL1 expression across all genotypes 20-fold (p=0.001). Compared to the homozygous G0 (ancestral) allele (0RV), high risk (2RV) HUVECs showed both depressed baseline and maximum mitochondrial oxygen consumption (p<0.01), and impaired mitochondrial networking on MitoTracker assays. These cells also demonstrated a contracted lysosome compartment (p<0.001), and an accumulation of autophagosomes suggesting a defect in autophagic flux. Treatment of 0RV HUVECs with a non-selective lysosome inhibitor, hydroxychloroquine, produced autophagosome accumulations similar to the 2RV cells, thus implicating lysosome dysfunction in blocking autophagy. Compared to 0RV and 2RV HUVECs, 1 RV cells demonstrated an intermediate autophagy defect which was exacerbated by IFNɣ. Our findings implicate dysfunction of mitochondrial respiration, lysosome, and autophagy in APOL1 RV-mediated endothelial cytotoxicity. IFNɣ amplified this phenotype even in variant heterozygous cells–a potential underpin of the APOL1/inflammation interaction. This is the first description of APOL1 pathobiology in variant heterozygous cell cultures.

How do species barriers decay? Concordance and local introgression in mosaic hybrid zones of mussels

The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study “replicated” instances of secondary contact between closely-related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry-informative panel of such SNPs. We then compared their frequencies in newly-sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bistable variants (Dobzhansky-Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact.